A bacterial glucanotransferase can replace the complex maltose metabolism required for starch to sucrose conversion in leaves at night

OsDPE2 Regulates Rice Panicle Morphogenesis by Modulating the Content of Starch

Starch is a carbon sink for most plants, and its biological role changes with response to the environment and during plant development. Disproportionating Enzyme 2 (DPE2) is a 4-α-glycosyltransferase involved in starch degradation in plants at night. LAX1 plays a vital role in axillary meristem initiation in rice. Herein, results showed that Oryza sativa Disproportionating Enzyme 2 (OsDPE2) could rescue the mutant phenotype of lax1-6, LAX1 mutant. OsDPE2 encodes rice DPE2 located in the cytoplasm. In this study, OsDPE2 affected the vegetative plant development of rice via DPE2 enzyme. Additionally, OsDPE2 regulated the reproductive plant development of rice by modulating starch content in young panicles. Furthermore, haplotype OsDPE2(AQ) with higher DPE2 enzyme activity increased the panicle yield of rice. In summary, OsDPE2 can regulate vegetative and reproductive plant development of rice by modulating starch content. Furthermore, DPE2 activities of OsDPE2 haplotypes are associated with the panicle yield of rice. This study provides guidance for rice breeding to improve panicle yield traits.

Cytosolic disproportionating enzyme2 is essential for pollen germination and pollen tube elongation in rice

Degradation of starch accumulated in pollen provides energy and cellular materials for pollen germination and pollen tube elongation. Little is known about the function of cytosolic disproportionating enzyme2 (DPE2) in rice (Oryza sativa). Here, we obtained several DPE2 knockout mutant (dpe2) lines via genomic editing and found that the mutants grew and developed normally but with greatly reduced seed-setting rates. Reciprocal crosses between dpe2 and wild-type plants demonstrated that the mutant was male sterile. In vitro and in vivo examinations revealed that the pollen of the dpe2 mutant developed and matured normally but was defective in germination and elongation. DPE2 deficiency increased maltose content in pollen, whereas it reduced the levels of starch, glucose, fructose, and adenosine triphosphate (ATP). Exogenous supply of glucose or ATP to the germination medium partially rescued the pollen germination defects of dpe2. The expression of cytosolic phosphorylase2 (Pho2) increased significantly in dpe2 pollen. Knockout of Pho2 resulted in a semi-sterile phenotype. We failed to obtain homozygous dpe2 pho2 double mutant lines. Our results demonstrate that maltose catalyzed by DPE2 to glucose is the main energy source for pollen germination and pollen tube elongation, while Pho2 might partially compensate for deficiency of DPE2.

Rising rates of starch degradation during daytime and trehalose 6-phosphate optimize carbon availability

Many plants, including Arabidopsis (Arabidopsis thaliana), accumulate starch in the light and remobilize it to support maintenance and growth at night. Starch synthesis and degradation are usually viewed as temporally separate processes. Recently, we reported that starch is also degraded in the light. Degradation rates are generally low early in the day but rise with time. Here, we show that the rate of degradation in the light depends on time relative to dawn rather than dusk. We also show that degradation in the light is inhibited by trehalose 6-phosphate, a signal for sucrose availability. The observed responses of degradation in the light can be simulated by a skeletal model in which the rate of degradation is a function of starch content divided by time remaining until dawn. The fit is improved by extension to include feedback inhibition of starch degradation by trehalose 6-phosphate. We also investigate possible functions of simultaneous starch synthesis and degradation in the light, using empirically parameterized models and experimental approaches. The idea that this cycle buffers growth against falling rates of photosynthesis at twilight is supported by data showing that rates of protein and cell wall synthesis remain high during a simulated dusk twilight. Degradation of starch in the light may also counter over-accumulation of starch in long photoperiods and stabilize signaling around dusk. We conclude that starch degradation in the light is regulated by mechanisms similar to those that operate at night and is important for stabilizing carbon availability and signaling, thus optimizing growth in natural light conditions.

