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);...

Parkinson’s disease is a slowly progressive neurodegenerative disease characterised by dysfunction and death of selectively vulnerable midbrain dopaminergic neurons and the development of human in vitro cellular models of the disease is a major challenge in Parkinson’s disease research. We constructed an automated cell culture platform optimised for long-term maintenance and monitoring of different cells in three dimensional microfluidic cell culture devices. The system can be flexibly adapted to various experimental protocols and features time-lapse imaging microscopy for quality control and electrophysiology monitoring to assess cellular activity. Using this system, we continuously monitored the differentiation of Parkinson’s disease patient derived human neuroepithelial stem cells into midbrain specific dopaminergic neurons. Calcium imaging confirmed the electrophysiological activity of differentiated neurons and immunostaining confirmed the efficiency of the differentiation protocol. This system is the first example of an automated Organ-on-a-Chip culture and has the potential to enable a versatile array of in vitro experiments for patient-specific disease modelling.

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.

A specific set of brainstem nuclei are susceptible to degeneration in Parkinson’s disease. We hypothesise that neuronal vulnerability reflects shared phenotypic characteristics that confer selective vulnerability to degeneration. Neuronal phenotypic specification is mainly the cumulative result of a transcriptional regulatory program that is active during the development. By manual curation of the developmental biology literature, we comprehensively reconstructed an anatomically resolved cellular developmental lineage for the adult neurons in five brainstem regions that are selectively vulnerable to degeneration in prodromal or early Parkinson’s disease. We synthesised the literature on transcription factors that are required to be active, or required to be inactive, in the development of each of these five brainstem regions, and at least two differentially vulnerable nuclei within each region. Certain transcription factors, e.g., Ascl1 and Lmx1b, seem to be required for specification of many brainstem regions that are susceptible to degeneration in early Parkinson’s disease. Some transcription factors can even distinguish between differentially vulnerable nuclei within the same brain region, e.g., Pitx3 is required for specification of the substantia nigra pars compacta, but not the ventral tegmental area. We do not suggest that Parkinson’s disease is a developmental disorder. In contrast, we consider identification of shared developmental trajectories as part of a broader effort to identify the molecular mechanisms that underlie the phenotypic features that are shared by selectively vulnerable neurons. Systematic in vivo assessment of fate determining transcription factors should be completed for all neuronal populations vulnerable to degeneration in early Parkinson’s disease.

One of the hallmarks of sporadic Parkinson's disease is degeneration of dopaminergic neurons in the pars compacta of the substantia nigra. The aetiopathogenesis of this degeneration is still not fully understood, with dysfunction of many biochemical pathways in different subsystems suggested to be involved. Recent advances in constraint-based modelling approaches hold great potential to systematically examine the relative contribution of dysfunction in disparate pathways to dopaminergic neuronal degeneration, but few studies have employed these methods in Parkinson's disease research. Therefore, this review outlines a framework for future constraint-based modelling of dopaminergic neuronal metabolism to decipher the multi-factorial mechanisms underlying the neuronal pathology of Parkinson's disease.