Flux balance analysis of metabolism during growth by osmotic cell expansion and its application to tomato fruits

Metabolic modeling reveals distinct roles of sugars and carboxylic acids in stomatal opening as well as unexpected carbon fluxes

Guard cell metabolism is crucial for stomatal dynamics, but a full understanding of its role is hampered by experimental limitations and the flexible nature of the metabolic network. To tackle this challenge, we constructed a time-resolved stoichiometric model of guard cell metabolism that accounts for energy and osmolyte requirements and which is integrated with the mesophyll. The model resolved distinct roles for starch, sugars, and malate in guard cell metabolism and revealed several unexpected flux patterns in central metabolism. During blue light-mediated stomatal opening, starch breakdown was the most efficient way to generate osmolytes with downregulation of glycolysis allowing starch-derived glucose to accumulate as a cytosolic osmolyte. Maltose could also accumulate as a cytosolic osmoticum, although this made the metabolic system marginally less efficient. The metabolic energy for stomatal opening was predicted to be derived independently of starch, using nocturnally accumulated citrate which was metabolized in the tricarboxylic acid cycle to malate to provide mitochondrial reducing power for ATP synthesis. In white light-mediated stomatal opening, malate transferred reducing equivalents from guard cell photosynthesis to mitochondria for ATP production. Depending on the capacity for guard cell photosynthesis, glycolysis showed little flux during the day but was crucial for energy metabolism at night. In summary, our analyses have corroborated recent findings in Arabidopsis guard cell research, resolved conflicting observations by highlighting the flexibility of guard cell metabolism, and proposed new metabolic flux modes for further experimental testing.

From Division to Death: Metabolomic Analysis of Nicotiana tabacum BY-2 Cells Reveals the Complexity of Life in Batch Culture

Tobacco BY-2 cell culture is one of the most widely used models in plant biology. The main advantage of BY-2 suspension cultures is the synchronization of cell development and the appearance of polar elongation. In batch culture, BY-2 cells passed through the lag, proliferation, elongation, and stationary phases. During this process, the composition of the growth medium changed dramatically. Sucrose was rapidly eliminated; hexose first accumulated and then depleted. The medium's pH initially decreased and then rose with aging. As a result of the crosstalk between the internal and external stimuli, cells pass through complicated systemic rearrangements, which cause metabolomic alterations. The early stages were characterized by high levels of amino acids and sterols, which could be interpreted as the result of synthetic activity. The most intense rearrangements occurred between the proliferation and active elongation stages, including repression of amino acid accumulation and up-regulation of sugar metabolism. Later stages were distinguished by higher levels of secondary metabolites, which may be a non-specific response to deteriorating conditions. Senescence was followed by some increase in fatty acids and sterols as well as amino acids, and probably led to self-destructive processes. A correlation analysis revealed relationships between metabolites' covariation, their biochemical ratio, and the growth phase.

Resource allocation modeling for autonomous prediction of plant cell phenotypes

Predicting the plant cell response in complex environmental conditions is a challenge in plant biology. Here we developed a resource allocation model of cellular and molecular scale for the leaf photosynthetic cell of Arabidopsis thaliana, based on the Resource Balance Analysis (RBA) constraint-based modeling framework. The RBA model contains the metabolic network and the major macromolecular processes involved in the plant cell growth and survival and localized in cellular compartments. We simulated the model for varying environmental conditions of temperature, irradiance, partial pressure of CO2 and O2, and compared RBA predictions to known resource distributions and quantitative phenotypic traits such as the relative growth rate, the C:N ratio, and finally to the empirical characteristics of CO2 fixation given by the well-established Farquhar model. In comparison to other standard constraint-based modeling methods like Flux Balance Analysis, the RBA model makes accurate quantitative predictions without the need for empirical constraints. Altogether, we show that RBA significantly improves the autonomous prediction of plant cell phenotypes in complex environmental conditions, and provides mechanistic links between the genotype and the phenotype of the plant cell.

The Arabidopsis leaf quantitative atlas: a cellular and subcellular mapping through unified data integration