Production of Large-Ring Cyclodextrins by Amylomaltases

Amylomaltase is a well-known glucan transferase that can produce large ring cyclodextrins (LR-CDs) or so-called cycloamyloses via cyclization reaction. Amylomaltases have been found in several microorganisms and their optimum temperatures are generally around 60–70 °C for thermostable amylomaltases and 30–45 °C for the enzymes from mesophilic bacteria and plants. The optimum pHs for mesophilic amylomaltases are around pH 6.0–7.0, while the thermostable amylomaltases are generally active at more acidic conditions. Size of LR-CDs depends on the source of amylomaltases and the reaction conditions including pH, temperature, incubation time, and substrate. For example, in the case of amylomaltase from Corynebacterium glutamicum, LR-CD productions at alkaline pH or at a long incubation time favored products with a low degree of polymerization. In this review, we explore the synthesis of LR-CDs by amylomaltases, structural information of amylomaltases, as well as current applications of LR-CDs and amylomaltases.

The challenges in production technology, health-associated functions, physico-chemical properties and food applications of isomaltooligosaccharides

Isomaltooligosaccharides (IMOs) are recognized as functional food ingredients with prebiotic potential that deliver health benefits. IMOs have attained commercial interest as they are produced from low-cost agricultural products that are widely available and have prospective applications in the food industry. The review examines the various production processes and the main challenges involved in deriving diverse structures of IMO with maximized yield and increased functionality. The different characterization and purification techniques employed for structural elucidation, the physico-chemical importance, technological properties, food-based applications and biological effects (in vitro and in vivo interventions) have been discussed in detail. The key finding is the need for research involving biotechnological and enzymology aspects to simplify the production technologies that meet the industrial and consumer requirements. The knowledge from this article delivers a clear insight to scientists, food technologists and the general public for the improved utilization of IMOs to support the emerging market for functional foods and nutraceuticals.

Starch granule initiation in Arabidopsis thaliana chloroplasts

The initiation of starch granule formation and the mechanism controlling the number of granules per plastid have been some of the most elusive aspects of starch metabolism. This review covers the advances made in the study of these processes. The analyses presented herein depict a scenario in which starch synthase isoform 4 (SS4) provides the elongating activity necessary for the initiation of starch granule formation. However, this protein does not act alone; other polypeptides are required for the initiation of an appropriate number of starch granules per chloroplast. The functions of this group of polypeptides include providing suitable substrates (maltooligosaccharides) to SS4, the localization of the starch initiation machinery to the thylakoid membranes, and facilitating the correct folding of SS4. The number of starch granules per chloroplast is tightly regulated and depends on the developmental stage of the leaves and their metabolic status. Plastidial phosphorylase (PHS1) and other enzymes play an essential role in this process since they are necessary for the synthesis of the substrates used by the initiation machinery. The mechanism of starch granule formation initiation in Arabidopsis seems to be generalizable to other plants and also to the synthesis of long-term storage starch. The latter, however, shows specific features due to the presence of more isoforms, the absence of constantly recurring starch synthesis and degradation, and the metabolic characteristics of the storage sink organs.

Evolution and function of carbohydrate reserve biosynthesis in parasitic protists

Nearly all eukaryotic cells synthesize carbohydrate reserves, such as glycogen, starch, or low-molecular-weight oligosaccharides. However, a number of parasitic protists have lost this capacity while others have lost, and subsequently evolved, entirely new pathways. Recent studies suggest that retention, loss, or acquisition of these pathways in different protists is intimately linked to their lifestyle. In particular, parasites with carbohydrate reserves often establish long-lived chronic infections and/or produce environmental cysts, whereas loss of these pathways is associated with parasites that have highly proliferative and metabolically active life-cycle stages. The evolution of mannogen biosynthesis in Leishmania and related parasites indicates that these pathways have played a role in defining the host range and niches occupied by some protists.

A putative novel starch-binding domain revealed by in silico analysis of the N-terminal domain in bacterial amylomaltases from the family GH77

