Systems-level modeling of cellular glycosylation reaction networks: O-linked glycan formation on natural selectin ligands

GlycoMME, a Markov modeling platform for studying N-glycosylation biosynthesis from glycomics data

Variations in N-glycosylation, which is crucial to glycoprotein functions, impact many diseases and the safety and efficacy of biotherapeutic drugs. Here, we present a protocol for using GlycoMME (Glycosylation Markov Model Evaluator) to study N-glycosylation biosynthesis from glycomics data. We describe steps for annotating glycomics data and quantifying perturbations to N-glycan biosynthesis with interpretable models. We then detail procedures to predict the impact of mutations in disease or potential glycoengineering strategies in drug development. For complete details on the use and execution of this protocol, please refer to Liang et al. (2020).¹.

GlycoEnzOnto: a GlycoEnzyme pathway and molecular function ontology

MOTIVATION

The 'glycoEnzymes' include a set of proteins having related enzymatic, metabolic, transport, structural and cofactor functions. Currently, there is no established ontology to describe glycoEnzyme properties and to relate them to glycan biosynthesis pathways.

RESULTS

We present GlycoEnzOnto, an ontology describing 403 human glycoEnzymes curated along 139 glycosylation pathways, 134 molecular functions and 22 cellular compartments. The pathways described regulate nucleotide-sugar metabolism, glycosyl-substrate/donor transport, glycan biosynthesis and degradation. The role of each enzyme in the glycosylation initiation, elongation/branching and capping/termination phases is described. IUPAC linear strings present systematic human/machine-readable descriptions of individual reaction steps and enable automated knowledge-based curation of biochemical networks. All GlycoEnzOnto knowledge is integrated with the Gene Ontology biological processes. GlycoEnzOnto enables improved transcript overrepresentation analyses and glycosylation pathway identification compared to other available schema, e.g. KEGG and Reactome. Overall, GlycoEnzOnto represents a holistic glycoinformatics resource for systems-level analyses.

AVAILABILITY AND IMPLEMENTATION

https://github.com/neel-lab/GlycoEnzOnto.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Computational Modeling of O-Linked Glycan Biosynthesis in CHO Cells

Glycan biosynthesis simulation research has progressed remarkably since 1997, when the first mathematical model for N-glycan biosynthesis was proposed. An O-glycan model has also been developed to predict O-glycan biosynthesis pathways in both forward and reverse directions. In this work, we started with a set of O-glycan profiles of CHO cells transiently transfected with various combinations of glycosyltransferases. The aim was to develop a model that encapsulated all the enzymes in the CHO transfected cell lines. Due to computational power restrictions, we were forced to focus on a smaller set of glycan profiles, where we were able to propose an optimized set of kinetics parameters for each enzyme in the model. Using this optimized model we showed that the abundance of more processed glycans could be simulated compared to observed abundance, while predicting the abundance of glycans earlier in the pathway was less accurate. The data generated show that for the accurate prediction of O-linked glycosylation, additional factors need to be incorporated into the model to better reflect the experimental conditions.

Big-Data Glycomics: Tools to Connect Glycan Biosynthesis to Extracellular Communication

Characteristically, cells must sense and respond to environmental cues. Despite the importance of cell-cell communication, our understanding remains limited and often lacks glycans. Glycans decorate proteins and cell membranes at the cell-environment interface, and modulate intercellular communication, from development to pathogenesis. Providing further challenges, glycan biosynthesis and cellular behavior are co-regulating systems. Here, we discuss how glycosylation contributes to extracellular responses and signaling. We further organize approaches for disentangling the roles of glycans in multicellular interactions using newly available datasets and tools, including glycan biosynthesis models, omics datasets, and systems-level analyses. Thus, emerging tools in big data analytics and systems biology are facilitating novel insights on glycans and their relationship with multicellular behavior.

Insights into Bioinformatic Applications for Glycosylation: Instigating an Awakening towards Applying Glycoinformatic Resources for Cancer Diagnosis and Therapy

Glycosylation plays a crucial role in various diseases and their etiology. This has led to a clear understanding on the functions of carbohydrates in cell communication, which eventually will result in novel therapeutic approaches for treatment of various disease. Glycomics has now become one among the top ten technologies that will change the future. The direct implication of glycosylation as a hallmark of cancer and for cancer therapy is well established. As in proteomics, where bioinformatics tools have led to revolutionary achievements, bioinformatics resources for glycosylation have improved its practical implication. Bioinformatics tools, algorithms and databases are a mandatory requirement to manage and successfully analyze large amount of glycobiological data generated from glycosylation studies. This review consolidates all the available tools and their applications in glycosylation research. The achievements made through the use of bioinformatics into glycosylation studies are also presented. The importance of glycosylation in cancer diagnosis and therapy is discussed and the gap in the application of widely available glyco-informatic tools for cancer research is highlighted. This review is expected to bring an awakening amongst glyco-informaticians as well as cancer biologists to bridge this gap, to exploit the available glyco-informatic tools for cancer.

