Metabolic network discovery by top-down and bottom-up approaches and paths for reconciliation

FastKnock: an efficient next-generation approach to identify all knockout strategies for strain optimization

Overproduction of desired native or nonnative biochemical(s) in (micro)organisms can be achieved through metabolic engineering. Appropriate rewiring of cell metabolism is performed by making rational changes such as insertion, up-/down-regulation and knockout of genes and consequently metabolic reactions. Finding appropriate targets (including proper sets of reactions to be knocked out) for metabolic engineering to design optimal production strains has been the goal of a number of computational algorithms. We developed FastKnock, an efficient next-generation algorithm for identifying all possible knockout strategies (with a predefined maximum number of reaction deletions) for the growth-coupled overproduction of biochemical(s) of interest. We achieve this by developing a special depth-first traversal algorithm that allows us to prune the search space significantly. This leads to a drastic reduction in execution time. We evaluate the performance of the FastKnock algorithm using various Escherichia coli genome-scale metabolic models in different conditions (minimal and rich mediums) for the overproduction of a number of desired metabolites. FastKnock efficiently prunes the search space to less than 0.2% for quadruple- and 0.02% for quintuple-reaction knockouts. Compared to the classic approaches such as OptKnock and the state-of-the-art techniques such as MCSEnumerator methods, FastKnock found many more beneficial and important practical solutions. The availability of all the solutions provides the opportunity to further characterize, rank and select the most appropriate intervention strategy based on any desired evaluation index. Our implementation of the FastKnock method in Python is publicly available at https://github.com/leilahsn/FastKnock .

FastKnock: An efficient next-generation approach to identify all knockout strategies for strain optimization

Overproduction of desired native or nonnative biochemical(s) in (micro)organisms can be achieved through metabolic engineering. Appropriate rewiring of cell metabolism is performed making rational changes such as insertion, up-/down-regulation and knockout of genes and consequently metabolic reactions. Finding appropriate targets (including proper sets of reactions to be knocked out) for metabolic engineering to design optimal production strains has been the goal of a number of computational algorithms. We developed FastKnock, an efficient next-generation algorithm for identifying all possible knockout strategies for the growth-coupled overproduction of biochemical(s) of interest. We achieve this by developing a special depth-first traversal algorithm that allows us to prune the search space significantly. This leads to a drastic reduction in execution time. We evaluate the performance of the FastKnock algorithm using three Escherichia coli genome-scale metabolic models in different conditions (minimal and rich mediums) for the overproduction of a number of desired metabolites. FastKnock efficiently prunes the search space to less than 0.2% for quadruple and 0.02% for quintuple-reaction knockouts. Compared to the classic approaches such as OptKnock and the state-of-the-art techniques such as MCSEnumerator methods, FastKnock found many more useful and important practical solutions. The availability of all the solutions provides the opportunity to further characterize and select the most appropriate intervention strategy based on any desired evaluation index. Our implementation of the FastKnock method in Python is publicly available at https://github.com/leilahsn/FastKnock.

Free Amino Acid Alterations in Patients with Gynecological and Breast Cancer: A Review

Gynecological and breast cancers still remain a significant health problem worldwide. Diagnostic methods are not sensitive and specific enough to detect the disease at an early stage. During carcinogenesis and tumor progression, the cellular need for DNA and protein synthesis increases leading to changes in the levels of amino acids. An important role of amino acids in many biological pathways, including biosynthesis of proteins, nucleic acids, enzymes, etc., which serve as an energy source and maintain redox balance, has been highlighted in many research articles. The aim of this review is a detailed analysis of the literature on metabolomic studies of gynecology and breast cancers with particular emphasis on alterations in free amino acid profiles. The work includes a brief overview of the metabolomic methodology and types of biological samples used in the studies. Special attention was paid to the possible role of selected amino acids in the carcinogenesis, especially proline and amino acids related to its metabolism. There is a clear need for further research and multiple external validation studies to establish the role of amino acid profiling in diagnosing gynecological and breast cancers.

Stoichiometric gene-to-reaction associations enhance model-driven analysis performance: Metabolic response to chronic exposure to Aldrin in prostate cancer

BACKGROUND