The family GH77 contains 4-α-glucanotransferase acting on α-1,4-glucans, known as amylomaltase in prokaryotes and disproportionating enzyme in plants. A group of bacterial GH77 members, represented by amylomaltases from Escherichia coli and Corynebacterium glutamicum, possesses an N-terminal extension that forms a distinct immunoglobulin-like fold domain, of which no function has been identified. Here, in silico analysis of 100 selected sequences of N-terminal domain homologues disclosed several well-conserved residues, among which Tyr108 (E. coli amylomaltase numbering) may be involved in α-glucan binding. These N-terminal domains, therefore, may represent a new type of starch-binding domain and define a new CBM family. This hypothesis is supported by docking of maltooligosaccharides to the N-terminal domain in amylomaltases, representing the four clusters of the phylogenetic tree.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02787-8.

Starch: A Flexible, Adaptable Carbon Store Coupled to Plant Growth

Research in the past decade has uncovered new and surprising information about the pathways of starch synthesis and degradation. This includes the discovery of previously unsuspected protein families required both for processes and for the long-sought mechanism of initiation of starch granules. There is also growing recognition of the central role of leaf starch turnover in making carbon available for growth across the day-night cycle. Sophisticated systems-level control mechanisms involving the circadian clock set rates of nighttime starch mobilization that maintain a steady supply of carbon until dawn and modulate partitioning of photosynthate into starch in the light, optimizing the fraction of assimilated carbon that can be used for growth. These discoveries also uncover complexities: Results from experiments with Arabidopsis leaves in conventional controlled environments are not necessarily applicable to other organs or species or to growth in natural, fluctuating environments.

Design starch: stochastic modeling of starch granule biogenesis

Starch is the most widespread and abundant storage carbohydrate in plants and the main source of carbohydrate in the human diet. Owing to its remarkable properties and commercial applications, starch is still of growing interest. Its unique granular structure made of intercalated layers of amylopectin and amylose has been unraveled thanks to recent progress in microscopic imaging, but the origin of such periodicity is still under debate. Both amylose and amylopectin are made of linear chains of α-1,4-bound glucose residues, with branch points formed by α-1,6 linkages. The net difference in the distribution of chain lengths and the branching pattern of amylose (mainly linear), compared with amylopectin (racemose structure), leads to different physico-chemical properties. Amylose is an amorphous and soluble polysaccharide, whereas amylopectin is insoluble and exhibits a highly organized structure of densely packed double helices formed between neighboring linear chains. Contrarily to starch degradation that has been investigated since the early 20th century, starch production is still poorly understood. Most enzymes involved in starch growth (elongation, branching, debranching, and partial hydrolysis) are now identified. However, their specific action, their interplay (cooperative or competitive), and their kinetic properties are still largely unknown. After reviewing recent results on starch structure and starch growth and degradation enzymatic activity, we discuss recent results and current challenges for growing polysaccharides on granular surface. Finally, we highlight the importance of novel stochastic models to support the analysis of recent and complex experimental results, and to address how macroscopic properties emerge from enzymatic activity and structural rearrangements.

Photometric assay of maltose and maltose-forming enzyme activity by using 4-alpha-glucanotransferase (DPE2) from higher plants

Maltose frequently occurs as intermediate of the central carbon metabolism of prokaryotic and eukaryotic cells. Various mutants possess elevated maltose levels. Maltose exists as two anomers, (α- and β-form) which are rapidly interconverted without requiring enzyme-mediated catalysis. As maltose is often abundant together with other oligoglucans, selective quantification is essential. In this communication, we present a photometric maltose assay using 4-alpha-glucanotransferase (AtDPE2) from Arabidopsis thaliana. Under in vitro conditions, AtDPE2 utilizes maltose as glucosyl donor and glycogen as acceptor releasing the other hexosyl unit as free glucose which is photometrically quantified following enzymatic phosphorylation and oxidation. Under the conditions used, DPE2 does not noticeably react with other di- or oligosaccharides. Selectivity compares favorably with that of maltase frequently used in maltose assays. Reducing end interconversion of the two maltose anomers is in rapid equilibrium and, therefore, the novel assay measures total maltose contents. Furthermore, an AtDPE2-based continuous photometric assay is presented which allows to quantify β-amylase activity and was found to be superior to a conventional test. Finally, the AtDPE2-based maltose assay was used to quantify leaf maltose contents of both Arabidopsis wild type and AtDPE2-deficient plants throughout the light-dark cycle. These data are presented together with assimilatory starch levels.