A consensus-based and readable extension of Linear Code for Reaction Rules (LiCoRR)

Systems glycobiology aims to provide models and analysis tools that account for the biosynthesis, regulation, and interactions with glycoconjugates. To facilitate these methods, there is a need for a clear glycan representation accessible to both computers and humans. Linear Code, a linearized and readily parsable glycan structure representation, is such a language. For this reason, Linear Code was adapted to represent reaction rules, but the syntax has drifted from its original description to accommodate new and originally unforeseen challenges. Here, we delineate the consensuses and inconsistencies that have arisen through this adaptation. We recommend options for a consensus-based extension of Linear Code that can be used for reaction rule specification going forward. Through this extension and specification of Linear Code to reaction rules, we aim to minimize inconsistent symbology thereby making glycan database queries easier. With a clear guide for generating reaction rule descriptions, glycan synthesis models will be more interoperable and reproducible thereby moving glycoinformatics closer to compliance with FAIR standards. Here, we present Linear Code for Reaction Rules (LiCoRR), version 1.0, an unambiguous representation for describing glycosylation reactions in both literature and code.

Golgi compartments enable controlled biomolecular assembly using promiscuous enzymes

The synthesis of eukaryotic glycans - branched sugar oligomers attached to cell-surface proteins and lipids - is organized like a factory assembly line. Specific enzymes within successive compartments of the Golgi apparatus determine where new monomer building blocks are linked to the growing oligomer. These enzymes act promiscuously and stochastically, causing microheterogeneity (molecule-to-molecule variability) in the final oligomer products. However, this variability is tightly controlled: a given eukaryotic protein type is typically associated with a narrow, specific glycan oligomer profile. Here, we use ideas from the mathematical theory of self-assembly to enumerate the enzymatic causes of oligomer variability and show how to eliminate each cause. We rigorously demonstrate that cells can specifically synthesize a larger repertoire of glycan oligomers by partitioning promiscuous enzymes across multiple Golgi compartments. This places limits on biomolecular assembly: glycan microheterogeneity becomes unavoidable when the number of compartments is limited, or enzymes are excessively promiscuous.

Modeling Glycan Processing Reveals Golgi-Enzyme Homeostasis upon Trafficking Defects and Cellular Differentiation

The decoration of proteins by carbohydrates is essential for eukaryotic life yet heterogeneous due to a lack of biosynthetic templates. This complex carbohydrate mixture—the glycan profile—is generated in the compartmentalized Golgi, in which level and localization of glycosylation enzymes are key determinants. Here, we develop and validate a computational model for glycan biosynthesis to probe how the biosynthetic machinery creates different glycan profiles. We combined stochastic modeling with Bayesian fitting that enables rigorous comparison to experimental data despite starting with uncertain initial parameters. This is an important development in the field of glycan modeling, which revealed biological insights about the glycosylation machinery in altered cellular states. We experimentally validated changes in N-linked glycan-modifying enzymes in cells with perturbed intra-Golgi-enzyme sorting and the predicted glycan-branching activity during osteogenesis. Our model can provide detailed information on altered biosynthetic paths, with potential for advancing treatments for glycosylation-related diseases and glyco-engineering of cells.

•

Developed a stochastic model of N-glycosylation coupled with Bayesian fitting

•

Validated predicted changes of Golgi organization in trafficking mutants

•

Model pinpointed functionally relevant glycan alterations in osteogenesis

A Systematic Review on the Implications of O-linked Glycan Branching and Truncating Enzymes on Cancer Progression and Metastasis

Glycosylation is the most commonly occurring post-translational modifications, and is believed to modify over 50% of all proteins. The process of glycan modification is directed by different glycosyltransferases, depending on the cell in which it is expressed. These small carbohydrate molecules consist of multiple glycan families that facilitate cell-cell interactions, protein interactions, and downstream signaling. An alteration of several types of O-glycan core structures have been implicated in multiple cancers, largely due to differential glycosyltransferase expression or activity. Consequently, aberrant O-linked glycosylation has been extensively demonstrated to affect biological function and protein integrity that directly result in cancer growth and progression of several diseases. Herein, we provide a comprehensive review of several initiating enzymes involved in the synthesis of O-linked glycosylation that significantly contribute to a number of different cancers.

A Markov model of glycosylation elucidates isozyme specificity and glycosyltransferase interactions for glycoengineering