Genome-scale metabolic models (GSMM) integrating transcriptomics have been widely used to study cancer metabolism. This integration is achieved through logical rules that describe the association between genes, proteins, and reactions (GPRs). However, current gene-to-reaction formulation lacks the stoichiometry describing the transcript copies necessary to generate an active catalytic unit, which limits our understanding of how genes modulate metabolism. The present work introduces a new state-of-the-art GPR formulation that considers the stoichiometry of the transcripts (S-GPR). As case of concept, this novel gene-to-reaction formulation was applied to investigate the metabolic effects of the chronic exposure to Aldrin, an endocrine disruptor, on DU145 prostate cancer cells. To this aim we integrated the transcriptomic data from Aldrin-exposed and non-exposed DU145 cells through S-GPR or GPR into a human GSMM by applying different constraint-based-methods.

RESULTS

Our study revealed a significant improvement of metabolite consumption/production predictions when S-GPRs are implemented. Furthermore, our computational analysis unveiled important alterations in carnitine shuttle and prostaglandine biosynthesis in Aldrin-exposed DU145 cells that is supported by bibliographic evidences of enhanced malignant phenotype.

CONCLUSIONS

The method developed in this work enables a more accurate integration of gene expression data into model-driven methods. Thus, the presented approach is conceptually new and paves the way for more in-depth studies of aberrant cancer metabolism and other diseases with strong metabolic component with important environmental and clinical implications.

JacLy: a Jacobian-based method for the inference of metabolic interactions from the covariance of steady-state metabolome data

Reverse engineering metabolome data to infer metabolic interactions is a challenging research topic. Here we introduce JacLy, a Jacobian-based method to infer metabolic interactions of small networks (<20 metabolites) from the covariance of steady-state metabolome data. The approach was applied to two different in silico small-scale metabolome datasets. The power of JacLy lies on the use of steady-state metabolome data to predict the Jacobian matrix of the system, which is a source of information on structure and dynamic characteristics of the system. Besides its advantage of inferring directed interactions, its superiority over correlation-based network inference was especially clear in terms of the required number of replicates and the effect of the use of priori knowledge in the inference. Additionally, we showed the use of standard deviation of the replicate data as a suitable approximation for the magnitudes of metabolite fluctuations inherent in the system.

The best models of metabolism

Biochemical systems are among of the oldest application areas of mathematical modeling. Spanning a time period of over one hundred years, the repertoire of options for structuring a model and for formulating reactions has been constantly growing, and yet, it is still unclear whether or to what degree some models are better than others and how the modeler is to choose among them. In fact, the variety of options has become overwhelming and difficult to maneuver for novices and experts alike. This review outlines the metabolic model design process and discusses the numerous choices for modeling frameworks and mathematical representations. It tries to be inclusive, even though it cannot be complete, and introduces the various modeling options in a manner that is as unbiased as that is feasible. However, the review does end with personal recommendations for the choices of default models. WIREs Syst Biol Med 2017, 9:e1391. doi: 10.1002/wsbm.1391 For further resources related to this article, please visit the WIREs website.

Current advances in biomarkers for targeted therapy in triple-negative breast cancer

Triple-negative breast cancer (TNBC) is a complex heterogeneous disease characterized by the absence of three hallmark receptors: human epidermal growth factor receptor 2, estrogen receptor, and progesterone receptor. Compared to other breast cancer subtypes, TNBC is more aggressive, has a higher prevalence in African-Americans, and more frequently affects younger patients. Currently, TNBC lacks clinically accepted targets for tailored therapy, warranting the need for candidate biomarkers. BiomarkerBase, an online platform used to find biomarkers reported in clinical trials, was utilized to screen all potential biomarkers for TNBC and select only the ones registered in completed TNBC trials through clinicaltrials.gov. The selected candidate biomarkers were classified as surrogate, prognostic, predictive, or pharmacodynamic (PD) and organized by location in the blood, on the cell surface, in the cytoplasm, or in the nucleus. Blood biomarkers include vascular endothelial growth factor/vascular endothelial growth factor receptor and interleukin-8 (IL-8); cell surface biomarkers include EGFR, insulin-like growth factor binding protein, c-Kit, c-Met, and PD-L1; cytoplasm biomarkers include PIK3CA, pAKT/S6/p4E-BP1, PTEN, ALDH1, and the PIK3CA/AKT/mTOR-related metabolites; and nucleus biomarkers include BRCA1, the gluco-corticoid receptor, TP53, and Ki67. Candidate biomarkers were further organized into a "cellular protein network" that demonstrates potential connectivity. This review provides an inventory and reference point for promising biomarkers for breakthrough targeted therapies in TNBC.