Functional and structural characterization of plastidic starch phosphorylase during barley endosperm development

The production of starch is essential for human nutrition and represents a major metabolic flux in the biosphere. The biosynthesis of starch in storage organs like barley endosperm operates via two main pathways using different substrates: starch synthases use ADP-glucose to produce amylose and amylopectin, the two major components of starch, whereas starch phosphorylase (Pho1) uses glucose-1-phosphate (G1P), a precursor for ADP-glucose production, to produce α-1,4 glucans. The significance of the Pho1 pathway in starch biosynthesis has remained unclear. To elucidate the importance of barley Pho1 (HvPho1) for starch biosynthesis in barley endosperm, we analyzed HvPho1 protein production and enzyme activity levels throughout barley endosperm development and characterized structure-function relationships of HvPho1. The molecular mechanisms underlying the initiation of starch granule biosynthesis, that is, the enzymes and substrates involved in the initial transition from simple sugars to polysaccharides, remain unclear. We found that HvPho1 is present as an active protein at the onset of barley endosperm development. Notably, purified recombinant protein can catalyze the de novo production of α-1,4-glucans using HvPho1 from G1P as the sole substrate. The structural properties of HvPho1 provide insights into the low affinity of HvPho1 for large polysaccharides like starch or amylopectin. Our results suggest that HvPho1 may play a role during the initiation of starch biosynthesis in barley.

Characterization of surface binding sites in glycoside hydrolases: A computational study

Structural properties of carbohydrate surface binding sites (SBSs) were investigated with computational methods. Eighty-five SBSs of 44 enzymes in 119 Protein Data Bank (PDB) files were collected as a dataset. On the basis of SBSs shape, they were divided into 3 categories: flat surfaces, clefts, and cavities (types A, B, and C, respectively). Ligand varieties showed the correlation between shape of SBSs and ligands size. To reduce cut-off differences in each SBSs with different ligand size, molecular docking were performed. Molecular docking results were used to refine SBSs classification and binding sites cut-off. Docking results predicted putative ligands positions and displayed dependence of the ligands binding mode to the structural type of SBSs. Physicochemical properties of SBSs were calculated for all docking results with YASARA Structure. The results showed that all SBSs are hydrophilic, while their charges could vary and depended to ligand size and defined cut-off. Surface binding sites type B had highest average values of solvent accessible surface area. Analysis of interactions showed that hydrophobic interactions occur more than hydrogen bonds, which is related to the presence of aromatic residues and carbohydrates interactions.

Using Carbohydrate Interaction Assays to Reveal Novel Binding Sites in Carbohydrate Active Enzymes

Carbohydrate active enzymes often contain auxiliary binding sites located either on independent domains termed carbohydrate binding modules (CBMs) or as so-called surface binding sites (SBSs) on the catalytic module at a certain distance from the active site. The SBSs are usually critical for the activity of their cognate enzyme, though they are not readily detected in the sequence of a protein, but normally require a crystal structure of a complex for their identification. A variety of methods, including affinity electrophoresis (AE), insoluble polysaccharide pulldown (IPP) and surface plasmon resonance (SPR) have been used to study auxiliary binding sites. These techniques are complementary as AE allows monitoring of binding to soluble polysaccharides, IPP to insoluble polysaccharides and SPR to oligosaccharides. Here we show that these methods are useful not only for analyzing known binding sites, but also for identifying new ones, even without structural data available. We further verify the chosen assays discriminate between known SBS/CBM containing enzymes and negative controls. Altogether 35 enzymes are screened for the presence of SBSs or CBMs and several novel binding sites are identified, including the first SBS ever reported in a cellulase. This work demonstrates that combinations of these methods can be used as a part of routine enzyme characterization to identify new binding sites and advance the study of SBSs and CBMs, allowing them to be detected in the absence of structural data.