Glycosylated biopharmaceuticals are important in the global pharmaceutical market. Despite the importance of their glycan structures, our limited knowledge of the glycosylation machinery still hinders controllability of this critical quality attribute. To facilitate discovery of glycosyltransferase specificity and predict glycoengineering efforts, here we extend the approach to model N-linked protein glycosylation as a Markov process. Our model leverages putative glycosyltransferase (GT) specificity to define the biosynthetic pathways for all measured glycans, and the Markov chain modelling is used to learn glycosyltransferase isoform activities and predict glycosylation following glycosyltransferase knock-in/knockout. We apply our methodology to four different glycoengineered therapeutics (i.e., Rituximab, erythropoietin, Enbrel, and alpha-1 antitrypsin) produced in CHO cells. Our model accurately predicted N-linked glycosylation following glycoengineering and further quantified the impact of glycosyltransferase mutations on reactions catalyzed by other glycosyltransferases. By applying these learned GT-GT interaction rules identified from single glycosyltransferase mutants, our model further predicts the outcome of multi-gene glycosyltransferase mutations on the diverse biotherapeutics. Thus, this modeling approach enables rational glycoengineering and the elucidation of relationships between glycosyltransferases, thereby facilitating biopharmaceutical research and aiding the broader study of glycosylation to elucidate the genetic basis of complex changes in glycosylation.

CHO-Omics Review: The Impact of Current and Emerging Technologies on Chinese Hamster Ovary Based Bioproduction

CHO cells are the most prevalent platform for modern bio-therapeutic production. Currently, there are several CHO cell lines used in bioproduction with distinct characteristics and unique genotypes and phenotypes. These differences limit advances in productivity and quality that can be achieved by the most common approaches to bioprocess optimization and cell line engineering. Incorporating omics-based approaches into current bioproduction processes will complement traditional methodologies to maximize gains from CHO engineering and bioprocess improvements. In order to highlight the utility of omics technologies in CHO bioproduction, the authors discuss current applications as well as limitations of genomics, transcriptomics, proteomics, metabolomics, lipidomics, fluxomics, glycomics, and multi-omics approaches and the potential they hold for the future of bioproduction. Multiple omics approaches are currently being used to improve CHO bioprocesses; however, the application of these technologies is still limited. As more CHO-omic datasets become available and integrated into systems models, the authors expect significant gains in product yield and quality. While individual omics technologies provide incremental improvements in bioproduction, the authors will likely see the most significant gains by applying multi-omics and systems biology approaches to individual CHO cell lines.

Model-based analysis of N-glycosylation in Chinese hamster ovary cells

The Chinese hamster ovary (CHO) cell is the gold standard for manufacturing of glycosylated recombinant proteins for production of biotherapeutics. The similarity of its glycosylation patterns to the human versions enable the products of this cell line favorable pharmacokinetic properties and lower likelihood of causing immunogenic responses. Because glycan structures are the product of the concerted action of intracellular enzymes, it is difficult to predict a priori how the effects of genetic manipulations alter glycan structures of cells and therapeutic properties. For that reason, quantitative models able to predict glycosylation have emerged as promising tools to deal with the complexity of glycosylation processing. For example, an earlier version of the same model used in this study was used by others to successfully predict changes in enzyme activities that could produce a desired change in glycan structure. In this study we utilize an updated version of this model to provide a comprehensive analysis of N-glycosylation in ten Chinese hamster ovary (CHO) cell lines that include a wild type parent and nine mutants of CHO, through interpretation of previously published mass spectrometry data. The updated N-glycosylation mathematical model contains up to 50,605 glycan structures. Adjusting the enzyme activities in this model to match N-glycan mass spectra produces detailed predictions of the glycosylation process, enzyme activity profiles and complete glycosylation profiles of each of the cell lines. These profiles are consistent with biochemical and genetic data reported previously. The model-based results also predict glycosylation features of the cell lines not previously published, indicating more complex changes in glycosylation enzyme activities than just those resulting directly from gene mutations. The model predicts that the CHO cell lines possess regulatory mechanisms that allow them to adjust glycosylation enzyme activities to mitigate side effects of the primary loss or gain of glycosylation function known to exist in these mutant cell lines. Quantitative models of CHO cell glycosylation have the potential for predicting how glycoengineering manipulations might affect glycoform distributions to improve the therapeutic performance of glycoprotein products.

What can mathematical modelling say about CHO metabolism and protein glycosylation?

Chinese hamster ovary cells have been in the spotlight for process optimization in recent years, due to being the major, long established cell factory for the production of recombinant proteins. A deep, quantitative understanding of CHO metabolism and mechanisms involved in protein glycosylation has proven to be attainable through the development of high throughput technologies. Here we review the most notable accomplishments in the field of modelling CHO metabolism and protein glycosylation.

Multi-level regulation of cellular glycosylation: from genes to transcript to enzyme to structure