Quantitative analyses and models are required to connect a plant's cellular organisation with its metabolism. However, quantitative data are often scattered over multiple studies, and finding such data and converting them into useful information is time-consuming. Consequently, there is a need to centralise the available data and to highlight the remaining knowledge gaps. Here, we present a step-by-step approach to manually extract quantitative data from various information sources, and to unify the data format. First, data from Arabidopsis leaf were collated, checked for consistency and correctness and curated by cross-checking sources. Second, quantitative data were combined by applying calculation rules. They were then integrated into a unique comprehensive, referenced, modifiable and reusable data compendium representing an Arabidopsis reference leaf. This atlas contains the metrics of the 15 cell types found in leaves at the cellular and subcellular levels.

Plant Heterotrophic Cultures: No Food, No Growth

Plant cells are capable of uptaking exogenous organic substances. This inherited trait allows the development of heterotrophic cell cultures in various plants. The most common of them are Nicotiana tabacum and Arabidopsis thaliana. Plant cells are widely used in academic studies and as factories for valuable substance production. The repertoire of compounds supporting the heterotrophic growth of plant cells is limited. The best growth of cultures is ensured by oligosaccharides and their cleavage products. Primarily, these are sucrose, raffinose, glucose and fructose. Other molecules such as glycerol, carbonic acids, starch, and mannitol have the ability to support growth occasionally, or in combination with another substrate. Culture growth is accompanied by processes of specialization, such as elongation growth. This determines the pattern of the carbon budget. Culture ageing is closely linked to substrate depletion, changes in medium composition, and cell physiological rearrangements. A lack of substrate leads to starvation, which results in a decrease in physiological activity and the mobilization of resources, and finally in the loss of viability. The cause of the instability of cultivated cells may be the non-optimal metabolism under cultural conditions or the insufficiency of internal regulation.

Dynamic shifts in primary metabolism across fruit development stages in Capsicum chinense (cv. Habanero)

The development of fleshy fruits involves changes in size and mass, followed by cell differentiation, which is associated with anatomical and histological changes. Parallel to these changes, metabolic alterations lead to the production of osmolytes and energy that modify cell turgor pressure, thereby promoting cell expansion and fruit growth. Detailed information is known about these processes in climacteric fruits (e.g. tomato); however, the regulation of metabolism and its association with anatomical changes in non-climacteric fruit development are poorly understood. In this study, we used detailed anatomical and histological analyses to define three developmental phases of chili pepper (Capsicum chinense cv. Habanero): cell division, cell expansion, and ripening. We showed that each was marked by distinct metabolic profiles, underpinning the switches in energy metabolism to support cellular processes. Interestingly, mitochondrial activity was high in the early stages of development and declined over time, with a modest increase in O2 consumption by pericarp tissues at the beginning of the ripening stage. This respiratory-like burst was associated with the degradation of starch and malate, which are the sources of energy and carbon required for other processes associated with fruit maturation.

A role for fermentation in aerobic conditions as revealed by computational analysis of maize root metabolism during growth by cell elongation

The root is a well-studied example of cell specialisation, yet little is known about the metabolism that supports the transport functions and growth of different root cell types. To address this, we used computational modelling to study metabolism in the elongation zone of a maize lateral root. A functional-structural model captured the cell-anatomical features of the root and modelled how they changed as the root elongated. From these data, we derived constraints for a flux balance analysis model that predicted metabolic fluxes of the 11 concentric rings of cells in the root. We discovered a distinct metabolic flux pattern in the cortical cell rings, endodermis and pericycle (but absent in the epidermis) that involved a high rate of glycolysis and production of the fermentation end-products lactate and ethanol. This aerobic fermentation was confirmed experimentally by metabolite analysis. The use of fermentation in the model was not obligatory but was the most efficient way to meet the specific demands for energy, reducing power and carbon skeletons of expanding cells. Cytosolic acidification was avoided in the fermentative mode due to the substantial consumption of protons by lipid synthesis. These results expand our understanding of fermentative metabolism beyond that of hypoxic niches and suggest that fermentation could play an important role in the metabolism of aerobic tissues.

Elucidating the impact of in vitro cultivation on Nicotiana tabacum metabolism through combined in silico modeling and multiomics analysis