Remarkable evolutionary relatedness among the enzymes and proteins from the α-amylase family

The α-amylase is a ubiquitous starch hydrolase catalyzing the cleavage of the α-1,4-glucosidic bonds in an endo-fashion. Various α-amylases originating from different taxonomic sources may differ from each other significantly in their exact substrate preference and product profile. Moreover, it also seems to be clear that at least two different amino acid sequences utilizing two different catalytic machineries have evolved to execute the same α-amylolytic specificity. The two have been classified in the Cabohydrate-Active enZyme database, the CAZy, in the glycoside hydrolase (GH) families GH13 and GH57. While the former and the larger α-amylase family GH13 evidently forms the clan GH-H with the families GH70 and GH77, the latter and the smaller α-amylase family GH57 has only been predicted to maybe define a future clan with the family GH119. Sequences and several tens of enzyme specificities found throughout all three kingdoms in many taxa provide an interesting material for evolutionarily oriented studies that have demonstrated remarkable observations. This review emphasizes just the three of them: (1) a close relatedness between the plant and archaeal α-amylases from the family GH13; (2) a common ancestry in the family GH13 of animal heavy chains of heteromeric amino acid transporter rBAT and 4F2 with the microbial α-glucosidases; and (3) the unique sequence features in the primary structures of amylomaltases from the genus Borrelia from the family GH77. Although the three examples cannot represent an exhaustive list of exceptional topics worth to be interested in, they may demonstrate the importance these enzymes possess in the overall scientific context.

Iminosugar inhibitors of carbohydrate-active enzymes that underpin cereal grain germination and endosperm metabolism

Starch is a major energy store in plants. It provides most of the calories in the human diet and, as a bulk commodity, it is used across broad industry sectors. Starch synthesis and degradation are not fully understood, owing to challenging biochemistry at the liquid/solid interface and relatively limited knowledge about the nature and control of starch degradation in plants. Increased societal and commercial demand for enhanced yield and quality in starch crops requires a better understanding of starch metabolism as a whole. Here we review recent advances in understanding the roles of carbohydrate-active enzymes in starch degradation in cereal grains through complementary chemical and molecular genetics. These approaches have allowed us to start dissecting aspects of starch degradation and the interplay with cell-wall polysaccharide hydrolysis during germination. With a view to improving and diversifying the properties and uses of cereal grains, it is possible that starch degradation may be amenable to manipulation through genetic or chemical intervention at the level of cell wall metabolism, rather than simply in the starch degradation pathway per se.

Underpinning Starch Biology with in vitro Studies on Carbohydrate-Active Enzymes and Biosynthetic Glycomaterials

Starch makes up more than half of the calories in the human diet and is also a valuable bulk commodity that is used across the food, brewing and distilling, medicines and renewable materials sectors. Despite its importance, our understanding of how plants make starch, and what controls the deposition of this insoluble, polymeric, liquid crystalline material, remains rather limited. Advances are hampered by the challenges inherent in analyzing enzymes that operate across the solid-liquid interface. Glyconanotechnology, in the form of glucan-coated sensor chips and metal nanoparticles, present novel opportunities to address this problem. Herein, we review recent developments aimed at the bottom-up generation and self-assembly of starch-like materials, in order to better understand which enzymes are required for starch granule biogenesis and metabolism.

Transcriptional modulation of squalene synthase genes in barley treated with PGPR

