Feeds:
Posts
Comments

Posts Tagged ‘metabolites’

The Reconstruction of Life Processes requires both Genomics and Metabolomics to explain Phenotypes and Phylogenetics

Writer and Curator: Larry H. Bernstein, MD, FCAP 

 

phylogenetics

phylogenetics

http://upload.wikimedia.org/wikipedia/commons/thumb/1/12/CollapsedtreeLabels-simplified.svg/200px-CollapsedtreeLabels-simplified.svg.png

 

This discussion that completes and is an epicrisis (summary and critical evaluation) of the series of discussions that preceded it.

  1. Innervation of Heart and Heart Rate
  2. Action of hormones on the circulation
  3. Allogeneic Transfusion Reactions
  4. Graft-versus Host reaction
  5. Unique problems of perinatal period
  6. High altitude sickness
  7. Deep water adaptation
  8. Heart-Lung-and Kidney
  9. Acute Lung Injury

The concept inherent in this series is that the genetic code is an imprint that is translated into a message.  It is much the same as a blueprint, or a darkroom photographic image that has to be converted to a print. It is biologically an innovation of evolutionary nature because it establishes a simple and reproducible standard for the transcription of the message through the transcription of the message using strings of nucleotides (oligonucleotides) that systematically transfer the message through ribonucleotides that communicate in the cytoplasm with the cytoskeleton based endoplasmic reticulum (ER), composing a primary amino acid sequence.  This process is a quite simple and convenient method of biological activity.  However, the simplicity ends at this step.  The metabolic components of the cell are organelles consisting of lipoprotein membranes and a cytosol which have particularly aligned active proteins, as in the inner membrane of the mitochondrion, or as in the liposome or phagosome, or the structure of the  ER, each of which is critical for energy transduction and respiration, in particular, for the mitochondria, cellular remodeling or cell death, with respect to the phagosome, and construction of proteins with respect to the ER, and anaerobic glycolysis and the hexose monophosphate shunt in the cytoplasmic domain.  All of this refers to structure and function, not to leave out the membrane assigned transport of inorganic, and organic ions (electrolytes and metabolites).

I have identified a specific role of the ER, the organelles, and cellular transactions within and between cells that is orchestrated.  But what I have outlined is a somewhat limited and rigid model that does not reach into the dynamics of cellular transactions.  The DNA has expression that may be old, no longer used messages, and this is perhaps only part of a significant portion of “dark matter”.  There is also nuclear DNA that is enmeshed with protein, mRNA that is a copy of DNA, and mDNA  is copied to ribosomal RNA (rRNA).  There is also rDNA. The classic model is DNA to RNA to protein.  However, there is also noncoding RNA, which plays an important role in regulation of transcription.

This has been discussed in other articles.  But the important point is that proteins have secondary structure through disulfide bonds, which is determined by position of sulfur amino acids, and by van der Waal forces, attraction and repulsion. They have tertiary structure, which is critical for 3-D structure.  When like subunits associate, or dissimilar oligomers, then you have heterodimers and oligomers.  These constructs that have emerged over time interact with metabolites within the cell, and also have an important interaction with the extracellular environment.

When you take this into consideration then a more complete picture emerges. The primitive cell or the multicellular organism lives in an environment that has the following characteristics – air composition, water and salinity, natural habitat, temperature, exposure to radiation, availability of nutrients, and exposure to chemical toxins or to predators.  In addition, there is a time dimension that proceeds from embryonic stage to birth in mammals, a rapid growth phase, a tapering, and a decline.  The time span is determined by body size, fluidity of adaptation, and environmental factors.  This is covered in great detail in this work.  The last two pieces are in the writing stage that completes the series. Much content has already be presented in previous articles.

The function of the heart, kidneys and metabolism of stressful conditions have already been extensively covered in http://pharmaceuticalintelligence.com  in the following and more:

The Amazing Structure and Adaptive Functioning of the Kidneys: Nitric Oxide – Part I

http://pharmaceuticalintelligence.com/2012/11/26/the-amazing-structure-and-adaptive-functioning-of-the-kidneys/

Nitric Oxide and iNOS have Key Roles in Kidney Diseases – Part II

http://pharmaceuticalintelligence.com/2012/11/26/nitric-oxide-and-inos-have-key-roles-in-kidney-diseases/

The pathological role of IL-18Rα in renal ischemia/reperfusion injury – Nature.com

http://pharmaceuticalintelligence.com/2014/10/24/the-pathological-role-of-il-18r%CE%B1-in-renal-ischemiareperfusion-injury-nature-com/

Summary, Metabolic Pathways

http://pharmaceuticalintelligence.com/2014/10/23/summary-metabolic-pathways/

 

Read Full Post »

Summary of Transcription, Translation ond Transcription Factors

Author and Curator:  Larry H. Bernstein, MD, FCAP  

 

Proteins are integral to the composition of the cytoskeleton, and also to the extracellular matrix.  Many proteins are actually enzymes, carrying out the transformation of some substrate, a derivative of the food we ingest.  They have a catalytic site, and they function with a cofactor – either a multivalent metal or a nucleotide. Proteins also are critically involved in the regulation of cell metabolism, and they are involved in translation of the DNA code, as they make up transcription factors (TFs). There are 20 essential amino acids that go into protein synthesis that are derived from animal or plant protein.   Protein synthesis is carried out by the transport of mRNA out of the nucleus to the ribosome, where tRNA is paired with a matching amino acid, and the primary sequence of a protein is constructed as a linear string of amino acids.

This is illustrated in the following three pictures:

protein synthesis

protein synthesis

mcell-transcription-translation

mcell-transcription-translation

transcription_translation

transcription_translation

Proteins synthesized at distal locations frequently contain intrinsically disordered segments. These regions are generally rich in assembly-promoting modules and are often regulated by post-translational modifications. Such proteins are tightly regulated but display distinct temporal dynamics upon stimulation with growth factors. Thus, proteins synthesized on-site may rapidly alter proteome composition and act as dynamically regulated scaffolds to promote the formation of reversible cellular assemblies.
RJ Weatheritt, et al. Nature Structural & Molecular Biology 24 Aug, 2014; 21: 833–839 http://dx.do.orgi:/10.1038/nsmb.2876

An overview of the potential advantages conferred by distal-site protein synthesis

An overview of the potential advantages conferred by distal-site protein synthesis

Turquoise and red filled circle represents off-target and correct interaction partners, respectively. Wavy lines represent a disordered region within a distal site synthesis protein. Grey and red line in graphs represents profiles of t…  http://www.nature.com/nsmb/journal/v21/n9/carousel/nsmb.2876-F5.jpg

In the the transcription process an RNA sequence is read.  This is essential for protein synthesis through the ordering of the amino acids in the primary structure. However, there are microRNAs and noncoding RNAs, and there are transcription factors.  The transcription factors bind to chromatin, and the RNAs also have some role in regulating the transcription process. (see picture above)

Transcription factors (TFs) interact dynamically in vivo with chromatin binding sites. Four different techniques are currently used to measure their kinetics in live cells,

  1. fluorescence recovery after photobleaching (FRAP),
  2. fluorescence correlation spectroscopy (FCS),
  3. single molecule tracking (SMT) and
  4. competition ChIP (CC).

A comparison of data from each of these techniques raises an important question:

  • do measured transcription kinetics reflect biologically functional interactions at specific sites (i.e. working TFs) or
  • do they reflect non-specific interactions (i.e. playing TFs)?

There are five key unresolved biological questions related to

  • the functionality of transient and prolonged binding events at both
  • specific promoter response elements as well as non-specific sites.

In support of functionality,

  • there are data suggesting that TF residence times are tightly regulated, and
  • that this regulation modulates transcriptional output at single genes.

In addition to this site-specific regulatory role, TF residence times

  • also determine the fraction of promoter targets occupied within a cell
  • thereby impacting the functional status of cellular gene networks.
  • TF residence times, then, are key parameters that could influence transcription in multiple ways.

Quantifying transcription factor kinetics: At work or at play? Mueller F., et al.  http://dx.doi.org:/10.3109/10409238.2013.833891

Dr. Virginie Mattot works in the team “Angiogenesis, endothelium activation and Cancer” directed by Dr. Fabrice Soncin at the Institut de Biologie de Lille in France where she studies the roles played by microRNAs in endothelial cells during physiological and pathological processes such as angiogenesis or endothelium activation. She has been using Target Site Blockers to investigate the role of microRNAs on putative targets.

A few years ago, the team identified

  • an endothelial cell-specific gene which
  • harbors a microRNA in its intronic sequence.

They have since been working on understanding the functions of

  • both this new gene and its intronic microRNA in endothelial cells.

While they were searching for the functions of the intronic microRNA,

  • theye identified an unknown gene as a putative target.

The aim of my project was to investigate if this unknown gene was actually a genuine target and

  • if regulation of this gene by the microRNA was involved in endothelial cell function.

They had already shown the endothelial cell phenotype is associated with the inhibition of the intronic microRNA.
They then used miRCURY LNA™ Target Site Blockers to demonstrate

  • the expression of this unknown gene is actually controlled by this microRNA.
  • the microRNA regulates specific endothelial cell properties through regulation of this unknown gene.

MicroRNA function in endothelial cells – Solving the mystery of an unknown target gene using Target Site Blockers to investigate the role of microRNAs on putative targets

We first verified that this TSB was functional by analyzing

  • the expression of the miRNA target against which the TSB was directed
  • we then showed the TSB induced similar phenotypes as those when we inhibited the microRNA in the same cells.

Target Site Blockers were shown to be efficient tools to demonstrate the specific involvement of

  • putative microRNA targets
  • in the function played by this microRNA.

Some genes are known to have several different alternatively spliced protein variants, but the Scripps Research Institute’s Paul Schimmel and his colleagues have uncovered almost 250 protein splice variants of an essential, evolutionarily conserved family of human genes. The results were published July 17 in Science.

Focusing on the 20-gene family of aminoacyl tRNA synthetases (AARSs),

  • the team captured AARS transcripts from human tissues—some fetal, some adult—and showed that
  • many of these messenger RNAs (mRNAs) were translated into proteins.

Previous studies have identified several splice variants of these enzymes that have novel functions, but uncovering so many more variants was unexpected, Schimmel said. Most of these new protein products

  • lack the catalytic domain but retain other AARS non-catalytic functional domains.

This study fundamentally effects how we view protein-synthesis, according to  Michael Ibba (who was not involved in the work), The Scientist reported. “The unexpected and potentially vast expanded functional networks that emerge from this study have the potential to influence virtually any aspect of cell growth.”

The team—comprehensively captured and sequenced the AARS mRNAs from six human tissue types using high-throughput deep sequencing. They next showed that a proportion of these transcripts, including those missing the catalytic domain, indeed resulted in stable protein products:

  • 48 of these splice variants associated with polysomes.

In vitro translation assays and the expression of more than 100 of these variants in cells confirmed that

  • many of these variants could be made into stable protein products.

The AARS enzymes—of which there’s one for each of the 20 amino acids—bring together an amino acid with its appropriate transfer RNA (tRNA) molecule. This reaction allows a ribosome to add the amino acid to a growing peptide chain during protein translation. AARS enzymes can be found in all living organisms and are thought to be among the first proteins to have originated on Earth.

One goal of human genetics is to understand how the information for precise and dynamic gene expression programs is encoded in the genome. The interactions of transcription factors (TFs) with DNA regulatory elements clearly

  • play an important role in determining gene expression outputs, yet
  • the regulatory logic underlying functional transcription factor binding is poorly understood.

An important question in genomics is to understand how a class of proteins called ‘‘transcription factors’’ controls the expression level of other genes in the genome in a cell type-specific manner – a process that is essential to human development. One major approach to this problem is to study where these transcription factors bind in the genome, but this does not tell us about the effect of that binding on gene expression levels and

  • it is generally accepted that much of the binding does not strongly influence gene expression.

DA Cusanovich et al. PLoS Genet 2014;10(3):e1004226.  http://dx.doi.org:/10.1371/journal.pgen.1004226

We knocked down 59 TFs and chromatin modifiers in one HapMap lymphoblastoid cell line

  • to evaluate the context of functional TF binding.

We then identified genes whose expression was affected by the knockdowns

  • by intersecting the gene expression data with transcription factor binding data
    (based on ChIP-seq and DNase-seq)
  • within 10 kb of the transcription start sites of expressed genes.

This combination of data allowed us to infer functional TF binding.
Only a small subset of genes bound by a factor were

  • differentially expressed following the knockdown of that factor,
  • suggesting that most interactions between TF and chromatin
  • do not result in measurable changes in gene expression levels
  • of putative target genes.

We found that functional TF binding is enriched

  • in regulatory elements that harbor a large number of TF binding sites,
  • at sites with predicted higher binding affinity, and
  • at sites that are enriched in genomic regions annotated as ‘‘active enhancers.’’

We aim to be able to predict the expression pattern of a gene based on its regulatory
sequence alone.

Combining a TF knockdown approach with TF binding data can help us to

  • distinguish functional binding from non-functional binding

This approach has previously been applied to the study of human TFs, although for the most part studies have only focused on

  • the regulatory relationship of a single factor with its downstream targets.

The FANTOM consortium knocked down 52 different transcription factors in

  • the THP-1 cell line, an acute monocytic leukemia-derived cell line, and
  • used a subset of these to validate certain regulatory predictions based on binding motif enrichments.

We and others previously studied the regulatory architecture of gene expression in

  • the model system of HapMap lymphoblastoid cell lines (LCLs) using both
  • binding map strategies and QTL mapping strategies.

We now sought to use knockdown experiments targeting transcription factors in a HapMap LCL

  • to refine our understanding of the gene regulatory circuitry of the human genome.

Therefore, We integrated the results of the knockdown experiments with previous data on TF binding to

  • better characterize the regulatory targets of 59 different factors and
  • to learn when a disruption in transcription factor binding
  • is most likely to be associated with variation in the expression level of a nearby gene.

Gene expression levels following the knockdown were compared to

  • expression data collected from six samples that were transfected with negative control siRNA.

Depending on the factor targeted, the knockdowns resulted in

  • between 39 and 3,892 differentially expressed genes at an FDR of 5%
    (Figure 1B; see Table S3 for a summary of the results).

The knockdown efficiency for the 59 factors ranged

  • from 50% to 90% (based on qPCR; Table S1).

The qPCR measurements of the knockdown level were significantly

  • correlated with estimates of the TF expression levels
  • based on the microarray data (P =0.001; Figure 1C).

 

Did the factors tended to have a consistent effect (either up- or down-regulation)

  • on the expression levels of genes they purportedly regulated?

All factors we tested are associated with both up- and down-regulation of downstream targets (Figure 6).

While there is compelling evidence for our inferences, the current chromatin functional annotations

  • do not fully explain the regulatory effects of the knockdown experiments.

For example, the enrichments for binding in ‘‘strong enhancer’’ regions of the genome range from 7.2% to 50.1% (median = 19.2%),

  • much beyond what is expected by chance alone, but far from accounting for all functional binding.

A slight majority of downstream target genes were expressed at higher levels

  • following the knockdown for 15 of the 29 factors for which we had binding information (Figure 6B).

The factor that is associated with the largest fraction (68.8%) of up-regulated target genes following the knockdown is EZH2,

  • the enzymatic component of the Polycomb group complex.

On the other end of the spectrum was JUND, a member of the AP-1 complex, for which

  • 66.7% of differentially expressed targets were down-regulated following the knockdown.

Our results, combined with the previous work from our group and others make for a complicated view

  • of the role of transcription factors in gene regulation as
  • it seems difficult to reconcile the inference from previous work that
  • many transcription factors should primarily act as activators with the results presented here.

One somewhat complicated hypothesis, which nevertheless can resolve the apparent discrepancy, is that

  • the ‘‘repressive’’ effects we observe for known activators may be
  • at sites in which the activator is acting as a weak enhancer of transcription and
  • that reducing the cellular concentration of the factor
  • releases the regulatory region to binding by an alternative, stronger activator.

Integrative study of Arabidopsis thaliana metabolomic and transcriptomic data
with the interactiveMarVis-Graph software

M Landesfeind, A Kaever, K Feussner, C Thurow, C Gatz, I Feussner and P Meinicke
PeerJ 2:e239;   http://dx.doi.org /10.7717/peerj.239

High-throughput technologies notoriously generate large datasets often including data from different omics platforms. Each dataset contains data for several thousand experimental markers, e.g., mass-to-charge ratios in mass spectrometry or spots in DNA microarray analysis. An experimental marker is associated with an intensity profile which may include several measurements according to different experimental conditions (Dettmer, Aronov & Hammock, 2007).

The combined analysis and visualization of data from different high-throughput technologies remains a key challenge in bioinformatics.We present here theMarVis-Graph software for integrative analysis of metabolic and transcriptomic data. All experimental data is investigated in terms of the full metabolic network obtained from a reference database. The reactions of the network are scored based on the associated data, and

  • sub-networks, according to connected high-scoring reactions, are identified.

Finally, MarVis-Graph scores the detected sub-networks,

  • evaluates them by means of a random permutation test and
  • presents them as a ranked list.

Furthermore, MarVis-Graph features an interactive network visualization that provides researchers with a convenient view on the results.

The key advantage ofMarVis-Graph is the analysis of reactions detached from their pathways so that

  • it is possible to identify new pathways or
  • to connect known pathways by previously unrelated reactions.

TheMarVis-Graph software is freely available for academic use and can be downloaded at: http://marvis.gobics.de/marvis-graph.

Significant differences or clusters may be explained by associated annotations, e.g., in terms of metabolic pathways or biological functions. During recent years, numerous specialized tools have been developed to aid biological researchers in automating all these steps (e.g., Medina et al., 2010; Kaever et al., 2009; Waegele et al., 2012). Comprehensive studies can be performed by combining technologies from different omics fields. The combination of transcriptomic and proteomic data sets revealed a strong
correlation between both kinds of data (Nie et al., 2007) and supported the detection of complex interactions, e.g., in RNA silencing (Haq et al., 2010). Moreover, correlations
were detected between RNA expression levels and metabolite abundances (Gibon et al., 2006). Therefore, tools that integrate, analyze and visualize experimental markers from different platforms are needed. To cope with the complexity of genome-wide studies, pathway models are utilized extensively as a simple abstraction of the underlying complex mechanisms. Set Enrichment Analysis (Subramanian et al., 2005) and Over-Representation Analysis (Huang, Sherman & Lempicki, 2009) have become state-of-the-art tools for analyzing large-scale datasets: both methods evaluate predefined sets of entities, e.g., the accumulation of differentially expressed genes in a pathway.

While manually curated pathways are convenient and easy to interpret, experimental studies have shown that all metabolic and signaling pathways are heavily interconnected (Kunkel & Brooks, 2002; Laule et al., 2003). Data from biomolecular databases support these studies: the metabolic network of Arabidopsis thaliana in the KEGG database (Kanehisa et al., 2012; Kanehisa & Goto, 2000) contains 1606 reactions from which 1464 are connected in a single sub-network (>91%), i.e., they
share a metabolite as product or substrate. In the AraCyc 10.0 database (Mueller, Zhang & Rhee, 2003; Rhee et al., 2006), more than 89% of the reactions are counted in a single sub-network. In both databases, most other reactions are completely disconnected. Additionally, Set Enrichment Analyses can not identify links between the predefined sets easily. This becomes even more important when analyzing smaller pathways as provided by the MetaCyc (Caspi et al., 2008; Caspi et al., 2012) database. Moreover, methods that utilize pathways as predefined sets ignore reactions and related biomolecular entities (e.g., metabolites, genes) which are not associated with a single pathway. For example, this affects 4000 reactions in MetaCyc and 2500 in KEGG, respectively (Altman et al., 2013). Therefore, it is desirable to develop additional methods

  • that do not require predefined sets but may detect enriched sub-networks in the full metabolic network.

While several tools support the statistical analysis of experimental markers from one or more omics technologies and then utilize variants of Set Enrichment Analysis (Xia et al., 2012; Chen et al., 2013; Howe et al., 2011),

  • no tool is able to explicitly search for connected reactions that include
  • most of the metabolites, genes, and enyzmes with experimental evidence.

However, the automatic identification of sub-networks has been proven useful in other contexts, e.g., in the analysis of protein–protein-interaction networks (Alcaraz et al., 2012; Baumbach et al., 2012; Maeyer et al., 2013).

MarVis-Graph imports experimental markers from different high-throughput experiments and

  • analyses them in the context of reaction-chains in full metabolic networks.

Then, MarVis-Graph scores the reactions in the metabolic network

  • according to the number of associated experimental markers and
  • identifies sub-networks consisting of subsequent, high-scoring reactions.

The resulting sub-networks are

  • ranked according to a scoring method and visualized interactively.

Hereby, sub-networks consisting of reactions from different pathways may be identified to be important

  • whereas the single pathways may not be found to be significantly enriched.

MarVis-Graph may also connect reactions without an assigned pathway

  • to reactions within a particular pathway.

TheMarVis-Graph tool was applied in a case-study investigating the wound response in Arabidopsis thaliana to analyze combined metabolomic and transcriptomic high-throughput data.

Figure 1 Schema of the metabolic network representation in MarVis-Graph. Metabolite markers are shown in gray, metabolites in red, reactions in blue, enzymes in green, genes in yellow, transcript markers in pink, and pathways in turquoise color. The edges are shown in black with labels that comply with the biological meaning. The orange arrows depict the flow of score for the initial scoring (described in section “Initial Scoring”). (not shown)

In MarVis-Graph, metabolite markers obtained from mass-spectrometry experiments additionally contain the experimental mass. The experimental mass has to be
calculated based on the mass-to-charge ratio (m/z-value) and specific isotope- or adduct-corrections (Draper et al., 2009) by means of specialized tools, e.g.,MarVis-Filter
(Kaever et al., 2012).

For each transcript marker the corresponding annotation has to be given. In DNA microarray experiments, each spot (transcript marker) is specific for a gene and can
therefore be used for annotation. For other technologies an annotation has to be provided by external tools.

In MarVis-Graph, each reaction is scored initially based on the associated experimental data (see “Initial scoring”). This initial scoring is refined (see “Refining the scoring”) and afterwards reactions with a score below a user-defined threshold are removed. The network is

  • decomposed into subsequent high-scoring reactions that constitute the sub-networks.

The weight of each experimental marker (see “Experimental markers”) is equally distributed over all metabolites and genes associated with the metabolite marker or
transcript marker, respectively. For all vertices, this is repeated as illustrated in Fig. 1 until the weights are accumulated by the reactions.

The initial reaction scores are used as input scoring for the random walk algorithm. The algorithm is performed as described by Glaab et al. (2012) with a user-defined
restart-probability r (default value 0.8). After convergence of the algorithm, reactions with a score lower than the user-defined threshold t (default value t = 1−r) are removed from the reaction network. During the removal process,

  • the network is decomposed into pairwise disconnected sub-networks containing only high-scoring reactions.

In the following, a resulting sub-network is denoted by a prime: G′ = (V′,L′) with V′ = M′ ∪C′ ∪R′ ∪E′ ∪G′ ∪T′ ∪P′.

The scores of the identified sub-networks can be assessed using a random permutation test, evaluating the marker annotations under the null hypothesis of being connected
randomly. Here, the assignments

  • from metabolite markers to metabolites and from transcript markers to genes are randomized.

For each association between a metabolite marker and a metabolite,

  • this connection is replaced by a connection between a randomly chosen metabolite marker and a randomly chosen metabolite.

The random metabolite marker is chosen from the pool of formerly connected metabolite markers. Each connected transcript marker

  • is associated with a randomly chosen gene.

Choosing from the list of already connected experimental markers ensures that

  • the sum of weights from the original and the permuted network are equal.

This method differs from the commonly utilized XSwap permutation (Hanhij¨arvi, Garriga & Puolam¨aki, 2009) that is based on swapping endpoints of two random edges. The main difference of our permutation method is that it results in a network with different topological structure, i.e., different degree of the metabolite and gene nodes.

Finally, the sub-networks are detected and scored with the same parameters applied for the original network. Based on the scores of the networks identified in the random
permutations, the family-wise-error-rate (FWER) and false-discovery-rate (FDR) are calculated for each originally identified sub-network.

MarVis-Graph was applied in a case study investigating the A. thaliana wound response. Data from a metabolite fingerprinting (Meinicke et al., 2008) and a DNA microarray
experiment (Yan et al., 2007) were imported into a metabolic network specific for A. thaliana created from the AraCyc 10.0 database (Lamesch et al., 2011). The metabolome
and transcriptome have been measured before wounding as control and at specific time points after wounding in wild-type and in the allene oxide synthase (AOS) knock-out
mutant dde-2-2 (Park et al., 2002) of A. thaliana Columbia (see Table 1). The AOS mutant was chosen, because AOS catalyzes the first specific step in the biosynthesis of the hormone jasmonic acid, which is the key regulator in wound response of plants (Wasternack & Hause, 2013).

Both datasets have been preprocessed with theMarVis-Filter tool (Kaever et al., 2012) utilizing the Kruskal–Wallis p-value calculation on the intensity profiles. Based on the ranking of ascending p-values,

  • the first 25% of the metabolite markers and 10% of the transcript markers have been selected for further investigation (Data S2).

The filtered metabolite and transcript markers were imported into the metabolic network. For metabolite markers, metabolites were associated

  • if the metabolite marker’s detected mass differs from the metabolites monoisotopic mass by a maximum of 0.005u.

Transcript markers were linked to the genes whose ID equaled the ID given in the CATMA database (Sclep et al., 2007) for that transcript marker.

Table 2 Vertices in the A. thaliana specific metabolic network after import of experimental markers. Number of objects in the metabolic network
in absolute counts and relative abundances. For experimental markers, the with annotation column gives the number of metabolite markers and
transcript markers that were annotated with a metabolite or gene, respectively. The direct evidence column contains the number of metabolites
and genes, that are associated with a metabolite marker or transcript marker. For enzymes, this is the number of enzymes encoded by a gene with
direct evidence. The number of vertices with an association to a reaction is given in the with reaction column. In the last column, this is given for
associations to metabolic pathways. (not shown)

MarVis-Graph detected a total of 133 sub-networks. The sub-networks were ranked according to size Ss, diameter Sd, and sum-of-weights Ssow
scores (Table S4). Interestingly, the different rankings show a high correlation with all pairwise correlations higher than 0.75 (Pearson correlation
coefficient) and 0.6 (Spearman rank correlation).

Allene-oxide cyclase sub-network
In all rankings, the sub-network allene-oxide cyclase (named after the reaction with the highest score in this sub-network) appeared as top candidate.

This sub-network is constituted of reactions from different pathways related to fatty acids. Figure 2 shows a visualization of the sub-network.
Jasmonic acid biosynthesis. The main part of the sub-network is formed by reactions from the “jasmonic acid biosynthesis” (PlantMetabolic Network, 2013)
resulting in jasmonic acid (jasmonate). The presence of this pathway is very well established because of its central role in mediating the plants wound response
(Reymond & Farmer, 1998; Creelman, Tierney & Mullet, 1992). Additionally, metabolites and transcripts from this pathway were expected to show prominent
expression profiles because AOS, a key enzyme in this pathway, is knocked-out in themutant plant. Jasmonic acid derivatives and hormones.

