COnstraint-Based Reconstruction and Analysis (COBRA) provides a molecular mechanistic framework for integrative analysis of experimental data and quantitative prediction of physicochemically and biochemically feasible phenotypic states. The COBRA Toolbox is a comprehensive software suite of interoperable COBRA methods. It has found widespread applications in biology, biomedicine, and biotechnology because its functions can be flexibly combined to implement tailored COBRA protocols for any biochemical network. Version 3.0 includes new methods for quality controlled reconstruction, modelling, topological analysis, strain and experimental design, network visualisation as well as network integration of chemoinformatic, metabolomic, transcriptomic, proteomic, and thermochemical data. New multi-lingual code integration also enables an expansion in COBRA application scope via high-precision, high-performance, and nonlinear numerical optimisation solvers for multi-scale, multi-cellular and reaction kinetic modelling, respectively. This protocol can be adapted for the generation and analysis of a constraint-based model in a wide variety of molecular systems biology scenarios. This protocol is an update to the COBRA Toolbox 1.0 and 2.0. The COBRA Toolbox 3.0 provides an unparalleled depth of constraint-based reconstruction and analysis methods. ]); 61 | The MUST sets are the sets of reactions that must increase or decrease their flux in order to achieve the desired phenotype in the mutant strain. As shown in Figure 6, the first order MUST sets are MustU and MustL while second order MUST sets are denoted as MustUU, MustLL, and MustUL. After parameters and constraints are defined, the functions findMustL and findMustU are run to determine the mustU and mustL sets, respectively. Define an ID of the run with:Each time the MUST sets are determined, folders are generated to read inputs and store outputs, i.e., reports. These folders are located in the directory defined by the uniquely defined runID.62 | In order to find the first order MUST sets, constraints should be defined: >> constrOpt = struct('rxnList', {{'EX_gluc', 'R75', 'EX_suc'}}, 'values', [-100; 0; 155.5]); 63 | The first order MUST set MustL is determined by running: >> [mustLSet, pos_mustL] = findMustL(model, minFluxesW, maxFluxesW, ... 'constrOpt', constrOpt, 'runID', runID);If runID is set to 'TestoptForceL', a folder TestoptForceL is created, in which two additional folders InputsMustL and OutputsMustL are created. The InputsMustL folder contains all the inputs required to run the function findMustL, while the OutputsMustL folder contains the mustL set found and a report that summarises all the inputs and outputs. In order to maintain a chronological order of computational experiments, the report is timestamped.64 | Display the reactions that belong to the mustL set using: >> disp(mustLSet) 65 | The first order MUST set MustU is determined by running: >> [mustUSet, pos_mustU] = findMustU(model, minFluxesW, maxFluxesW, ... 'constrOpt', constrOpt, 'runID', runID);...

Patient-derived cellular models are a powerful approach to study human disease, especially neurodegenerative diseases, such as Parkinson's disease, where affected primary neurons, e.g., substantia nigra dopaminergic neurons, are almost inaccessible. Starting with a comprehensive generic reconstruction of human metabolism, Recon3D, we generated a high-quality, constraint-based, genome-scale, in silico model of human dopaminergic neuronal metabolism (iDopaNeuro1). It is a synthesis of extensive manual curation of the biochemical literature on neuronal metabolism, together with novel, quantitative, transcriptomic and targeted exometabolomic data from human stem cell-derived, midbrain-specific, dopaminergic neurons in vitro. Thermodynamic constraint-based modelling with iDopaNeuro1 is qualitatively accurate (92% correct) and quantitatively accurate (Spearman rank 0.7) at predicting metabolite secretion or uptake, given quantitative exometabolomic constraints on uptakes, or secretions, respectively. iDopaNeuro1 is also qualitatively accurate at predicting the consequences of metabolic perturbations, e.g., complex I inhibition (Spearman rank 0.69) in a manner consistent with literature on monogenic mitochondrial Parkinson's disease. The iDopaNeuro1 model provides a foundation for a quantitative systems biochemistry approach to metabolic dysfunction in Parkinson's disease. Moreover, the plethora of novel mathematical and computational approaches required to develop it are generalisable to study any other disease associated with metabolic dysfunction.

The endoplasmic reticulum (ER) is a central eukaryotic organelle with a tubular network made of hairpin proteins linked by hydrolysis of guanosine triphosphate nucleotides. Among posttranslational modifications initiated at the ER level, glycosylation is the most common reaction. However, our understanding of the impact of glycosylation on the ER structure remains unclear. Here, we show that exostosin-1 (EXT1) glycosyltransferase, an enzyme involved in N-glycosylation, is a key regulator of ER morphology and dynamics. We have integrated multiomics and superresolution imaging to characterize the broad effect of EXT1 inactivation, including the ER shape-dynamics-function relationships in mammalian cells. We have observed that inactivating EXT1 induces cell enlargement and enhances metabolic switches such as protein secretion. In particular, suppressing EXT1 in mouse thymocytes causes developmental dysfunctions associated with the ER network extension. Last, our data illuminate the physical and functional aspects of the ER proteome-glycome-lipidome structure axis, with implications in biotechnology and medicine.

SUMMARYThe endoplasmic reticulum (ER) is a central eukaryotic organelle with a tubular network made of hairpin proteins linked by hydrolysis of GTP nucleotides. Among post-translational modifications initiated at the ER level, glycosylation is the most common reaction. However, our understanding of the impact of glycosylation on ER structure remains unclear. Here, we show that Exostosin-1 (EXT1) glycosyltransferase, an enzyme involved in N-glycosylation, is a key regulator of ER morphology and dynamics. We have integrated multi-omics data and super-resolution imaging to characterize the broad effect of EXT1 inactivation, including ER shape-dynamics-function relationships in mammalian cells. We have observed that, inactivating EXT1 induces cell enlargement and enhances metabolic switches such as protein secretion. In particular, suppressing EXT1 in mouse thymocytes causes developmental dysfunctions associated to ER network extension. Our findings suggest that EXT1 drives glycosylation reactions involving ER structural proteins and high-energy nucleotide sugars, which might also apply to other organelles.

Metabolic models are reconstructed based on an organism's available genome annotation and provide predictive tools to study metabolic processes at a systemslevel. Genome-scale metabolic models may include gaps as well as reactions that are unverified experimentally. Reconstructed models of newly isolated microalgal species will result in weaknesses due to these gaps, as there is usually sparse biochemical evidence available for the metabolism of such isolates. The phenotype microarray (PM) technology is an effective, high-throughput method that functionally determines cellular metabolic activities in response to a wide array of entry metabolites. Combining the high throughput phenotypic assays with metabolic modeling can allow existing metabolic network models to be rapidly reconstructed or optimized by providing biochemical evidence to support and expand genomic evidence. This work will show the use of PM assays for the study of microalgae by using the green microalgal model species Chlamydomonas reinhardtii as an example. Experimental evidence for over 254 reactions obtained by PM was used in this study to expand and refine a genome-scale C. reinhardtii metabolic network model, iRC1080, by approximately 25 percent. The protocol created here can be used as a basis for functionally profiling the metabolism of other microalgae, including known microalgae mutants and new isolates.network models can guide the rational designs for the rapid development of optimization strategies 1 , 2 , 3 , 4 .Although approximately 160 microalgal species have been sequenced 5 , there are, to our knowledge, only 44 algal