Phytosterol contents and food quality of plant produce is directly associated with transcription of gene squalene synthase (SS). In current study, barley plants were treated with different rhizobacterial strains under semi controlled (27 ± 3°C) greenhouse conditions in order to modulate expression of SS gene. Plant samples were analyzed through semi-quantitative PCR to evaluate effect of rhizobacterial application on transcriptional status of SS. Results revealed that among four SS genes (i.e., SSA, SS1, SS2, and SS3), the most expressive gene was SSA; while, SS2 was screened out as the second best induced gene due to Acetobacter aceti. The most efficient bacterial strain which recorded maximum gene expression was A. aceti AC8. Moreover, AC7 was reported as the least efficient bacterial species for inducing SS gene expression. AC8 enhanced the share of SSA and SS2 up to 43 and 31%, respectively. The study also described ribosomal sequence of the most efficient bacterial strain AC8, which was used to determine its phylogenetic relationships with other microbial strains. The study would be helpful to improve quality of plant produce by modulating transcription of SS genes.

Intracellular and cell-to-apoplast compartmentation of carbohydrate metabolism

In most plants, carbohydrates represent the major energy store as well as providing the building blocks for essential structural polymers. Although the major pathways for carbohydrate biosynthesis, degradation, and transport are well characterized, several key steps have only recently been discovered. In addition, several novel minor metabolic routes have been uncovered in the past few years. Here we review current studies of plant carbohydrate metabolism detailing the expanding compendium of functionally characterized transport proteins as well as our deeper comprehension of more minor and conditionally activated metabolic pathways. We additionally explore the pertinent questions that will allow us to enhance our understanding of the response of both major and minor carbohydrate fluxes to changing cellular circumstances.

Structural Basis for the Interconversion of Maltodextrins by MalQ, the Amylomaltase of Escherichia coli

Amylomaltase MalQ is essential for the metabolism of maltose and maltodextrins in Escherichia coli. It catalyzes transglycosylation/disproportionation reactions in which glycosyl or dextrinyl units are transferred among linear maltodextrins of various lengths. To elucidate the molecular basis of transglycosylation by MalQ, we have determined three crystal structures of this enzyme, i.e. the apo-form, its complex with maltose, and an inhibitor complex with the transition state analog acarviosine-glucose-acarbose, at resolutions down to 2.1 Å. MalQ represents the first example of a mesophilic bacterial amylomaltase with known structure and exhibits an N-terminal extension of about 140 residues, in contrast with previously described thermophilic enzymes. This moiety seems unique to amylomaltases from Enterobacteriaceae and folds into two distinct subdomains that associate with different parts of the catalytic core. Intriguingly, the three MalQ crystal structures appear to correspond to distinct states of this enzyme, revealing considerable conformational changes during the catalytic cycle. In particular, the inhibitor complex highlights the requirement of both a 3-OH group and a 4-OH group (or α1-4-glycosidic bond) at the acceptor subsite +1 for the catalytically competent orientation of the acid/base catalyst Glu-496. Using an HPLC-based MalQ enzyme assay, we could demonstrate that the equilibrium concentration of maltodextrin products depends on the length of the initial substrate; with increasing numbers of glycosidic bonds, less glucose is formed. Thus, both structural and enzymatic data are consistent with the extremely low hydrolysis rates observed for amylomaltases and underline the importance of MalQ for the metabolism of maltodextrins in E. coli.

In silico analysis of family GH77 with focus on amylomaltases from borreliae and disproportionating enzymes DPE2 from plants and bacteria