Jasmonic acid derivatives and hormones. Jasmonate is a precursor for a broad variety of plant hormones (Wasternack & Hause, 2013), e.g., the derivative (-)-
jasmonic acid methyl ester (also Methyl Jasmonic Acid; MeJA) is a volatile, airborne signal mediating wound response between plants (Farmer&Ryan, 1990).
Reactions from the jasmonoyl-amino acid conjugates biosynthesis I (PMN, 2013a) pathway connect jasmonate to different amino acids, including L-valine,
L-leucine, and L-isoleucine. Via these amino acids, this sub-network is connected to the indole-3-acetylamino acid biosynthesis (PMN, 2013b) (IAA biosynthesis).
Again, this pathway produces a well known plant hormone: Auxine (Woodward & Bartel, 2005). Even though, jasmonate and auxin are both plant hormones, their
connection in this subnetwork is of minor relevance because amino acid conjugates are often utilized as active or storage forms of signaling molecules.While
jasmonoyl-amino acid conjugates represent the active signaling form of jasmonates, IAA amino acid conjugates are the storage form of this hormone (Staswick et al.,
2005).

polyhydroxy fatty acids synthesis

polyhydroxy fatty acids synthesis

 

Figure 2 Schema of the allene-oxide cyclase sub-network. Metabolites are shown in red, reactions in blue, and enzymes in green color. Metabolites and reactions without direct experimental evidence are marked by a dashed outline and a brighter color while enzymes without experimental evidence are hidden. The metabolic pathways described in section “Resulting sub-networks” are highlighted with different colors. The orange and green parts indicate the reaction chains required to build jasmonate and its amino acid conjugates. The coloring of pathways was done manually after export from MarVis-Graph.

The ω-3-fatty acid desaturase should catalyze a reaction from linoleate to α-linolenate. Metabolite markers that match the mass of crepenynic acid do also match α-linolenate
because both molecules have the same sum-formula and monoisotopic mass. As mentioned above, MarVis-Graph compiled the metabolic network for this study
from the AraCyc database version 10.0. On June 4th, a curator changed the database to remove theΔ12-fatty acid dehydrogenase prior to the release of AraCyc version 11.0.

The presented new software tool MarVis-Graph supports the investigation and visualization of omics data from different fields of study. The introduced algorithm for
identification of sub-networks is able to identify reaction-chains across different pathways and includes reactions that are not associated with a single pathway. The application of MarVis-Graph in the case study on A. thaliana wound response resulted in a convenient graphical representation of high-throughput data which allows the analysis of the complex dynamics in a metabolic network.

 

Read Full Post »

Preface to Metabolomics as a Discipline in Medicine

Author: Larry H. Bernstein, MD, FCAP

 

The family of ‘omics fields has rapidly outpaced its siblings over the decade since
the completion of the Human Genome Project.  It has derived much benefit from
the development of Proteomics, which has recently completed a first draft of the
human proteome.  Since genomics, transcriptomics, and proteomics, have matured
considerably, it has become apparent that the search for a driver or drivers of cellular signaling and metabolic pathways could not depend on a full clarity of the genome. There have been unresolved issues, that are not solely comprehended from assumptions about mutations.

The most common diseases affecting mankind are derangements in metabolic
pathways, develop at specific ages periods, and often in adulthood or in the
geriatric period, and are at the intersection of signaling pathways.  Moreover,
the organs involved and systemic features are heavily influenced by physical
activity, and by the air we breathe and the water we drink.

The emergence of the new science is also driven by a large body of work
on protein structure, mechanisms of enzyme action, the modulation of gene
expression, the pH dependent effects on protein binding and conformation.
Beyond what has just been said, a significant portion of DNA has been
designated as “dark matter”. It turns out to have enormous importance in
gene regulation, even though it is not transcriptional, effected in a
modulatory way by “noncoding RNAs.  Metabolomics is the comprehensive
analysis of small molecule metabolites. These might be substrates of
sequenced enzyme reactions, or they might be “inhibiting” RNAs just
mentioned.  In either case, they occur in the substructures of the cell
called organelles, the cytoplasm, and in the cytoskeleton.

The reactions are orchestrated, and they can be modified with respect to
the flow of metabolites based on pH, temperature, membrane structural
modifications, and modulators.  Since most metabolites are generated by
enzymatic proteins that result from gene expression, and metabolites give
organisms their biochemical characteristics, the metabolome links
genotype with phenotype.

Metabolomics is still developing, and the continued development has
relied on two major events. The first is chromatographic separation and
mass  spectroscopy (MS), MS/MS, as well as advances in fluorescence
ultrasensitive optical photonic methods, and the second, as crucial,
is the developments in computational biology. The continuation of
this trend brings expectations of an impact on pharmaceutical and
on neutraceutical developments, which will have an impact on medical
practice. What has lagged behind, and may continue to contribute to the
lag is the failure to develop a suitable electronic medical record to
assist the physician in decisions confronted with so much as yet,
hidden data, the ready availability of which could guide more effective
diagnosis and management of the patient. Put all of this together, and
we can meet series challenges as the research community
interprets and integrates the complex data they are acquiring.

.

Read Full Post »

Metabolomics Summary and Perspective

Metabolomics Summary and Perspective

Author and Curator: Larry H Bernstein, MD, FCAP 

 

This is the final article in a robust series on metabolism, metabolomics, and  the “-OMICS-“ biological synthesis that is creating a more holistic and interoperable view of natural sciences, including the biological disciplines, climate science, physics, chemistry, toxicology, pharmacology, and pathophysiology with as yet unforeseen consequences.

There have been impressive advances already in the research into developmental biology, plant sciences, microbiology, mycology, and human diseases, most notably, cancer, metabolic , and infectious, as well as neurodegenerative diseases.

Acknowledgements:

I write this article in honor of my first mentor, Harry Maisel, Professor and Emeritus Chairman of Anatomy, Wayne State University, Detroit, MI and to my stimulating mentors, students, fellows, and associates over many years:

Masahiro Chiga, MD, PhD, Averill A Liebow, MD, Nathan O Kaplan, PhD, Johannes Everse, PhD, Norio Shioura, PhD, Abraham Braude, MD, Percy J Russell, PhD, Debby Peters, Walter D Foster, PhD, Herschel Sidransky, MD, Sherman Bloom, MD, Matthew Grisham, PhD, Christos Tsokos, PhD,  IJ Good, PhD, Distinguished Professor, Raool Banagale, MD, Gustavo Reynoso, MD,Gustave Davis, MD, Marguerite M Pinto, MD, Walter Pleban, MD, Marion Feietelson-Winkler, RD, PhD,  John Adan,MD, Joseph Babb, MD, Stuart Zarich, MD,  Inder Mayall, MD, A Qamar, MD, Yves Ingenbleek, MD, PhD, Emeritus Professor, Bette Seamonds, PhD, Larry Kaplan, PhD, Pauline Y Lau, PhD, Gil David, PhD, Ronald Coifman, PhD, Emeritus Professor, Linda Brugler, RD, MBA, James Rucinski, MD, Gitta Pancer, Ester Engelman, Farhana Hoque, Mohammed Alam, Michael Zions, William Fleischman, MD, Salman Haq, MD, Jerard Kneifati-Hayek, Madeleine Schleffer, John F Heitner, MD, Arun Devakonda,MD, Liziamma George,MD, Suhail Raoof, MD, Charles Oribabor,MD, Anthony Tortolani, MD, Prof and Chairman, JRDS Rosalino, PhD, Aviva Lev Ari, PhD, RN, Rosser Rudolph, MD, PhD, Eugene Rypka, PhD, Jay Magidson, PhD, Izaak Mayzlin, PhD, Maurice Bernstein, PhD, Richard Bing, Eli Kaplan, PhD, Maurice Bernstein, PhD.

This article has EIGHT parts, as follows:

Part 1

Metabolomics Continues Auspicious Climb

Part 2

Biologists Find ‘Missing Link’ in the Production of Protein Factories in Cells

Part 3

Neuroscience

Part 4

Cancer Research

Part 5

Metabolic Syndrome

Part 6

Biomarkers

Part 7

Epigenetics and Drug Metabolism

Part 8

Pictorial

genome cartoon

genome cartoon

 iron metabolism

iron metabolism

personalized reference range within population range

personalized reference range within population range

Part 1.  MetabolomicsSurge

metagraph  _OMICS

metagraph _OMICS

Metabolomics Continues Auspicious Climb

Jeffery Herman, Ph.D.
GEN May 1, 2012 (Vol. 32, No. 9)

Aberrant biochemical and metabolite signaling plays an important role in

  • the development and progression of diseased tissue.

This concept has been studied by the science community for decades. However, with relatively

  1. recent advances in analytical technology and bioinformatics as well as
  2. the development of the Human Metabolome Database (HMDB),

metabolomics has become an invaluable field of research.

At the “International Conference and Exhibition on Metabolomics & Systems Biology” held recently in San Francisco, researchers and industry leaders discussed how

  • the underlying cellular biochemical/metabolite fingerprint in response to
  1. a specific disease state,
  2. toxin exposure, or
  3. pharmaceutical compound
  • is useful in clinical diagnosis and biomarker discovery and
  • in understanding disease development and progression.

Developed by BASF, MetaMap® Tox is

  • a database that helps identify in vivo systemic effects of a tested compound, including
  1. targeted organs,
  2. mechanism of action, and
  3. adverse events.

Based on 28-day systemic rat toxicity studies, MetaMap Tox is composed of

  • differential plasma metabolite profiles of rats
  • after exposure to a large variety of chemical toxins and pharmaceutical compounds.

“Using the reference data,

  • we have developed more than 110 patterns of metabolite changes, which are
  • specific and predictive for certain toxicological modes of action,”

said Hennicke Kamp, Ph.D., group leader, department of experimental toxicology and ecology at BASF.

With MetaMap Tox, a potential drug candidate

  • can be compared to a similar reference compound
  • using statistical correlation algorithms,
  • which allow for the creation of a toxicity and mechanism of action profile.

“MetaMap Tox, in the context of early pre-clinical safety enablement in pharmaceutical development,” continued Dr. Kamp,

  • has been independently validated “
  • by an industry consortium (Drug Safety Executive Council) of 12 leading biopharmaceutical companies.”

Dr. Kamp added that this technology may prove invaluable

  • allowing for quick and accurate decisions and
  • for high-throughput drug candidate screening, in evaluation
  1. on the safety and efficacy of compounds
  2. during early and preclinical toxicological studies,
  3. by comparing a lead compound to a variety of molecular derivatives, and
  • the rapid identification of the most optimal molecular structure
  • with the best efficacy and safety profiles might be streamlined.
Dynamic Construct of the –Omics

Dynamic Construct of the –Omics

Targeted Tandem Mass Spectrometry

Biocrates Life Sciences focuses on targeted metabolomics, an important approach for

  • the accurate quantification of known metabolites within a biological sample.

Originally used for the clinical screening of inherent metabolic disorders from dried blood-spots of newborn children, Biocrates has developed

  • a tandem mass spectrometry (MS/MS) platform, which allows for
  1. the identification,
  2. quantification, and
  3. mapping of more than 800 metabolites to specific cellular pathways.

It is based on flow injection analysis and high-performance liquid chromatography MS/MS.

Clarification of Pathway-Specific Inhibition by Fourier Transform Ion Cyclotron Resonance.Mass Spectrometry-Based Metabolic Phenotyping Studies F5.large

common drug targets

common drug targets

The MetaDisIDQ® Kit is a

  • “multiparamatic” diagnostic assay designed for the “comprehensive assessment of a person’s metabolic state” and
  • the early determination of pathophysiological events with regards to a specific disease.

MetaDisIDQ is designed to quantify

  • a diverse range of 181 metabolites involved in major metabolic pathways
  • from a small amount of human serum (10 µL) using isotopically labeled internal standards,

This kit has been demonstrated to detect changes in metabolites that are commonly associated with the development of

  • metabolic syndrome, type 2 diabetes, and diabetic nephropathy,

Dr. Dallman reports that data generated with the MetaDisIDQ kit correlates strongly with

  • routine chemical analyses of common metabolites including glucose and creatinine

Biocrates has also developed the MS/MS-based AbsoluteIDQ® kits, which are

  • an “easy-to-use” biomarker analysis tool for laboratory research.

The kit functions on MS machines from a variety of vendors, and allows for the quantification of 150-180 metabolites.

The SteroIDQ® kit is a high-throughput standardized MS/MS diagnostic assay,

  • validated in human serum, for the rapid and accurate clinical determination of 16 known steroids.

Initially focusing on the analysis of steroid ranges for use in hormone replacement therapy, the SteroIDQ Kit is expected to have a wide clinical application.

Hormone-Resistant Breast Cancer

Scientists at Georgetown University have shown that

  • breast cancer cells can functionally coordinate cell-survival and cell-proliferation mechanisms,
  • while maintaining a certain degree of cellular metabolism.

To grow, cells need energy, and energy is a product of cellular metabolism. For nearly a century, it was thought that

  1. the uncoupling of glycolysis from the mitochondria,
  2. leading to the inefficient but rapid metabolism of glucose and
  3. the formation of lactic acid (the Warburg effect), was

the major and only metabolism driving force for unchecked proliferation and tumorigenesis of cancer cells.

Other aspects of metabolism were often overlooked.

“.. we understand now that

  • cellular metabolism is a lot more than just metabolizing glucose,”

said Robert Clarke, Ph.D., professor of oncology and physiology and biophysics at Georgetown University. Dr. Clarke, in collaboration with the Waters Center for Innovation at Georgetown University (led by Albert J. Fornace, Jr., M.D.), obtained

  • the metabolomic profile of hormone-sensitive and -resistant breast cancer cells through the use of UPLC-MS.

They demonstrated that breast cancer cells, through a rather complex and not yet completely understood process,

  1. can functionally coordinate cell-survival and cell-proliferation mechanisms,
  2. while maintaining a certain degree of cellular metabolism.

This is at least partly accomplished through the upregulation of important pro-survival mechanisms; including

  • the unfolded protein response;
  • a regulator of endoplasmic reticulum stress and
  • initiator of autophagy.

Normally, during a stressful situation, a cell may

  • enter a state of quiescence and undergo autophagy,
  • a process by which a cell can recycle organelles
  • in order to maintain enough energy to survive during a stressful situation or,

if the stress is too great,

  • undergo apoptosis.

By integrating cell-survival mechanisms and cellular metabolism

  • advanced ER+ hormone-resistant breast cancer cells
  • can maintain a low level of autophagy
  • to adapt and resist hormone/chemotherapy treatment.

This adaptation allows cells

  • to reallocate important metabolites recovered from organelle degradation and
  • provide enough energy to also promote proliferation.

With further research, we can gain a better understanding of the underlying causes of hormone-resistant breast cancer, with

  • the overall goal of developing effective diagnostic, prognostic, and therapeutic tools.

NMR

Over the last two decades, NMR has established itself as a major tool for metabolomics analysis. It is especially adept at testing biological fluids. [Bruker BioSpin]

Historically, nuclear magnetic resonance spectroscopy (NMR) has been used for structural elucidation of pure molecular compounds. However, in the last two decades, NMR has established itself as a major tool for metabolomics analysis. Since

  • the integral of an NMR signal is directly proportional to
  • the molar concentration throughout the dynamic range of a sample,

“the simultaneous quantification of compounds is possible

  • without the need for specific reference standards or calibration curves,” according to Lea Heintz of Bruker BioSpin.

NMR is adept at testing biological fluids because of

  1.  high reproducibility,
  2. standardized protocols,
  3. low sample manipulation, and
  4. the production of a large subset of data,

Bruker BioSpin is presently involved in a project for the screening of inborn errors of metabolism in newborn children from Turkey, based on their urine NMR profiles. More than 20 clinics are participating to the project that is coordinated by INFAI, a specialist in the transfer of advanced analytical technology into medical diagnostics. The construction of statistical models are being developed

  • for the detection of deviations from normality, as well as
  • automatic quantification methods for indicative metabolites

Bruker BioSpin recently installed high-resolution magic angle spinning NMR (HRMAS-NMR) systems that can rapidly analyze tissue biopsies. The main objective for HRMAS-NMR is to establish a rapid and effective clinical method to assess tumor grade and other important aspects of cancer during surgery.

Combined NMR and Mass Spec

There is increasing interest in combining NMR and MS, two of the main analytical assays in metabolomic research, as a means

  • to improve data sensitivity and to
  • fully elucidate the complex metabolome within a given biological sample.
  •  to realize a potential for cancer biomarker discovery in the realms of diagnosis, prognosis, and treatment.

.

Using combined NMR and MS to measure the levels of nearly 250 separate metabolites in the patient’s blood, Dr. Weljie and other researchers at the University of Calgary were able to rapidly determine the malignancy of a  pancreatic lesion (in 10–15% of the cases, it is difficult to discern between benign and malignant), while avoiding unnecessary surgery in patients with benign lesions.

When performing NMR and MS on a single biological fluid, ultimately “we are,” noted Dr. Weljie,

  1. “splitting up information content, processing, and introducing a lot of background noise and error and
  2. then trying to reintegrate the data…
    It’s like taking a complex item, with multiple pieces, out of an IKEA box and trying to repackage it perfectly into another box.”

By improving the workflow between the initial splitting of the sample, they improved endpoint data integration, proving that

  • a streamlined approach to combined NMR/MS can be achieved,
  • leading to a very strong, robust and precise metabolomics toolset.

Metabolomics Research Picks Up Speed

Field Advances in Quest to Improve Disease Diagnosis and Predict Drug Response

John Morrow Jr., Ph.D.
GEN May 1, 2011 (Vol. 31, No. 9)

As an important discipline within systems biology, metabolomics is being explored by a number of laboratories for

  • its potential in pharmaceutical development.

Studying metabolites can offer insights into the relationships between genotype and phenotype, as well as between genotype and environment. In addition, there is plenty to work with—there are estimated to be some 2,900 detectable metabolites in the human body, of which

  1. 309 have been identified in cerebrospinal fluid,
  2. 1,122 in serum,
  3. 458 in urine, and
  4. roughly 300 in other compartments.

Guowang Xu, Ph.D., a researcher at the Dalian Institute of Chemical Physics.  is investigating the causes of death in China,

  • and how they have been changing over the years as the country has become a more industrialized nation.
  •  the increase in the incidence of metabolic disorders such as diabetes has grown to affect 9.7% of the Chinese population.

Dr. Xu,  collaborating with Rainer Lehman, Ph.D., of the University of Tübingen, Germany, compared urinary metabolites in samples from healthy individuals with samples taken from prediabetic, insulin-resistant subjects. Using mass spectrometry coupled with electrospray ionization in the positive mode, they observed striking dissimilarities in levels of various metabolites in the two groups.

“When we performed a comprehensive two-dimensional gas chromatography, time-of-flight mass spectrometry analysis of our samples, we observed several metabolites, including

  • 2-hydroxybutyric acid in plasma,
  •  as potential diabetes biomarkers,” Dr. Xu explains.

In other, unrelated studies, Dr. Xu and the German researchers used a metabolomics approach to investigate the changes in plasma metabolite profiles immediately after exercise and following a 3-hour and 24-hour period of recovery. They found that

  • medium-chain acylcarnitines were the most distinctive exercise biomarkers, and
  • they are released as intermediates of partial beta oxidation in human myotubes and mouse muscle tissue.

Dr. Xu says. “The traditional approach of assessment based on a singular biomarker is being superseded by the introduction of multiple marker profiles.”

Typical of the studies under way by Dr. Kaddurah-Daouk and her colleaguesat Duke University

  • is a recently published investigation highlighting the role of an SNP variant in
  • the glycine dehydrogenase gene on individual response to antidepressants.
  •  patients who do not respond to the selective serotonin uptake inhibitors citalopram and escitalopram
  • carried a particular single nucleotide polymorphism in the GD gene.

“These results allow us to pinpoint a possible

  • role for glycine in selective serotonin reuptake inhibitor response and
  • illustrate the use of pharmacometabolomics to inform pharmacogenomics.

These discoveries give us the tools for prognostics and diagnostics so that

  • we can predict what conditions will respond to treatment.

“This approach to defining health or disease in terms of metabolic states opens a whole new paradigm.

By screening hundreds of thousands of molecules, we can understand

  • the relationship between human genetic variability and the metabolome.”

Dr. Kaddurah-Daouk talks about statins as a current

  • model of metabolomics investigations.

It is now known that the statins  have widespread effects, altering a range of metabolites. To sort out these changes and develop recommendations for which individuals should be receiving statins will require substantial investments of energy and resources into defining the complex web of biochemical changes that these drugs initiate.
Furthermore, Dr. Kaddurah-Daouk asserts that,

  • “genetics only encodes part of the phenotypic response.

One needs to take into account the

  • net environment contribution in order to determine
  • how both factors guide the changes in our metabolic state that determine the phenotype.”

Interactive Metabolomics

Researchers at the University of Nottingham use diffusion-edited nuclear magnetic resonance spectroscopy to assess the effects of a biological matrix on metabolites. Diffusion-edited NMR experiments provide a way to

  • separate the different compounds in a mixture
  • based on the differing translational diffusion coefficients (which reflect the size and shape of the molecule).

The measurements are carried out by observing

  • the attenuation of the NMR signals during a pulsed field gradient experiment.

Clare Daykin, Ph.D., is a lecturer at the University of Nottingham, U.K. Her field of investigation encompasses “interactive metabolomics,”which she defines as

“the study of the interactions between low molecular weight biochemicals and macromolecules in biological samples ..

  • without preselection of the components of interest.

“Blood plasma is a heterogeneous mixture of molecules that

  1. undergo a variety of interactions including metal complexation,
  2. chemical exchange processes,
  3. micellar compartmentation,
  4. enzyme-mediated biotransformations, and
  5. small molecule–macromolecular binding.”

Many low molecular weight compounds can exist

  • freely in solution,
  • bound to proteins, or
  • within organized aggregates such as lipoprotein complexes.

Therefore, quantitative comparison of plasma composition from

  • diseased individuals compared to matched controls provides an incomplete insight to plasma metabolism.

“It is not simply the concentrations of metabolites that must be investigated,

  • but their interactions with the proteins and lipoproteins within this complex web.

Rather than targeting specific metabolites of interest, Dr. Daykin’s metabolite–protein binding studies aim to study

  • the interactions of all detectable metabolites within the macromolecular sample.

Such activities can be studied through the use of diffusion-edited nuclear magnetic resonance (NMR) spectroscopy, in which one can assess

  • the effects of the biological matrix on the metabolites.

“This can lead to a more relevant and exact interpretation

  • for systems where metabolite–macromolecule interactions occur.”

Diffusion-edited NMR experiments provide a way to separate the different compounds in a mixture based on

  • the differing translational diffusion coefficients (which reflect the size and shape of the molecule).

The measurements are carried out by observing

  • the attenuation of the NMR signals during a pulsed field gradient experiment.

Pushing the Limits

It is widely recognized that many drug candidates fail during development due to ancillary toxicity. Uwe Sauer, Ph.D., professor, and Nicola Zamboni, Ph.D., researcher, both at the Eidgenössische Technische Hochschule, Zürich (ETH Zürich), are applying

  • high-throughput intracellular metabolomics to understand
  • the basis of these unfortunate events and
  • head them off early in the course of drug discovery.

“Since metabolism is at the core of drug toxicity, we developed a platform for

  • measurement of 50–100 targeted metabolites by
  • a high-throughput system consisting of flow injection
  • coupled to tandem mass spectrometry.”

Using this approach, Dr. Sauer’s team focused on

  • the central metabolism of the yeast Saccharomyces cerevisiae, reasoning that
  • this core network would be most susceptible to potential drug toxicity.

Screening approximately 41 drugs that were administered at seven concentrations over three orders of magnitude, they observed changes in metabolome patterns at much lower drug concentrations without attendant physiological toxicity.

The group carried out statistical modeling of about

  • 60 metabolite profiles for each drug they evaluated.

This data allowed the construction of a “profile effect map” in which

  • the influence of each drug on metabolite levels can be followed, including off-target effects, which
  • provide an indirect measure of the possible side effects of the various drugs.

Dr. Sauer says.“We have found that this approach is

  • at least 100 times as fast as other omics screening platforms,”

“Some drugs, including many anticancer agents,

  • disrupt metabolism long before affecting growth.”
killing cancer cells

killing cancer cells

Furthermore, they used the principle of 13C-based flux analysis, in which

  • metabolites labeled with 13C are used to follow the utilization of metabolic pathways in the cell.

These 13C-determined intracellular responses of metabolic fluxes to drug treatment demonstrate

  • the functional performance of the network to be rather robust,
conformational changes leading to substrate efflux.

conformational changes leading to substrate efflux.

leading Dr. Sauer to the conclusion that

  • the phenotypic vigor he observes to drug challenges
  • is achieved by a flexible make up of the metabolome.

Dr. Sauer is confident that it will be possible to expand the scope of these investigations to hundreds of thousands of samples per study. This will allow answers to the questions of

  • how cells establish a stable functioning network in the face of inevitable concentration fluctuations.

Is Now the Hour?

There is great enthusiasm and agitation within the biotech community for

  • metabolomics approaches as a means of reversing the dismal record of drug discovery

that has accumulated in the last decade.

While the concept clearly makes sense and is being widely applied today, there are many reasons why drugs fail in development, and metabolomics will not be a panacea for resolving all of these questions. It is too early at this point to recognize a trend or a track record, and it will take some time to see how this approach can aid in drug discovery and shorten the timeline for the introduction of new pharmaceutical agents.

Degree of binding correlated with function

Degree of binding correlated with function

Diagram_of_a_two-photon_excitation_microscope_

Diagram_of_a_two-photon_excitation_microscope_

Part 2.  Biologists Find ‘Missing Link’ in the Production of Protein Factories in Cells

Biologists at UC San Diego have found

  • the “missing link” in the chemical system that
  • enables animal cells to produce ribosomes

—the thousands of protein “factories” contained within each cell that

  • manufacture all of the proteins needed to build tissue and sustain life.
‘Missing Link’

‘Missing Link’

Their discovery, detailed in the June 23 issue of the journal Genes & Development, will not only force

  • a revision of basic textbooks on molecular biology, but also
  • provide scientists with a better understanding of
  • how to limit uncontrolled cell growth, such as cancer,
  • that might be regulated by controlling the output of ribosomes.

Ribosomes are responsible for the production of the wide variety of proteins that include

  1. enzymes;
  2. structural molecules, such as hair,
  3. skin and bones;
  4. hormones like insulin; and
  5. components of our immune system such as antibodies.

Regarded as life’s most important molecular machine, ribosomes have been intensively studied by scientists (the 2009 Nobel Prize in Chemistry, for example, was awarded for studies of its structure and function). But until now researchers had not uncovered all of the details of how the proteins that are used to construct ribosomes are themselves produced.

In multicellular animals such as humans,

  • ribosomes are made up of about 80 different proteins
    (humans have 79 while some other animals have a slightly different number) as well as
  • four different kinds of RNA molecules.

In 1969, scientists discovered that

  • the synthesis of the ribosomal RNAs is carried out by specialized systems using two key enzymes:
  • RNA polymerase I and RNA polymerase III.

But until now, scientists were unsure if a complementary system was also responsible for

  • the production of the 80 proteins that make up the ribosome.