Glycosylation is a ubiquitous mammalian post-translational modification that both decorates a majority of expressed proteins and regulates their function. Cellular glycan biosynthesis is facilitated by a few hundred enzymes that are collectively termed 'glycoenzymes'. The expression and activity of these enzymes is controlled at the transcription, translation and post-translation levels. New wet-lab advances are providing analytical methods to collect large-scale data at these multiple levels, relational databases are starting to collate these results, and computer models are beginning to integrate this information across scales in order to gain new knowledge. These activities are likely to enable the qualitative and quantitative mapping of pathways regulating glycan production and function in proteins, cells and tissue.

A systematic framework to derive N-glycan biosynthesis process and the automated construction of glycosylation networks

Background

Abnormalities in glycan biosynthesis have been conclusively related to various diseases, whereas the complexity of the glycosylation process has impeded the quantitative analysis of biochemical experimental data for the identification of glycoforms contributing to disease. To overcome this limitation, the automatic construction of glycosylation reaction networks in silico is a critical step.

Results

In this paper, a framework K2014 is developed to automatically construct N-glycosylation networks in MATLAB with the involvement of the 27 most-known enzyme reaction rules of 22 enzymes, as an extension of previous model KB2005. A toolbox named Glycosylation Network Analysis Toolbox (GNAT) is applied to define network properties systematically, including linkages, stereochemical specificity and reaction conditions of enzymes. Our network shows a strong ability to predict a wider range of glycans produced by the enzymes encountered in the Golgi Apparatus in human cell expression systems.

Conclusions

Our results demonstrate a better understanding of the underlying glycosylation process and the potential of systems glycobiology tools for analyzing conventional biochemical or mass spectrometry-based experimental data quantitatively in a more realistic and practical way.

A Knowledge-Based System for Display and Prediction of O-Glycosylation Network Behaviour in Response to Enzyme Knockouts

O-linked glycosylation is an important post-translational modification of mucin-type protein, changes to which are important biomarkers of cancer. For this study of the enzymes of O-glycosylation, we developed a shorthand notation for representing GalNAc-linked oligosaccharides, a method for their graphical interpretation, and a pattern-matching algorithm that generates networks of enzyme-catalysed reactions. Software for generating glycans from the enzyme activities is presented, and is also available online. The degree distributions of the resulting enzyme-reaction networks were found to be Poisson in nature. Simple graph-theoretic measures were used to characterise the resulting reaction networks. From a study of in-silico single-enzyme knockouts of each of 25 enzymes known to be involved in mucin O-glycan biosynthesis, six of them, β-1,4-galactosyltransferase (β4Gal-T4), four glycosyltransferases and one sulfotransferase, play the dominant role in determining O-glycan heterogeneity. In the absence of β4Gal-T4, all Lewis X, sialyl-Lewis X, Lewis Y and Sda/Cad glycoforms were eliminated, in contrast to knockouts of the N-acetylglucosaminyltransferases, which did not affect the relative abundances of O-glycans expressing these epitopes. A set of 244 experimentally determined mucin-type O-glycans obtained from the literature was used to validate the method, which was able to predict up to 98% of the most common structures obtained from human and engineered CHO cell glycoforms.

A systematic analysis of acceptor specificity and reaction kinetics of five human α(2,3)sialyltransferases: Product inhibition studies illustrate reaction mechanism for ST3Gal-I

Sialyltransferases (STs) catalyze the addition of sialic acids to the non-reducing ends of glycoproteins and glycolipids. In this work, we examined the acceptor specificity of five human α(2,3)sialyltransferases, namely ST3Gal -I, -II, -III, -IV and -VI. KM values for each of these enzymes is presented using radioactivity for acceptors containing Type-I (Galβ1,3GlcNAc), Type-II (Galβ1,4GlcNAc), Type-III (Galβ1,3GalNAc) and Core-2 (Galβ1,3(GlcNAcβ1,6)GalNAc) reactive groups. Several variants of acceptors inhibited ST3Gal activity emphasizing structural role of acceptor in enzyme-catalyzed reactions. In some cases, mass spectrometry was performed for structural verification. The results demonstrate human ST3Gal-I catalysis towards Type-III and Core-2 acceptors with KM = 5-50 μM and high VMax values. The KM for ST3Gal-I and ST3Gal-II was 100 and 30-fold lower, respectively, for Type-III compared to Type-I acceptors. Variants of Type-I and Type-II structures characterized ST3Gal-III, -IV and -VI for their catalytic specificity. This manuscript also estimates KM for human ST3Gal-VI using Type-I and Type-II substrates. Together, these findings built a platform for designing inhibitors of STs having therapeutic potential.

Integration of systems glycobiology with bioinformatics toolboxes, glycoinformatics resources, and glycoproteomics data