The CAZy glycoside hydrolase (GH) family GH77 is a monospecific family containing 4-α-glucanotransferases that if from prokaryotes are known as amylomaltases and if from plants including algae are known as disproportionating enzymes (DPE). The family GH77 is a member of the α-amylase clan GH-H. The main difference discriminating a GH77 4-α-glucanotransferase from the main GH13 α-amylase family members is the lack of domain C succeeding the catalytic (β/α)8-barrel. Of more than 2400 GH77 members, bacterial amylomaltases clearly dominate with more than 2300 sequences; the rest being approximately equally represented by Archaea and Eucarya. The main goal of the present study was to deliver a detailed bioinformatics study of family GH77 (416 collected sequences) focused on amylomaltases from borreliae (containing unique sequence substitutions in functionally important positions) and plant DPE2 representatives (possessing an insert of ~140 residues between catalytic nucleophile and proton donor). The in silico analysis reveals that within the genus of Borrelia a gradual evolutionary transition from typical bacterial Thermus-like amylomaltases may exist to family-GH77 amylomaltase versions that currently possess progressively mutated the most important and otherwise invariantly conserved positions. With regard to plant DPE2, a large group of bacterial amylomaltases represented by the amylomaltase from Escherichia coli with a longer N-terminus was identified as a probable intermediary connection between Thermus-like and DPE2-like (existing also among bacteria) family GH77 members. The presented results concerning both groups, i.e. amylomaltases from borreliae and plant DPE2 representatives (with their bacterial counterpart), may thus indicate the direction for future experimental studies.

Pathway thermodynamics highlights kinetic obstacles in central metabolism

In metabolism research, thermodynamics is usually used to determine the directionality of a reaction or the feasibility of a pathway. However, the relationship between thermodynamic potentials and fluxes is not limited to questions of directionality: thermodynamics also affects the kinetics of reactions through the flux-force relationship, which states that the logarithm of the ratio between the forward and reverse fluxes is directly proportional to the change in Gibbs energy due to a reaction (Δ_rG′). Accordingly, if an enzyme catalyzes a reaction with a Δ_rG′ of -5.7 kJ/mol then the forward flux will be roughly ten times the reverse flux. As Δ_rG′ approaches equilibrium (Δ_rG′ = 0 kJ/mol), exponentially more enzyme counterproductively catalyzes the reverse reaction, reducing the net rate at which the reaction proceeds. Thus, the enzyme level required to achieve a given flux increases dramatically near equilibrium. Here, we develop a framework for quantifying the degree to which pathways suffer these thermodynamic limitations on flux. For each pathway, we calculate a single thermodynamically-derived metric (the Max-min Driving Force, MDF), which enables objective ranking of pathways by the degree to which their flux is constrained by low thermodynamic driving force. Our framework accounts for the effect of pH, ionic strength and metabolite concentration ranges and allows us to quantify how alterations to the pathway structure affect the pathway's thermodynamics. Applying this methodology to pathways of central metabolism sheds light on some of their features, including metabolic bypasses (e.g., fermentation pathways bypassing substrate-level phosphorylation), substrate channeling (e.g., of oxaloacetate from malate dehydrogenase to citrate synthase), and use of alternative cofactors (e.g., quinone as an electron acceptor instead of NAD). The methods presented here place another arrow in metabolic engineers' quiver, providing a simple means of evaluating the thermodynamic and kinetic quality of different pathway chemistries that produce the same molecules.

Given data about enzyme kinetics and reaction thermodynamics, traditional metabolic control analysis (MCA) can pinpoint the enzymes whose expression will have the largest effect on steady-state flux through the pathway. These analyses can aid experimentalists in tuning enzyme expression levels along a metabolic pathway. In this work, we offer a framework that is complementary to MCA. Rather than focusing on the relationship between enzyme levels and pathway flux, we examine a pathway's stoichiometry and thermodynamics and ask whether it is likely to support high flux in cellular conditions. Our framework calculates a single thermodynamically-derived metric (the MDF) for each pathway, which is convenient for selecting the promising pathways from a large collection. This approach has several advantages. First, enzyme kinetic properties are laborious to measure and differ between organisms and isozymes, but no kinetic data is required to calculate the MDF. Second, as our framework accounts for pH, ionic strength and allowed concentration ranges, it is simple to model the effect of these parameters on the MDF. Finally, as it can be difficult to control the exact expression level of enzymes within cells, the MDF helps identify alternative pathways that are less sensitive to the levels of their constituent enzymes.