That’s essentially what the UC San Diego researchers headed by Jim Kadonaga, a professor of biology, set out to examine. What they found was the missing link—the specialized

  • system that allows ribosomal proteins themselves to be synthesized by the cell.

Kadonaga says that he and coworkers found that ribosomal proteins are synthesized via

  • a novel regulatory system with the enzyme RNA polymerase II and
  • a factor termed TRF2,”

“For the production of most proteins,

  1. RNA polymerase II functions with
  2. a factor termed TBP,
  3. but for the synthesis of ribosomal proteins, it uses TRF2.”
  •  this specialized TRF2-based system for ribosome biogenesis
  • provides a new avenue for the study of ribosomes and
  • its control of cell growth, and

“it should lead to a better understanding and potential treatment of diseases such as cancer.”

Coordination of the transcriptome and metabolome

Coordination of the transcriptome and metabolome

the potential advantages conferred by distal-site protein synthesis

the potential advantages conferred by distal-site protein synthesis

Other authors of the paper were UC San Diego biologists Yuan-Liang Wang, Sascha Duttke and George Kassavetis, and Kai Chen, Jeff Johnston, and Julia Zeitlinger of the Stowers Institute for Medical Research in Kansas City, Missouri. Their research was supported by two grants from the National Institutes of Health (1DP2OD004561-01 and R01 GM041249).

Turning Off a Powerful Cancer Protein

Scientists have discovered how to shut down a master regulatory transcription factor that is

  • key to the survival of a majority of aggressive lymphomas,
  • which arise from the B cells of the immune system.

The protein, Bcl6, has long been considered too complex to target with a drug since it is also crucial

  • to the healthy functioning of many immune cells in the body, not just B cells gone bad.

The researchers at Weill Cornell Medical College report that it is possible

  • to shut down Bcl6 in diffuse large B-cell lymphoma (DLBCL)
  • while not affecting its vital function in T cells and macrophages
  • that are needed to support a healthy immune system.

If Bcl6 is completely inhibited, patients might suffer from systemic inflammation and atherosclerosis. The team conducted this new study to help clarify possible risks, as well as to understand

  • how Bcl6 controls the various aspects of the immune system.

The findings in this study were inspired from

  • preclinical testing of two Bcl6-targeting agents that Dr. Melnick and his Weill Cornell colleagues have developed
  • to treat DLBCLs.

These experimental drugs are

  • RI-BPI, a peptide mimic, and
  • the small molecule agent 79-6.

“This means the drugs we have developed against Bcl6 are more likely to be

  • significantly less toxic and safer for patients with this cancer than we realized,”

says Ari Melnick, M.D., professor of hematology/oncology and a hematologist-oncologist at NewYork-Presbyterian Hospital/Weill Cornell Medical Center.

Dr. Melnick says the discovery that

  • a master regulatory transcription factor can be targeted
  • offers implications beyond just treating DLBCL.

Recent studies from Dr. Melnick and others have revealed that

  • Bcl6 plays a key role in the most aggressive forms of acute leukemia, as well as certain solid tumors.

Bcl6 can control the type of immune cell that develops in the bone marrow—playing many roles

  • in the development of B cells, T cells, macrophages, and other cells—including a primary and essential role in
  • enabling B-cells to generate specific antibodies against pathogens.

According to Dr. Melnick, “When cells lose control of Bcl6,

  • lymphomas develop in the immune system.

Lymphomas are ‘addicted’ to Bcl6, and therefore

  • Bcl6 inhibitors powerfully and quickly destroy lymphoma cells,” .

The big surprise in the current study is that rather than functioning as a single molecular machine,

  • Bcl6 functions like a Swiss Army knife,
  • using different tools to control different cell types.

This multifunction paradigm could represent a general model for the functioning of other master regulatory transcription factors.

“In this analogy, the Swiss Army knife, or transcription factor, keeps most of its tools folded,

  • opening only the one it needs in any given cell type,”

He makes the following analogy:

  • “For B cells, it might open and use the knife tool;
  • for T cells, the cork screw;
  • for macrophages, the scissors.”

“this means that you only need to prevent the master regulator from using certain tools to treat cancer. You don’t need to eliminate the whole knife,” . “In fact, we show that taking out the whole knife is harmful since

  • the transcription factor has many other vital functions that other cells in the body need.”

Prior to these study results, it was not known that a master regulator could separate its functions so precisely. Researchers hope this will be a major benefit to the treatment of DLBCL and perhaps other disorders that are influenced by Bcl6 and other master regulatory transcription factors.

The study is published in the journal Nature Immunology, in a paper titled “Lineage-specific functions of Bcl-6 in immunity and inflammation are mediated by distinct biochemical mechanisms”.

Part 3. Neuroscience

Vesicles influence function of nerve cells 
Oct, 06 2014        source: http://feeds.sciencedaily.com

Neurons (blue) which have absorbed exosomes (green) have increased levels of the enzyme catalase (red), which helps protect them against peroxides.

Neurons (blue) which have absorbed exosomes (green) have increased levels of the enzyme catalase (red), which helps protect them against peroxides.

Neurons (blue) which have absorbed exosomes (green) have increased levels of the enzyme catalase (red), which helps protect them against peroxides.

Tiny vesicles containing protective substances

  • which they transmit to nerve cells apparently
  • play an important role in the functioning of neurons.

As cell biologists at Johannes Gutenberg University Mainz (JGU) have discovered,

  • nerve cells can enlist the aid of mini-vesicles of neighboring glial cells
  • to defend themselves against stress and other potentially detrimental factors.

These vesicles, called exosomes, appear to stimulate the neurons on various levels:

  • they influence electrical stimulus conduction,
  • biochemical signal transfer, and
  • gene regulation.

Exosomes are thus multifunctional signal emitters

  • that can have a significant effect in the brain.
Exosome

Exosome

The researchers in Mainz already observed in a previous study that

  • oligodendrocytes release exosomes on exposure to neuronal stimuli.
  • these are absorbed by the neurons and improve neuronal stress tolerance.

Oligodendrocytes, a type of glial cell, form an

  • insulating myelin sheath around the axons of neurons.

The exosomes transport protective proteins such as

  • heat shock proteins,
  • glycolytic enzymes, and
  • enzymes that reduce oxidative stress from one cell type to another,
  • but also transmit genetic information in the form of ribonucleic acids.

“As we have now discovered in cell cultures, exosomes seem to have a whole range of functions,” explained Dr. Eva-Maria Krmer-Albers. By means of their transmission activity, the small bubbles that are the vesicles

  • not only promote electrical activity in the nerve cells, but also
  • influence them on the biochemical and gene regulatory level.

“The extent of activities of the exosomes is impressive,” added Krmer-Albers. The researchers hope that the understanding of these processes will contribute to the development of new strategies for the treatment of neuronal diseases. Their next aim is to uncover how vesicles actually function in the brains of living organisms.

http://labroots.com/user/news/article/id/217438/title/vesicles-influence-function-of-nerve-cells

The above story is based on materials provided by Universitt Mainz.

Universitt Mainz. “Vesicles influence function of nerve cells.” ScienceDaily. ScienceDaily, 6 October 2014. www.sciencedaily.com/releases/2014/10/141006174214.htm

Neuroscientists use snail research to help explain “chemo brain”

10/08/2014
It is estimated that as many as half of patients taking cancer drugs experience a decrease in mental sharpness. While there have been many theories, what causes “chemo brain” has eluded scientists.

In an effort to solve this mystery, neuroscientists at The University of Texas Health Science Center at Houston (UTHealth) conducted an experiment in an animal memory model and their results point to a possible explanation. Findings appeared in The Journal of Neuroscience.

In the study involving a sea snail that shares many of the same memory mechanisms as humans and a drug used to treat a variety of cancers, the scientists identified

  • memory mechanisms blocked by the drug.

Then, they were able to counteract or

  • unblock the mechanisms by administering another agent.

“Our research has implications in the care of people given to cognitive deficits following drug treatment for cancer,” said John H. “Jack” Byrne, Ph.D., senior author, holder of the June and Virgil Waggoner Chair and Chairman of the Department of Neurobiology and Anatomy at the UTHealth Medical School. “There is no satisfactory treatment at this time.”

Byrne’s laboratory is known for its use of a large snail called Aplysia californica to further the understanding of the biochemical signaling among nerve cells (neurons).  The snails have large neurons that relay information much like those in humans.

When Byrne’s team compared cell cultures taken from normal snails to

  • those administered a dose of a cancer drug called doxorubicin,

the investigators pinpointed a neuronal pathway

  • that was no longer passing along information properly.

With the aid of an experimental drug,

  • the scientists were able to reopen the pathway.

Unfortunately, this drug would not be appropriate for humans, Byrne said. “We want to identify other drugs that can rescue these memory mechanisms,” he added.

According the American Cancer Society, some of the distressing mental changes cancer patients experience may last a short time or go on for years.

Byrne’s UT Health research team includes co-lead authors Rong-Yu Liu, Ph.D., and Yili Zhang, Ph.D., as well as Brittany Coughlin and Leonard J. Cleary, Ph.D. All are affiliated with the W.M. Keck Center for the Neurobiology of Learning and Memory.

Byrne and Cleary also are on the faculty of The University of Texas Graduate School of Biomedical Sciences at Houston. Coughlin is a student at the school, which is jointly operated by UT Health and The University of Texas MD Anderson Cancer Center.

The study titled “Doxorubicin Attenuates Serotonin-Induced Long-Term Synaptic Facilitation by Phosphorylation of p38 Mitogen-Activated Protein Kinase” received support from National Institutes of Health grant (NS019895) and the Zilkha Family Discovery Fellowship.

Doxorubicin Attenuates Serotonin-Induced Long-Term Synaptic Facilitation by Phosphorylation of p38 Mitogen-Activated Protein Kinase

Source: Univ. of Texas Health Science Center at Houston

http://www.rdmag.com/news/2014/10/neuroscientists-use-snail-research-help-explain-E2_9_Cchemo-brain

Doxorubicin Attenuates Serotonin-Induced Long-Term Synaptic Facilitation by Phosphorylation of p38 Mitogen-Activated Protein Kinase

Rong-Yu Liu*,  Yili Zhang*,  Brittany L. Coughlin,  Leonard J. Cleary, and  John H. Byrne   +Show Affiliations
The Journal of Neuroscience, 1 Oct 2014, 34(40): 13289-13300;
http://dx.doi.org:/10.1523/JNEUROSCI.0538-14.2014

Doxorubicin (DOX) is an anthracycline used widely for cancer chemotherapy. Its primary mode of action appears to be

  • topoisomerase II inhibition, DNA cleavage, and free radical generation.

However, in non-neuronal cells, DOX also inhibits the expression of

  • dual-specificity phosphatases (also referred to as MAPK phosphatases) and thereby
  1. inhibits the dephosphorylation of extracellular signal-regulated kinase (ERK) and
  2. p38 mitogen-activated protein kinase (p38 MAPK),
  3. two MAPK isoforms important for long-term memory (LTM) formation.

Activation of these kinases by DOX in neurons, if present,

  • could have secondary effects on cognitive functions, such as learning and memory.

The present study used cultures of rat cortical neurons and sensory neurons (SNs) of Aplysia

  • to examine the effects of DOX on levels of phosphorylated ERK (pERK) and
  • phosphorylated p38 (p-p38) MAPK.

In addition, Aplysia neurons were used to examine the effects of DOX on

  • long-term enhanced excitability, long-term synaptic facilitation (LTF), and
  • long-term synaptic depression (LTD).

DOX treatment led to elevated levels of

  • pERK and p-p38 MAPK in SNs and cortical neurons.

In addition, it increased phosphorylation of

  • the downstream transcriptional repressor cAMP response element-binding protein 2 in SNs.

DOX treatment blocked serotonin-induced LTF and enhanced LTD induced by the neuropeptide Phe-Met-Arg-Phe-NH2. The block of LTF appeared to be attributable to

  • overriding inhibitory effects of p-p38 MAPK, because
  • LTF was rescued in the presence of an inhibitor of p38 MAPK
    (SB203580 [4-(4-fluorophenyl)-2-(4-methylsulfinylphenyl)-5-(4-pyridyl)-1H-imidazole]) .

These results suggest that acute application of DOX might impair the formation of LTM via the p38 MAPK pathway.
Terms: Aplysia chemotherapy ERK  p38 MAPK serotonin synaptic plasticity

Technology that controls brain cells with radio waves earns early BRAIN grant

10/08/2014

bright spots = cells with increased calcium after treatment with radio waves,  allows neurons to fire

bright spots = cells with increased calcium after treatment with radio waves, allows neurons to fire

BRAIN control: The new technology uses radio waves to activate or silence cells remotely. The bright spots above represent cells with increased calcium after treatment with radio waves, a change that would allow neurons to fire.

A proposal to develop a new way to

  • remotely control brain cells

from Sarah Stanley, a research associate in Rockefeller University’s Laboratory of Molecular Genetics, headed by Jeffrey M. Friedman, is

  • among the first to receive funding from U.S. President Barack Obama’s BRAIN initiative.

The project will make use of a technique called

  • radiogenetics that combines the use of radio waves or magnetic fields with
  • nanoparticles to turn neurons on or off.

The National Institutes of Health is one of four federal agencies involved in the BRAIN (Brain Research through Advancing Innovative Neurotechnologies) initiative. Following in the ambitious footsteps of the Human Genome Project, the BRAIN initiative seeks

  • to create a dynamic map of the brain in action,

a goal that requires the development of new technologies. The BRAIN initiative working group, which outlined the broad scope of the ambitious project, was co-chaired by Rockefeller’s Cori Bargmann, head of the Laboratory of Neural Circuits and Behavior.

Stanley’s grant, for $1.26 million over three years, is one of 58 projects to get BRAIN grants, the NIH announced. The NIH’s plan for its part of this national project, which has been pitched as “America’s next moonshot,” calls for $4.5 billion in federal funds over 12 years.

The technology Stanley is developing would

  • enable researchers to manipulate the activity of neurons, as well as other cell types,
  • in freely moving animals in order to better understand what these cells do.

Other techniques for controlling selected groups of neurons exist, but her new nanoparticle-based technique has a

  • unique combination of features that may enable new types of experimentation.
  • it would allow researchers to rapidly activate or silence neurons within a small area of the brain or
  • dispersed across a larger region, including those in difficult-to-access locations.

Stanley also plans to explore the potential this method has for use treating patients.

“Francis Collins, director of the NIH, has discussed

  • the need for studying the circuitry of the brain,
  • which is formed by interconnected neurons.

Our remote-control technology may provide a tool with which researchers can ask new questions about the roles of complex circuits in regulating behavior,” Stanley says.
Rockefeller University’s Laboratory of Molecular Genetics
Source: Rockefeller Univ.

Part 4.  Cancer

Two Proteins Found to Block Cancer Metastasis

Why do some cancers spread while others don’t? Scientists have now demonstrated that

  • metastatic incompetent cancers actually “poison the soil”
  • by generating a micro-environment that blocks cancer cells
  • from settling and growing in distant organs.

The “seed and the soil” hypothesis proposed by Stephen Paget in 1889 is now widely accepted to explain how

  • cancer cells (seeds) are able to generate fertile soil (the micro-environment)
  • in distant organs that promotes cancer’s spread.

However, this concept had not explained why some tumors do not spread or metastasize.

The researchers, from Weill Cornell Medical College, found that

  • two key proteins involved in this process work by
  • dramatically suppressing cancer’s spread.

The study offers hope that a drug based on these

  • potentially therapeutic proteins, prosaposin and Thrombospondin 1 (Tsp-1),

might help keep human cancer at bay and from metastasizing.

Scientists don’t understand why some tumors wouldn’t “want” to spread. It goes against their “job description,” says the study’s senior investigator, Vivek Mittal, Ph.D., an associate professor of cell and developmental biology in cardiothoracic surgery and director of the Neuberger Berman Foundation Lung Cancer Laboratory at Weill Cornell Medical College. He theorizes that metastasis occurs when

  • the barriers that the body throws up to protect itself against cancer fail.

But there are some tumors in which some of the barriers may still be intact. “So that suggests

  • those primary tumors will continue to grow, but that
  • an innate protective barrier still exists that prevents them from spreading and invading other organs,”

The researchers found that, like typical tumors,

  • metastasis-incompetent tumors also send out signaling molecules
  • that establish what is known as the “premetastatic niche” in distant organs.

These niches composed of bone marrow cells and various growth factors have been described previously by others including Dr. Mittal as the fertile “soil” that the disseminated cancer cell “seeds” grow in.

Weill Cornell’s Raúl Catena, Ph.D., a postdoctoral fellow in Dr. Mittal’s laboratory, found an important difference between the tumor types. Metastatic-incompetent tumors

  • systemically increased expression of Tsp-1, a molecule known to fight cancer growth.
  • increased Tsp-1 production was found specifically in the bone marrow myeloid cells
  • that comprise the metastatic niche.

These results were striking, because for the first time Dr. Mittal says

  • the bone marrow-derived myeloid cells were implicated as
  • the main producers of Tsp-1,.

In addition, Weill Cornell and Harvard researchers found that

  • prosaposin secreted predominantly by the metastatic-incompetent tumors
  • increased expression of Tsp-1 in the premetastatic lungs.

Thus, Dr. Mittal posits that prosaposin works in combination with Tsp-1

  • to convert pro-metastatic bone marrow myeloid cells in the niche
  • into cells that are not hospitable to cancer cells that spread from a primary tumor.
  • “The very same myeloid cells in the niche that we know can promote metastasis
  • can also be induced under the command of the metastatic incompetent primary tumor to inhibit metastasis,”

The research team found that

  • the Tsp-1–inducing activity of prosaposin
  • was contained in only a 5-amino acid peptide region of the protein, and
  • this peptide alone induced Tsp-1 in the bone marrow cells and
  • effectively suppressed metastatic spread in the lungs
  • in mouse models of breast and prostate cancer.

This 5-amino acid peptide with Tsp-1–inducing activity

  • has the potential to be used as a therapeutic agent against metastatic cancer,

The scientists have begun to test prosaposin in other tumor types or metastatic sites.

Dr. Mittal says that “The clinical implications of the study are:

  • “Not only is it theoretically possible to design a prosaposin-based drug or drugs
  • that induce Tsp-1 to block cancer spread, but
  • you could potentially create noninvasive prognostic tests
  • to predict whether a cancer will metastasize.”

The study was reported in the April 30 issue of Cancer Discovery, in a paper titled “Bone Marrow-Derived Gr1+ Cells Can Generate a Metastasis-Resistant Microenvironment Via Induced Secretion of Thrombospondin-1”.

Disabling Enzyme Cripples Tumors, Cancer Cells

First Step of Metastasis

First Step of Metastasis

Published: Sep 05, 2013  http://www.technologynetworks.com/Metabolomics/news.aspx?id=157138

Knocking out a single enzyme dramatically cripples the ability of aggressive cancer cells to spread and grow tumors.

The paper, published in the journal Proceedings of the National Academy of Sciences, sheds new light on the importance of lipids, a group of molecules that includes fatty acids and cholesterol, in the development of cancer.

Researchers have long known that cancer cells metabolize lipids differently than normal cells. Levels of ether lipids – a class of lipids that are harder to break down – are particularly elevated in highly malignant tumors.

“Cancer cells make and use a lot of fat and lipids, and that makes sense because cancer cells divide and proliferate at an accelerated rate, and to do that,

  • they need lipids, which make up the membranes of the cell,”

said study principal investigator Daniel Nomura, assistant professor in UC Berkeley’s Department of Nutritional Sciences and Toxicology. “Lipids have a variety of uses for cellular structure, but what we’re showing with our study is that

  • lipids can send signals that fuel cancer growth.”

In the study, Nomura and his team tested the effects of reducing ether lipids on human skin cancer cells and primary breast tumors. They targeted an enzyme,

  • alkylglycerone phosphate synthase, or AGPS,
  • known to be critical to the formation of ether lipids.

The researchers confirmed that

  1. AGPS expression increased when normal cells turned cancerous.
  2. inactivating AGPS substantially reduced the aggressiveness of the cancer cells.

“The cancer cells were less able to move and invade,” said Nomura.

The researchers also compared the impact of

  • disabling the AGPS enzyme in mice that had been injected with cancer cells.

Nomura. observes -“Among the mice that had the AGPS enzyme inactivated,

  • the tumors were nonexistent,”

“The mice that did not have this enzyme

  • disabled rapidly developed tumors.”

The researchers determined that

  • inhibiting AGPS expression depleted the cancer cells of ether lipids.
  • AGPS altered levels of other types of lipids important to the ability of the cancer cells to survive and spread, including
    • prostaglandins and acyl phospholipids.

“What makes AGPS stand out as a treatment target is that the enzyme seems to simultaneously

  • regulate multiple aspects of lipid metabolism
  • important for tumor growth and malignancy.”

Future steps include the

  • development of AGPS inhibitors for use in cancer therapy,

“This study sheds considerable light on the important role that AGPS plays in ether lipid metabolism in cancer cells, and it suggests that

  • inhibitors of this enzyme could impair tumor formation,”

said Benjamin Cravatt, Professor and Chair of Chemical Physiology at The Scripps Research Institute, who is not part of the UC.

Agilent Technologies Thought Leader Award Supports Translational Research Program
Published: Mon, March 04, 2013

The award will support Dr DePinho’s research into

  • metabolic reprogramming in the earliest stages of cancer.

Agilent Technologies Inc. announces that Dr. Ronald A. DePinho, a world-renowned oncologist and researcher, has received an Agilent Thought Leader Award.

DePinho is president of the University of Texas MD Anderson Cancer Center. DePinho and his team hope to discover and characterize

  • alterations in metabolic flux during tumor initiation and maintenance, and to identify biomarkers for early detection of pancreatic cancer together with
  • novel therapeutic targets.

Researchers on his team will work with scientists from the university’s newly formed Institute of Applied Cancer Sciences.

The Agilent Thought Leader Award provides funds to support personnel as well as a state-of-the-art Agilent 6550 iFunnel Q-TOF LC/MS system.

“I am extremely pleased to receive this award for metabolomics research, as the survival rates for pancreatic cancer have not significantly improved over the past 20 years,” DePinho said. “This technology will allow us to

  • rapidly identify new targets that drive the formation, progression and maintenance of pancreatic cancer.

Discoveries from this research will also lead to

  • the development of effective early detection biomarkers and novel therapeutic interventions.”

“We are proud to support Dr. DePinho’s exciting translational research program, which will make use of

  • metabolomics and integrated biology workflows and solutions in biomarker discovery,”

said Patrick Kaltenbach, Agilent vice president, general manager of the Liquid Phase Division, and the executive sponsor of this award.

The Agilent Thought Leader Program promotes fundamental scientific advances by support of influential thought leaders in the life sciences and chemical analysis fields.

The covalent modifier Nedd8 is critical for the activation of Smurf1 ubiquitin ligase in tumorigenesis

Ping Xie, Minghua Zhang, Shan He, Kefeng Lu, Yuhan Chen, Guichun Xing, et al.
Nature Communications
  2014; 5(3733).  http://dx.doi.org:/10.1038/ncomms4733

Neddylation, the covalent attachment of ubiquitin-like protein Nedd8, of the Cullin-RING E3 ligase family

  • regulates their ubiquitylation activity.

However, regulation of HECT ligases by neddylation has not been reported to date. Here we show that

  • the C2-WW-HECT ligase Smurf1 is activated by neddylation.

Smurf1 physically interacts with

  1. Nedd8 and Ubc12,
  2. forms a Nedd8-thioester intermediate, and then
  3. catalyses its own neddylation on multiple lysine residues.

Intriguingly, this autoneddylation needs

  • an active site at C426 in the HECT N-lobe.

Neddylation of Smurf1 potently enhances

  • ubiquitin E2 recruitment and
  • augments the ubiquitin ligase activity of Smurf1.

The regulatory role of neddylation

  • is conserved in human Smurf1 and yeast Rsp5.

Furthermore, in human colorectal cancers,

  • the elevated expression of Smurf1, Nedd8, NAE1 and Ubc12
  • correlates with cancer progression and poor prognosis.

These findings provide evidence that

  • neddylation is important in HECT ubiquitin ligase activation and
  • shed new light on the tumour-promoting role of Smurf1.
 Swinging domains in HECT E3

Swinging domains in HECT E3

Subject terms: Biological sciences Cancer Cell biology

Figure 1: Smurf1 expression is elevated in colorectal cancer tissues.

Smurf1 expression is elevated in colorectal cancer tissues.

Smurf1 expression is elevated in colorectal cancer tissues.

(a) Smurf1 expression scores are shown as box plots, with the horizontal lines representing the median; the bottom and top of the boxes representing the 25th and 75th percentiles, respectively; and the vertical bars representing the ra

Figure 2: Positive correlation of Smurf1 expression with Nedd8 and its interacting enzymes in colorectal cancer.

Positive correlation of Smurf1 expression with Nedd8 and its interacting enzymes in colorectal cancer

Positive correlation of Smurf1 expression with Nedd8 and its interacting enzymes in colorectal cancer