A kidney-specific genome-scale metabolic network model for analyzing focal segmental glomerulosclerosis

Focal Segmental Glomerulosclerosis (FSGS) is a type of nephrotic syndrome which accounts for 20 and 40 % of such cases in children and adults, respectively. The high prevalence of FSGS makes it the most common primary glomerular disorder causing end-stage renal disease. Although the pathogenesis of this disorder has been widely investigated, the exact mechanism underlying this disease is still to be discovered. Current therapies seek to stop the progression of FSGS and often fail to cure the patients since progression to end-stage renal failure is usually inevitable. In the present work, we use a kidney-specific metabolic network model to study FSGS. The model was obtained by merging two previously published kidney-specific metabolic network models. The validity of the new model was checked by comparing the inactivating reaction genes identified in silico to the list of kidney disease implicated genes. To model the disease state, we used a complete list of FSGS metabolic biomarkers extracted from transcriptome and proteome profiling of patients as well as genetic deficiencies known to cause FSGS. We observed that some specific pathways including chondroitin sulfate degradation, eicosanoid metabolism, keratan sulfate biosynthesis, vitamin B6 metabolism, and amino acid metabolism tend to show variations in FSGS model compared to healthy kidney. Furthermore, we computationally searched for the potential drug targets that can revert the diseased metabolic state to the healthy state. Interestingly, only one drug target, N-acetylgalactosaminidase, was found whose inhibition could alter cellular metabolism towards healthy state.

Systems Medicine: The Application of Systems Biology Approaches for Modern Medical Research and Drug Development

The exponential development of highly advanced scientific and medical research technologies throughout the past 30 years has arrived to the point where the high number of characterized molecular agents related to pathogenesis cannot be readily integrated or processed by conventional analytical approaches. Indeed, the realization that several moieties are signatures of disease has partly led to the increment of complex diseases being characterized. Scientists and clinicians can now investigate and analyse any individual dysregulations occurring within the genomic, transcriptomic, miRnomic, proteomic, and metabolomic levels thanks to currently available advanced technologies. However, there are drawbacks within this scientific brave new age in that only isolated molecular levels are individually investigated for their influence in affecting any particular health condition. Since their conception in 1992, systems biology/medicine focuses mainly on the perturbations of overall pathway kinetics for the consequent onset and/or deterioration of the investigated condition/s. Systems medicine approaches can therefore be employed for shedding light in multiple research scenarios, ultimately leading to the practical result of uncovering novel dynamic interaction networks that are critical for influencing the course of medical conditions. Consequently, systems medicine also serves to identify clinically important molecular targets for diagnostic and therapeutic measures against such a condition.

A review on computational systems biology of pathogen-host interactions

Pathogens manipulate the cellular mechanisms of host organisms via pathogen-host interactions (PHIs) in order to take advantage of the capabilities of host cells, leading to infections. The crucial role of these interspecies molecular interactions in initiating and sustaining infections necessitates a thorough understanding of the corresponding mechanisms. Unlike the traditional approach of considering the host or pathogen separately, a systems-level approach, considering the PHI system as a whole is indispensable to elucidate the mechanisms of infection. Following the technological advances in the post-genomic era, PHI data have been produced in large-scale within the last decade. Systems biology-based methods for the inference and analysis of PHI regulatory, metabolic, and protein-protein networks to shed light on infection mechanisms are gaining increasing demand thanks to the availability of omics data. The knowledge derived from the PHIs may largely contribute to the identification of new and more efficient therapeutics to prevent or cure infections. There are recent efforts for the detailed documentation of these experimentally verified PHI data through Web-based databases. Despite these advances in data archiving, there are still large amounts of PHI data in the biomedical literature yet to be discovered, and novel text mining methods are in development to unearth such hidden data. Here, we review a collection of recent studies on computational systems biology of PHIs with a special focus on the methods for the inference and analysis of PHI networks, covering also the Web-based databases and text-mining efforts to unravel the data hidden in the literature.