Funding, Deals & Partnerships: BIOLOGICS & MEDICAL DEVICES; BioMed e-Series; Medicine and Life Sciences Scientific Journal – http://PharmaceuticalIntelligence.com
This discussion that completes and is an epicrisis (summary and critical evaluation) of the series of discussions that preceded it.
Innervation of Heart and Heart Rate
Action of hormones on the circulation
Allogeneic Transfusion Reactions
Graft-versus Host reaction
Unique problems of perinatal period
High altitude sickness
Deep water adaptation
Heart-Lung-and Kidney
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
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
mcell-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
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,
fluorescence recovery after photobleaching (FRAP),
fluorescence correlation spectroscopy (FCS),
single molecule tracking (SMT) and
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.
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.
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.
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
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.
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.
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
iron metabolism
personalized reference range within population range
Part 1. MetabolomicsSurge
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
recent advances in analytical technology and bioinformatics as well as
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
a specific disease state,
toxin exposure, or
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
targeted organs,
mechanism of action, and
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
on the safety and efficacy of compounds
during early and preclinical toxicological studies,
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
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
the identification,
quantification, and
mapping of more than 800 metabolites to specific cellular pathways.
It is based on flow injection analysis and high-performance liquid chromatography MS/MS.
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
the uncoupling of glycolysis from the mitochondria,
leading to the inefficient but rapid metabolism of glucose and
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,
can functionally coordinate cell-survival and cell-proliferation mechanisms,
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
high reproducibility,
standardized protocols,
low sample manipulation, and
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,
“splitting up information content, processing, and introducing a lot of background noise and error and
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
309 have been identified in cerebrospinal fluid,
1,122 in serum,
458 in urine, and
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
undergo a variety of interactions including metal complexation,
chemical exchange processes,
micellar compartmentation,
enzyme-mediated biotransformations, and
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
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.
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
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’
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
enzymes;
structural molecules, such as hair,
skin and bones;
hormones like insulin; and
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,
RNA polymerase II functions with
a factor termed TBP,
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
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”.
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
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.
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
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
inhibits the dephosphorylation of extracellular signal-regulated kinase (ERK) and
p38 mitogen-activated protein kinase (p38 MAPK),
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
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
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”.
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
AGPS expression increased when normal cells turned cancerous.
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
Figure 1: 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
(a) Representative images from immunohistochemical staining of Smurf1, Ubc12, NAE1 and Nedd8 in the same colorectal cancer tumour. Scale bars, 100 μm. (b–d) 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
(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
(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.
(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…
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.
HETEs harm the mitochondria, which then
fail to produce sufficient energy to enable
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.”
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
mad cow disease,
Type 1 diabetes and
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.
While work with biomarkers continues to grow, scientists are also grappling with research-related bottlenecks, such as
affinity reagent development,
platform reproducibility, and
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:
where to find them,
when they may appear,
what form they may take, and
how they can be used to diagnose a condition or
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
to garner information about the molecular characteristics of the cancer
that will help with cancer detection and
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
mass spectrometry and
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
build rigorous, formal mathematical models that will allow something measured in the blood
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
a biologically inert, biophysically stable protein scaffold
containing three variable regions into which
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
rheumatoid arthritis,
psoriatic arthritis,
SLE, or
giant cell arteritis.
Epigenetic Biomarkers
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
tumor type,
histology, or
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
somatic mutations,
gene copy number aberrations,
gene expression abnormalities,
protein and miRNA expression, and
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
prognostic drivers
predictive markers for taxanes and
monoclonal antibodies and
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 system 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
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
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.
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
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.
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.
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
enterocytes (Sahoo and Thiele2013),
macrophages (Bordbar et al. 2010),
adipocytes (Mardinoglu et al. 2013),
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
cancer metabolism (Jerby and Ruppin, 2012).
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
Such analyses have also been used to reveal the effects of
enzymopathies on red blood cells (Price et al. 2004),
to study effects of diet on diabetes (Thiele et al. 2005) and
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.
metabolic differences between two lymphoblastic leukemia cell lines (Fig. 1A).
Fig. 1
metabol leukem cell lines11306_2014_721_Fig1_HTML
A Combined experimental and computational pipeline to study human metabolism.
Experimental work and omics data analysis steps precede computational modeling.
Model predictions are validated based on targeted experimental data.
Metabolomic and transcriptomic data are used for model refinement and submodel extraction.
Functional analysis methods are used to characterize the metabolism of the cell-line models and compare it to additional experimental data.
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.
A number of metabolite exchanges mapped to the model were unique to one cell line.
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:
data acquisition,
data analysis,
metabolic modeling and
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
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
^ 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
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
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
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).
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
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
metabolic pathways 1476-4598-10-70-1
Metabolic Systems Research Team fig2
Metabolic control analysis of respiration in human cancer tissue. fphys-04-00151-g001
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
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
how changes in the extracellular metabolome can be used
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 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
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.
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.
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
metabolic states relating potassium limitation and
ammonium excess conditions to one another.
The model-based analysis of both
separately published extracellular metabolome datasets
suggests a relationship between
glutamate,
threonine and
folate metabolism,
which are collectively perturbed when ammonium assimilation processes are broadly disrupted
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.
Minimization of Metabolic Adjustment (MoMA) [39], and
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).
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.
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
A 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
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 j 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
where Nj is the number of reactions in Rj and mmet,j is calculated as
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 Rk
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
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
904 genes,
1,228 individual metabolites, and
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),
which are required in the assembly of glycogen and
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
changes to the compartmental location of a gene and
its corresponding reaction(s),
changes in reaction reversibility and cofactor specificity, and
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],
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 (p < 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 (p < 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.
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.
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.
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 p < 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
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],
which reported significant reduction in the pentose phosphate pathway flux
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.
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
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.
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.
Thiele I, Vo TD, Price ND, Palsson B: An Expanded Metabolic Reconstruction of Helicobacter pylori (i IT341 GSM/GPR): An in silico genome-scale characterization of single and double deletion mutants.
Vo TD, Greenberg HJ, Palsson BO: Reconstruction and functional characterization of the human mitochondrial metabolic network based on proteomic and biochemical data.
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.
Becker SA, Feist AM, Mo ML, Hannum G, Palsson BO, Herrgard MJ: Quantitative Prediction of Cellular Metabolism with Constraint-based Models: The COBRA Toolbox.
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.
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.
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.
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.
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.
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.
Oliveira AP, Patil KR, Nielsen J: Architecture of transcriptional regulatory circuits is knitted over the topology of bio-molecular interaction networks.
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.
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.
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.
Price ND, Thiele I, Palsson BO: Candidate states of Helicobacter pylori’s genome-scale metabolic network upon application of “loop law” thermodynamic constraints.
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.
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.
Nissen TL, Kielland-Brandt MC, Nielsen J, Villadsen J: Optimization of ethanol production in Saccharomyces cerevisiae by metabolic engineering of the ammonium assimilation.
Denis V, Daignan-Fornier B: Synthesis of glutamine, glycine and 10-formyl tetrahydrofolate is coregulated with purine biosynthesis in Saccharomyces cerevisiae.
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.
Harnessing Personalized Medicine for Cancer Management, Prospects of Prevention and Cure: Opinions of Cancer Scientific Leaders @ http://pharmaceuticalintelligence.com