(a) Representative images from immunohistochemical staining of Smurf1, Ubc12, NAE1 and Nedd8 in the same colorectal cancer tumour. Scale bars, 100 μm. (bd) The expression scores of Nedd8 (b, n=283 ), NAE1 (c, n=281) and Ubc12 (d, n=19…

Figure 3: Smurf1 interacts with Ubc12.

Smurf1 interacts with Ubc12

Smurf1 interacts with Ubc12

(a) GST pull-down assay of Smurf1 with Ubc12. Both input and pull-down samples were subjected to immunoblotting with anti-His and anti-GST antibodies. Smurf1 interacted with Ubc12 and UbcH5c, but not with Ubc9. (b) Mapping the regions…

Figure 4: Nedd8 is attached to Smurf1through C426-catalysed autoneddylation.

Nedd8 is attached to Smurf1through C426-catalysed autoneddylation

Nedd8 is attached to Smurf1through C426-catalysed autoneddylation

(a) Covalent neddylation of Smurf1 in vitro.Purified His-Smurf1-WT or C699A proteins were incubated with Nedd8 and Nedd8-E1/E2. Reactions were performed as described in the Methods section. Samples were analysed by western blotting wi…

Figure 5: Neddylation of Smurf1 activates its ubiquitin ligase activity.

Neddylation of Smurf1 activates its ubiquitin ligase activity.

Neddylation of Smurf1 activates its ubiquitin ligase activity.

(a) In vivo Smurf1 ubiquitylation assay. Nedd8 was co-expressed with Smurf1 WT or C699A in HCT116 cells (left panels). Twenty-four hours post transfection, cells were treated with MG132 (20 μM, 8 h). HCT116 cells were transfected with…

http://www.nature.com/ncomms/2014/140513/ncomms4733/carousel/ncomms4733-f1.jpg

http://www.nature.com/ncomms/2014/140513/ncomms4733/carousel/ncomms4733-f2.jpg

http://www.nature.com/ncomms/2014/140513/ncomms4733/carousel/ncomms4733-f3.jpg

http://www.nature.com/ncomms/2014/140513/ncomms4733/carousel/ncomms4733-f4.jpg

http://www.nature.com/ncomms/2014/140513/ncomms4733/carousel/ncomms4733-f5.jpg

http://www.nature.com/ncomms/2014/140513/ncomms4733/carousel/ncomms4733-f6.jpg

http://www.nature.com/ncomms/2014/140513/ncomms4733/carousel/ncomms4733-f7.jpg

http://www.nature.com/ncomms/2014/140513/ncomms4733/carousel/ncomms4733-f8.jpg

The deubiquitylase USP33 discriminates between RALB functions in autophagy and innate immune response

M Simicek, S Lievens, M Laga, D Guzenko, VN. Aushev, et al.
Nature Cell Biology 2013; 15, 1220–1230    http://dx.doi.org:/10.1038/ncb2847

The RAS-like GTPase RALB mediates cellular responses to nutrient availability or viral infection by respectively

  • engaging two components of the exocyst complex, EXO84 and SEC5.
  1. RALB employs SEC5 to trigger innate immunity signalling, whereas
  2. RALB–EXO84 interaction induces autophagocytosis.

How this differential interaction is achieved molecularly by the RAL GTPase remains unknown.

We found that whereas GTP binding

  • turns on RALB activity,

ubiquitylation of RALB at Lys 47

  • tunes its activity towards a particular effector.

Specifically, ubiquitylation at Lys 47

  • sterically inhibits RALB binding to EXO84, while
  • facilitating its interaction with SEC5.

Double-stranded RNA promotes

  • RALB ubiquitylation and
  • SEC5–TBK1 complex formation.

In contrast, nutrient starvation

  • induces RALB deubiquitylation
  • by accumulation and relocalization of the deubiquitylase USP33
  • to RALB-positive vesicles.

Deubiquitylated RALB

  • promotes the assembly of the RALB–EXO84–beclin-1 complexes
  • driving autophagosome formation. Thus,
  • ubiquitylation within the effector-binding domain
  • provides the switch for the dual functions of RALB in
    • autophagy and innate immune responses.

Part 5. Metabolic Syndrome

Single Enzyme is Necessary for Development of Diabetes

Published: Aug 20, 2014 http://www.technologynetworks.com/Metabolomics/news.aspx?ID=169416

12-LO enzyme promotes the obesity-induced oxidative stress in the pancreatic cells.

An enzyme called 12-LO promotes the obesity-induced oxidative stress in the pancreatic cells that leads

  • to pre-diabetes, and diabetes.

12-LO’s enzymatic action is the last step in

  • the production of certain small molecules that harm the cell,

according to a team from Indiana University School of Medicine, Indianapolis.

The findings will enable the development of drugs that can interfere with this enzyme, preventing or even reversing diabetes. The research is published ahead of print in the journal Molecular and Cellular Biology.

In earlier studies, these researchers and their collaborators at Eastern Virginia Medical School showed that

  • 12-LO (which stands for 12-lipoxygenase) is present in these cells
  • only in people who become overweight.

The harmful small molecules resulting from 12-LO’s enzymatic action are known as HETEs, short for hydroxyeicosatetraenoic acid.

  1. HETEs harm the mitochondria, which then
  2. fail to produce sufficient energy to enable
  3. the pancreatic cells to manufacture the necessary quantities of insulin.

For the study, the investigators genetically engineered mice that

  • lacked the gene for 12-LO exclusively in their pancreas cells.

Mice were either fed a low-fat or high-fat diet.

Both the control mice and the knockout mice on the high fat diet

  • developed obesity and insulin resistance.

The investigators also examined the pancreatic beta cells of both knockout and control mice, using both microscopic studies and molecular analysis. Those from the knockout mice were intact and healthy, while

  • those from the control mice showed oxidative damage,
  • demonstrating that 12-LO and the resulting HETEs
  • caused the beta cell failure.

Mirmira notes that fatty diet used in the study was the Western Diet, which comprises mostly saturated-“bad”-fats. Based partly on a recent study of related metabolic pathways, he says that

  • the unsaturated and mono-unsaturated fats-which comprise most fats in the healthy,
  • relatively high fat Mediterranean diet-are unlikely to have the same effects.

“Our research is the first to show that 12-LO in the beta cell

  • is the culprit in the development of pre-diabetes, following high fat diets,” says Mirmira.

“Our work also lends important credence to the notion that

  • the beta cell is the primary defective cell in virtually all forms of diabetes and pre-diabetes.”

A New Player in Lipid Metabolism Discovered

Published: Aug18, 2014  http://www.technologynetworks.com/Metabolomics/news.aspx?ID=169356

Specially engineered mice gained no weight, and normal counterparts became obese

  • on the same high-fat, obesity-inducing Western diet.

Specially engineered mice that lacked a particular gene did not gain weight

  • when fed a typical high-fat, obesity-inducing Western diet.

Yet, these mice ate the same amount as their normal counterparts that became obese.

The mice were engineered with fat cells that lacked a gene called SEL1L,

  • known to be involved in the clearance of mis-folded proteins
  • in the cell’s protein making machinery called the endoplasmic reticulum (ER).

When mis-folded proteins are not cleared but accumulate,

  • they destroy the cell and contribute to such diseases as
  1. mad cow disease,
  2. Type 1 diabetes and
  3. cystic fibrosis.

“The million-dollar question is why don’t these mice gain weight? Is this related to its inability to clear mis-folded proteins in the ER?” said Ling Qi, associate professor of molecular and biochemical nutrition and senior author of the study published online July 24 in Cell Metabolism. Haibo Sha, a research associate in Qi’s lab, is the paper’s lead author.

Interestingly, the experimental mice developed a host of other problems, including

  • postprandial hypertriglyceridemia,
  • and fatty livers.

“Although we are yet to find out whether these conditions contribute to the lean phenotype, we found that

  • there was a lipid partitioning defect in the mice lacking SEL1L in fat cells,
  • where fat cells cannot store fat [lipids], and consequently
  • fat goes to the liver.

During the investigation of possible underlying mechanisms, we discovered

  • a novel function for SEL1L as a regulator of lipid metabolism,” said Qi.

Sha said “We were very excited to find that

  • SEL1L is required for the intracellular trafficking of
  • lipoprotein lipase (LPL), acting as a chaperone,” .

and added that “Using several tissue-specific knockout mouse models,

  • we showed that this is a general phenomenon,”

Without LPL, lipids remain in the circulation;

  • fat and muscle cells cannot absorb fat molecules for storage and energy combustion,

People with LPL mutations develop

  • postprandial hypertriglyceridemia similar to
  • conditions found in fat cell-specific SEL1L-deficient mice, said Qi.

Future work will investigate the

  • role of SEL1L in human patients carrying LPL mutations and
  • determine why fat cell-specific SEL1L-deficient mice remain lean under Western diets, said Sha.

Co-authors include researchers from Cedars-Sinai Medical Center in Los Angeles; Wageningen University in the Netherlands; Georgia State University; University of California, Los Angeles; and the Medical College of Soochow University in China.

The study was funded by the U.S. National Institutes of Health, the Netherlands Organization for Health Research and Development National Institutes of Health, the Cedars-Sinai Medical Center, Chinese National Science Foundation, the American Diabetes Association, Cornell’s Center for Vertebrate Genomics and the Howard Hughes Medical Institute.

Part 6. Biomarkers

Biomarkers Take Center Stage

Josh P. Roberts
GEN May 1, 2013 (Vol. 33, No. 9)  http://www.genengnews.com/

While work with biomarkers continues to grow, scientists are also grappling with research-related bottlenecks, such as

  1. affinity reagent development,
  2. platform reproducibility, and
  3. sensitivity.

Biomarkers by definition indicate some state or process that generally occurs

  • at a spatial or temporal distance from the marker itself, and

it would not be an exaggeration to say that biomedicine has become infatuated with them:

  1. where to find them,
  2. when they may appear,
  3. what form they may take, and
  4. how they can be used to diagnose a condition or
  5. predict whether a therapy may be successful.

Biomarkers are on the agenda of many if not most industry gatherings, and in cases such as Oxford Global’s recent “Biomarker Congress” and the GTC “Biomarker Summit”, they hold the naming rights. There, some basic principles were built upon, amended, and sometimes challenged.

In oncology, for example, biomarker discovery is often predicated on the premise that

  • proteins shed from a tumor will traverse to and persist in, and be detectable in, the circulation.

By quantifying these proteins—singularly or as part of a larger “signature”—the hope is

  1. to garner information about the molecular characteristics of the cancer
  2. that will help with cancer detection and
  3. personalization of the treatment strategy.

Yet this approach has not yet turned into the panacea that was hoped for. Bottlenecks exist in

  • affinity reagent development,
  • platform reproducibility, and
  • sensitivity.

There is also a dearth of understanding of some of the

  • fundamental principles of biomarker biology that we need to know the answers to,

said Parag Mallick, Ph.D., whose lab at Stanford University is “working on trying to understand where biomarkers come from.”

There are dogmas saying that

  • circulating biomarkers come solely from secreted proteins.

But Dr. Mallick’s studies indicate that fully

  • 50% of circulating proteins may come from intracellular sources or
  • proteins that are annotated as such.

“We don’t understand the processes governing

  • which tumor-derived proteins end up in the blood.”

Other questions include “how does the size of a tumor affect how much of a given protein will be in the blood?”—perhaps

  • the tumor is necrotic at the center, or
  • it’s hypervascular or hypovascular.

He points out “The problem is that these are highly nonlinear processes at work, and

  • there is a large number of factors that might affect the answer to that question,” .

Their research focuses on using

  1. mass spectrometry and
  2. computational analysis
  • to characterize the biophysical properties of the circulating proteome, and
  • relate these to measurements made of the tumor itself.

Furthermore, he said – “We’ve observed that the proteins that are likely to

  • first show up and persist in the circulation, ..
  • are more stable than proteins that don’t,”
  • “we can quantify how significant the effect is.”

The goal is ultimately to be able to

  1. build rigorous, formal mathematical models that will allow something measured in the blood
  2. to be tied back to the molecular biology taking place in the tumor.

And conversely, to use those models

  • to predict from a tumor what will be found in the circulation.

“Ultimately, the models will allow you to connect the dots between

  • what you measure in the blood and the biology of the tumor.”

Bound for Affinity Arrays

Affinity reagents are the main tools for large-scale protein biomarker discovery. And while this has tended to mean antibodies (or their derivatives), other affinity reagents are demanding a place in the toolbox.

Affimers, a type of affinity reagent being developed by Avacta, consist of

  1. a biologically inert, biophysically stable protein scaffold
  2. containing three variable regions into which
  3. distinct peptides are inserted.

The resulting three-dimensional surface formed by these peptides

  • interacts and binds to proteins and other molecules in solution,
  • much like the antigen-binding site of antibodies.

Unlike antibodies, Affimers are relatively small (13 KDa),

  • non-post-translationally modified proteins
  • that can readily be expressed in bacterial culture.

They may be made to bind surfaces through unique residues

  • engineered onto the opposite face of the Affimer,
  • allowing the binding site to be exposed to the target in solution.

“We don’t seem to see in what we’ve done so far

  • any real loss of activity or functionality of Affimers when bound to surfaces—

they’re very robust,” said CEO Alastair Smith, Ph.D.

Avacta is taking advantage of this stability and its large libraries of Affimers to develop

  • very large affinity microarrays for
  • drug and biomarker discovery.

To date they have printed arrays with around 20–25,000 features, and Dr. Smith is “sure that we can get toward about 50,000 on a slide,” he said. “There’s no real impediment to us doing that other than us expressing the proteins and getting on with it.”

Customers will be provided with these large, complex “naïve” discovery arrays, readable with standard equipment. The plan is for the company to then “support our customers by providing smaller arrays with

  • the Affimers that are binding targets of interest to them,” Dr. Smith foretold.

And since the intellectual property rights are unencumbered,

  • Affimers in those arrays can be licensed to the end users
  • to develop diagnostics that can be validated as time goes on.

Around 20,000-Affimer discovery arrays were recently tested by collaborator Professor Ann Morgan of the University of Leeds with pools of unfractionated serum from patients with symptoms of inflammatory disease. The arrays

  • “rediscovered” elevated C-reactive protein (CRP, the clinical gold standard marker)
  • as well as uncovered an additional 22 candidate biomarkers.
  • other candidates combined with CRP, appear able to distinguish between different diseases such as
  1. rheumatoid arthritis,
  2. psoriatic arthritis,
  3. SLE, or
  4. giant cell arteritis.

Epigenetic Biomarkers

Methylation of adenine

Sometimes biomarkers are used not to find disease but

  • to distinguish healthy human cell types, with
  •  examples being found in flow cytometry and immunohistochemistry.

These widespread applications, however, are difficult to standardize, being

  • subject to arbitrary or subjective gating protocols and other imprecise criteria.

Epiontis instead uses an epigenetic approach. “What we need is a unique marker that is

  • demethylated only in one cell type and
  • methylated in all the other cell types,”

Each cell of the right cell type will have

  • two demethylated copies of a certain gene locus,
  • allowing them to be enumerated by quantitative PCR.

The biggest challenge is finding that unique epigenetic marker. To do so they look through the literature for proteins and genes described as playing a role in the cell type’s biology, and then

  • look at the methylation patterns to see if one can be used as a marker,

They also “use customized Affymetrix chips to look at the

  • differential epigenetic status of different cell types on a genomewide scale.”

explained CBO and founder Ulrich Hoffmueller, Ph.D.

The company currently has a panel of 12 assays for 12 immune cell types. Among these is an assay for

  • regulatory T (Treg) cells that queries the Foxp3 gene—which is uniquely demethylated in Treg
  • even though it is transiently expressed in activated T cells of other subtypes.

Also assayed are Th17 cells, difficult to detect by flow cytometry because

  • “the cells have to be stimulated in vitro,” he pointed out.

Developing New Assays for Cancer Biomarkers

Researchers at Myriad RBM and the Cancer Prevention Research Institute of Texas are collaborating to develop

  • new assays for cancer biomarkers on the Myriad RBM Multi-Analyte Profile (MAP) platform.

The release of OncologyMAP 2.0 expanded Myriad RBM’s biomarker menu to over 250 analytes, which can be measured from a small single sample, according to the company. Using this menu, L. Stephen et al., published a poster, “Analysis of Protein Biomarkers in Prostate and Colorectal Tumor Lysates,” which showed the results of

  • a survey of proteins relevant to colorectal (CRC) and prostate (PC) tumors
  • to identify potential proteins of interest for cancer research.

The study looked at CRC and PC tumor lysates and found that 102 of the 115 proteins showed levels above the lower limit of quantification.

  • Four markers were significantly higher in PC and 10 were greater in CRC.

For most of the analytes, duplicate sections of the tumor were similar, although some analytes did show differences. In four of the CRC analytes, tumor number four showed differences for CEA and tumor number 2 for uPA.

Thirty analytes were shown to be

  • different in CRC tumor compared to its adjacent tissue.
  • Ten of the analytes were higher in adjacent tissue compared to CRC.
  • Eighteen of the markers examined demonstrated  —-

significant correlations of CRC tumor concentration to serum levels.

“This suggests.. that the Oncology MAP 2.0 platform “provides a good method for studying changes in tumor levels because many proteins can be assessed with a very small sample.”

Clinical Test Development with MALDI-ToF

While there have been many attempts to translate results from early discovery work on the serum proteome into clinical practice, few of these efforts have progressed past the discovery phase.

Matrix-assisted laser desorption/ionization-time of flight (MALDI-ToF) mass spectrometry on unfractionated serum/plasma samples offers many practical advantages over alternative techniques, and does not require

  • a shift from discovery to development and commercialization platforms.

Biodesix claims it has been able to develop the technology into

  • a reproducible, high-throughput tool to
  • routinely measure protein abundance from serum/plasma samples.

“.. we improved data-analysis algorithms to

  • reproducibly obtain quantitative measurements of relative protein abundance from MALDI-ToF mass spectra.

Heinrich Röder, CTO points out that the MALDI-ToF measurements

  • are combined with clinical outcome data using
  • modern learning theory techniques
  • to define specific disease states
  • based on a patient’s serum protein content,”

The clinical utility of the identification of these disease states can be investigated through a retrospective analysis of differing sample sets. For example, Biodesix clinically validated its first commercialized serum proteomic test, VeriStrat®, in 85 different retrospective sample sets.

Röder adds that “It is becoming increasingly clear that

  • the patients whose serum is characterized as VeriStrat Poor show
  • consistently poor outcomes irrespective of
  1. tumor type,
  2. histology, or
  3. molecular tumor characteristics,”

MALDI-ToF mass spectrometry, in its standard implementation,

  • allows for the observation of around 100 mostly high-abundant serum proteins.

Further, “while this does not limit the usefulness of tests developed from differential expression of these proteins,

  • the discovery potential would be greatly enhanced
  • if we could probe deeper into the proteome
  • while not giving up the advantages of the MALDI-ToF approach,”

Biodesix reports that its new MALDI approach, Deep MALDI™, can perform

  • simultaneous quantitative measurement of more than 1,000 serum protein features (or peaks) from 10 µL of serum in a high-throughput manner.
  • it increases the observable signal noise ratio from a few hundred to over 50,000,
  • resulting in the observation of many lower-abundance serum proteins.

Breast cancer, a disease now considered to be a collection of many complexes of symptoms and signatures—the dominant ones are labeled Luminal A, Luminal B, Her2, and Basal— which suggests different prognose, and

  • these labels are considered too simplistic for understanding and managing a woman’s cancer.

Studies published in the past year have looked at

  1. somatic mutations,
  2. gene copy number aberrations,
  3. gene expression abnormalities,
  4. protein and miRNA expression, and
  5. DNA methylation,

coming up with a list of significantly mutated genes—hot spots—in different categories of breast cancers. Targeting these will inevitably be the focus of much coming research.

“We’ve been taking these large trials and profiling these on a variety of array or sequence platforms. We think we’ll get

  1. prognostic drivers
  2. predictive markers for taxanes and
  3. monoclonal antibodies and
  4. tamoxifen and aromatase inhibitors,”
    explained Brian Leyland-Jones, Ph.D., director of Edith Sanford Breast Cancer Research. “We will end up with 20–40 different diseases, maybe more.”

Edith Sanford Breast Cancer Research is undertaking a pilot study in collaboration with The Scripps Research Institute, using a variety of tests on 25 patients to see how the information they provide complements each other, the overall flow, and the time required to get and compile results.

Laser-captured tumor samples will be subjected to low passage whole-genome, exome, and RNA sequencing (with targeted resequencing done in parallel), and reverse-phase protein and phosphorylation arrays, with circulating nucleic acids and circulating tumor cells being queried as well. “After that we hope to do a 100- or 150-patient trial when we have some idea of the best techniques,” he said.

Dr. Leyland-Jones predicted that ultimately most tumors will be found

  • to have multiple drivers,
  • with most patients receiving a combination of two, three, or perhaps four different targeted therapies.

Reduce to Practice

According to Randox, the evidence Investigator is a sophisticated semi-automated biochip sys­tem designed for research, clinical, forensic, and veterinary applications.

Once biomarkers that may have an impact on therapy are discovered, it is not always routine to get them into clinical practice. Leaving regulatory and financial, intellectual property and cultural issues aside, developing a diagnostic based on a biomarker often requires expertise or patience that its discoverer may not possess.

Andrew Gribben is a clinical assay and development scientist at Randox Laboratories, based in Northern Ireland, U.K. The company utilizes academic and industrial collaborators together with in-house discovery platforms to identify biomarkers that are

  • augmented or diminished in a particular pathology
  • relative to appropriate control populations.

Biomarkers can be developed to be run individually or

  • combined into panels of immunoassays on its multiplex biochip array technology.

Specificity can also be gained—or lost—by the affinity of reagents in an assay. The diagnostic potential of Heart-type fatty acid binding protein (H-FABP) abundantly expressed in human myocardial cells was recognized by Jan Glatz of Maastricht University, The Netherlands, back in 1988. Levels rise quickly within 30 minutes after a myocardial infarction, peaking at 6–8 hours and return to normal within 24–30 hours. Yet at the time it was not known that H-FABP was a member of a multiprotein family, with which the polyclonal antibodies being used in development of an assay were cross-reacting, Gribben related.

Randox developed monoclonal antibodies specific to H-FABP, funded trials investigating its use alone, and multiplexed with cardiac biomarker assays, and, more than 30 years after the biomarker was identified, in 2011, released a validated assay for H-FABP as a biomarker for early detection of acute myocardial infarction.

Ultrasensitive Immunoassays for Biomarker Development

Research has shown that detection and monitoring of biomarker concentrations can provide

  • insights into disease risk and progression.

Cytokines have become attractive biomarkers and candidates

  • for targeted therapies for a number of autoimmune diseases, including rheumatoid arthritis (RA), Crohn’s disease, and psoriasis, among others.

However, due to the low-abundance of circulating cytokines, such as IL-17A, obtaining robust measurements in clinical samples has been difficult.

Singulex reports that its digital single-molecule counting technology provides

  • increased precision and detection sensitivity over traditional ELISA techniques,
  • helping to shed light on biomarker verification and validation programs.

The company’s Erenna® immunoassay system, which includes optimized immunoassays, offers LLoQ to femtogram levels per mL resolution—even in healthy populations, at an improvement of 1-3 fold over standard ELISAs or any conventional technology and with a dynamic range of up to 4-logs, according to a Singulex official, who adds that

  • this sensitivity improvement helps minimize undetectable samples that
  • could otherwise delay or derail clinical studies.

The official also explains that the Singulex solution includes an array of products and services that are being applied to a number of programs and have enabled the development of clinically relevant biomarkers, allowing translation from discovery to the clinic.

In a poster entitled “Advanced Single Molecule Detection: Accelerating Biomarker Development Utilizing Cytokines through Ultrasensitive Immunoassays,” a case study was presented of work performed by Jeff Greenberg of NYU to show how the use of the Erenna system can provide insights toward

  • improving the clinical utility of biomarkers and
  • accelerating the development of novel therapies for treating inflammatory diseases.

A panel of inflammatory biomarkers was examined in DMARD (disease modifying antirheumatic drugs)-naïve RA (rheumatoid arthritis) vs. knee OA (osteoarthritis) patient cohorts. Markers that exhibited significant differences in plasma concentrations between the two cohorts included

  • CRP, IL-6R alpha, IL-6, IL-1 RA, VEGF, TNF-RII, and IL-17A, IL-17F, and IL-17A/F.

Among the three tested isoforms of IL-17,

  • the magnitude of elevation for IL-17F in RA patients was the highest.

“Singulex provides high-resolution monitoring of baseline IL-17A concentrations that are present at low levels,” concluded the researchers. “The technology also enabled quantification of other IL-17 isoforms in RA patients, which have not been well characterized before.”

The Singulex Erenna System has also been applied to cardiovascular disease research, for which its

  • cardiac troponin I (cTnI) digital assay can be used to measure circulating
  • levels of cTnI undetectable by other commercial assays.

Recently presented data from Brigham and Women’s Hospital and the TIMI-22 study showed that

  • using the Singulex test to serially monitor cTnI helps
  • stratify risk in post-acute coronary syndrome patients and
  • can identify patients with elevated cTnI
  • who have the most to gain from intensive vs. moderate-dose statin therapy,

according to the scientists involved in the research.

The study poster, “Prognostic Performance of Serial High Sensitivity Cardiac Troponin Determination in Stable Ischemic Heart Disease: Analysis From PROVE IT-TIMI 22,” was presented at the 2013 American College of Cardiology (ACC) Annual Scientific Session & Expo by R. O’Malley et al.

Biomarkers Changing Clinical Medicine

Better Diagnosis, Prognosis, and Drug Targeting Are among Potential Benefits

  1. John Morrow Jr., Ph.D.

Researchers at EMD Chemicals are developing biomarker immunoassays

  • to monitor drug-induced toxicity including kidney damage.

The pace of biomarker development is accelerating as investigators report new studies on cancer, diabetes, Alzheimer disease, and other conditions in which the evaluation and isolation of workable markers is prominently featured.

Wei Zheng, Ph.D., leader of the R&D immunoassay group at EMD Chemicals, is overseeing a program to develop biomarker immunoassays to

  • monitor drug-induced toxicity, including kidney damage.

“One of the principle reasons for drugs failing during development is because of organ toxicity,” says Dr. Zheng.
“proteins liberated into the serum and urine can serve as biomarkers of adverse response to drugs, as well as disease states.”

Through collaborative programs with Rules-Based Medicine (RBM), the EMD group has released panels for the profiling of human renal impairment and renal toxicity. These urinary biomarker based products fit the FDA and EMEA guidelines for assessment of drug-induced kidney damage in rats.

The group recently performed a screen for potential protein biomarkers in relation to

  • kidney toxicity/damage on a set of urine and plasma samples
  • from patients with documented renal damage.

Additionally, Dr. Zheng is directing efforts to move forward with the multiplexed analysis of

  • organ and cellular toxicity.

Diseases thought to involve compromised oxidative phosphorylation include

  • diabetes, Parkinson and Alzheimer diseases, cancer, and the aging process itself.

Good biomarkers allow Dr. Zheng to follow the mantra, “fail early, fail fast.” With robust, multiplexible biomarkers, EMD can detect bad drugs early and kill them before they move into costly large animal studies and clinical trials. “Recognizing the severe liability that toxicity presents, we can modify the structure of the candidate molecule and then rapidly reassess its performance.”

Scientists at Oncogene Science a division of Siemens Healthcare Diagnostics, are also focused on biomarkers. “We are working on a number of antibody-based tests for various cancers, including a test for the Ca-9 CAIX protein, also referred to as carbonic anhydrase,” Walter Carney, Ph.D., head of the division, states.

CAIX is a transmembrane protein that is

  • overexpressed in a number of cancers, and, like Herceptin and the Her-2 gene,
  • can serve as an effective and specific marker for both diagnostic and therapeutic purposes.
  • It is liberated into the circulation in proportion to the tumor burden.

Dr. Carney and his colleagues are evaluating patients after tumor removal for the presence of the Ca-9 CAIX protein. If

  • the levels of the protein in serum increase over time,
  • this suggests that not all the tumor cells were removed and the tumor has metastasized.

Dr. Carney and his team have developed both an immuno-histochemistry and an ELISA test that could be used as companion diagnostics in clinical trials of CAIX-targeted drugs.

The ELISA for the Ca-9 CAIX protein will be used in conjunction with Wilex’ Rencarex®, which is currently in a

  • Phase III trial as an adjuvant therapy for non-metastatic clear cell renal cancer.

Additionally, Oncogene Science has in its portfolio an FDA-approved test for the Her-2 marker. Originally approved for Her-2/Neu-positive breast cancer, its indications have been expanded over time, and was approved

  • for the treatment of gastric cancer last year.

It is normally present on breast cancer epithelia but

  • overexpressed in some breast cancer tumors.

“Our products are designed to be used in conjunction with targeted therapies,” says Dr. Carney. “We are working with companies that are developing technology around proteins that are

  • overexpressed in cancerous tissues and can be both diagnostic and therapeutic targets.”

The long-term goal of these studies is to develop individualized therapies, tailored for the patient. Since the therapies are expensive, accurate diagnostics are critical to avoid wasting resources on patients who clearly will not respond (or could be harmed) by the particular drug.

“At this time the rate of response to antibody-based therapies may be very poor, as

  • they are often employed late in the course of the disease, and patients are in such a debilitated state
  • that they lack the capacity to react positively to the treatment,” Dr. Carney explains.

Nanoscale Real-Time Proteomics

Stanford University School of Medicine researchers, working with Cell BioSciences, have developed a

  • nanofluidic proteomic immunoassay that measures protein charge,
  • similar to immunoblots, mass spectrometry, or flow cytometry.
  • unlike these platforms, this approach can measure the amount of individual isoforms,
  • specifically, phosphorylated molecules.

“We have developed a nanoscale device for protein measurement, which I believe could be useful for clinical analysis,” says Dean W. Felsher, M.D., Ph.D., associate professor at Stanford University School of Medicine.

Critical oncogenic transformations involving

  • the activation of the signal-related kinases ERK-1 and ERK-2 can now be followed with ease.

“The fact that we measure nanoquantities with accuracy means that

  • we can interrogate proteomic profiles in clinical patients,

by drawing tiny needle aspirates from tumors over the course of time,” he explains.

“This allows us to observe the evolution of tumor cells and

  • their response to therapy
  • from a baseline of the normal tissue as a standard of comparison.”

According to Dr. Felsher, 20 cells is a large enough sample to obtain a detailed description. The technology is easy to automate, which allows

  • the inclusion of hundreds of assays.

Contrasting this technology platform with proteomic analysis using microarrays, Dr. Felsher notes that the latter is not yet workable for revealing reliable markers.

Dr. Felsher and his group published a description of this technology in Nature Medicine. “We demonstrated that we could take a set of human lymphomas and distinguish them from both normal tissue and other tumor types. We can

  • quantify changes in total protein, protein activation, and relative abundance of specific phospho-isoforms
  • from leukemia and lymphoma patients receiving targeted therapy.

Even with very small numbers of cells, we are able to show that the results are consistent, and

  • our sample is a random profile of the tumor.”

Splice Variant Peptides

“Aberrations in alternative splicing may generate

  • much of the variation we see in cancer cells,”

says Gilbert Omenn, Ph.D., director of the center for computational medicine and bioinformatics at the University of Michigan School of Medicine. Dr. Omenn and his colleague, Rajasree Menon, are

  • using this variability as a key to new biomarker identification.

It is becoming evident that splice variants play a significant role in the properties of cancer cells, including

  • initiation, progression, cell motility, invasiveness, and metastasis.

Alternative splicing occurs through multiple mechanisms

  • when the exons or coding regions of the DNA transcribe mRNA,
  • generating initiation sites and connecting exons in protein products.

Their translation into protein can result in numerous protein isoforms, and

  • these isoforms may reflect a diseased or cancerous state.

Regulatory elements within the DNA are responsible for selecting different alternatives; thus

  • the splice variants are tempting targets for exploitation as biomarkers.
Analyses of the splice-site mutation

Analyses of the splice-site mutation

Despite the many questions raised by these observations, splice variation in tumor material has not been widely studied. Cancer cells are known for their tremendous variability, which allows them to

  • grow rapidly, metastasize, and develop resistance to anticancer drugs.

Dr. Omenn and his collaborators used

  • mass spec data to interrogate a custom-built database of all potential mRNA sequences
  • to find alternative splice variants.

When they compared normal and malignant mammary gland tissue from a mouse model of Her2/Neu human breast cancers, they identified a vast number (608) of splice variant proteins, of which

  • peptides from 216 were found only in the tumor sample.

“These novel and known alternative splice isoforms

  • are detectable both in tumor specimens and in plasma and
  • represent potential biomarker candidates,” Dr. Omenn adds.

Dr. Omenn’s observations and those of his colleague Lewis Cantley, Ph.D., have also

  • shed light on the origins of the classic Warburg effect,
  • the shift to anaerobic glycolysis in tumor cells.

The novel splice variant M2, of muscle pyruvate kinase,

  • is observed in embryonic and tumor tissue.

It is associated with this shift, the result of

  • the expression of a peptide splice variant sequence.

It is remarkable how many different areas of the life sciences are tied into the phenomenon of splice variation. The changes in the genetic material can be much greater than point mutations, which have been traditionally considered to be the prime source of genetic variability.

“We now have powerful methods available to uncover a whole new category of variation,” Dr. Omenn says. “High-throughput RNA sequencing and proteomics will be complementary in discovery studies of splice variants.”

Splice variation may play an important role in rapid evolutionary changes, of the sort discussed by Susumu Ohno and Stephen J. Gould decades ago. They, and other evolutionary biologists, argued that

  • gene duplication, combined with rapid variability, could fuel major evolutionary jumps.

At the time, the molecular mechanisms of variation were poorly understood, but today

  • the tools are available to rigorously evaluate the role of
  • splice variation and other contributors to evolutionary change.

“Biomarkers derived from studies of splice variants, could, in the future, be exploited

  • both for diagnosis and prognosis and
  • for drug targeting of biological networks,
  • in situations such as the Her-2/Neu breast cancers,” Dr. Omenn says.

Aminopeptidase Activities

“By correlating the proteolytic patterns with disease groups and controls, we have shown that

  • exopeptidase activities contribute to the generation of not only cancer-specific
  • but also cancer type specific serum peptides.

according to Paul Tempst, Ph.D., professor and director of the Protein Center at the Memorial Sloan-Kettering Cancer Center.

So there is a direct link between peptide marker profiles of disease and differential protease activity.” For this reason Dr. Tempst argues that “the patterns we describe may have value as surrogate markers for detection and classification of cancer.”

To investigate this avenue, Dr. Tempst and his colleagues have followed

  • the relationship between exopeptidase activities and metastatic disease.

“We monitored controlled, de novo peptide breakdown in large numbers of biological samples using mass spectrometry, with relative quantitation of the metabolites,” Dr. Tempst explains. This entailed the use of magnetic, reverse-phase beads for analyte capture and a MALDI-TOF MS read-out.

“In biomarker discovery programs, functional proteomics is usually not pursued,” says Dr. Tempst. “For putative biomarkers, one may observe no difference in quantitative levels of proteins, while at the same time, there may be substantial differences in enzymatic activity.”

In a preliminary prostate cancer study, the team found a significant difference

  • in activity levels of exopeptidases in serum from patients with metastatic prostate cancer
  • as compared to primary tumor-bearing individuals and normal healthy controls.

However, there were no differences in amounts of the target protein, and this potential biomarker would have been missed if quantitative levels of protein had been the only criterion of selection.

It is frequently stated that “practical fusion energy is 30 years in the future and always will be.” The same might be said of functional, practical biomarkers that can pass muster with the FDA. But splice variation represents a new handle on this vexing problem. It appears that we are seeing the emergence of a new approach that may finally yield definitive diagnostic tests, detectable in serum and urine samples.

Part 7. Epigenetics and Drug Metabolism

DNA Methylation Rules: Studying Epigenetics with New Tools

The tools to unravel the epigenetic control mechanisms that influence how cells control access of transcriptional proteins to DNA are just beginning to emerge.

Patricia Fitzpatrick Dimond, Ph.D.

http://www.genengnews.com/media/images/AnalysisAndInsight/Feb7_2013_24454248_GreenPurpleDNA_EpigeneticsToolsII3576166141.jpg

New tools may help move the field of epigenetic analysis forward and potentially unveil novel biomarkers for cellular development, differentiation, and disease.

DNA sequencing has had the power of technology behind it as novel platforms to produce more sequencing faster and at lower cost have been introduced. But the tools to unravel the epigenetic control mechanisms that influence how cells control access of transcriptional proteins to DNA are just beginning to emerge.

Among these mechanisms, DNA methylation, or the enzymatically mediated addition of a methyl group to cytosine or adenine dinucleotides,

  • serves as an inherited epigenetic modification that
  • stably modifies gene expression in dividing cells.

The unique methylomes are largely maintained in differentiated cell types, making them critical to understanding the differentiation potential of the cell.

In the DNA methylation process, cytosine residues in the genome are enzymatically modified to 5-methylcytosine,

  • which participates in transcriptional repression of genes during development and disease progression.

5-methylcytosine can be further enzymatically modified to 5-hydroxymethylcytosine by the TET family of methylcytosine dioxygenases. DNA methylation affects gene transcription by physically

  • interfering with the binding of proteins involved in gene transcription.

Methylated DNA may be bound by methyl-CpG-binding domain proteins (MBDs) that can

  • then recruit additional proteins. Some of these include histone deacetylases and other chromatin remodeling proteins that modify histones, thereby
  • forming compact, inactive chromatin, or heterochromatin.

While DNA methylation doesn’t change the genetic code,

  • it influences chromosomal stability and gene expression.

Epigenetics and Cancer Biomarkers

multistage chemical carcinogenesis

multistage chemical carcinogenesis

And because of the increasing recognition that DNA methylation changes are involved in human cancers, scientists have suggested that these epigenetic markers may provide biological markers for cancer cells, and eventually point toward new diagnostic and therapeutic targets. Cancer cell genomes display genome-wide abnormalities in DNA methylation patterns,

  • some of which are oncogenic and contribute to genome instability.

In particular, de novo methylation of tumor suppressor gene promoters

  • occurs frequently in cancers, thereby silencing them and promoting transformation.

Cytosine hydroxymethylation (5-hydroxymethylcytosine, or 5hmC), the aforementioned DNA modification resulting from the enzymatic conversion of 5mC into 5-hydroxymethylcytosine by the TET family of oxygenases, has been identified

  • as another key epigenetic modification marking genes important for
  • pluripotency in embryonic stem cells (ES), as well as in cancer cells.

The base 5-hydroxymethylcytosine was recently identified as an oxidation product of 5-methylcytosine in mammalian DNA. In 2011, using sensitive and quantitative methods to assess levels of 5-hydroxymethyl-2′-deoxycytidine (5hmdC) and 5-methyl-2′-deoxycytidine (5mdC) in genomic DNA, scientists at the Department of Cancer Biology, Beckman Research Institute of the City of Hope, Duarte, California investigated

  • whether levels of 5hmC can distinguish normal tissue from tumor tissue.

They showed that in squamous cell lung cancers, levels of 5hmdC showed

  • up to five-fold reduction compared with normal lung tissue.

In brain tumors,5hmdC showed an even more drastic reduction

  • with levels up to more than 30-fold lower than in normal brain,
  • but 5hmdC levels were independent of mutations in isocitrate dehydrogenase-1, the enzyme that converts 5hmC to 5hmdC.

Immunohistochemical analysis indicated that 5hmC is “remarkably depleted” in many types of human cancer.

  • there was an inverse relationship between 5hmC levels and cell proliferation with lack of 5hmC in proliferating cells.

Their data suggest that 5hmdC is strongly depleted in human malignant tumors,

  • a finding that adds another layer of complexity to the aberrant epigenome found in cancer tissue.

In addition, a lack of 5hmC may become a useful biomarker for cancer diagnosis.

Enzymatic Mapping

But according to New England Biolabs’ Sriharsa Pradhan, Ph.D., methods for distinguishing 5mC from 5hmC and analyzing and quantitating the cell’s entire “methylome” and “hydroxymethylome” remain less than optimal.

The protocol for bisulphite conversion to detect methylation remains the “gold standard” for DNA methylation analysis. This method is generally followed by PCR analysis for single nucleotide resolution to determine methylation across the DNA molecule. According to Dr. Pradhan, “.. bisulphite conversion does not distinguish 5mC and 5hmC,”

Recently we found an enzyme, a unique DNA modification-dependent restriction endonuclease, AbaSI, which can

  • decode the hydryoxmethylome of the mammalian genome.

You easily can find out where the hydroxymethyl regions are.”

AbaSI, recognizes 5-glucosylatedmethylcytosine (5gmC) with high specificity when compared to 5mC and 5hmC, and

  • cleaves at narrow range of distances away from the recognized modified cytosine.

By mapping the cleaved ends, the exact 5hmC location can, the investigators reported, be determined.

Dr. Pradhan and his colleagues at NEB; the Department of Biochemistry, Emory University School of Medicine, Atlanta; and the New England Biolabs Shanghai R&D Center described use of this technique in a paper published in Cell Reports this month, in which they described high-resolution enzymatic mapping of genomic hydroxymethylcytosine in mouse ES cells.

In the current report, the authors used the enzyme technology for the genome-wide high-resolution hydroxymethylome, describing simple library construction even with a low amount of input DNA (50 ng) and the ability to readily detect 5hmC sites with low occupancy.

As a result of their studies, they propose that

factors affecting the local 5mC accessibility to TET enzymes play important roles in the 5hmC deposition

  • including include chromatin compaction, nucleosome positioning, or TF binding.
  •  the regularly oscillating 5hmC profile around the CTCF-binding sites, suggests 5hmC ‘‘writers’’ may be sensitive to the nucleosomal environment.
  • some transiently stable 5hmCs may indicate a poised epigenetic state or demethylation intermediate, whereas others may suggest a locally accessible chromosomal environment for the TET enzymatic apparatus.

“We were able to do complete mapping in mouse embryonic cells and are pleased about what this enzyme can do and how it works,” Dr. Pradhan said.

And the availability of novel tools that make analysis of the methylome and hypomethylome more accessible will move the field of epigenetic analysis forward and potentially novel biomarkers for cellular development, differentiation, and disease.

Patricia Fitzpatrick Dimond, Ph.D. (pdimond@genengnews.com), is technical editor at Genetic Engineering & Biotechnology News.

Epigenetic Regulation of ADME-Related Genes: Focus on Drug Metabolism and Transport

Published: Sep 23, 2013

Epigenetic regulation of gene expression refers to heritable factors that are functionally relevant genomic modifications but that do not involve changes in DNA sequence.

Examples of such modifications include

  • DNA methylation, histone modifications, noncoding RNAs, and chromatin architecture.

Epigenetic modifications are crucial for

packaging and interpreting the genome, and they have fundamental functions in regulating gene expression and activity under the influence of physiologic and environmental factors.

In this issue of Drug Metabolism and Disposition, a series of articles is presented to demonstrate the role of epigenetic factors in regulating

  • the expression of genes involved in drug absorption, distribution, metabolism, and excretion in organ development, tissue-specific gene expression, sexual dimorphism, and in the adaptive response to xenobiotic exposure, both therapeutic and toxic.

The articles also demonstrate that, in addition to genetic polymorphisms, epigenetics may also contribute to wide inter-individual variations in drug metabolism and transport. Identification of functionally relevant epigenetic biomarkers in human specimens has the potential to improve prediction of drug responses based on patient’s epigenetic profiles.

http://www.technologynetworks.com/Metabolomics/news.aspx?ID=157804

This study is published online in Drug Metabolism and Disposition

Part 8.  Pictorial Maps

 Prediction of intracellular metabolic states from extracellular metabolomic data

MK Aurich, G Paglia, Ottar Rolfsson, S Hrafnsdottir, M Magnusdottir, MM Stefaniak, BØ Palsson, RMT Fleming &

Ines Thiele

Metabolomics Aug 14, 2014;

http://dx.doi.org:/10.1007/s11306-014-0721-3

http://link.springer.com/article/10.1007/s11306-014-0721-3/fulltext.html#Sec1

http://link.springer.com/static-content/images/404/art%253A10.1007%252Fs11306-014-0721-3/MediaObjects/11306_2014_721_Fig1_HTML.gif

Metabolic models can provide a mechanistic framework

  • to analyze information-rich omics data sets, and are
  • increasingly being used to investigate metabolic alternations in human diseases.

An expression of the altered metabolic pathway utilization is the selection of metabolites consumed and released by cells. However, methods for the

  • inference of intracellular metabolic states from extracellular measurements in the context of metabolic models remain underdeveloped compared to methods for other omics data.

Herein, we describe a workflow for such an integrative analysis

  • emphasizing on extracellular metabolomics data.

We demonstrate,

  • using the lymphoblastic leukemia cell lines Molt-4 and CCRF-CEM,

how our methods can reveal differences in cell metabolism. Our models explain metabolite uptake and secretion by predicting

  • a more glycolytic phenotype for the CCRF-CEM model and
  • a more oxidative phenotype for the Molt-4 model,
  • which was supported by our experimental data.

Gene expression analysis revealed altered expression of gene products at

  • key regulatory steps in those central metabolic pathways, and

literature query emphasized the role of these genes in cancer metabolism.

Moreover, in silico gene knock-outs identified unique

  •  control points for each cell line model, e.g., phosphoglycerate dehydrogenase for the Molt-4 model.

Thus, our workflow is well suited to the characterization of cellular metabolic traits based on

  • -extracellular metabolomic data, and it allows the integration of multiple omics data sets
  • into a cohesive picture based on a defined model context.

Keywords Constraint-based modeling _ Metabolomics _ Multi-omics _ Metabolic network _ Transcriptomics

1 Introduction

Modern high-throughput techniques have increased the pace of biological data generation. Also referred to as the ‘‘omics avalanche’’, this wealth of data provides great opportunities for metabolic discovery. Omics data sets

  • contain a snapshot of almost the entire repertoire of mRNA, protein, or metabolites at a given time point or

under a particular set of experimental conditions. Because of the high complexity of the data sets,

  • computational modeling is essential for their integrative analysis.

Currently, such data analysis is a bottleneck in the research process and methods are needed to facilitate the use of these data sets, e.g., through meta-analysis of data available in public databases [e.g., the human protein atlas (Uhlen et al. 2010) or the gene expression omnibus (Barrett et al.  2011)], and to increase the accessibility of valuable information for the biomedical research community.

Constraint-based modeling and analysis (COBRA) is

  • a computational approach that has been successfully used to
  • investigate and engineer microbial metabolism through the prediction of steady-states (Durot et al.2009).

The basis of COBRA is network reconstruction: networks are assembled in a bottom-up fashion based on

  • genomic data and extensive
  • organism-specific information from the literature.

Metabolic reconstructions capture information on the

  • known biochemical transformations taking place in a target organism
  • to generate a biochemical, genetic and genomic knowledge base (Reed et al. 2006).

Once assembled, a

  • metabolic reconstruction can be converted into a mathematical model (Thiele and Palsson 2010), and
  • model properties can be interrogated using a great variety of methods (Schellenberger et al. 2011).

The ability of COBRA models

  • to represent genotype–phenotype and environment–phenotype relationships arises
  • through the imposition of constraints, which
  • limit the system to a subset of possible network states (Lewis et al. 2012).

Currently, COBRA models exist for more than 100 organisms, including humans (Duarte et al. 2007; Thiele et al. 2013).

Since the first human metabolic reconstruction was described [Recon 1 (Duarte et al. 2007)],

  • biomedical applications of COBRA have increased (Bordbar and Palsson 2012).

One way to contextualize networks is to

  • define their system boundaries according to the metabolic states of the system, e.g., disease or dietary regimes.

The consequences of the applied constraints can

  • then be assessed for the entire network (Sahoo and Thiele 2013).

Additionally, omics data sets have frequently been used

  • to generate cell-type or condition-specific metabolic models.

Models exist for specific cell types, such as

  1. enterocytes (Sahoo and Thiele2013),
  2. macrophages (Bordbar et al. 2010),
  3. adipocytes (Mardinoglu et al. 2013),
  4. even multi-cell assemblies that represent the interactions of brain cells (Lewis et al. 2010).

All of these cell type specific models, except the enterocyte reconstruction

  • were generated based on omics data sets.

Cell-type-specific models have been used to study

  • diverse human disease conditions.

For example, an adipocyte model was generated using

  • transcriptomic, proteomic, and metabolomics data.

This model was subsequently used to investigate metabolic alternations in adipocytes

  • that would allow for the stratification of obese patients (Mardinoglu et al. 2013).

The biomedical applications of COBRA have been

  1. cancer metabolism (Jerby and Ruppin, 2012).
  2. predicting drug targets (Folger et al. 2011; Jerby et al. 2012).

A cancer model was generated using

  • multiple gene expression data sets and subsequently used
  • to predict synthetic lethal gene pairs as potential drug targets
  • selective for the cancer model, but non-toxic to the global model (Recon 1),

a consequence of the reduced redundancy in the cancer specific model (Folger et al. 2011).

In a follow up study, lethal synergy between FH and enzymes of the heme metabolic pathway

  • were experimentally validated and resolved the mechanism by which FH deficient cells,
    e.g., in renal-cell cancer cells survive a non-functional TCA cycle (Frezza et al. 2011).

Contextualized models, which contain only the subset of reactions active in a particular tissue (or cell-) type,

  • can be generated in different ways (Becker and Palsson, 2008; Jerby et al. 2010).

However, the existing algorithms mainly consider

  • gene expression and proteomic data
  • to define the reaction sets that comprise the contextualized metabolic models.

These subset of reactions are usually defined

  • based on the expression or absence of expression of the genes or proteins (present and absent calls),
  • or inferred from expression values or differential gene expression.

Comprehensive reviews of the methods are available (Blazier and Papin, 2012; Hyduke et al. 2013). Only the compilation of a large set of omics data sets

  • can result in a tissue (or cell-type) specific metabolic model, whereas

the representation of one particular experimental condition is achieved

  • through the integration of omics data set generated from one experiment only (condition-specific cell line model).

Recently, metabolomic data sets have become more comprehensive and

  • using these data sets allow direct determination of the metabolic network components (the metabolites).

Additionally, metabolomics has proven to be stable, relatively inexpensive, and highly reproducible (Antonucci et al. 2012). These factors make metabolomic data sets particularly valuable for

  • interrogation of metabolic phenotypes.

Thus, the integration of these data sets is now an active field of research (Li et al. 2013; Mo et al. 2009; Paglia et al. 2012b; Schmidt et al. 2013).

Generally, metabolomic data can be incorporated into metabolic networks as

  • qualitative, quantitative, and thermodynamic constraints (Fleming et al. 2009; Mo et al. 2009).

Mo et al. used metabolites detected in the

  • spent medium of yeast cells to determine intracellular flux states through a sampling analysis (Mo et al. 2009),
  • which allowed unbiased interrogation of the possible network states (Schellenberger and Palsson 2009) and
  • prediction of internal pathway use.
Modes of transcriptional regulation during the YMC

Modes of transcriptional regulation during the YMC

Such analyses have also been used to reveal the effects of

  1. enzymopathies on red blood cells (Price et al. 2004),
  2. to study effects of diet on diabetes (Thiele et al. 2005) and
  3. to define macrophage metabolic states (Bordbar et al. 2010).

This type of analysis is available as a function in the COBRA toolbox (Schellenberger et al. 2011).

In this study, we established a workflow

  • for the generation and analysis of condition-specific metabolic cell line models
  • that can facilitate the interpretation of metabolomic data.

Our modeling yields meaningful predictions regarding

  • metabolic differences between two lymphoblastic leukemia cell lines (Fig. 1A).

Fig. 1

metabol leukem cell lines11306_2014_721_Fig1_HTML

metabol leukem cell lines11306_2014_721_Fig1_HTML

A Combined experimental and computational pipeline to study human metabolism.

  1. Experimental work and omics data analysis steps precede computational modeling.
  2. Model predictions are validated based on targeted experimental data.
  3. Metabolomic and transcriptomic data are used for model refinement and submodel extraction.
  4. Functional analysis methods are used to characterize the metabolism of the cell-line models and compare it to additional experimental data.
  5. The validated models are subsequently used for the prediction of drug targets.

B Uptake and secretion pattern of model metabolites. All metabolite uptakes and secretions that were mapped during model generation are shown.

  • Metabolite uptakes are depicted on the left, and
  • secreted metabolites are shown on the right.
  1. A number of metabolite exchanges mapped to the model were unique to one cell line.
  2. Differences between cell lines were used to set quantitative constraints for the sampling analysis.

C Statistics about the cell line-specific network generation.

D Quantitative constraints.

For the sampling analysis, an additional set of constraints was imposed on the cell line specific models,

  • emphasizing the differences in metabolite uptake and secretion between cell lines.

Higher uptake of a metabolite was allowed

  • in the model of the cell line that consumed more of the metabolite in vitro, whereas
  • the supply was restricted for the model with lower in vitro uptake.

This was done by establishing the same ratio between the models bounds as detected in vitro.

X denotes the factor (slope ratio) that distinguishes the bounds, and

  • which was individual for each metabolite.

(a) The uptake of a metabolite could be x times higher in CCRF-CEM cells,

(b) the metabolite uptake could be x times higher in Molt-4,

(c) metabolite secretion could be x times higher in CCRF-CEM, or

(d) metabolite secretion could be x times higher in Molt-4 cells.LOD limit of detection.

The consequence of the adjustment was, in case of uptake, that one model was constrained to a lower metabolite uptake (A, B), and the difference depended on the ratio detected in vitro. In case of secretion, one model

  • had to secrete more of the metabolite, and again
  • the difference depended on the experimental difference detected between the cell lines

2 Results

We set up a pipeline that could be used to infer intracellular metabolic states

  • from semi-quantitative data regarding metabolites exchanged between cells and their environment.

Our pipeline combined the following four steps:

  1. data acquisition,
  2. data analysis,
  3. metabolic modeling and
  4. experimental validation of the model predictions (Fig. 1A).

We demonstrated the pipeline and the predictive potential to predict metabolic alternations in diseases such as cancer based on

^two lymphoblastic leukemia cell lines.

The resulting Molt-4 and CCRF-CEM condition-specific cell line models could explain

^  metabolite uptake and secretion
^  by predicting the distinct utilization of central metabolic pathways by the two cell lines.
^  the CCRF-CEM model resembled more a glycolytic, commonly referred to as ‘Warburg’ phenotype,
^  our model predicted a more respiratory phenotype for the Molt-4 model.

We found these predictions to be in agreement with measured gene expression differences

  • at key regulatory steps in the central metabolic pathways, and they were also
  • consistent with additional experimental data regarding the energy and redox states of the cells.

After a brief discussion of the data generation and analysis steps, the results derived from model generation and analysis will be described in detail.

2.1 Pipeline for generation of condition-specific metabolic cell line models

integration of exometabolomic (EM) data

integration of exometabolomic (EM) data

2.1.1 Generation of experimental data

We monitored the growth and viability of lymphoblastic leukemia cell lines in serum-free medium (File S2, Fig. S1). Multiple omics data sets were derived from these cells.Extracellular metabolomics (exo-metabolomic) data,

integration of exometabolomic (EM) data

integration of exometabolomic (EM) data

^  comprising measurements of the metabolites in the spent medium of the cell cultures (Paglia et al. 2012a),
^ were collected along with transcriptomic data, and these data sets were used to construct the models.

2.1.4 Condition-specific models for CCRF-CEM and Molt-4 cells

To determine whether we had obtained two distinct models, we evaluated the reactions, metabolites, and genes of the two models. Both the Molt-4 and CCRF-CEM models contained approximately half of the reactions and metabolites present in the global model (Fig. 1C). They were very similar to each other in terms of their reactions, metabolites, and genes (File S1, Table S5A–C).

(1) The Molt-4 model contained seven reactions that were not present in the CCRF-CEM model (Co-A biosynthesis pathway and exchange reactions).
(2) The CCRF-CEM contained 31 unique reactions (arginine and proline metabolism, vitamin B6 metabolism, fatty acid activation, transport, and exchange reactions).
(3) There were 2 and 15 unique metabolites in the Molt-4 and CCRF-CEM models, respectively (File S1, Table S5B).
(4) Approximately three quarters of the global model genes remained in the condition-specific cell line models (Fig. 1C).
(5) The Molt-4 model contained 15 unique genes, and the CCRF-CEM model had 4 unique genes (File S1, Table S5C).
(6) Both models lacked NADH dehydrogenase (complex I of the electron transport chain—ETC), which was determined by the absence of expression of a mandatory subunit (NDUFB3, Entrez gene ID 4709).

Rather, the ETC was fueled by FADH2 originating from succinate dehydrogenase and from fatty acid oxidation, which through flavoprotein electron transfer

FADH2

FADH2

  • could contribute to the same ubiquinone pool as complex I and complex II (succinate dehydrogenase).

Despite their different in vitro growth rates (which differed by 11 %, see File S2, Fig. S1) and
^^^ differences in exo-metabolomic data (Fig. 1B) and transcriptomic data,
^^^ the internal networks were largely conserved in the two condition-specific cell line models.

2.1.5 Condition-specific cell line models predict distinct metabolic strategies

Despite the overall similarity of the metabolic models, differences in their cellular uptake and secretion patterns suggested distinct metabolic states in the two cell lines (Fig. 1B and see “Materials and methods” section for more detail). To interrogate the metabolic differences, we sampled the solution space of each model using an Artificial Centering Hit-and-Run (ACHR) sampler (Thiele et al. 2005). For this analysis, additional constraints were applied, emphasizing the quantitative differences in commonly uptaken and secreted metabolites. The maximum possible uptake and maximum possible secretion flux rates were reduced
^^^ according to the measured relative differences between the cell lines (Fig. 1D, see “Materials and methods” section).

We plotted the number of sample points containing a particular flux rate for each reaction. The resulting binned histograms can be understood as representing the probability that a particular reaction can have a certain flux value.

A comparison of the sample points obtained for the Molt-4 and CCRF-CEM models revealed

  • a considerable shift in the distributions, suggesting a higher utilization of glycolysis by the CCRF-CEM model
    (File S2, Fig. S2).

This result was further supported by differences in medians calculated from sampling points (File S1, Table S6).
The shift persisted throughout all reactions of the pathway and was induced by the higher glucose uptake (34 %) from the extracellular medium in CCRF-CEM cells.

The sampling median for glucose uptake was 34 % higher in the CCRF-CEM model than in Molt-4 model (File S2, Fig. S2).

The usage of the TCA cycle was also distinct in the two condition-specific cell-line models (Fig. 2). Interestingly,
the models used succinate dehydrogenase differently (Figs. 2, 3).

TCA_reactions

TCA_reactions

The Molt-4 model utilized an associated reaction to generate FADH2, whereas

  • in the CCRF-CEM model, the histogram was shifted in the opposite direction,
  • toward the generation of succinate.

Additionally, there was a higher efflux of citrate toward amino acid and lipid metabolism in the CCRF-CEM model (Fig. 2). There was higher flux through anaplerotic and cataplerotic reactions in the CCRF-CEM model than in the Molt-4 model (Fig. 2); these reactions include

(1) the efflux of citrate through ATP-citrate lyase,
(2) uptake of glutamine,
(3) generation of glutamate from glutamine,
(4) transamination of pyruvate and glutamate to alanine and to 2-oxoglutarate,
(5) secretion of nitrogen, and
(6) secretion of alanine.

energetics-of-cellular-respiration

energetics-of-cellular-respiration

The Molt-4 model showed higher utilization of oxidative phosphorylation (Fig. 3), again supported by
elevated median flux through ATP synthase (36 %) and other enzymes, which contributed to higher oxidative metabolism. The sampling analysis therefore revealed different usage of central metabolic pathways by the condition-specific models.

Fig. 2

Differences in the use of  the TCA cycle by the CCRF-CEM model (red) and the Molt-4 model (blue).

Differences in the use of the TCA cycle by the CCRF-CEM model (red) and the Molt-4 model (blue).

Differences in the use of the TCA cycle by the CCRF-CEM model (red) and the Molt-4 model (blue).

The table provides the median values of the sampling results. Negative values in histograms and in the table describe reversible reactions with flux in the reverse direction. There are multiple reversible reactions for the transformation of isocitrate and α-ketoglutarate, malate and fumarate, and succinyl-CoA and succinate. These reactions are unbounded, and therefore histograms are not shown. The details of participating cofactors have been removed.

Figure 3.

Molt-4 has higher median flux through ETC reactions II–IV 11306_2014_721_Fig3_HTML

Molt-4 has higher median flux through ETC reactions II–IV 11306_2014_721_Fig3_HTML

Atp ATP, cit citrate, adp ADP, pi phosphate, oaa oxaloacetate, accoa acetyl-CoA, coa coenzyme-A, icit isocitrate, αkg α-ketoglutarate, succ-coa succinyl-CoA, succ succinate, fumfumarate, mal malate, oxa oxaloacetate,
pyr pyruvate, lac lactate, ala alanine, gln glutamine, ETC electron transport chain

Ingenuity network analysis showing up (red) and downregulation (green) of miRNAs involved in PC and their target genes

Ingenuity network analysis showing up (red) and downregulation (green) of miRNAs involved in PC and their target genes

metabolic pathways 1476-4598-10-70-1

metabolic pathways 1476-4598-10-70-1

Metabolic Systems Research Team fig2

Metabolic Systems Research Team fig2

Metabolic control analysis of respiration in human cancer tissue. fphys-04-00151-g001

Metabolic control analysis of respiration in human cancer tissue. fphys-04-00151-g001

Metabolome Informatics Research fig1

Metabolome Informatics Research fig1

Modelling of Central Metabolism network3

Modelling of Central Metabolism network3

N. gaditana metabolic pathway map ncomms1688-f4

N. gaditana metabolic pathway map ncomms1688-f4

protein changes in biological mechanisms

protein changes in biological mechanisms

Read Full Post »

Extracellular evaluation of intracellular flux in yeast cells

Larry H. Bernstein, MD, FCAP, Reviewer and Curator

Leaders in Pharmaceutical Intelligence

This is the fourth article in a series on metabolomics, which is a major development in -omics, integrating transcriptomics, proteomics,  genomics, metabolic pathways analysis, metabolic and genomic regulatory control using computational mapping.  In the previous two part presentation, flux analysis was not a topic for evaluation, but here it is the major focus.  It is a study of yeast cells, and bears some relationship to the comparison of glycemia, oxidative phosphorylation, TCA cycle, and ETC in leukemia cell lines.  In the previous study – system flux was beyond the scope of analysis, and explicitly stated.  The inferences made in comparing the two lymphocytic leukemia cells was of intracellular metabolism from extracellular measurements.  The study of yeast cells is aimed at looking at cellular effluxes, which is also an important method for studying pharmacological effects and drug resistance.

Metabolomic series

1.  Metabolomics, Metabonomics and Functional Nutrition: the next step in nutritional metabolism and biotherapeutics

http://pharmaceuticalintelligence.com/2014/08/22/metabolomics-metabonomics-and-functional-nutrition-the-next-step-in-nutritional-metabolism-and-biotherapeutics/

2.  Metabolomic analysis of two leukemia cell lines. I

http://pharmaceuticalintelligence.com/2014/08/23/metabolomic-analysis-of-two-leukemia-cell-lines-_i/

3.  Metabolomic analysis of two leukemia cell lines. II.

 http://pharmaceuticalintelligence.com/2014/08/24/metabolomic-analysis-of-two-leukemia-cell-lines-ii/

4.  Extracellular evaluation of intracellular flux in yeast cells

Q1. What is efflux?

Q2. What measurements were excluded from the previous study that would not allow inference about fluxes?

Q3. Would this study bear any relationship to the Pasteur effect?

Q4 What is a genome scale network reconstruction?

Q5 What type of information is required for a network prediction model?

Q6. Is there a difference between the metabolites profiles for yeast grown under aerobic and anaerobuc conditions – under the constrainsts?

Q7.  If there is a difference in the S metabolism, would there be an effect on ATP production?

 

 

Connecting extracellular metabolomic measurements to intracellular flux
states in yeast

Monica L Mo1Bernhard Ø Palsson1 and Markus J Herrgård12*

Author Affiliations

1 Department of Bioengineering, University of California, San Diego, La Jolla, CA 92093, USA

2 Current address: Synthetic Genomics, Inc, 11149 N Torrey Pines Rd, La Jolla, CA 92037, USA

For all author emails, please log on.

BMC Systems Biology 2009, 3:37  doi:10.1186/1752-0509-3-37

 

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1752-0509/3/37

 

Received: 15 December 2008
Accepted: 25 March 2009
Published: 25 March 2009

© 2009 Mo et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Metabolomics has emerged as a powerful tool in the

  • quantitative identification of physiological and disease-induced biological states.

Extracellular metabolome or metabolic profiling data, in particular,

  • can provide an insightful view of intracellular physiological states in a noninvasive manner.

Results

We used an updated genome-scale

  • metabolic network model of Saccharomyces cerevisiae, iMM904, to investigate
  1. how changes in the extracellular metabolome can be used
  2. to study systemic changes in intracellular metabolic states.

The iMM904 metabolic network was reconstructed based on

  • an existing genome-scale network, iND750,
  • and includes 904 genes and 1,412 reactions.

The network model was first validated by

  • comparing 2,888 in silico single-gene deletion strain growth phenotype predictions
  • to published experimental data.

Extracellular metabolome data measured

  • of ammonium assimilation pathways 
  • in response to environmental and genetic perturbations

was then integrated with the iMM904 network

  • in the form of relative overflow secretion constraints and
  • a flux sampling approach was used to characterize candidate flux distributions allowed by these constraints.

Predicted intracellular flux changes were

  • consistent with published measurements
  • on intracellular metabolite levels and fluxes.

Patterns of predicted intracellular flux changes

  • could also be used to correctly identify the regions of
  • the metabolic network that were perturbed.

Conclusion

Our results indicate that

  • integrating quantitative extracellular metabolomic profiles
  • in a constraint-based framework
  • enables inferring changes in intracellular metabolic flux states.

Similar methods could potentially be applied

  • towards analyzing biofluid metabolome variations
  • related to human physiological and disease states.

Background

“Omics” technologies are rapidly generating high amounts of data

  • at varying levels of biological detail.

In addition, there is a rapidly growing literature and

  • accompanying databases that compile this information.

This has provided the basis for the assembly of

  • genome-scale metabolic networks for various microbial and eukaryotic organisms [111].

These network reconstructions serve

  • as manually curated knowledge bases of
  • biological information as well as
  • mathematical representations of biochemical components and
  • interactions specific to each organism.

genome-scale network reconstruction is

  • structured collection of genes, proteins, biochemical reactions, and metabolites
  • determined to exist and operate within a particular organism.

This network can be converted into a predictive model

  • that enables in silico simulations of allowable network states based on
  • governing physico-chemical and genetic constraints [12,13].

A wide range of constraint-based methods have been developed and applied

  • to analyze network metabolic capabilities under
  • different environmental and genetic conditions [13].

These methods have been extensively used to

  • study genome-scale metabolic networks and have successfully predicted, for example,
  1. optimal metabolic states,
  2. gene deletion lethality, and
  3. adaptive evolutionary endpoints [1416].

Most of these applications utilize

  • optimization-based methods such as flux balance analysis (FBA)
  • to explore the metabolic flux space.

However, the behavior of genome-scale metabolic networks can also be studied

  • using unbiased approaches such as
  • uniform random sampling of steady-state flux distributions [17].

Instead of identifying a single optimal flux distribution based on

  • a given optimization criterion (e.g. biomass production),

these methods allow statistical analysis of

  • a large range of possible alternative flux solutions determined by
  • constraints imposed on the network.

Sampling methods have been previously used to study

  1. global organization of E. coli metabolism [18] as well as
  2. to identify candidate disease states in the cardiomyocyte mitochondria [19].

Network reconstructions provide a structured framework

  • to systematically integrate and analyze disparate datasets
  • including transcriptomic, proteomic, metabolomic, and fluxomic data.

Metabolomic data is one of the more relevant data types for this type of analysis as

  1. network reconstructions define the biochemical links between metabolites, and
  2. recent advancements in analytical technologies have allowed increasingly comprehensive
  • intracellular and extracellular metabolite level measurements [20,21].

The metabolome is

  1. the set of metabolites present under a given physiological condition
  2. at a particular time and is the culminating phenotype resulting from
  • various “upstream” control mechanisms of metabolic processes.

Of particular interest to this present study are

  • the quantitative profiles of metabolites that are secreted into the extracellular environment
  • by cells under different conditions.

Recent advances in profiling the extracellular metabolome (EM) have allowed

  • obtaining insightful biological information on cellular metabolism
  • without disrupting the cell itself.

This information can be obtained through various

  • analytical detection,
  • identification, and
  • quantization techniques

for a variety of systems ranging from

  • unicellular model organisms to human biofluids [2023].

Metabolite secretion by a cell reflects its internal metabolic state, and

  • its composition varies in response to
  • genetic or experimental perturbations
  • due to changes in intracellular pathway activities
  • involved in the production and utilization of extracellular metabolites [21].

Variations in metabolic fluxes can be reflected in EM changes which can

  • provide insight into the intracellular pathway activities related to metabolite secretion.

The extracellular metabolomic approach has already shown promise

  • in a variety of applications, including
  1. capturing detailed metabolite biomarker variations related to disease and
  2. drug-induced states and
  3. characterizing gene functions in yeast [2427].

However, interpreting changes in the extracellular metabolome can be challenging

  • due to the indirect relationship between the proximal cause of the change
    (e.g. a mutation)
  • and metabolite secretion.

Since metabolic networks describe

  • mechanistic,
  • biochemical links between metabolites,

integrating such data can allow a systematic approach

  • to identifying altered pathways linked to
  • quantitative changes in secretion profiles.

Measured secretion rates of major byproduct metabolites

  • can be applied as additional exchange flux constraints
  • that define observed metabolic behavior.

For example, a recent study integrating small-scale EM data

  • with a genome-scale yeast model
  • correctly predicted oxygen consumption and ethanol production capacities
  • in mutant strains with respiratory deficiencies [28].

The respiratory deficient mutant study

  • used high accuracy measurements for a small number of
  • major byproduct secretion rates
  • together with an optimization-based method well suited for such data.

Here, we expand the application range of the model-based method used in [28]

  • to extracellular metabolome profiles,
  • which represent a temporal snapshot of the relative abundance
  • for a larger number of secreted metabolites.

Our approach is complementary to

  • statistical (i.e. “top-down”) approaches to metabolome analysis [29]
  • and can potentially be used in applications such as biofluid-based diagnostics or
  • large-scale characterization of mutants strains using metabolite profiles.

This study implements a constraint-based sampling approach on

  • an updated genome-scale network of yeast metabolism
  • to systematically determine how EM level variations

are linked to global changes in intracellular metabolic flux states.

By using a sampling-based network approach and statistical methods (Figure 1),

  • EM changes were linked to systemic intracellular flux perturbations
    in an unbiased manner
  • without relying on defining single optimal flux distributions
  • used in the previously mentioned study [28].

The inferred perturbations in intracellular reaction fluxes were further analyzed

  • using reporter metabolite and subsystem (i.e., metabolic pathway) approaches [30]
  • in order to identify dominant metabolic features that are collectively perturbed (Figure 2).

The sampling-based approach also has the additional benefit of

  • being less sensitive to inaccuracies in metabolite secretion profiles than
  • optimization-based methods and can effectively be used – in biofluid metabolome analysis.
integration of exometabolomic (EM) data

integration of exometabolomic (EM) data

Figure 1. Schematic illustrating the integration of exometabolomic (EM) data with the constraint-based framework.

(A) Cells are subjected to genetic and/or environmental perturbations to secrete metabolite patterns unique to that condition.
(B) EM is detected, identified, and quantified.
(C) EM data is integrated as required secretion flux constraints to define allowable solution space.
(D) Random sampling of solution space yields the range of feasible flux distributions for intracellular reactions.
(E) Sampled fluxes were compared to sampled fluxes of another condition to determine

  • which metabolic regions were altered between the two conditions (see Figure 2).

(F) Significantly altered metabolic regions were identified.

http://www.biomedcentral.com/content/figures/1752-0509-3-37-1.jpg

 

sampling and scoring analysis to determine intracellular flux changes

sampling and scoring analysis to determine intracellular flux changes

Figure 2. Schematic of sampling and scoring analysis to determine intracellular flux changes.

(A) Reaction fluxes are sampled for two conditions.
(B & C) Sample of flux differences is calculated by selecting random flux values from each condition

  • to obtain a distribution of flux differences for each reaction.

(D) Standardized reaction Z-scores are determined, which represent

  • how far the sampled flux differences deviates from a zero flux change.

Reaction scores can be used in

  1. visualizing perturbation subnetworks and
  2. analyzing reporter metabolites and subsystems.

http://www.biomedcentral.com/content/figures/1752-0509-3-37-2.jpg

This study was divided into two parts and describes:

(i) the reconstruction and validation of an expanded S. cerevisiae metabolic network, iMM904; and
(ii) the systematic inference of intracellular metabolic states from

  • two yeast EM data sets using a constraint-based sampling approach.

The first EM data set compares wild type yeast to the gdh1/GDH2 (glutamate dehydrogenase) strain [31],

  • which indicated good agreement between predicted metabolic changes
  • of intracellular metabolite levels and fluxes [31,32].

The second EM data set focused on secreted amino acid measurements

  • from a separate study of yeast cultured in different
    ammonium and potassium concentrations [33].

We analyzed the EM data to gain further insight into

  • perturbed ammonium assimilation processes as well as
  1. metabolic states relating potassium limitation and
  2. ammonium excess conditions to one another.

The model-based analysis of both

  • separately published extracellular metabolome datasets
  • suggests a relationship between
  1. glutamate,
  2. threonine and
  3. folate metabolism,
  • which are collectively perturbed when
    ammonium assimilation processes are broadly disrupted
  1. either by environmental (excess ammonia) or
  2. genetic (gene deletion/overexpression) perturbations.

The methods herein present an approach to

  • interpreting extracellular metabolome data and
  • associating these measured secreted metabolite variations
  • to changes in intracellular metabolic network states.

Additional file 1. iMM904 network content.

The data provided represent the content description of the iMM904 metabolic network and
detailed information on the expanded content.

Format: XLS Size: 2.7MB Download file

This file can be viewed with: Microsoft Excel Viewer

Additional file 2. iMM904 model files.

The data provided are the model text files of the iMM904 metabolic network
that is compatible with the available COBRA Toolbox [13]. The model structure
can be loaded into Matlab using the ‘SimPhenyPlus’ format with GPR and compound information.

Format: ZIP Size: 163KB Download file

Conversion of the network to a predictive model

The network reconstruction was converted to a constraint-based model using established procedures [13].

Network reactions and metabolites were assembled into a stoichiometric matrix 

  • containing the stoichiometric coefficients of the reactions in the network.

The steady-state solution space containing possible flux distributions

  • is determined by calculating the null space of S= 0,

where is the reaction flux vector.

Minimal media conditions were set through constraints on exchange fluxes

  • corresponding to the experimental measured substrate uptake rates.

All the model-based calculations were done using the Matlab COBRA Toolbox [13]

  • utilizing the glpk or Tomlab/CPLEX (Tomopt, Inc.) optimization solvers.

Chemostat growth simulations

The iMM904 model was initially validated by

  1. simulating wild type yeast growth in aerobic and anaerobic
    carbon-limited chemostat conditions
  2. and comparing the simulation results to published experimental data

on substrate uptake and byproduct secretion in these conditions [34].

The study was performed following the approach taken to validate the iFF708 model in a previous study [35].

The predicted glucose uptake rates were determined

  1. by setting the in silico growth rate to the measured dilution rate,
    – equivalent under continuous culture growth,
  2. and minimizing the glucose uptake rate.

The accuracy of in silico predictions of

  • substrate uptake and byproduct secretion by the iMM904 model
  • was similar to the accuracy obtained using the iFF708 model
  • and results are shown in Figure S1 [see Additional file 3].

Additional file 3. Supplemental figures. 

The file provides the supplemental figures and descriptions of S1, S2, S3, and S4.

Format: PDF Size: 513KB Download file

This file can be viewed with: Adobe Acrobat Reader

Genome-scale gene deletion phenotype predictions

The iMM904 network was further validated by

  • performing genome-scale gene lethality computations
  • following established procedures to determine growth phenotypes
  1. under minimal medium conditions and
  2. compared to published data.

A modified version of the biomass function used in previous iND750 studies

  1. was set as the objective to be maximized and
  2. gene deletions were simulated by

setting the flux through the corresponding reaction(s) to zero.

The biomass function was based on the experimentally measured

  1. composition of major cellular constituents
  2. during exponential growth of yeast cells and
  3. was reformulated to include trace amounts of
  4. additional cofactors and metabolites
  5. with the assumed fractional contribution of 10-.

These additional biomass compounds were included

according to the biomass formulation used in the iLL672 study

  • to improve lethality predictions through
  • the inclusion of additional essential biomass components [3].

The model was constrained by limiting

  1. the carbon source uptake to 10 mmol/h/gDW
  2. and oxygen uptake to 2 mmol/h/gDW.

Ammonia, phosphate, and sulfate were assumed to be non-limiting.

The experimental phenotyping data was obtained

  • using strains that were auxotrophic for
  1. methionine,
  2. leucine,
  3. histidine, and
  4. uracil.

These auxotrophies were simulated

  1. by deleting the appropriate genes from the model and
  2. supplementing the in silico strain with the appropriate supplements
  3. at non-limiting, but low levels.

Furthermore, trace amounts of essential nutrients that are present

  • in the experimental minimal media formulation
  1. 4-aminobenzoate,
  2. biotin,
  3. inositol,
  4. nicotinate,
  5. panthothenate,
  6. thiamin)
  • were supplied in the simulations [3].