The glycome constitutes the entire complement of free carbohydrates and glycoconjugates expressed on whole cells or tissues. 'Systems Glycobiology' is an emerging discipline that aims to quantitatively describe and analyse the glycome. Here, instead of developing a detailed understanding of single biochemical processes, a combination of computational and experimental tools are used to seek an integrated or 'systems-level' view. This can explain how multiple biochemical reactions and transport processes interact with each other to control glycome biosynthesis and function. Computational methods in this field commonly build in silico reaction network models to describe experimental data derived from structural studies that measure cell-surface glycan distribution. While considerable progress has been made, several challenges remain due to the complex and heterogeneous nature of this post-translational modification. First, for the in silico models to be standardized and shared among laboratories, it is necessary to integrate glycan structure information and glycosylation-related enzyme definitions into the mathematical models. Second, as glycoinformatics resources grow, it would be attractive to utilize 'Big Data' stored in these repositories for model construction and validation. Third, while the technology for profiling the glycome at the whole-cell level has been standardized, there is a need to integrate mass spectrometry derived site-specific glycosylation data into the models. The current review discusses progress that is being made to resolve the above bottlenecks. The focus is on how computational models can bridge the gap between 'data' generated in wet-laboratory studies with 'knowledge' that can enhance our understanding of the glycome.

Systems glycobiology for glycoengineering

Glycosylation serves essential functions on many proteins produced in biopharmaceutical manufacturing, making it mandatory to thoroughly consider its biogenesis during the production process. Glycoengineering efforts involve the rational design of glycosylation through adjustments in culturing conditions or genetic modifications. Computational models have been developed to aid this process, aiming to offer cheaper and faster alternatives to costly screening strategies. Recently, these models have been successfully utilized to predict glycosylation of products of industrial relevance. Furthermore, systems-level analyses of glycan diversity are elucidating deeper insights into the mechanisms underlying glycosylation. As computational models of glycosylation continue to be expanded, refined, and leveraged for detailed analysis of glycomics data, they will become invaluable resources for cell line development and glycoengineering.

Enhancing MS(n) mass spectrometry strategy for carbohydrate analysis: A b2 ion spectral library

Searchable mass spectral libraries for glycans may be enhanced using a B2 ion library. Using a quadrupole ion-trap mass spectrometer, successive fragmentations of sodiated oligosaccharides were carried out in the positive ion mode. In B,Y-type fragmentation, disaccharide B2 ions are generated which correspond to specific glycosidic linkages using progressive MS stages. Fragmentation of "B2 ions" corresponding to glycosidic linkages such as Hex-Fuc, Hex-Hex, Hex-HexNAc, HexNAc-Hex and HexNAc-HexNAc, were systematically studied in low energy CID and collected to form a "B2 library". Linkages produce characteristic fragmentation patterns in the absence of cross-ring fragmentation. Patterns of "B2 ions" rely on relative stability of glycosidic bonds and carbohydrate-metal complexes in the gas phase. MS(n) studies of linear, branched trisaccharides and tetrasaccharides show that isomers for which B2 ion information is not available are rarely a problem in practice by their absence in an isomeric sequence or by their scarcity in nature. This MS strategy for linkage determination of carbohydrates aided by a "B2 library" was developed with a scope for expansion, providing an improved tool for glycomics. We validated this method examining levels of expressed activities of two glycosyl transferases in cancer cell lines: β3(B3GALNT2) and β4GalNAcT(B4GALNT3&4) that generate GalNAcβ3GlcNAcβ and GalNAcβ4GlcNAcβ.

BIOLOGICAL SIGNIFICANCE