The systematical characterization and understanding of the metabolic behaviors are the basis of the efficient plant metabolic engineering and synthetic biology. Genome-scale metabolic networks (GSMNs) are indispensable tools for the comprehensive characterization of overall metabolic profile. Here we first constructed a GSMN of tobacco, which is one of the most widely used plant chassis, and then combined the tobacco GSMN and multiomics analysis to systematically elucidate the impact of in-vitro cultivation on the tobacco metabolic network. In-vitro cultivation is a widely used technique for plant cultivation, not only in the field of basic research but also for the rapid propagation of valuable horticultural and pharmaceutical plants. However, the systemic effects of in-vitro cultivation on overall plant metabolism could easily be overlooked and are still poorly understood. We found that in-vitro tobacco showed slower growth, less biomass and suppressed photosynthesis than soil-grown tobacco. Many changes of metabolites and metabolic pathways between in-vitro and soil-grown tobacco plants were identified, which notably revealed a significant increase of the amino acids content under in-vitro condition. The in silico investigation showed that in-vitro tobacco downregulated photosynthesis and primary carbon metabolism, while significantly upregulated the GS/GOGAT cycle, as well as producing more energy and less NADH/NADPH to acclimate in-vitro growth demands. Altogether, the combination of experimental and in silico analyses offers an unprecedented view of tobacco metabolism, with valuable insights into the impact of in-vitro cultivation, enabling more efficient utilization of in-vitro techniques for plant propagation and metabolic engineering.

Root phenotypes for improved nitrogen capture

Background

Suboptimal nitrogen availability is a primary constraint for crop production in low-input agroecosystems, while nitrogen fertilization is a primary contributor to the energy, economic, and environmental costs of crop production in high-input agroecosystems. In this article we consider avenues to develop crops with improved nitrogen capture and reduced requirement for nitrogen fertilizer.

Scope

Intraspecific variation for an array of root phenotypes has been associated with improved nitrogen capture in cereal crops, including architectural phenotypes that colocalize root foraging with nitrogen availability in the soil; anatomical phenotypes that reduce the metabolic costs of soil exploration, improve penetration of hard soil, and exploit the rhizosphere; subcellular phenotypes that reduce the nitrogen requirement of plant tissue; molecular phenotypes exhibiting optimized nitrate uptake kinetics; and rhizosphere phenotypes that optimize associations with the rhizosphere microbiome. For each of these topics we provide examples of root phenotypes which merit attention as potential selection targets for crop improvement. Several cross-cutting issues are addressed including the importance of soil hydrology and impedance, phenotypic plasticity, integrated phenotypes, in silico modeling, and breeding strategies using high throughput phenotyping for co-optimization of multiple phenes.

Conclusions

Substantial phenotypic variation exists in crop germplasm for an array of root phenotypes that improve nitrogen capture. Although this topic merits greater research attention than it currently receives, we have adequate understanding and tools to develop crops with improved nitrogen capture. Root phenotypes are underutilized yet attractive breeding targets for the development of the nitrogen efficient crops urgently needed in global agriculture.

Cell-type-specific metabolism in plants

Every plant organ contains tens of different cell types, each with a specialized function. These functions are intrinsically associated with specific metabolic flux distributions that permit the synthesis of the ATP, reducing equivalents and biosynthetic precursors demanded by the cell. Investigating such cell-type-specific metabolism is complicated by the mosaic of different cells within each tissue combined with the relative scarcity of certain types. However, techniques for the isolation of specific cells, their analysis in situ by microscopy, or modeling of their function in silico have permitted insight into cell-type-specific metabolism. In this review we present some of the methods used in the analysis of cell-type-specific metabolism before describing what we know about metabolism in several cell types that have been studied in depth; (i) leaf source and sink cells; (ii) glandular trichomes that are capable of rapid synthesis of specialized metabolites; (iii) guard cells that must accumulate large quantities of the osmolytes needed for stomatal opening; (iv) cells of seeds involved in storage of reserves; and (v) the mesophyll and bundle sheath cells of C4 plants that participate in a CO₂ concentrating cycle. Metabolism is discussed in terms of its principal features, connection to cell function and what factors affect the flux distribution. Demand for precursors and energy, availability of substrates and suppression of deleterious processes are identified as key factors in shaping cell-type-specific metabolism.

Spatial resolution of an integrated C4+CAM photosynthetic metabolism