Three distinct methods to simulate the outcome of gene deletions were utilized:

  1. Flux-balance analysis (FBA) [36-38],
  2. Minimization of Metabolic Adjustment (MoMA) [39], and
  3. a linear version of MoMA (linearMoMA).

In the linearMoMA method, minimization of the quadratic objective function
of the original MoMA algorithm

  • was replaced by minimization of the corresponding 1-norm objective function
    (i.e. sum of the absolute values of the differences of wild type FBA solution
    and the knockout strain flux solution).

The computed results were then compared to growth phenotype data
(viable/lethal) from a previously published experimental gene deletion study [3].

The comparison between experimental and in silico deletion phenotypes involved

  • choosing a threshold for the predicted relative growth rate of
  • a deletion strain that is considered to be viable.

We used standard ROC curve analysis

  • to assess the accuracy of different prediction methods and models
  • across the full range of the viability threshold parameter,
    results shown in Figure S2 [see Additional file 3].

The ROC curve plots the true viable rate against the false viable rate

  • allowing comparison of different models and methods
  • without requiring arbitrarily choosing this parameter a priori [40].

The optimal prediction performance corresponds to

  • the point closest to the top left corner of the ROC plot
    (i.e. 100% true viable rate, 0% false viable rate).

Table 1

Table 1 Comparison of iMM904 and iLL672 gene deletion predictions and experimental data under minimal media conditions
Media Model Method True viable False viable False lethal True lethal True viable % False viable % MCC
Glucose iMM904 full FBA 647 10 32 33 95.29 23.26 0.6
iMM904 full linMOMA 644 10 35 33 94.85 23.26 0.58
iMM904 full MOMA 644 10 35 33 94.85 23.26 0.58
iMM904 red FBA 440 9 28 33 94.02 21.43 0.61
iMM904 red linMOMA 437 9 31 33 93.38 21.43 0.6
iMM904 red MOMA 437 9 31 33 93.38 21.43 0.6
iLL672 full MOMA 433 9 35 33 92.52 21.43 0.57
Galactose iMM904 full FBA 595 32 36 59 94.29 35.16 0.58
iMM904 full linMOMA 595 32 36 59 94.29 35.16 0.58
iMM904 full MOMA 595 32 36 59 94.29 35.16 0.58
iMM904 red FBA 409 12 33 56 92.53 17.65 0.67
iMM904 red linMOMA 409 12 33 56 92.53 17.65 0.67
iMM904 red MOMA 409 12 33 56 92.53 17.65 0.67
iLL672 full MOMA 411 19 31 49 92.99 27.94 0.61
Glycerol iMM904 full FBA 596 43 36 47 94.3 47.78 0.48
iMM904 full linMOMA 595 44 37 46 94.15 48.89 0.47
iMM904 full MOMA 598 44 34 46 94.62 48.89 0.48
iMM904 red FBA 410 20 34 46 92.34 30.3 0.57
iMM904 red linMOMA 409 21 35 45 92.12 31.82 0.56
iMM904 red MOMA 412 21 32 45 92.79 31.82 0.57
iLL672 full MOMA 406 20 38 46 91.44 30.3 0.55
Ethanol iMM904 full FBA 593 45 29 55 95.34 45 0.54
iMM904 full linMOMA 592 45 30 55 95.18 45 0.54
iMM904 full MOMA 592 44 30 56 95.18 44 0.55
iMM904 red FBA 408 21 27 54 93.79 28 0.64
iMM904 red linMOMA 407 21 28 54 93.56 28 0.63
iMM904 red MOMA 407 20 28 55 93.56 26.67 0.64
iLL672 full MOMA 401 13 34 62 92.18 17.33 0.68
MCC, Matthews correlation coefficient (see Methods). Note that the iLL672 predictions were obtained directly from [3] and thus the viability threshold was not optimized using the maximum MCC approach.
Mo et al. BMC Systems Biology 2009 3:37  http://dx.doi.org:/10.1186/1752-0509-3-37

 

