Refine
Document Type
- Article (9)
- Doctoral Thesis (1)
- Postprint (1)
Language
- English (11)
Is part of the Bibliography
- yes (11) (remove)
Keywords
- constraint-based modeling (2)
- 13C labeling (1)
- Arabidopsis thaliana (1)
- B12-dependent 1,2-propanediol degradation (1)
- Genome-wide (1)
- Genomic selection (1)
- INST-MFA (1)
- Metabolic models (1)
- Single-nucleotide polymorphisms (1)
- algae (1)
Institute
The metabolic state of an organism reflects the entire phenotype that is jointly affected by genetic and environmental changes. Due to the complexity of metabolism, system-level modelling approaches have become indispensable tools to obtain new insights into biological functions. In particular, simulation and analysis of metabolic networks using constraint-based modelling approaches have helped the analysis of metabolic fluxes. However, despite ongoing improvements in prediction of reaction flux through a system, approaches to directly predict metabolite concentrations from large-scale metabolic networks remain elusive. In this thesis, we present a computational approach for inferring concentration ranges from genome-scale metabolic models endowed with mass action kinetics. The findings specify a molecular mechanism underling facile control of concentration ranges for components in large-scale metabolic networks. Most importantly, an extended version of the approach can be used to predict concentration ranges without knowledge of kinetic parameters, provided measurements of concentrations in a reference state. We show that the approach is applicable with large-scale kinetic and stoichiometric metabolic models of organisms from different kingdoms of life. By challenging the predictions of concentration ranges in the genome-scale metabolic network of Escherichia coli with real-world data sets, we further demonstrate the prediction power and limitations of the approach. To predict concentration ranges in other species, e.g. model plant species Arabidopsis thaliana, we would rely on estimates of kinetic parameters (i.e. enzyme catalytic rates) since plant-specific enzyme catalytic rates are poorly documented. Using the constraint-based approach of Davidi et al. for estimation of enzyme catalytic rates, we obtain values for 168 plant enzymes. The approach depends on quantitative proteomics data and flux estimates obtained from constraint-based model of plant leaf metabolism integrating maximal rates of selected enzymes, plant-specific constraints on fluxes through canonical pathways, and growth measurements from Arabidopsis thaliana rosette under ten conditions. We demonstrate a low degree of plant enzyme saturation, supported by the agreement between concentrations of nicotinamide adenine dinucleotide, adenosine triphosphate, and glyceraldehyde 3-phosphate, based on our maximal in vivo catalytic rates, and available quantitative metabolomics data. Hence, our results show genome-wide estimation for plant-specific enzyme catalytic rates is feasible. These can now be readily employed to study resource allocation, to predict enzyme and metabolite concentrations using recent constrained-based modelling approaches. Constraint-based methods do not directly account for kinetic mechanisms and corresponding parameters. Therefore, a number of workflows have already been proposed to approximate reaction kinetics and to parameterize genome-scale kinetic models. We present a systems biology strategy to build a fully parameterized large-scale model of Chlamydomonas reinhardtii accounting for microcompartmentalization in the chloroplast stroma. Eukaryotic algae comprise a microcompartment, the pyrenoid, essential for the carbon concentrating mechanism (CCM) that improves their photosynthetic performance. Since the experimental study of the effects of microcompartmentation on metabolic pathways is challenging, we employ our model to investigate compartmentation of fluxes through the Calvin-Benson cycle between pyrenoid and stroma. Our model predicts that ribulose-1,5-bisphosphate, the substrate of Rubisco, and 3-phosphoglycerate, its product, diffuse in and out of the pyrenoid. We also find that there is no major diffusional barrier to metabolic flux between the pyrenoid and stroma. Therefore, our computational approach represents a stepping stone towards understanding of microcompartmentalized CCM in other organisms. This thesis provides novel strategies to use genome-scale metabolic networks to predict and integrate metabolite concentrations. Therefore, the presented approaches represent an important step in broadening the applicability of large-scale metabolic models to a range of biotechnological and medical applications.
Characterization of maximal enzyme catalytic rates in central metabolism of Arabidopsis thaliana
(2020)
Availability of plant-specific enzyme kinetic data is scarce, limiting the predictive power of metabolic models and precluding identification of genetic factors of enzyme properties. Enzyme kinetic data are measuredin vitro, often under non-physiological conditions, and conclusions elicited from modeling warrant caution. Here we estimate maximalin vivocatalytic rates for 168 plant enzymes, including photosystems I and II, cytochrome-b6f complex, ATP-citrate synthase, sucrose-phosphate synthase as well as enzymes from amino acid synthesis with previously undocumented enzyme kinetic data in BRENDA. The estimations are obtained by integrating condition-specific quantitative proteomics data, maximal rates of selected enzymes, growth measurements fromArabidopsis thalianarosette with and fluxes through canonical pathways in a constraint-based model of leaf metabolism. In comparison to findings inEscherichia coli, we demonstrate weaker concordance between the plant-specificin vitroandin vivoenzyme catalytic rates due to a low degree of enzyme saturation. This is supported by the finding that concentrations of nicotinamide adenine dinucleotide (phosphate), adenosine triphosphate and uridine triphosphate, calculated based on our maximalin vivocatalytic rates, and available quantitative metabolomics data are below reportedKMvalues and, therefore, indicate undersaturation of respective enzymes. Our findings show that genome-wide profiling of enzyme kinetic properties is feasible in plants, paving the way for understanding resource allocation.
Understanding the complexity of metabolic networks has implications for manipulation of their functions. The complexity of metabolic networks can be characterized by identifying multireaction dependencies that are challenging to determine due to the sheer number of combinations to consider. Here, we propose the concept of concordant complexes that captures multireaction dependencies and can be efficiently determined from the algebraic structure and operational constraints of metabolic networks. The concordant complexes imply the existence of concordance modules based on which the apparent complexity of 12 metabolic networks of organisms from all kingdoms of life can be reduced by at least 78%. A comparative analysis against an ensemble of randomized metabolic networks shows that the metabolic network of Escherichia coli contains fewer concordance modules and is, therefore, more tightly coordinated than expected by chance. Together, our findings demonstrate that metabolic networks are considerably simpler than what can be perceived from their structure alone.
Successfully designed and implemented plant-specific synthetic metabolic pathways hold promise to increase crop yield and nutritional value. Advances in synthetic biology have already demonstrated the capacity to design artificial biological pathways whose behavior can be predicted and controlled in microbial systems. However, the transfer of these advances to model plants and crops faces the lack of characterization of plant cellular pathways and increased complexity due to compartmentalization and multicellularity. Modern computational developments provide the means to test the feasibility of plant synthetic metabolic pathways despite gaps in the accumulated knowledge of plant metabolism. Here, we provide a succinct systematic review of optimization-based and retrobiosynthesis approaches that can be used to design and in silico test synthetic metabolic pathways in large-scale plant context-specific metabolic models. In addition, by surveying the existing case studies, we highlight the challenges that these approaches face when applied to plants. Emphasis is placed on understanding the effect that metabolic designs can have on native metabolism, particularly with respect to metabolite concentrations and thermodynamics of biochemical reactions. In addition, we discuss the computational developments that may help to transform the identified challenges into opportunities for plant synthetic biology.
Cells and organelles are not homogeneous but include microcompartments that alter the spatiotemporal characteristics of cellular processes. The effects of microcompartmentation on metabolic pathways are however difficult to study experimentally. The pyrenoid is a microcompartment that is essential for a carbon concentrating mechanism (CCM) that improves the photosynthetic performance of eukaryotic algae. Using Chlamydomonas reinhardtii, we obtained experimental data on photosynthesis, metabolites, and proteins in CCM-induced and CCM-suppressed cells. We then employed a computational strategy to estimate how fluxes through the Calvin-Benson cycle are compartmented between the pyrenoid and the stroma. Our model predicts that ribulose-1,5-bisphosphate (RuBP), the substrate of Rubisco, and 3-phosphoglycerate (3PGA), its product, diffuse in and out of the pyrenoid, respectively, with higher fluxes in CCM-induced cells. It also indicates that there is no major diffusional barrier to metabolic flux between the pyrenoid and stroma. Our computational approach represents a stepping stone to understanding microcompartmentalized CCM in other organisms.
Cells and organelles are not homogeneous but include microcompartments that alter the spatiotemporal characteristics of cellular processes. The effects of microcompartmentation on metabolic pathways are however difficult to study experimentally. The pyrenoid is a microcompartment that is essential for a carbon concentrating mechanism (CCM) that improves the photosynthetic performance of eukaryotic algae. Using Chlamydomonas reinhardtii, we obtained experimental data on photosynthesis, metabolites, and proteins in CCM-induced and CCM-suppressed cells. We then employed a computational strategy to estimate how fluxes through the Calvin-Benson cycle are compartmented between the pyrenoid and the stroma. Our model predicts that ribulose-1,5-bisphosphate (RuBP), the substrate of Rubisco, and 3-phosphoglycerate (3PGA), its product, diffuse in and out of the pyrenoid, respectively, with higher fluxes in CCM-induced cells. It also indicates that there is no major diffusional barrier to metabolic flux between the pyrenoid and stroma. Our computational approach represents a stepping stone to understanding microcompartmentalized CCM in other organisms.
Introduction Flux phenotypes from different organisms and growth conditions allow better understanding of differential metabolic networks functions. Fluxes of metabolic reactions represent the integrated outcome of transcription, translation, and post-translational modifications, and directly affect growth and fitness. However, fluxes of intracellular metabolic reactions cannot be directly measured, but are estimated via metabolic flux analysis (MFA) that integrates data on isotope labeling patterns of metabolites with metabolic models. While the application of metabolomics technologies in photosynthetic organisms have resulted in unprecedented data from 13CO2-labeling experiments, the bottleneck in flux estimation remains the application of isotopically nonstationary MFA (INST-MFA). INST-MFA entails fitting a (large) system of coupled ordinary differential equations, with metabolite pools and reaction fluxes as parameters. Here, we focus on the Calvin-Benson cycle (CBC) as a key pathway for carbon fixation in photosynthesizing organisms and ask if approaches other than classical INST-MFA can provide reliable estimation of fluxes for reactions comprising this pathway.
Methods First, we show that flux estimation with the labeling patterns of all CBC intermediates can be formulated as a single constrained regression problem, avoiding the need for repeated simulation of time-resolved labeling patterns.
Results We then compare the flux estimates of the simulation-free constrained regression approach with those obtained from the classical INST-MFA based on labeling patterns of metabolites from the microalgae Chlamydomonas reinhardtii, Chlorella sorokiniana and Chlorella ohadii under different growth conditions.
Discussion Our findings indicate that, in data-rich scenarios, simulation-free regression-based approaches provide a suitable alternative for flux estimation from classical INST-MFA since we observe a high qualitative agreement (rs=0.89) to predictions obtained from INCA, a state-of-the-art tool for INST-MFA.
Large-scale biochemical models are of increasing sizes due to the consideration of interacting organisms and tissues. Model reduction approaches that preserve the flux phenotypes can simplify the analysis and predictions of steady-state metabolic phenotypes. However, existing approaches either restrict functionality of reduced models or do not lead to significant decreases in the number of modelled metabolites. Here, we introduce an approach for model reduction based on the structural property of balancing of complexes that preserves the steady-state fluxes supported by the network and can be efficiently determined at genome scale. Using two large-scale mass-action kinetic models of Escherichia coli, we show that our approach results in a substantial reduction of 99% of metabolites. Applications to genome-scale metabolic models across kingdoms of life result in up to 55% and 85% reduction in the number of metabolites when arbitrary and mass-action kinetics is assumed, respectively. We also show that predictions of the specific growth rate from the reduced models match those based on the original models. Since steady-state flux phenotypes from the original model are preserved in the reduced, the approach paves the way for analysing other metabolic phenotypes in large-scale biochemical networks.
The deficiency of a (bio)chemical reaction network can be conceptually interpreted as a measure of its ability to support exotic dynamical behavior and/or multistationarity. The classical definition of deficiency relates to the capacity of a network to permit variations of the complex formation rate vector at steady state, irrespective of the network kinetics. However, the deficiency is by definition completely insensitive to the fine details of the directionality of reactions as well as bounds on reaction fluxes. While the classical definition of deficiency can be readily applied in the analysis of unconstrained, weakly reversible networks, it only provides an upper bound in the cases where relevant constraints on reaction fluxes are imposed. Here we propose the concept of effective deficiency, which provides a more accurate assessment of the network’s capacity to permit steady state variations at the complex level for constrained networks of any reversibility patterns. The effective deficiency relies on the concept of nonstoichiometric balanced complexes, which we have already shown to be present in real-world biochemical networks operating under flux constraints. Our results demonstrate that the effective deficiency of real-world biochemical networks is smaller than the classical deficiency, indicating the effects of reaction directionality and flux bounds on the variation of the complex formation rate vector at steady state.
The current trends of crop yield improvements are not expected to meet the projected rise in demand. Genomic selection uses molecular markers and machine learning to identify superior genotypes with improved traits, such as growth. Plant growth directly depends on rates of metabolic reactions which transform nutrients into the building blocks of biomass. Here, we predict growth of Arabidopsis thaliana accessions by employing genomic prediction of reaction rates estimated from accession-specific metabolic models. We demonstrate that, comparing to classical genomic selection on the available data sets for 67 accessions, our approach improves the prediction accuracy for growth within and across nitrogen environments by 32.6% and 51.4%, respectively, and from optimal nitrogen to low carbon environment by 50.4%. Therefore, integration of molecular markers into metabolic models offers an approach to predict traits directly related to metabolism, and its usefulness in breeding can be examined by gathering matching datasets in crops. An increase in genomic selection (GS) accuracy can accelerate genetic gain by shortening the breeding cycles. Here, the authors introduce a network-based GS method that uses metabolic models and improves the prediction accuracy of Arabidopsis growth within and across environments.