Glycosylation is an important class of the "postranslationome", which includes manifold aspects of post-translational protein modification, affecting protein conformation, providing ligands for protein receptors [1-5], and encoding unique haptenic [6,7] or antigenic markers for oncology [8-11] and other applications. Identification of individual monomeric units, linkages, ring size, branching and anomerity has posed significant challenges to mass spectrometrists. MS(n) is a growing key instrumental method to differentiate among isomers [12]. While the potential isomers in oligosaccharides are impossibly large [12], likely possibilities can be limited by the biological system, including the expressed glycosyl transferases [13-20]. Mass spectra from sequential stages of collision activation (MS(n)) can supply structural details for precise characterization of linkage, monomer ID, substitutions, anomerity and branching [21-25]. There is a fundamental need for high throughput tools in glycomics to complement proteome studies. In that regard, nothing could be more important than searchable spectral library files for structural confirmation. The National Academy of Science (NAS) report (http://glyco.nas.edu) recommends the need of more than 10,000 synthetic structures of carbohydrates to advance the field of glycomics. This study demonstrates that the general reproducibility of ion trap spectra, and energy independence from modes of ionization and collisional activation, make compiling an MS(n) library for carbohydrate identification an achievable research target [26]. We intend to use the new B2 library for carbohydrate differences found on cancers, where we profile the glycosyltransferases to predict classes of potential structures, and use the library for MS identification of the expected cohort of altered structures.

A computational framework for the automated construction of glycosylation reaction networks

Glycosylation is among the most common and complex post-translational modifications identified to date. It proceeds through the catalytic action of multiple enzyme families that include the glycosyltransferases that add monosaccharides to growing glycans, and glycosidases which remove sugar residues to trim glycans. The expression level and specificity of these enzymes, in part, regulate the glycan distribution or glycome of specific cell/tissue systems. Currently, there is no systematic method to describe the enzymes and cellular reaction networks that catalyze glycosylation. To address this limitation, we present a streamlined machine-readable definition for the glycosylating enzymes and additional methodologies to construct and analyze glycosylation reaction networks. In this computational framework, the enzyme class is systematically designed to store detailed specificity data such as enzymatic functional group, linkage and substrate specificity. The new classes and their associated functions enable both single-reaction inference and automated full network reconstruction, when given a list of reactants and/or products along with the enzymes present in the system. In addition, graph theory is used to support functions that map the connectivity between two or more species in a network, and that generate subset models to identify rate-limiting steps regulating glycan biosynthesis. Finally, this framework allows the synthesis of biochemical reaction networks using mass spectrometry (MS) data. The features described above are illustrated using three case studies that examine: i) O-linked glycan biosynthesis during the construction of functional selectin-ligands; ii) automated N-linked glycosylation pathway construction; and iii) the handling and analysis of glycomics based MS data. Overall, the new computational framework enables automated glycosylation network model construction and analysis by integrating knowledge of glycan structure and enzyme biochemistry. All the implemented features are provided as part of the Glycosylation Network Analysis Toolbox (GNAT), an open-source, platform-independent, MATLAB based toolbox for studies of Systems Glycobiology.

Glycosylation Network Analysis Toolbox: a MATLAB-based environment for systems glycobiology

Systems glycobiology studies the interaction of various pathways that regulate glycan biosynthesis and function. Software tools for the construction and analysis of such pathways are not yet available. We present GNAT, a platform-independent, user-extensible MATLAB-based toolbox that provides an integrated computational environment to construct, manipulate and simulate glycans and their networks. It enables integration of XML-based glycan structure data into SBML (Systems Biology Markup Language) files that describe glycosylation reaction networks. Curation and manipulation of networks is facilitated using class definitions and glycomics database query tools. High quality visualization of networks and their steady-state and dynamic simulation are also supported.

AVAILABILITY

The software package including source code, help documentation and demonstrations are available at http://sourceforge.net/projects/gnatmatlab/files/.

Understanding glycomechanics using mathematical modeling: a review of current approaches to simulate cellular glycosylation reaction networks

Following the footsteps of genomics and proteomics, recent years have witnessed the growth of large-scale experimental methods in the field of glycomics. In parallel, there has also been growing interest in developing Systems Biology based methods to study the glycome. The combined goals of these endeavors is to identify glycosylation-dependent mechanisms regulating human physiology, check points that can control the progression of pathophysiology, and modifications to reaction pathways that can result in more uniform biopharmaceutical processes. In these efforts, mathematical models of N- and O-linked glycosylation have emerged as paradigms for the field. While these are relatively few in number, nevertheless, the existing models provide a basic framework that can be used to develop more sophisticated analysis strategies for glycosylation in the future. The current review surveys these computational models with focus on the underlying mathematics and assumptions, and with respect to their ability to generate experimentally testable hypotheses.

The engagement of selectins and their ligands in colorectal cancer liver metastases

The colonization of the liver by colorectal cancer (CRC) cells is a complicated process which includes many stages, until macrometastases occur. The entrapment of malignant cells within the hepatic sinusoids and their interactions with resident non-parenchymal cells are considered very important for the whole metastatic sequence. In the sinusoids, cell connection and signalling is mediated by multiple cell adhesion molecules, such as the selectins. The three members of the selectin family, E-, P- and L-selectin, in conjunction with sialylated Lewis ligands and CD44 variants, regulate colorectal cell communication and adhesion with platelets, leucocytes, sinusoidal endothelial cells and stellate cells. Their role in CRC liver metastases has been investigated in animal models and human tissue, in vivo and in vitro, in static and shear flow conditions, and their key-function in several molecular pathways has been displayed. Therefore, trials have already commenced aiming to exploit selectins and their ligands in the treatment of benign and malignant diseases. Multiple pharmacological agents have been developed that are being tested for potential therapeutic applications.

Module-based multiscale simulation of angiogenesis in skeletal muscle

Background