The values reported in Table 1 correspond to selecting

  • the optimal viability threshold based on this criterion.

We summarized the overall prediction accuracy of a model and method

  • using the Matthews Correlation Coefficient (MCC) [40].

The MCC ranges from -1 (all predictions incorrect) to +1 (all predictions correct) and

  • is suitable for summarizing overall prediction performance

in our case where there are substantially more viable than lethal gene deletions.

ROC plots were produced in Matlab (Mathworks, Inc.).

 

Table 1. Comparison of iMM904 and iLL672

  • gene deletion predictions and
  • experimental data

Inferring perturbed metabolic regions based on EM profiles

The method implemented in this study is shown schematically in Figures 1 and 2

Constraining the iMM904 network 

Relative levels of quantitative EM data were incorporated into the constraint-based framework

  • as overflow secretion exchange fluxes to simulate the required low-level production of
  • experimentally observed excreted metabolites.

The primary objective of this study is to associate

  • relative metabolite levels that are generally measured for metabonomic or biofluid analyses
  • to the quantitative ranges of intracellular reaction fluxes required to produce them.

However, without detailed kinetic information or dynamic metabolite measurements available,

  • we approximated EM datasets of relative quantitative metabolite levels
  • to be proportional to the rate in which they are secreted and detected
  • (at a steady state) – into the extracellular media.

This approach is analogous to approximating uptake rates based

  • on metabolite concentrations from a previous study performing sampling analysis
  • on a cardiomyocyte mitochondrial network
  • to identify differential flux distribution ranges

for various environmental (i.e. substrate uptake) conditions [19].