C4 and CAM photosynthesis have repeatedly evolved in plants over the past 30 million years. Because both repurpose the same set of enzymes but differ in their spatial and temporal deployment, they have long been considered as distinct and incompatible adaptations. Portulaca contains multiple C4 species that perform CAM when droughted. Spatially explicit analyses of gene expression reveal that C4 and CAM systems are completely integrated in Portulaca oleracea, with CAM and C4 carbon fixation occurring in the same cells and CAM-generated metabolites likely incorporated directly into the C4 cycle. Flux balance analysis corroborates the gene expression findings and predicts an integrated C4+CAM system under drought. This first spatially explicit description of a C4+CAM photosynthetic metabolism presents a potential new blueprint for crop improvement.

Bruceine D Identified as a Drug Candidate against Breast Cancer by a Novel Drug Selection Pipeline and Cell Viability Assay

The multi-target effects of natural products allow us to fight complex diseases like cancer on multiple fronts. Unlike docking techniques, network-based approaches such as genome-scale metabolic modelling can capture multi-target effects. However, the incompleteness of natural product target information reduces the prediction accuracy of in silico gene knockout strategies. Here, we present a drug selection workflow based on context-specific genome-scale metabolic models, built from the expression data of cancer cells treated with natural products, to predict cell viability. The workflow comprises four steps: first, in silico single-drug and drug combination predictions; second, the assessment of the effects of natural products on cancer metabolism via the computation of a dissimilarity score between the treated and control models; third, the identification of natural products with similar effects to the approved drugs; and fourth, the identification of drugs with the predicted effects in pathways of interest, such as the androgen and estrogen pathway. Out of the initial 101 natural products, nine candidates were tested in a 2D cell viability assay. Bruceine D, emodin, and scutellarein showed a dose-dependent inhibition of MCF-7 and Hs 578T cell proliferation with IC50 values between 0.7 to 65 μM, depending on the drug and cell line. Bruceine D, extracted from Brucea javanica seeds, showed the highest potency.

A hybrid kinetic and constraint-based model of leaf metabolism allows predictions of metabolic fluxes in different environments

While flux balance analysis (FBA) provides a framework for predicting steady-state leaf metabolic network fluxes, it does not readily capture the response to environmental variables without being coupled to other modelling formulations. To address this, we coupled an FBA model of 903 reactions of soybean (Glycine max) leaf metabolism with e-photosynthesis, a dynamic model that captures the kinetics of 126 reactions of photosynthesis and associated chloroplast carbon metabolism. Successful coupling was achieved in an iterative formulation in which fluxes from e-photosynthesis were used to constrain the FBA model and then, in turn, fluxes computed from the FBA model used to update parameters in e-photosynthesis. This process was repeated until common fluxes in the two models converged. Coupling did not hamper the ability of the kinetic module to accurately predict the carbon assimilation rate, photosystem II electron flux, and starch accumulation of field-grown soybean at two CO2 concentrations. The coupled model also allowed accurate predictions of additional parameters such as nocturnal respiration, as well as analysis of the effect of light intensity and elevated CO2 on leaf metabolism. Predictions included an unexpected decrease in the rate of export of sucrose from the leaf at high light, due to altered starch-sucrose partitioning, and altered daytime flux modes in the tricarboxylic acid cycle at elevated CO2 . Mitochondrial fluxes were notably different between growing and mature leaves, with greater anaplerotic, tricarboxylic acid cycle and mitochondrial ATP synthase fluxes predicted in the former, primarily to provide carbon skeletons and energy for protein synthesis.

Reproductive water supply is prioritized during drought in tomato

Reproductive success largely defines the fitness of plant species. Understanding how heat and drought affect plant reproduction is thus key to predicting future plant fitness under rising global temperatures. Recent work suggests reproductive tissues are highly vulnerable to water stress in perennial plants where reproductive sacrifice could preserve plant survival. However, most crop species are annuals where such a strategy would theoretically reduce fitness. We examined the reproductive strategy of tomato (Solanum lycopersicum var. Rheinlands Ruhm) to determine whether water supply to fruits is prioritized above vegetative tissues during drought. Using optical methods, we mapped xylem cavitation and tissue shrinkage in vegetative and reproductive organs during dehydration to determine the priority of water flow under acute water stress. Stems and peduncles of tomato showed significantly greater xylem cavitation resistance than leaves. This maintenance of intact water supply enabled tomato fruit to continue to expand during acute water stress, utilizing xylem water made available by tissue collapse and early cavitation of leaves. Here, tomato plants prioritize water supply to reproductive tissues, maintaining fruit development under drought conditions. These results emphasize the critical role of water transport in shaping life history and suggest a broad relevance of hydraulic prioritization in plant ecology.