Mathematical modeling of angiogenesis has been gaining momentum as a means to shed new light on the biological complexity underlying blood vessel growth. A variety of computational models have been developed, each focusing on different aspects of the angiogenesis process and occurring at different biological scales, ranging from the molecular to the tissue levels. Integration of models at different scales is a challenging and currently unsolved problem.

Results

We present an object-oriented module-based computational integration strategy to build a multiscale model of angiogenesis that links currently available models. As an example case, we use this approach to integrate modules representing microvascular blood flow, oxygen transport, vascular endothelial growth factor transport and endothelial cell behavior (sensing, migration and proliferation). Modeling methodologies in these modules include algebraic equations, partial differential equations and agent-based models with complex logical rules. We apply this integrated model to simulate exercise-induced angiogenesis in skeletal muscle. The simulation results compare capillary growth patterns between different exercise conditions for a single bout of exercise. Results demonstrate how the computational infrastructure can effectively integrate multiple modules by coordinating their connectivity and data exchange. Model parameterization offers simulation flexibility and a platform for performing sensitivity analysis.

Conclusions

This systems biology strategy can be applied to larger scale integration of computational models of angiogenesis in skeletal muscle, or other complex processes in other tissues under physiological and pathological conditions.

Systems glycobiology: biochemical reaction networks regulating glycan structure and function

There is a growing use of bioinformatics based methods in the field of Glycobiology. These have been used largely to curate glycan structures, organize array-based experimental data and display existing knowledge of glycosylation-related pathways in silico. Although the cataloging of vast amounts of data is beneficial, it is often a challenge to gain meaningful mechanistic insight from this exercise alone. The development of specific analysis tools to query the database is necessary. If these queries can integrate existing knowledge of glycobiology, new insights may be gained. Such queries that couple biochemical knowledge and mathematics have been developed in the field of Systems Biology. The current review summarizes the current state of the art in the application of computational modeling in the field of Glycobiology. It provides (i) an overview of experimental and online resources that can be used to construct glycosylation reaction networks, (ii) mathematical methods to formulate the problem including a description of ordinary differential equation and logic-based reaction networks, (iii) optimization techniques that can be applied to fit experimental data for the purpose of model reconstruction and for evaluating unknown model parameters, (iv) post-simulation analysis methods that yield experimentally testable hypotheses and (v) a summary of available software tools that can be used by non-specialists to perform many of the above functions.

Protein glycosylation control in mammalian cell culture: past precedents and contemporary prospects

Protein glycosylation is a post-translational modification of paramount importance for the function, immunogenicity, and efficacy of recombinant glycoprotein therapeutics. Within the repertoire of post-translational modifications, glycosylation stands out as having the most significant proven role towards affecting pharmacokinetics and protein physiochemical characteristics. In mammalian cell culture, the understanding and controllability of the glycosylation metabolic pathway has achieved numerous successes. However, there is still much that we do not know about the regulation of the pathway. One of the frequent conclusions regarding protein glycosylation control is that it needs to be studied on a case-by-case basis since there are often conflicting results with respect to a control variable and the resulting glycosylation. In attempts to obtain a more multivariate interpretation of these potentially controlling variables, gene expression analysis and systems biology have been used to study protein glycosylation in mammalian cell culture. Gene expression analysis has provided information on how glycosylation pathway genes both respond to culture environmental cues, and potentially facilitate changes in the final glycoform profile. Systems biology has allowed researchers to model the pathway as well-defined, inter-connected systems, allowing for the in silico testing of pathway parameters that would be difficult to test experimentally. Both approaches have facilitated a macroscopic and microscopic perspective on protein glycosylation control. These tools have and will continue to enhance our understanding and capability of producing optimal glycoform profiles on a consistent basis.

Genome-scale metabolic reconstruction and in silico analysis of methylotrophic yeast Pichia pastoris for strain improvement

Background

Pichia pastoris has been recognized as an effective host for recombinant protein production. A number of studies have been reported for improving this expression system. However, its physiology and cellular metabolism still remained largely uncharacterized. Thus, it is highly desirable to establish a systems biotechnological framework, in which a comprehensive in silico model of P. pastoris can be employed together with high throughput experimental data analysis, for better understanding of the methylotrophic yeast's metabolism.

Results

A fully compartmentalized metabolic model of P. pastoris (iPP668), composed of 1,361 reactions and 1,177 metabolites, was reconstructed based on its genome annotation and biochemical information. The constraints-based flux analysis was then used to predict achievable growth rate which is consistent with the cellular phenotype of P. pastoris observed during chemostat experiments. Subsequent in silico analysis further explored the effect of various carbon sources on cell growth, revealing sorbitol as a promising candidate for culturing recombinant P. pastoris strains producing heterologous proteins. Interestingly, methanol consumption yields a high regeneration rate of reducing equivalents which is substantial for the synthesis of valuable pharmaceutical precursors. Hence, as a case study, we examined the applicability of P. pastoris system to whole-cell biotransformation and also identified relevant metabolic engineering targets that have been experimentally verified.