The raw data was normalized by the raw maximum value of the dataset
(thus the maximum secretion flux was 1 mmol/hr/gDW) with

  • an assumed error of 10%
  • to set the lower and upper bounds and thus
  • inherently accounting for sampling calculation sensitivity.

The gdh1/GDH2 strains were flask cultured under minimal glucose media conditions; thus,

  • glucose and oxygen uptake rates were set at 15 and 2 mmol/hr/gDW, respectively,
  • for the gdh1/GDH2 strain study.

In the anaerobic case the oxygen uptake rate was set to zero, and

  • sterols and fatty acids were provided as in silico supplements as described in [35].

For the potassium limitation/ammonium toxicity study

  • the growth rate was set at 0.17 1/h, and
  • the glucose uptake rate was minimized
  • to mimic experimental chemostat cultivation conditions.

These input constraints were constant for each perturbation and comparative wild-type condition

  • such that the calculated solution spaces between the conditions
  • differed based only on variations in the output secretion constraints.

FBA optimization of EM-constrained networks

A modified FBA method with minimization of the 1-norm objective function

  • between two optimal flux distributions was used
  • to determine optimal intracellular fluxes
  • based on the EM-constrained metabolic models.

This method determines two optimal flux distributions simultaneously

  • for two differently constrained models (e.g. wild type vs. mutant) –
  • these flux distributions maximize biomass production in each case and
  • the 1-norm distance between the distributions is as small as possible
  • given the two sets of constraints.

This approach avoids problems with

  • alternative optimal solutions when comparing two FBA-computed flux distributions
  • by assuming minimal rerouting of flux distibution between a perturbed network and its reference network.

Reaction flux changes from the FBA optimization results were determined

  • by computing the relative percentage fold change for each reaction
  • between the mutant and wild-type flux distributions.

Random sampling of the steady-state solution space

We utilized artificial centering hit-and-run (ACHR) Monte Carlo sampling [19,41]

  • to uniformly sample the metabolic flux solution space
  • defined by the constraints described above.

Reactions, and their participating metabolites, found to participate in intracellular loops [42]

  • were discarded from further analysis as these reactions can have arbitrary flux values.

The following sections describe the approaches used for the analysis of the different datasets.

Sampling approach used in the gdh1/GDH2 study

Due to the overall shape of the metabolic flux solution space,

  • most of the sampled flux distributions resided close to the minimally allowed growth rate
    (i.e. biomass production) and
  • corresponded to various futile cycles that utilized substrates but
  • did not produce significant biomass.

In order to study more physiologically relevant portions of the flux space

  • we restricted the sampling to the part of the solution space
  • where the growth rate was at least 50% of the maximum growth rate
  • for the condition as determined by FBA.

This assumes that cellular growth remains an important overall objective by the yeast cells

  • even in batch cultivation conditions, but
  • that the intracellular flux distributions
  • may not correspond to maximum biomass production [43].

To test the sensitivity of the results to the minimum growth rate threshold,

  • separate Monte Carlo samples were created for each minimum threshold
  • ranging from 50% to 100% at 5% increments.

We also tested the sensitivity of the results

  • to the relative magnitude of the extracellular metabolite secretion rates
  • by performing the sampling at three different relative levels

(0 corresponding to no extracellular metabolite secretion, maximum rate of 0.5 mmol/hr/gDW,
and maximum rate of 1.0 mmol/hr/gDW).

For each minimum growth rate threshold and extracellular metabolite secretion rate,

  • the ACHR sampler was run for 5 million steps and
  • a flux distribution was stored every 5000 steps.

The sensitivity analysis results are presented in Figures S3 and S4 [see Additional File 3], and

  • the results indicate that the reaction Z-scores (see below) are not significantly affected by
  1. either the portion of the solution space sampled or
  2. the exact scaling of secretion rates.

The final overall sample used was created by combining the samples for all minimum growth rate thresholds

  • for the highest extracellular metabolite secretion rate (maximum 1 mmol/hr/gDW).

This approach allowed biasing the sampling towards

  • physiologically relevant parts of the solution space
  • without imposing the requirement of strictly maximizing a predetermined objective function.

The samples obtained with no EM data were used as control samples

  • to filter reporter metabolites/subsystems whose scores were significantly high
  • due to only random differences between sampling runs.

Sampling approach used in the potassium limitation/ammonium toxicity study

Since the experimental data used in this study was generated in chemostat conditions, and

  • previous studies have indicated that chemostat flux patterns predicted by FBA are
  • close to the experimentally measured ones [43],
  • we assumed that sampling of the optimal solution space was appropriate for this study.

In order to sample a physiologically reasonable range of flux distributions,

  • samples for four different oxygen uptake rates
    (1, 2, 3, and 4 mmol/hr/gDW with 5 million steps each)
  • were combined in the final analysis.

Standardized scoring of flux differences between perturbation and control conditions

Z-score based approach was implemented to quantify differences in flux samples between two conditions (Figure 2).
First, two flux vectors were chosen randomly,

  • one from each of the two samples to be compared and
  • the difference between the flux vectors was computed.

This approach was repeated to create a sample of 10,000 (n) flux difference vectors

  • for each pair of conditions considered (e.g. mutant or perturbed environment vs. wild type).

Based on this flux difference sample, the sample mean (μdiff,i) and standard deviation (σdiff,i)

  • between the two conditions was calculated for each reaction i. The reaction Z-score was calculated as:

 

reaction Z-score

reaction Z-score

which describes the sampled mean difference deviation

  • from a population mean change of zero (i.e. no flux difference
    between perturbation and wild type).

Note that this approach allows accounting for uncertainty in the

  • flux distributions inferred based on the extracellular metabolite secretion constraints.

This is in contrast to approaches such as FBA or MoMA that would predict

  • a single flux distribution for each condition and thus potentially
  • overestimate differences between conditions.

The reaction Z-scores can then be further used in analysis

  • to identify significantly perturbed regions of the metabolic network
  • based on reporter metabolite [44] or subsystem [30] Z-scores.

These reporter regions indicate, or “report”, dominant perturbation features

  • at the metabolite and pathway levels for a particular condition.

The reporter metabolite Z-score for any metabolite can be derived from the reaction Z-scores

  • of the reactions consuming or producing j (set of reactions denoted as Rj) as:

 

reporter z-score for any metabolite j

reporter z-score for any metabolite j

where Nis the number of reactions in Rand mmet,is calculated as

 

distributional correction for m_met,j SQRT

distributional correction for m_met,j SQRT

To account and correct for background distribution, the metabolite Z-score was normalized

  • by computing μmet,Nj and σmet,,Nj corresponding to the mean mmet and
  • its standard deviation for 1,000 randomly generated reaction sets of size Nj.

Z-scores for subsystems were calculated similarly by considering the set of reactions R

  • that belongs to each subsystem k.

Hence, positive metabolite and subsystem scores indicate a significantly perturbed metabolic region

  • relative to other regions, whereas
  • a negative score indicate regions that are not perturbed
  • more significantly than what is expected by random chance.

Perturbation subnetworks of reactions and connecting metabolites were visualized using Cytoscape [45].

Results and discussion

  1. Reconstruction and validation of iMM904 network iMM904 network content 

A previously reconstructed S. cerevisiae network, iND750,

  • was used as the basis for the construction of the expanded iMM904 network.
  • Prior to its presentation here, the
    iMM904 network content was the basis for a consensus jamboree network that was recently published
  • but has not yet been adapted for FBA calculations [46].

The majority of iND750 content was carried over and

  • further expanded on to construct iMM904, which accounts for
  1. 904 genes,
  2. 1,228 individual metabolites, and
  3. 1,412 reactions of which
  •                       395 are transport reactions.

Both the number of gene-associated reactions and the number of metabolites

  • increased in iMM904 compared with the iND750 network.

Additional genes and reactions included in the network primarily expanded the

  • lipid,
  • transport, and
  • carbohydrate subsystems.

The lipid subsystem includes

  • new genes and
  • reactions involving the degradation of sphingolipids and glycerolipids.

Sterol metabolism was also expanded to include

  • the formation and degradation of steryl esters, the
  •                      storage form of sterols.

The majority of the new transport reactions were added

  • to connect network gaps between intracellular compartments
  • to enable the completion of known physiological functions.

We also added a number of new secretion pathways

  • based on experimentally observed secreted metabolites [31].

A number of gene-protein-reaction (GPR) relationships were modified

  • to include additional gene products that are required to catalyze a reaction.

For example, the protein compounds

  • thioredoxin and
  • ferricytochrome C

were explicitly represented as compounds in iND750 reactions, but

  • the genes encoding these proteins were not associated with their corresponding GPRs.

Other examples include glycogenin and NADPH cytochrome p450 reductases (CPRs),

  1. which are required in the assembly of glycogen and
  2. to sustain catalytic activity in cytochromes p450, respectively.

These additional proteins were included in iMM904 as

  • part of protein complexes to provide a more complete
  • representation of the genes and
  • their corresponding products necessary for a catalytic activity to occur.

Major modifications to existing reactions were in cofactor biosynthesis, namely in

  • quinone,
  • beta-alanine, and
  • riboflavin biosynthetic pathways.

Reactions from previous S. cerevisiae networks associated with

  • quinone,
  • beta-alanine, and
  • riboflavin biosynthetic pathways

were essentially inferred from known reaction mechanisms based on

  • reactions in previous network reconstructions of E. coli [2,47].

These pathways were manually reviewed

  • based on current literature and subsequently replaced by
  • reactions and metabolites specific to yeast.

Additional changes in other subsystems were also made, such as

  1. changes to the compartmental location of a gene and
  2. its corresponding reaction(s),
  3. changes in reaction reversibility and cofactor specificity, and
  4. the elucidation of particular transport mechanisms.

A comprehensive listing of iMM904 network contents as well as

  • a detailed list of changes between iND750 and iMM904 is included
    [see Additional file 1].

Predicting deletion growth phenotypes

The updated genome-scale iMM904 metabolic network was validated

  • by comparing in silico single-gene deletion predictions to
  • in vivo results from a previous study used
  • to analyze another S. cerevisiae metabolic model, iLL672 [3].

This network was constructed based on the iFF708 network [22],

  • which was also the starting point for
  • reconstructing the iND750 network [2].

The experimental data used to validate the iLL672 model consisted of

3,360 single-gene knockout strain phenotypes evaluated

  • under minimal media growth conditions with
  1. glucose,
  2. galactose,
  3. glycerol, and
  4. ethanol

as sole carbon sources. Growth phenotypes for the iMM904 network were predictedusing

  1. FBA [3234],
  2. MoMA [35], and
  3. linear MoMA methods

as described in Methods and subsequently compared to the experimental data (Table 1).

Each deleted gene growth prediction comparison was classified as

  1. true lethal,
  2. true viable,
  3. false lethal, or
  4. false viable.

The growth rate threshold for considering a prediction viable was chosen

  • for each condition and method separately
  • to optimize the tradeoff between true viable and false viable predictions
    (maximum Matthews correlation coefficient, see Methods).

Since iMM904 has 212 more genes than iLL672 with experimental data, we also present results

  • for the subset of iMM904 predictions with genes included in iLL672 (reduced iMM904 set).

When the same gene sets are compared, iMM904 improves gene lethality predictions under

  • glucose,
  • galactose, and
  • glycerol conditions

over iLL672 somewhat, but is less accurate

  • at predicting growth phenotypes under the ethanol condition.

It should be noted that the iLL672 predictions were obtained directly from [3]

  • thus the growth rate threshold was not optimized similarly to iMM904 predictions.

Overall, when viability cutoff is chosen

  • as indicated above for each method separately,
  • the three prediction methods perform similarly
  1. FBA,
  2. MOMA, and
  3. linear MOMA) .

While the full gene complement in iMM904 greatly increased

  • the number of true viable predictions,
  • the full model also made significantly more false viable predictions
  • compared with reduced iMM904 and iLL672 predictions.

However, it is important to note that 143 reactions involved in dead-end biosynthetic pathways were actually

  • removed from iFF708 to build the iLL672 reconstruction [3].

These dead-ends are considered “knowledge gaps” in pathways

  • that have not been fully characterized and, as a result,
  • lead to false viable predictions when determining gene essentiality
  • if the pathway is in fact required for growth under a certain condition [2,26].

As more of these pathways are elucidated and

  • included in the model to
  • fill in existing network gaps,
  • we can expect false viable prediction rates to consequently decrease.

Thus, while a larger network has a temporarily reduced capacity to accurately predict gene deletion phenotypes,

  • it captures a more complete picture of currently known metabolic functions and
  • provides a framework for network expansion as new pathways are elucidated [48].

 

Inferring intracellular perturbation states from metabolic profiles – Aerobic and anaerobic gdh1/GDH2 mutant behavior

The gdh1/GDH2 mutant strain was previously developed [49,50]

  • to lower NADPH consumption in ammonia assimilation, which would
  • favor the NADPH-dependent fermentation of xylose.

In this strain, the NADPH-dependent glutamate dehydrogenase, Gdh1, was

  • deleted and the NADH-dependent form of the enzyme, Gdh2,
  •                     was overexpressed.

The net effect is to allow efficient assimilation of ammonia

  • into glutamate using NADH instead of NADPH as a cofactor.

While growth characteristics remained unaffected,

  • relative quantities of secreted metabolites differed between the wild-type and mutant strain
  • under aerobic and anaerobic conditions.

We analyzed EM data for the gdh1/GDH2 and wild-type strains reported

  • in [31] under aerobic and anaerobic conditions separately using
  • both FBA optimization and
  • sampling-based approaches as described in Methods.

43 measured extracellular and intracellular metabolites from the original dataset [31],

  • primarily of central carbon and amino acid metabolism,
  • were explicitly represented in the iMM904 network [see Additional file 4].

Extracellular metabolite levels were used

  • to formulate secretion constraints and
  • differential intracellular metabolites were used
  • to compare and validate the intracellular flux predictions.

Perturbed reactions from the FBA results were

  • determined by calculating relative flux changes, and
  • reaction Z-scores were calculated from the sampling analysis
  • to quantify flux changes between the mutant and wild-type strains,
  • with Z reaction > 1.96 corresponding to a two-tailed p-value < 0.05 and
  • considered to be significantly perturbed [see Additional file 4].

Additional file 4. Gdh mutant aerobic and anaerobic analysis results. 

The data provided are the full results for the exometabolomic analysis of aerobic and anerobic gdh1/GDH2 mutant.

Format: XLS Size: 669KB Download file

This file can be viewed with: Microsoft Excel Viewer

To validate the predicted results, reaction flux changes from both FBA and sampling methods were compared to differential intracellular metabolite level data measured from the same study. Intracellular metabolites involved in highly perturbed reactions (i.e. reactants and products) predicted from FBA and sampling analyses were identified and
compared to metabolites that were experimentally identified as significantly changed (< 0.05) between mutant and wild-type. Statistical measures of recall, accuracy, and
precision were calculated and represent the predictive sensitivity, exactness, and reproducibility respectively. From the sampling analysis, a considerably larger number of
significantly perturbed reactions are predicted in the anaerobic case (505 reactions, or 70.7% of active reactions) than in aerobic (394 reactions, or 49.8% of active reactions). The top percentile of FBA flux changes equivalent to the percentage of significantly perturbed sampling reactions were compared to the intracellular data. Results from both analyses are summarized in Table 2. Sampling predictions were considerably higher in recall than FBA predictions for both conditions, with respective ranges of 0.83–1
compared to 0.48–0.96. Accuracy was also higher in sampling predictions; however, precision was slightly better in the FBA predictions as expected due to the smaller
number of predicted changes. Overall, the sampling predictions of perturbed intracellular metabolites are strongly consistent with the experimental data and significantly
outperforms that of FBA optimization predictions in accurately predicting differential metabolites involved in perturbed intracellular fluxes.

Table 2. Statistical comparison of the differential intracellular metabolite data set (< 0.05) with metabolites involved in perturbed reactions predicted by FBA optimization and sampling analyses for aerobic and anaerobic gdh1/GDH2 mutant.

 

Table 2 Statistical comparison of the differential intracellular metabolite data set (p < 0.05)
with metabolites involved in perturbed reactions predicted by FBA optimization and
sampling analyses for aerobic and anaerobic gdh1/GDH2 mutant.
                           Aerobic                         Anaerobic                             Overall
FBA Sampling FBA Sampling FBA
Recall 0.48 0.83 0.96 1 0.71 0.91
Accuracy 0.55 0.62 0.64 0.64 0.6 0.63
Precision 0.78 0.69 0.64 0.63 0.68 0.66
Overall statistics indicate combined results of both conditions.
Mo et al. BMC Systems Biology 2009 3:37   http://dx.doi.org:/10.1186/1752-0509-3-37


Figure 3.
 Perturbation reaction subnetwork of gdh1/GDH2 mutant under aerobic conditions.

The network illustrates a simplified subset of highly perturbedPerturbation subnetworks can be drawn to visualize predicted significantly perturbed intracellular reactions and illustrate their connection to the observed secreted metabolites in the aerobic and anaerobic gdh1/GDH2 mutants.

Perturbation reaction subnetwork of gdh1.GDH2 mutant under aerobic conditions.

Perturbation reaction subnetwork of gdh1.GDH2 mutant under aerobic conditions.

Figure 3 shows an example of a simplified aerobic perturbation subnetwork consisting primarily of proximal pathways connected directly to a subset of major secreted
metabolites

  • glutamate,
  • proline,
  • D-lactate, and
  • 2-hydroxybuturate.

Figure 4 displays anaerobic reactions with Z-scores of similar magnitude to the perturbed reactions in Figure 3. The same subset of metabolites is also present in the
larger anaerobic perturbation network and indicates that the NADPH/NADH balance perturbation induced by the gdh1/GDH2 manipulation has widespread effects
beyond just altering glutamate metabolism anaerobically.

Interestingly, it is clear that the majority of the secreted metabolite pathways involve connected perturbed reactions that broadly converge on glutamate.

Note that Figures 3 and 4 only show the subnetworks that consisted of two or more connected reactions  for a number of secreted metabolites no contiguous perturbed pathway could be identified by the sampling approach. This indicates that the secreted metabolite pattern alone is not sufficient to determine which specific
production and secretion pathways are used by the cell for these metabolites.

Reactions connected to aerobically-secreted metabolites predicted from the sampling analysis of the gdh1/GDH2 mutant strain.
The major secreted metabolites

  • glutamate,
  • proline,
  • D-lactate, and
  • 2-hydroxybuturate

were also detected in the anaerobic condition. Metabolite abbreviations are found in Additional file 1.

Figure 4.

Perturbation reaction subnetwork of gdh1/GDH2 mutant under anaerobic conditions.

Perturbation reaction subnetwork of gdh1.GDH2 mutant under anaerobic conditions

Perturbation reaction subnetwork of gdh1.GDH2 mutant under anaerobic conditions

Subnetwork illustrates the highly perturbed anaerobic reactions of similar Z-reaction magnitude to the reactions in Figure 3.

A significantly larger number of reactions indicates mutant metabolic effects are more widespread in the anaerobic environment.
The network shows that perturbed pathways converge on glutamate, the main site in which the gdh1/GDH2 modification was introduced, which
suggests that the direct genetic perturbation effects are amplified under this environment. Metabolite abbreviations are found in Additional file 1.

To further highlight metabolic regions that have been systemically affected by the gdh1/GDH2 modification, reporter metabolite and subsystem methods [30] were used to
summarize reaction scores around specific metabolites and in specific metabolic subsystems. The top ten significant scores for metabolites/subsystems associated with more
than three reactions are summarized in Tables 3 (aerobic) and 4 (anaerobic), with Z > 1.64 corresponding to < 0.05 for a one-tailed distribution. Full data for all reactions,
reporter metabolites, and reporter subsystems is included [see Additional file 4].

Table 3. List of the top ten significant reporter metabolite and subsystem scores for the gdh1/GDH2 vs. wild type comparison in aerobic conditions.

Table 3
List of the top ten significant reporter metabolite and subsystem scores for the gdh1/GDH2 vs. wild type comparison in aerobic conditions.
Reporter metabolite Z-score No of reactions*
L-proline [c] 2.71 4
Carbon dioxide [m] 2.51 15
Proton [m] 2.19 51
Glyceraldehyde 3-phosphate [c] 1.93 7
Ubiquinone-6 [m] 1.82 5
Ubiquinol-6 [m] 1.82 5
Ribulose-5-phosphate [c] 1.8 4
Uracil [c] 1.74 4
L-homoserine [c] 1.72 4
Alpha-ketoglutarate [m] 1.71 8
Reporter subsystem Z-score No of reactions
Citric Acid Cycle 4.58 7
Pentose Phosphate Pathway 3.29 12
Glycine and Serine Metabolism 2.69 17
Alanine and Aspartate Metabolism 2.65 6
Oxidative Phosphorylation 1.79 8
Thiamine Metabolism 1.54 8
Arginine and Proline Metabolism 1.44 20
Other Amino Acid Metabolism 1.28 5
Glycolysis/Gluconeogenesis 0.58 14
Anaplerotic reactions 0.19 9
*Number of reactions categorized in a subsystem or found to be neighboring each metabolite
Mo et al. BMC Systems Biology 2009 3:37   http://dx.doi.org:/10.1186/1752-0509-3-37

Table 4. List of top ten significant reporter metabolite and subsystem scores for the gdh1/GDH2 vs. wild type comparison in anaerobic conditions.

 

Table 4
List of top ten significant reporter metabolite and subsystem scores for the gdh1/GDH2 vs. wild type comparison in anaerobic conditions.
Reporter metabolite Z-score No of reactions
Glutamate [c] 4.52 35
Aspartate [c] 3.21 11
Alpha-ketoglutarate [c] 2.66 17
Glycine [c] 2.65 7
Pyruvate [m] 2.56 7
Ribulose-5-phosphate [c] 2.43 4
Threonine [c] 2.28 6
10-formyltetrahydrofolate [c] 2.27 5
Fumarate [c] 2.27 5
L-proline [c] 2.04 4
Reporter subsystem Z-score No of reactions
Valine, Leucine, and Isoleucine Metabolism 3.97 15
Tyrosine, Tryptophan, and Phenylalanine Metabolism 3.39 23
Pentose Phosphate Pathway 3.29 11
Purine and Pyrimidine Biosynthesis 3.08 40
Arginine and Proline Metabolism 2.96 19
Threonine and Lysine Metabolism 2.74 14
NAD Biosynthesis 2.66 7
Alanine and Aspartate Metabolism 2.65 6
Histidine Metabolism 2.24 10
Cysteine Metabolism 1.85 10
Mo et al. BMC Systems Biology 2009 3:37   http://dx.doi.org:/10.1186/1752-0509-3-37
Open Data

Perturbations under aerobic conditions largely consisted of pathways involved in mediating the NADH and NADPH balance. Among the highest scoring aerobic subsystems
are TCA cycle and pentose phosphate pathway – key pathways directly involved in the generation of NADH and NADPH. Reporter metabolites involved in these
subsystems –

  • glyceraldehyde-3-phosphate,
  • ribulose-5-phosphate, and
  • alpha-ketoglutarate – were also identified.

These results are consistent with flux and enzyme activity measurements

  • of the gdh1/GDH2 strain under aerobic conditions [32],
  1. which reported significant reduction in the pentose phosphate pathway flux
  2. with concomitant changes in other central metabolic pathways.

Levels of several TCA cycle intermediates (e.g. fumarate, succinate, malate) were also elevated

  • in the gdh1/GDH2 mutant according to the differential intracellular metabolite data.

Altered energy metabolism, as indicated by

  • reporter metabolites (i.e. ubiquinone- , ubiquinol, mitochondrial proton)
  • and subsystem (oxidative phosphorylation),

is certainly feasible as NADH is a primary reducing agent for ATP production.

Pentose phosphate pathway and NAD biosynthesis also appears

  • among the most perturbed anaerobic subsystems, further suggesting
  • perturbed cofactor balance as a common, dominant effect under both conditions.

Glutamate dehydrogenase is a critical enzyme of amino acid biosynthesis as it acts as

  • the entry point for ammonium assimilation via glutamate.

Consequently, metabolic subsystems involved in amino acid biosynthesis were broadly perturbed

  • as a result of the gdh1/GDH2 modification in both aerobic and anaerobic conditions.

For example, the proline biosynthesis pathway that uses glutamate as a precursor

  • was significantly perturbed in both conditions,
  • with significantly changed intracellular and extracellular levels.

There were differences, however, in that more amino acid related subsystems were

  • significantly affected in the anaerobic case (Table 4),
  • further highlighting that altered ammonium assimilation in the mutant
  • has a more widespread effect under anaerobic conditions.

This effect is especially pronounced for

  • threonine and nucleotide metabolism,
  • which were predicted to be significantly perturbed only in anaerobic conditions.

Intracellular threonine levels were amongst the most significantly reduced

  • relative to other differential intracellular metabolites in the anaerobically grown gdh1/GDH2 strain
    (see [31] and Additional file 4), and
  • the relationship between threonine and nucleotide biosynthesis is further supported

by threonine’s recently discovered role as a key precursor in yeast nucleotide biosynthesis [51].

Other key anaerobic reporter metabolites are

  • glycine and 10-formyltetrahydrofolate,
  • both of which are involved in the cytosolic folate cycle (one-carbon metabolism).

Folate is intimately linked to biosynthetic pathways of

  • glycine (with threonine as its precursor) and purines
  • by mediating one-carbon reaction transfers necessary in their metabolism and
  • is a key cofactor in cellular growth [52].

Thus, the anaerobic perturbations identified in the analysis emphasize the close relationship

  • between threonine, folate, and nucleotide metabolic pathways as well as
  • their potential connection to perturbed ammonium assimilation processes.

Interestingly, this association has been previously demonstrated at the transcriptional level

  • as yeast ammonium assimilation (via glutamine synthesis) was found to be
  • co-regulated with genes involved in glycine, folate, and purine synthesis [53].

In summary, the overall differences in predicted gdh1/GDH2 mutant behavior

  • under aerobic and anaerobic conditions show that changes in flux states
  • directly related to modified ammonium assimilation pathway
  1. are amplified anaerobically whereas the
  2. indirect effects through NADH/NADPH balance are more significant aerobically.

Perturbed metabolic regions under aerobic conditions were predominantly

  • in central metabolic pathways involved in responding to the changed NADH/NADPH demand
  • and did not necessarily emphasize that glutamate dehydrogenase was the site of the genetic modification.

The majority of affected anaerobic pathways were involved directly

  • in modified ammonium assimilation as evidenced by

1) significantly perturbed amino acid subsystems,

2) a broad perturbation subnetwork converging on glutamate (Figure 4), and

3) glutamate as the most significant reporter metabolite (Table 4).

Potassium-limited and excess ammonium environments

A recent study reported that potassium limitation resulted in significant

  • growth retardation effect in yeast due to excess ammonium uptake
  • when ammonium was provided as the sole nitrogen source [33].

The proposed mechanism for this effect was that ammonium

  • could to be freely transported through potassium channels
  • when potassium concentrations were low in the media environment, thereby
  • resulting in excess ammonium uptake [33].

As a result, yeast incurred a significant metabolic cost

  • in assimilating ammonia to glutamate and
  • secreting significant amounts of glutamate and other amino acids
  • in potassium-limited conditions as a means to detoxify the excess ammonium.