Plant cell cultures as heterologous bio-factories for secondary metabolite production

Synthetic biology has been developing rapidly in the last decade and is attracting increasing attention from many plant biologists. The production of high-value plant-specific secondary metabolites is, however, limited mostly to microbes. This is potentially problematic because of incorrect post-translational modification of proteins and differences in protein micro-compartmentalization, substrate availability, chaperone availability, product toxicity, and cytochrome p450 reductase enzymes. Unlike other heterologous systems, plant cells may be a promising alternative for the production of high-value metabolites. Several commercial plant suspension cell cultures from different plant species have been used successfully to produce valuable metabolites in a safe, low cost, and environmentally friendly manner. However, few metabolites are currently being biosynthesized using plant platforms, with the exception of the natural pigment anthocyanin. Both Arabidopsis thaliana and Nicotiana tabacum cell cultures can be developed by multiple gene transformations and CRISPR-Cas9 genome editing. Given that the introduction of heterologous biosynthetic pathways into Arabidopsis and N. tabacum is not widely used, the biosynthesis of foreign metabolites is currently limited; however, therein lies great potential. Here, we discuss the exemplary use of plant cell cultures and prospects for using A. thaliana and N. tabacum cell cultures to produce valuable plant-specific metabolites.

Advances in flux balance analysis by integrating machine learning and mechanism-based models

The availability of multi-omics data sets and genome-scale metabolic models for various organisms provide a platform for modeling and analyzing genotype-to-phenotype relationships. Flux balance analysis is the main tool for predicting flux distributions in genome-scale metabolic models and various data-integrative approaches enable modeling context-specific network behavior. Due to its linear nature, this optimization framework is readily scalable to multi-tissue or -organ and even multi-organism models. However, both data and model size can hamper a straightforward biological interpretation of the estimated fluxes. Moreover, flux balance analysis simulates metabolism at steady-state and thus, in its most basic form, does not consider kinetics or regulatory events. The integration of flux balance analysis with complementary data analysis and modeling techniques offers the potential to overcome these challenges. In particular machine learning approaches have emerged as the tool of choice for data reduction and selection of most important variables in big data sets. Kinetic models and formal languages can be used to simulate dynamic behavior. This review article provides an overview of integrative studies that combine flux balance analysis with machine learning approaches, kinetic models, such as physiology-based pharmacokinetic models, and formal graphical modeling languages, such as Petri nets. We discuss the mathematical aspects and biological applications of these integrated approaches and outline challenges and future perspectives.

Modelling predicts tomatoes can be bigger and sweeter if biophysical factors and transmembrane transports are fine-tuned during fruit development

The trade-off between yield and quality, a major problem for the production of fleshy fruits, involves fruit expansive growth and sugar metabolism. Here we developed an integrative model by coupling a biophysical model of fleshy fruit growth processes, including water and carbon fluxes and organ expansion, with an enzyme-based kinetic model of sugar metabolism to better understand the interactions between these two processes. The integrative model was initially tested on tomato fruit, a model system for fleshy fruit. The integrative model closely simulated the biomass and major carbon metabolites of tomato fruit developing under optimal or stress conditions. The model also performed robustly when simulating the fruit size and sugar concentrations of different tomato genotypes including wild species. The validated model was used to explore ways of uncoupling the size-sweetness trade-off in fruit. Model-based virtual experiments suggested that larger sweeter tomatoes could be obtained by simultaneously manipulating certain biophysical factors and transmembrane transports. The integrative fleshy fruit model provides a promising tool to facilitate the targeted bioengineering and breeding of tomatoes and other fruits.

Recent advances in constraint and machine learning-based metabolic modeling by leveraging stoichiometric balances, thermodynamic feasibility and kinetic law formalisms

