Photorespiration causes a substantial decrease in crop yield because of mitochondrial decarboxylation. Alternative pathways (APs) have been designed to relocate the decarboxylation step or even fix additional carbon. To improve the success of transferring those engineered APs from model species to crops, we must understand how they will interact with metabolism and how plant physiology affects their performance. Here, we used multiple mathematical modeling techniques to analyze and compare existing AP designs. We show that carbon-fixing APs are the most promising candidates to replace native photorespiration in major crop species. Our results demonstrate the different metabolic routes that APs use to increase yield and which plant physiology can profit the most from them. We anticipate our results to guide the design of new APs and to help improve existing ones.

SUMMARY: Several commercially important secondary metabolites are produced and accumulated in high amounts by glandular trichomes, giving the prospect of using them as metabolic cell factories. Due to extremely high metabolic fluxes through glandular trichomes, previous research focused on how such flows are achieved. The question regarding their bioenergetics became even more interesting with the discovery of photosynthetic activity in some glandular trichomes. Despite recent advances, how primary metabolism contributes to the high metabolic fluxes in glandular trichomes is still not fully elucidated. Using computational methods and available multi‐omics data, we first developed a quantitative framework to investigate the possible role of photosynthetic energy supply in terpenoid production and next tested experimentally the simulation‐driven hypothesis. With this work, we provide the first reconstruction of specialised metabolism in Type‐VI photosynthetic glandular trichomes of Solanum lycopersicum. Our model predicted that increasing light intensities results in a shift of carbon partitioning from catabolic to anabolic reactions driven by the energy availability of the cell. Moreover, we show the benefit of shifting between isoprenoid pathways under different light regimes, leading to a production of different classes of terpenes. Our computational predictions were confirmed in vivo, demonstrating a significant increase in production of monoterpenoids while the sesquiterpenes remained unchanged under higher light intensities. The outcomes of this research provide quantitative measures to assess the beneficial role of chloroplast in glandular trichomes for enhanced production of secondary metabolites and can guide the design of new experiments that aim at modulating terpenoid production.

Mathematical modeling of metabolic networks is a powerful approach to investigate the underlying principles of metabolism and growth. Such approaches include, among others, differential-equation-based modeling of metabolic systems, constraint-based modeling and metabolic network expansion of metabolic networks. Most of these methods are well established and are implemented in numerous software packages, but these are scattered between different programming languages, packages and syntaxes. This complicates establishing straight forward pipelines integrating model construction and simulation. We present a Python package moped that serves as an integrative hub for reproducible construction, modification, curation and analysis of metabolic models. moped supports draft reconstruction of models directly from genome/proteome sequences and pathway/genome databases utilizing GPR annotations, providing a completely reproducible model construction and curation process within executable Python scripts. Alternatively, existing models published in SBML format can be easily imported. Models are represented as Python objects, for which a wide spectrum of easy-to-use modification and analysis methods exist. The model structure can be manually altered by adding, removing or modifying reactions, and gap-filling reactions can be found and inspected. This greatly supports the development of draft models, as well as the curation and testing of models. Moreover, moped provides several analysis methods, in particular including the calculation of biosynthetic capacities using metabolic network expansion. The integration with other Python-based tools is facilitated through various model export options. For example, a model can be directly converted into a CobraPy object for constraint-based analyses. moped is a fully documented and expandable Python package. We demonstrate the capability to serve as a hub for integrating reproducible model construction and curation, database import, metabolic network expansion and export for constraint-based analyses.

Photosynthetic organisms use photosynthesis to harvest sunlight and convert the solar energy into chemical energy, which is then used to reduce atmospheric carbon dioxide into organic molecules. This process forms the basis of all life on Earth, and stands at the beginning of the food chain which feeds the world population. Not surprisingly, many research efforts are currently ongoing aiming at improving growth and product yield of photosynthetic organisms, and several of these activities directly target the photosynthetic pathways. Metabolic Control Analysis (MCA) shows that, in general, the control over a metabolic flux, such as carbon fixation, is distributed among several steps and highly dependent on the external conditions. Therefore, the concept of a single 'rate-limiting' step is hardly ever applicable, and as a consequence, any strategy relying on improving a single molecular process in a complex metabolic system is bound to fail to yield the expected results. In photosynthesis, reports on which processes exert the highest control over carbon fixation are contradictory. This refers to both the photosynthetic 'light' reactions harvesting photons and the 'dark' reactions of the Calvin-Benson-Bassham Cycle (CBB cycle). Here, we employ a recently developed mathematical model, which describes photosynthesis as an interacting supply-demand system, to systematically study how external conditions affect the control over carbon fixation fluxes.

The modelbase package is a free expandable Python package for building and analysing dynamic mathematical models of biological systems. Originally it was designed for the simulation of metabolic systems, but it can be used for virtually any deterministic chemical processes. modelbase provides easy construction methods to define reactions and their rates. Based on the rates and stoichiometries, the system of differential equations is assembled automatically. modelbase minimises the constraints imposed on the user, allowing for easy and dynamic access to all variables, including derived ones, in a convenient manner. A simple incorporation of algebraic equations is, for example, convenient to study systems with rapid equilibrium or quasi steady-state approximations. Moreover, modelbase provides construction methods that automatically build all isotope-specific versions of a particular reaction, making it a convenient tool to analyse non-steady state isotope-labelling experiments.

Terpenoids make up the largest class of specialized metabolites with over 180,000 reported compounds currently across all kingdoms of life. Their synthesis accentuates one of natures most choreographed enzymatic and non-reversible chemistries, leading to an extensive range of structural functionality and diversity. Current terpenoid repositories provide a seemingly endless landscape to systematically survey for information regarding structure, sourcing, and synthesis. Efforts here investigate entries for the 20-carbon diterpenoid variants and deconstruct the complex patterns into simple, categorical groups. This deconstruction approach reduces over 60,000 unique diterpenoid structures to less than 1,000 categorical structures. Furthermore, the majority of diterpene entries (over 75%) can be represented by less than 25 core skeletons. Natural diterpenoid abundance was mapped throughout the tree of life and structural diversity was correlated at an atom-and-bond resolution. Additionally, all identified core structures provide guidelines for predicting how diterpene diversity originates via the mechanisms catalyzed by diterpene synthases. Over 95% of diterpenoid structures rely on cyclization. Here a reconstructive approach is reapplied based on known biochemical rules to model the birth of compound diversity. Reconstruction enabled prediction of highly probable synthesis mechanisms for bioactive taxane-relatives, which were discovered over three decades ago. This computational synthesis validates previously identified reaction products and pathways, as well as enables predicting trajectories for synthesizing real and theoretical compounds. This deconstructive and reconstructive approach applied to the diterpene landscape provides modular, flexible, and an easy-to-use toolset for categorically simplifying otherwise complex or hidden patterns. Significance Statement: We take a deconstructive and reconstructive approach to explore the origins of the diterpene landscape. Introduction of a navigational toolset enables users to survey compound libraries in ways formerly uncharted. Their utility demonstrated here, maps out diterpene cyclization routes, critical intermediate waypoints, and guidance for how to arrive at compounds previously off-the-map. Information acquired from these tools may imply the diterpene landscape is vastly unexplored, with the plateau for discovery potentially still out of sight.