Conclusion

The genome-scale metabolic model characterizes the cellular physiology of P. pastoris, thus allowing us to gain valuable insights into the metabolism of methylotrophic yeast and devise possible strategies for strain improvement through in silico simulations. This computational approach, combined with synthetic biology techniques, potentially forms a basis for rational analysis and design of P. pastoris metabolic network to enhance humanized glycoprotein production.

Functional network of glycan-related molecules: glyco-net in glycoconjugate data bank

Background

Glycans are involved in a wide range of biological process, and they play an essential role in functions such as cell differentiation, cell adhesion, pathogen-host recognition, toxin-receptor interactions, signal transduction, cancer metastasis, and immune responses. Elucidating pathways related to post-translational modifications (PTMs) such as glycosylation are of growing importance in post-genome science and technology. Graphical networks describing the relationships among glycan-related molecules, including genes, proteins, lipids and various biological events are considered extremely valuable and convenient tools for the systematic investigation of PTMs. However, there is no database which dynamically draws functional networks related to glycans.

Description

We have created a database called Glyco-Net http://www.glycoconjugate.jp/functions/, with many binary relationships among glycan-related molecules. Using search results, we can dynamically draw figures of the functional relationships among these components with nodes and arrows. A certain molecule or event corresponds to a node in the network figures, and the relationship between the molecule and the event are indicated by arrows. Since all components are treated equally, an arrow is also a node.

Conclusions

In this paper, we describe our new database, Glyco-Net, which is the first database to dynamically show networks of the functional profiles of glycan related molecules. The graphical networks will assist in the understanding of the role of the PTMs. In addition, since various kinds of bio-objects such as genes, proteins, and inhibitors are equally treated in Glyco-Net, we can obtain a large amount of information on the PTMs.

Understanding inhibition of viral proteins on type I IFN signaling pathways with modeling and optimization

The interferon system provides a powerful and universal intracellular defense mechanism against viruses. As one part of their survival strategies, many viruses have evolved mechanisms to counteract the host type I interferon (IFN-alpha/beta) responses. In this study, we attempt to investigate virus- and double-strand RNA (dsRNA)-triggered type I IFN signaling pathways and understand the inhibition of IFN-alpha/beta induction by viral proteins using mathematical modeling and quantitative analysis. Based on available literature and our experimental data, we develop a mathematical model of virus- and dsRNA-triggered signaling pathways leading to type I IFN gene expression during the primary response, and use the genetic algorithm to optimize all rate constants in the model. The consistency between numerical simulation results and biological experimental data demonstrates that our model is reasonable. Further, we use the model to predict the following phenomena: (1) the dose-dependent inhibition by classical swine fever virus (CSFV) N(pro) or E(rns) protein is observed at a low dose and can reach a saturation above a certain dose, not an increase; (2) E(rns) and N(pro) have no synergic inhibitory effects on IFN-beta induction; (3) the different characters in an important transcription factor, phosphorylated IRF3 (IRF3p), are exhibited because N(pro) or E(rns) counteracted dsRNA- and virus-triggered IFN-beta induction by targeting the different molecules in the signaling pathways and (4) N(pro) inhibits the IFN-beta expression not only by interacting with IFR3 but also by affecting its complex with MITA. Our approaches help to gain insight into system properties and rational therapy design, as well as to generate hypotheses for further research.

Genome-scale modeling and in silico analysis of mouse cell metabolic network

Genome-scale metabolic modeling has been successfully applied to a multitude of microbial systems, thus improving our understanding of their cellular metabolisms. Nevertheless, only a handful of works have been done for describing mammalian cells, particularly mouse, which is one of the important model organisms, providing various opportunities for both biomedical research and biotechnological applications. Presented herein is a genome-scale mouse metabolic model that was systematically reconstructed by improving and expanding the previous generic model based on integrated biochemical and genomic data of Mus musculus. The key features of the updated model include additional information on gene-protein-reaction association, and improved network connectivity through lipid, amino acid, carbohydrate and nucleotide biosynthetic pathways. After examining the model predictability both quantitatively and qualitatively using constraints-based flux analysis, the structural and functional characteristics of the mouse metabolism were investigated by evaluating network statistics/centrality, gene/metabolite essentiality and their correlation. The results revealed that overall mouse metabolic network is topologically dominated by highly connected and bridging metabolites, and functionally by lipid metabolism that most of essential genes and metabolites are from. The current in silico mouse model can be exploited for understanding and characterizing the cellular physiology, identifying potential cell engineering targets for the enhanced production of recombinant proteins and developing diseased state models for drug targeting.