Understanding the governing principles behind organisms' metabolism and growth underpins their effective deployment as bioproduction chassis. A central objective of metabolic modeling is predicting how metabolism and growth are affected by both external environmental factors and internal genotypic perturbations. The fundamental concepts of reaction stoichiometry, thermodynamics, and mass action kinetics have emerged as the foundational principles of many modeling frameworks designed to describe how and why organisms allocate resources towards both growth and bioproduction. This review focuses on the latest algorithmic advancements that have integrated these foundational principles into increasingly sophisticated quantitative frameworks.

An insight to flux-balance analysis for biochemical networks

Systems biology is one of the integrated ways to study biological systems and is more favourable than the earlier used approaches. It includes metabolic pathway analysis, modelling, and regulatory as well as signal transduction for getting insights into cellular behaviour. Among the various techniques of modelling, simulation, analysis of networks and pathways, flux-based analysis (FBA) has been recognised because of its extensibility as well as simplicity. It is widely accepted because it is not like a mechanistic simulation which depends on accurate kinetic data. The study of fluxes through the network is informative and can give insights even in the absence of kinetic data. FBA is one of the widely used tools to study biochemical networks and needs information of reaction stoichiometry, growth requirements, specific measurement parameters of the biological system, in particular the reconstruction of the metabolic network for the genome-scale, many of which have already been built previously. It defines the boundaries of flux distributions which are possible and achievable with a defined set of genes. This review article gives an insight into FBA, from the extension of flux balancing to mathematical representation followed by a discussion about the formulation of flux-balance analysis problems, defining constraints for the stoichiometry of the pathways and the tools that can be used in FBA such as FASIMA, COBRA toolbox, and OptFlux. It also includes broader areas in terms of applications which can be covered by FBA as well as the queries which can be addressed through FBA.

Model-assisted comparison of sugar accumulation patterns in ten fleshy fruits highlights differences between herbaceous and woody species

BACKGROUND AND AIMS

Sugar concentration is a key determinant of fruit quality. Soluble sugars and starch concentrations in fruits vary greatly from one species to another. The aim of this study was to investigate similarities and differences in sugar accumulation strategies across ten contrasting fruit species using a modelling approach.

METHODS

We developed a coarse-grained model of primary metabolism based on the description of the main metabolic and hydraulic processes (synthesis of compounds other than sugar and starch, synthesis and hydrolysis of starch, and water dilution) involved in the accumulation of soluble sugars during fruit development.

KEY RESULTS

Statistical analyses based on metabolic rates separated the species into six groups according to the rate of synthesis of compounds other than sugar and starch. Herbaceous species (cucumber, tomato, eggplant, pepper and strawberry) were characterized by a higher synthesis rate than woody species (apple, nectarine, clementine, grape and kiwifruit). Inspection of the dynamics of the processes involved in sugar accumulation revealed that net sugar importation, metabolism and dilution processes were remarkably synchronous in most herbaceous plants, whereas in kiwifruit, apple and nectarine, processes related to starch metabolism were temporally separated from other processes. Strawberry, clementine and grape showed a distinct dynamic compared with all other species.

CONCLUSIONS

Overall, these results provide fresh insights into species-specific regulatory strategies and into the role of starch metabolism in the accumulation of soluble sugars in fleshy fruits. In particular, inter-specific differences in development period shape the co-ordination of metabolic processes and affect priorities for carbon allocation across species. The six metabolic groups identified by our analysis do not show a clear separation into climacteric and non-climacteric species, possibly suggesting that the metabolic processes related to sugar concentration are not greatly affected by ethylene-associated events.

Modeling the metabolic interplay between a parasitic worm and its bacterial endosymbiont allows the identification of novel drug targets

The filarial nematode Brugia malayi represents a leading cause of disability in the developing world, causing lymphatic filariasis in nearly 40 million people. Currently available drugs are not well-suited to mass drug administration efforts, so new treatments are urgently required. One potential vulnerability is the endosymbiotic bacteria Wolbachia-present in many filariae-which is vital to the worm. Genome scale metabolic networks have been used to study prokaryotes and protists and have proven valuable in identifying therapeutic targets, but have only been applied to multicellular eukaryotic organisms more recently. Here, we present iDC625, the first compartmentalized metabolic model of a parasitic worm. We used this model to show how metabolic pathway usage allows the worm to adapt to different environments, and predict a set of 102 reactions essential to the survival of B. malayi. We validated three of those reactions with drug tests and demonstrated novel antifilarial properties for all three compounds.