A similar effect was observed when yeast was grown

  • with no potassium limitation,
  • but with excess ammonia in the environment.

While the observed effect of both environments (low potassium or excess ammonia) was similar,

  • quantitatively unique amino acid secretion profiles suggested that
  • internal metabolic states in these conditions are potentially different.

In order to elucidate the differences in internal metabolic states, we utilized

  • the iMM904 model and the EM profile analysis method to analyze amino acid secretion profiles
  • for a range of low potassium and high ammonia conditions reported in [33].

As before, we utilized amino acid secretion patterns as constraints to the iMM904 model,

  1. sampled the allowable solution space,
  2. computed reaction Z-scores for changes from a reference condition (normal potassium and ammonia), and
  3. finally summarized the resulting changes using reporter metabolites.

Figure 5 shows a clustering of the most significant reporter metabolites (Z ≥ 1.96 in any of the four conditions studied)

  • obtained from this analysis across the four conditions studied.

Interestingly, the potassium-limited environment perturbed only a subset of

  • the significant reporter metabolites identified in the high ammonia environments.

Both low potassium environments shared a consistent pattern of

  • highly perturbed amino acids and related precursor biosynthesis metabolites
    (e.g. pyruvate, PRPP, alpha-ketoglutarate)
  • with high ammonium environments.

The amino acid perturbation pattern (indicated by red labels in Figure 5) was present in

  • the ammonium-toxic environments, although the pattern was
  • slightly weaker for the lower ammonium concentration.

Nevertheless, the results clearly indicate that a similar

  • ammonium detoxifying mechanism that primarily perturbs pathways
  • directly related to amino acid metabolism
  • exists under both types of media conditions.

Figure 5.

Clustergram of top reporter metabolites - y in ammonium-toxic and potassium-limited conditions

Clustergram of top reporter metabolites – y in ammonium-toxic and potassium-limited conditions

Clustergram of top reporter metabolites (i.e. in yellow) in ammonium-toxic and potassium-limited conditions.

Amino acid perturbation patterns (shown in red) were shown to be consistently scored across conditions, indicating that potassium-limited environments K1 (lowest
concentration) and K2 (low concentration) elicited a similar ammonium detoxification response as ammonium-toxic environments N1 (high concentration) and N2
(highest concentration). Metabolites associated with folate metabolism (highlighted in green) are also highly perturbed in ammonium-toxic conditions. Metabolite
abbreviations are found in Additional file 1.

In addition to perturbed amino acids, a secondary effect notably appears at high ammonia levels in which metabolic regions related to folate metabolism are significantly affected. As highlighted in green in Figure 3, we predicted significantly perturbed key metabolites involved in the cytosolic folate cycle. These include tetrahydrofolate derivatives and other metabolites connected to the folate pathway, namely glycine and the methionine-derived methylation cofactors S-adenosylmethionine and S-adenosyl-homocysteine. Additionally, threonine was identified to be a key perturbed metabolite in excess ammonium conditions. These results further illustrate the close
connection between threonine biosynthesis, folate metabolism involving glycine derived from its threonine precursor, and nucleotide biosynthesis [51] that was discussed in
conjunction with the gdh1/GDH2 strain data. Taken together with the anaerobic gdh1/GDH2 data, the results consistently suggest highly perturbed threonine and folate
metabolism when amino acid-related pathways are broadly affected.

In both ammonium-toxic and potassium-limited environments, impaired cellular growth was observed, which can be attributed to high energetic costs of increased
ammonium assimilation to synthesize and excrete amino acids. However, under high ammonium environments, reporter metabolites related to threonine and folate
metabolism indicated that their perturbation, and thus purine supply, may be an additional factor in decreasing cellular viability as there is a direct relationship between
intracellular folate levels and growth rate [54]. Based on these results, we concluded that while potassiumlimited growth in yeast indeed shares physiological features with
growth in ammonium excess, its effects are not as detrimental as actual ammonium excess. The effects on proximal amino acid metabolic pathways are similar in both
environments as indicated by the secretion of the majority of amino acids. However, when our method was applied to analyze the physiological basis behind differences in
secretion profiles between low potassium and high ammonium conditions, ammonium excess was predicted to likely disrupt physiological ammonium assimilation processes,
which in turn potentially impacts folate metabolism and associated cellular growth.

Conclusion

The method presented in this study presents an approach to connecting intracellular flux states to metabolites that are excreted under various physiological conditions. We
showed that well-curated genome-scale metabolic networks can be used to integrate and analyze quantitative EM data by systematically identifying altered intracellular
pathways related to measured changes in the extracellular metabolome. We were able to identify statistically significant metabolic regions that were altered as a result of
genetic (gdh1/GD2 mutant) and environmental (excess ammonium and limited potassium) perturbations, and the predicted intracellular metabolic changes were consistent
with previously published experimental data including measurements of intracellular metabolite levels and metabolic fluxes. Our reanalysis of previously published EM data
on ammonium assimilation-related genetic and environmental perturbations also resulted in testable hypotheses about the role of threonine and folate pathways in mediating
broad responses to changes in ammonium utilization. These studies also demonstrated that the samplingbased method can be readily applied when only partial secreted
metabolite profiles (e.g. only amino acids) are available.

With the emergence of metabolite biofluid biomarkers as a diagnostic tool in human disease [55,56] and the availability of genome-scale human metabolic networks [1],
extensions of the present method would allow identifying potential pathway changes linked to these biomarkers. Employing such a method for studying yeast metabolism was possible as the metabolomic data was measured under controllable environmental conditions where the inputs and outputs of the system were defined. Measured metabolite biomarkers in a clinical setting, however, is far from a controlled environment with significant variations in genetic, nutritional, and environmental factors between different
patients. While there are certainly limitations for clinical applications, the method introduced here is a progressive step towards applying genome-scale metabolic networks
towards analyzing biofluid metabolome data as it 1) avoids the need to only study optimal metabolic states based on a predetermined objective function, 2) allows dealing with noisy experimental data through the sampling approach, and 3) enables analysis even with limited identification of metabolites in the data. The ability to establish potential
connections between extracellular markers and intracellular pathways would be valuable in delineating the genetic and environmental factors associated with a particular
disease.

Authors’ contributions

Conceived and designed the experiments: MLM MJH BOP. Performed experiments: MLM MJH. Analyzed the data: MLM MJH. Wrote the paper: MLM MJH BOP. All authors have read and approved the final manuscript.

Acknowledgements

We thank Jens Nielsen for providing the raw metabolome data for the mutant strain, and Jan Schellenberger and Ines Thiele for valuable discussions. This work was supported by NIH grant R01 GM071808. BOP serves on the scientific advisory board of Genomatica Inc.

 

References

  1. Duarte NC, Becker SA, Jamshidi N, Thiele I, Mo ML, Vo TD, Srivas R, Palsson BO: Global reconstruction of the human metabolic network based on genomic and bibliomic data. 

Proc Natl Acad Sci USA 2007, 104(6):1777-1782. PubMed Abstract | Publisher Full Text | PubMed Central Full Text

  1. Duarte NC, Herrgard MJ, Palsson B: Reconstruction and Validation of Saccharomyces cerevisiae iND750, a Fully Compartmentalized Genome-Scale Metabolic Model. 

Genome Res 2004, 14(7):1298-1309. PubMed Abstract | Publisher Full Text | PubMed Central Full Text

  1. Kuepfer L, Sauer U, Blank LM: Metabolic functions of duplicate genes in Saccharomyces cerevisiae. 

Genome Res 2005, 15(10):1421-1430. PubMed Abstract | Publisher Full Text | PubMed Central Full Text

  1. Nookaew I, Jewett MC, Meechai A, Thammarongtham C, Laoteng K, Cheevadhanarak S, Nielsen J, Bhumiratana S: The genome-scale metabolic model iIN800 of Saccharomyces cerevisiae and its validation: a scaffold to query lipid metabolism. 

BMC Syst Biol 2008, 2:71. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text

  1. Edwards JS, Palsson BO: Systems properties of the Haemophilus influenzae Rd metabolic genotype. 

J biol chem 1999, 274(25):17410-17416. PubMed Abstract | Publisher Full Text

  1. Edwards JS, Palsson BO: The Escherichia coli MG1655 in silico metabolic genotype: Its definition, characteristics, and capabilities. 

Proc Natl Acad Sci USA 2000, 97(10):5528-5533. PubMed Abstract | Publisher Full Text | PubMed Central Full Text

  1. Thiele I, Vo TD, Price ND, Palsson B: An Expanded Metabolic Reconstruction of Helicobacter pylori (IT341 GSM/GPR): An in silico genome-scale characterization of single and double deletion mutants. 

J Bacteriol 2005, 187(16):5818-5830. PubMed Abstract | Publisher Full Text | PubMed Central Full Text

  1. Vo TD, Greenberg HJ, Palsson BO: Reconstruction and functional characterization of the human mitochondrial metabolic network based on proteomic and biochemical data. 

J Biol Chem 2004, 279(38):39532-39540. PubMed Abstract | Publisher Full Text

  1. Reed JL, Vo TD, Schilling CH, Palsson BO: An expanded genome-scale model of Escherichia coli K-12 (iJR904 GSM/GPR). 

Genome Biology 2003, 4(9):R54.51-R54.12. BioMed Central Full Text

  1. Feist AM, Henry CS, Reed JL, Krummenacker M, Joyce AR, Karp PD, V H, Palsson BO: A genome-scale metabolic reconstruction for Escherichia coli K-12 MG1655 that accounts for 1261 ORFs and thermodynamic information. 

Molecular Systems Biology 2007, 3:121. PubMed Abstract | Publisher Full Text | PubMed Central Full Text

  1. Suthers PF, Dasika MS, Kumar VS, Denisov G, Glass JI, Maranas CD: A genome-scale metabolic reconstruction of Mycoplasma genitalium, iPS189. 

PLos Comp Biol 2009, 5(2):e1000285. Publisher Full Text

  1. Price ND, Reed JL, Palsson BO: Genome-scale models of microbial cells: evaluating the consequences of constraints. 

Nat Rev Microbiol 2004, 2(11):886-897. PubMed Abstract | Publisher Full Text

  1. Becker SA, Feist AM, Mo ML, Hannum G, Palsson BO, Herrgard MJ: Quantitative Prediction of Cellular Metabolism with Constraint-based Models: The COBRA Toolbox. 

Nature protocols 2007, 2(3):727-738. PubMed Abstract | Publisher Full Text

  1. Reed JL, Palsson BO: Genome-Scale In Silico Models of E. coli Have Multiple Equivalent Phenotypic States: Assessment of Correlated Reaction Subsets That Comprise Network States. 

Genome Res 2004, 14(9):1797-1805. PubMed Abstract | Publisher Full Text | PubMed Central Full Text

  1. Fong SS, Palsson BO: Metabolic gene deletion strains of Escherichia coli evolve to computationally predicted growth phenotypes. 

Nature Genetics 2004, 36(10):1056-1058. PubMed Abstract | Publisher Full Text

  1. Ibarra RU, Edwards JS, Palsson BO: Escherichia coli K-12 undergoes adaptive evolution to achieve in silico predicted optimal growth. 

Nature 2002, 420(6912):186-189. PubMed Abstract | Publisher Full Text

  1. Schellenberger J, Palsson BØ: Use of randomized sampling for analysis of metabolic networks. 

J Biol Chem 2009, 284(9):5457-5461. PubMed Abstract | Publisher Full Text

  1. Almaas E, Kovács B, Vicsek T, Oltvai ZN, Barabási AL: Global organization of metabolic fluxes in the bacterium Escherichia coli. 

Nature 2004, 427(6977):839-843. PubMed Abstract | Publisher Full Text

  1. Thiele I, Price ND, Vo TD, Palsson BO: Candidate metabolic network states in human mitochondria: Impact of diabetes, ischemia, and diet. 

J Biol Chem 2005, 280(12):11683-11695. PubMed Abstract | Publisher Full Text

  1. Kell DB: Metabolomics and systems biology: making sense of the soup. 

Curr Opin Microbiol 2004, 7(3):296-307. PubMed Abstract | Publisher Full Text

  1. Kell DB, Brown M, Davey HM, Dunn WB, Spasic I, Oliver SG: Metabolic footprinting and systems biology: the medium is the message. 

Nat Rev Microbiol 2005, 3(7):557-565. PubMed Abstract | Publisher Full Text

  1. Goodacre R, Vaidyanathan S, Dunn WB, Harrigan GG, Kell DB: Metabolomics by numbers: acquiring and understanding global metabolite data. 

Trends Biotechnol 2004, 22(5):245-252. PubMed Abstract | Publisher Full Text

  1. Lenz EM, Bright J, Wilson ID, Morgan SR, Nash AF: A 1H NMR-based metabonomic study of urine and plasma samples obtained from healthy human subjects. 

J Pharm Biomed Anal 2003, 33(5):1103-1115. PubMed Abstract | Publisher Full Text

  1. Allen J, Davey HM, Broadhurst D, Heald JK, Rowland JJ, Oliver SG, Kell DB: High-throughput classification of yeast mutants for functional genomics using metabolic footprinting. 

Nat Biotech 2003, 21(6):692-696. Publisher Full Text

  1. Nicholson JK, Connelly J, Lindon JC, Holmes E: Metabonomics: a platform for studying drug toxicity and gene function. 

Nat Rev Drug Discov 2002, 1(2):153-161. PubMed Abstract | Publisher Full Text

  1. Mortishire-Smith RJ, Skiles GL, Lawrence JW, Spence S, Nicholls AW, Johnson BA, Nicholson JK: Use of metabonomics to identify impaired fatty acid metabolism as the mechanism of a drug-induced toxicity. 

Chem Res Toxicol 2004, 17(2):165-173. PubMed Abstract | Publisher Full Text

  1. Sabatine MS, Liu E, Morrow DA, Heller E, McCarroll R, Wiegand R, Berriz GF, Roth FP, Gerszten RE: Metabolomic identification of novel biomarkers of myocardial ischemia. 

Circulation 2005, 112(25):3868-3875. PubMed Abstract | Publisher Full Text

  1. Cakir T, Efe C, Dikicioglu D, Hortaçsu AKB, Oliver SG: Flux balance analysis of a genome-scale yeast model constrained by exometabolomic data allows metabolic system identification of genetically different strains. 

Biotechnol Prog 2007, 23(2):320-326. PubMed Abstract | Publisher Full Text

  1. Bang JW, Crockford DJ, Holmes E, Pazos F, Sternberg MJ, Muggleton SH, Nicholson JK: Integrative top-down system metabolic modeling in experimental disease states via data-driven Bayesian methods. 

J Proteome Res 2008, 7(2):497-503. PubMed Abstract | Publisher Full Text

  1. Oliveira AP, Patil KR, Nielsen J: Architecture of transcriptional regulatory circuits is knitted over the topology of bio-molecular interaction networks. 

BMC Syst Biol 2008, 2:17. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text

  1. Villas-Boas SG, Moxley JF, Akesson M, Stephanopoulos G, Nielsen J: High-throughput metabolic state analysis: the missing link in integrated functional genomics of yeasts. 

Biochem J 2005, 388(Pt 2):669-677. PubMed Abstract | Publisher Full Text | PubMed Central Full Text

  1. Moreira dos Santos M, Thygesen G, Kötter P, Olsson L, Nielsen J: Aerobic physiology of redox-engineered Saccharomyces cerevisiae strains modified in the ammonium assimilation for increased NADPH availability. 

FEMS Yeast Res 2003, 4(1):59-68. PubMed Abstract | Publisher Full Text

  1. Hess DC, Lu W, Rabinowitz JD, Botstein D: Ammonium toxicity and potassium limitation in yeast. 

PLoS Biol 2006, 4(11):e351. PubMed Abstract | Publisher Full Text | PubMed Central Full Text

  1. Nissen TL, Schulze U, Nielsen J, Villadsen J: Flux distributions in anaerobic, glucose-limited continuous cultures of Saccharomyces cerevisiae. 

Microbiology 1997, 143(Pt 1):203-218. PubMed Abstract | Publisher Full Text

  1. Famili I, Forster J, Nielsen J, Palsson BO: Saccharomyces cerevisiae phenotypes can be predicted by using constraint-based analysis of a genome-scale reconstructed metabolic network. 

Proc Natl Acad Sci USA 2003, 100(23):13134-13139. PubMed Abstract | Publisher Full Text | PubMed Central Full Text

  1. Bonarius HPJ, Schmid G, Tramper J: Flux analysis of underdetermined metabolic networks: The quest for the missing constraints. 

Trends in Biotechnology 1997, 15(8):308-314. Publisher Full Text

  1. Edwards JS, Palsson BO: Metabolic flux balance analysis and the in silico analysis of Escherichia coli K-12 gene deletions. 

BMC Bioinformatics 2000, 1:1. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text

  1. Varma A, Palsson BO: Metabolic Flux Balancing: Basic concepts, Scientific and Practical Use. 

Nat Biotechnol 1994, 12:994-998. Publisher Full Text

  1. Segre D, Vitkup D, Church GM: Analysis of optimality in natural and perturbed metabolic networks. 

Proc Natl Acad Sci USA 2002, 99(23):15112-15117. PubMed Abstract | Publisher Full Text | PubMed Central Full Text

  1. Baldi P, Brunak S, Chauvin Y, Andersen CA, Nielsen H: Assessing the accuracy of prediction algorithms for classification: an overview. 

Bioinformatics 2000, 16(5):412-424. PubMed Abstract | Publisher Full Text

  1. Price ND, Schellenberger J, Palsson BO: Uniform Sampling of Steady State Flux Spaces: Means to Design Experiments and to Interpret Enzymopathies. 

Biophysical Journal 2004, 87(4):2172-2186. PubMed Abstract | Publisher Full Text | PubMed Central Full Text

  1. Price ND, Thiele I, Palsson BO: Candidate states of Helicobacter pylori’s genome-scale metabolic network upon application of “loop law” thermodynamic constraints. 

Biophysical J 2006, 90(11):3919-3928. Publisher Full Text

  1. Schuetz R, Kuepfer L, Sauer U: Systematic evaluation of objective functions for predicting intracellular fluxes in Escherichia coli. 

Mol Syst Biol 2007, 3:119. PubMed Abstract | Publisher Full Text | PubMed Central Full Text

  1. Patil KR, Nielsen J: Uncovering transcriptional regulation of metabolism by using metabolic network topology. 

Proc Natl Acad Sci USA 2005, 102(8):2685-2689. PubMed Abstract | Publisher Full Text | PubMed Central Full Text

  1. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: a software environment for integrated models of biomolecular interaction networks. 

Genome Res 2003, 13(11):2498-2504. PubMed Abstract | Publisher Full Text | PubMed Central Full Text

  1. Herrgard MJ, Swainston N, Dobson P, Dunn WB, Arga KY, Arvas M, Bluthgen N, Borger S, Costenoble R, Heinemann M, et al.: A consensus yeast metabolic network reconstruction obtained from a community approach to systems biology. 

Nat Biotech 2008, 26:1155-1160. Publisher Full Text

  1. Forster J, Famili I, Fu PC, Palsson BO, Nielsen J: Genome-Scale Reconstruction of the Saccharomyces cerevisiae Metabolic Network. 

Genome Research 2003, 13(2):244-253. PubMed Abstract | Publisher Full Text | PubMed Central Full Text

  1. Reed JL, Patel TR, Chen KH, Joyce AR, Applebee MK, Herring CD, Bui OT, Knight EM, Fong SS, Palsson BO: Systems Approach to Genome Annotation: Prediction and Validation of Metabolic Functions. 

Proc Natl Acad Sci USA 2006, 103(46):17480-17484. PubMed Abstract | Publisher Full Text | PubMed Central Full Text

  1. Nissen TL, Kielland-Brandt MC, Nielsen J, Villadsen J: Optimization of ethanol production in Saccharomyces cerevisiae by metabolic engineering of the ammonium assimilation. 

Metab Eng 2000, 2(1):69-77. PubMed Abstract | Publisher Full Text

  1. Roca C, Nielsen J, Olsson L: Metabolic engineering of ammonium assimilation in xylose-fermenting Saccharomyces cerevisiae improves ethanol production. 

Appl Environ Microbiol 2003, 69(8):4732-4736. PubMed Abstract | Publisher Full Text | PubMed Central Full Text

  1. Hartman JL IV: Buffering of deoxyribonucleotide pool homeostasis by threonine metabolism. 

Proc Natl Acad Sci USA 2007, 104(28):11700-11705. PubMed Abstract | Publisher Full Text | PubMed Central Full Text

  1. Gelling CL, Piper MD, Hong SP, Kornfeld GD, Dawes IW: Identification of a novel one-carbon metabolism regulon in Saccharomyces cerevisiae. 

J Biol Chem 2004, 279(8):7072-7081. PubMed Abstract | Publisher Full Text

  1. Denis V, Daignan-Fornier B: Synthesis of glutamine, glycine and 10-formyl tetrahydrofolate is coregulated with purine biosynthesis in Saccharomyces cerevisiae. 

Mol Gen Genet 1998, 259(3):246-255. PubMed Abstract | Publisher Full Text

  1. Hjortmo S, Patring J, Andlid T: Growth rate and medium composition strongly affect folate content in Saccharomyces cerevisiae. 

Int J Food Microbiol 2008, 123(1–2):93-100. PubMed Abstract | Publisher Full Text

  1. Kussmann MRF, Affolter M: OMICS-driven biomarker discovery in nutrition and health. 

J Biotechnol 2006, 124(4):758-787. PubMed Abstract | Publisher Full Text

  1. Serkova NJ, Niemann CU: Pattern recognition and biomarker validation using quantitative 1H-NMR-based metabolomics. 

Expert Rev Mol Diagn 2006, 6(5):717-731. PubMed Abstract | Publisher Full Text

 

Read Full Post »

Metabolomics: its Applications in Food and Nutrition Research

Reporter and Curator: Sudipta Saha, Ph.D.

 

Metabolomics is a relatively new field of “omics” research concerned with the high-throughput identification and quantification of small molecule (<1500 Da) metabolites in the metabolome. The metabolome is formally defined as the collection of all small molecule metabolites or chemicals that can be found in a cell, organ or organism. These small molecules can include a range of endogenous and exogenous chemical entities such as peptides, amino acids, nucleic acids, carbohydrates, organic acids, vitamins, polyphenols, alkaloids, minerals and just about any other chemical that can be used, ingested or synthesized by a given cell or organism.

Metabolomics is ideally positioned to be used in many areas of food science and nutrition research including food component analysis, food quality/authenticity assessment, food consumption monitoring and physiological monitoring in food intervention studies. However, the potential impact of metabolomics is still limited by two factors: (1) technology and (2) databases. In terms of instrumentation, it is clear that significant improvements need to be made to make metabolite detection and quantification technology more robust, automated and comprehensive. While promising advances have been made, current techniques are only capable of detecting perhaps 1/10th of the relevant metabolome. This expanded breadth and depth of coverage is particularly important in food and nutrition studies.

Many more reference spectral or chromatographic databases on metabolites, food components and phytochemicals need to be developed and made public. It is only through these databases that nutritionally relevant compounds can be routinely identified or quantified. Indeed a comprehensive effort, similar to that undertaken to annotate the human metabolome, needs to be made to complete and annotate the “food metabolome”. Similar efforts also need to be directed towards creating publicly accessible, comprehensive nutritional phenotype databases that include quantitative metabolomic (and other omic) data collected from diet-challenge or food intervention experiments. While these kinds of endeavours may take years to complete and cost millions of dollars, hopefully the food science community (and its funding agencies) will find a way of coordinating its activities to complete these efforts. Indeed, having public resource like a food metabolome database or a nutritional phenotype database could be as valuable to food scientists as GenBank has been to molecular biologists.

Source References:

http://www.sciencedirect.com/science/article/pii/S0924224408000770

http://www.sciencedirect.com/science/article/pii/B9780123945983000010

http://www.sciencedirect.com/science/article/pii/S092422440900226X

http://www.sciencedirect.com/science/article/pii/S1359644605036093

http://www.sciencedirect.com/science/article/pii/B9780080885049000520

http://www.sciencedirect.com/science/article/pii/B9780123744135000051

Other articles related to this topic were published on this Open Access Online Scientific Journal, including the following:

Ca2+ signaling: transcriptional control

Larry H. Bernstein, MD, FCAP, Reporter, RN 03/06/2013

http://pharmaceuticalintelligence.com/2013/03/06/ca2-signaling-transcriptional-control/

Harnessing Personalized Medicine for Cancer Management, Prospects of Prevention and Cure: Opinions of Cancer Scientific Leaders @ http://pharmaceuticalintelligence.com

Aviva Lev-Ari, PhD, RN 01/12/2013

http://pharmaceuticalintelligence.com/2013/01/12/harnessing-personalized-medicine-for-cancer-management-prospects-of-prevention-and-cure-opinions-of-cancer-scientific-leaders-httppharmaceuticalintelligence-com/

Breakthrough Digestive Disorders Research: Conditions affecting the Gastrointestinal Tract.

Aviva Lev-Ari, PhD, RN 12/12/2012

http://pharmaceuticalintelligence.com/2012/12/12/breakthrough-digestive-disorders-research-conditions-affecting-the-gastrointestinal-tract/

A Second Look at the Transthyretin Nutrition Inflammatory Conundrum

Larry H. Bernstein, MD, FCAP, Reporter, RN 12/03/2012

http://pharmaceuticalintelligence.com/2012/12/03/a-second-look-at-the-transthyretin-nutrition-inflammatory-conundrum/

Metabolic drivers in aggressive brain tumors

Prabodh Kandala, PhD, RN 11/11/2012

http://pharmaceuticalintelligence.com/2012/11/11/metabolic-drivers-in-aggressive-brain-tumors/

Metabolite Identification Combining Genetic and Metabolic Information: Genetic association links unknown metabolites to functionally related genes

Aviva Lev-Ari, PhD, RN 10/22/2012

http://pharmaceuticalintelligence.com/2012/10/22/metabolite-identification-combining-genetic-and-metabolic-information-genetic-association-links-unknown-metabolites-to-functionally-related-genes/

Advances in Separations Technology for the “OMICs” and Clarification of Therapeutic Targets

Larry H. Bernstein, MD, FCAP, Reporter, RN 10/22/2012

http://pharmaceuticalintelligence.com/2012/10/22/advances-in-separations-technology-for-the-omics-and-clarification-of-therapeutic-targets/

Expanding the Genetic Alphabet and linking the genome to the metabolome

Larry H. Bernstein, MD, FCAP, Reporter, RN 09/24/2012

http://pharmaceuticalintelligence.com/2012/09/24/expanding-the-genetic-alphabet-and-linking-the-genome-to-the-metabolome/

Therapeutic Targets for Diabetes and Related Metabolic Disorders

Aviva Lev-Ari, PhD, RN 08/20/2012

http://pharmaceuticalintelligence.com/2012/08/20/therapeutic-targets-for-diabetes-and-related-metabolic-disorders/

The Automated Second Opinion Generator

Larry H. Bernstein, MD, FCAP, Reporter, RN 08/13/2012

http://pharmaceuticalintelligence.com/2012/08/13/the-automated-second-opinion-generator/

 

Read Full Post »

%d bloggers like this: