BayFlux: A Bayesian method to quantify metabolic Fluxes and their uncertainty at the genome scale.

Metabolic fluxes, the number of metabolites traversing each biochemical reaction in a cell per unit time, are crucial for assessing and understanding cell function. 13C Metabolic Flux Analysis (13C MFA) is considered to be the gold standard for measuring metabolic fluxes. 13C MFA typically works by leveraging extracellular exchange fluxes as well as data from 13C labeling experiments to calculate the flux profile which best fit the data for a small, central carbon, metabolic model. However, the nonlinear nature of the 13C MFA fitting procedure means that several flux profiles fit the experimental data within the experimental error, and traditional optimization methods offer only a partial or skewed picture, especially in "non-gaussian" situations where multiple very distinct flux regions fit the data equally well. Here, we present a method for flux space sampling through Bayesian inference (BayFlux), that identifies the full distribution of fluxes compatible with experimental data for a comprehensive genome-scale model. This Bayesian approach allows us to accurately quantify uncertainty in calculated fluxes. We also find that, surprisingly, the genome-scale model of metabolism produces narrower flux distributions (reduced uncertainty) than the small core metabolic models traditionally used in 13C MFA. The different results for some reactions when using genome-scale models vs core metabolic models advise caution in assuming strong inferences from 13C MFA since the results may depend significantly on the completeness of the model used. Based on BayFlux, we developed and evaluated novel methods (P-13C MOMA and P-13C ROOM) to predict the biological results of a gene knockout, that improve on the traditional MOMA and ROOM methods by quantifying prediction uncertainty.

Scalable and automated CRISPR-based strain engineering using droplet microfluidics.

We present a droplet-based microfluidic system that enables CRISPR-based gene editing and high-throughput screening on a chip. The microfluidic device contains a 10 × 10 element array, and each element contains sets of electrodes for two electric field-actuated operations: electrowetting for merging droplets to mix reagents and electroporation for transformation. This device can perform up to 100 genetic modification reactions in parallel, providing a scalable platform for generating the large number of engineered strains required for the combinatorial optimization of genetic pathways and predictable bioengineering. We demonstrate the system's capabilities through the CRISPR-based engineering of two test cases: (1) disruption of the function of the enzyme galactokinase (galK) in E. coli and (2) targeted engineering of the glutamine synthetase gene (glnA) and the blue-pigment synthetase gene (bpsA) to improve indigoidine production in E. coli.

Genome-Scale 13C Fluxomics Modeling for Metabolic Engineering of Saccharomyces cerevisiae.

Synthetic biology is a rapidly developing field that pursues the application of engineering principles and development approaches to biological engineering. Synthetic biology is poised to change the way biology is practiced, and has important practical applications: for example, building genetically engineered organisms to produce biofuels, medicines, and other chemicals. Traditionally, synthetic biology has focused on manipulating a few genes (e.g., in a single pathway or genetic circuit), but its combination with systems biology holds the promise of creating new cellular architectures and constructing complex biological systems from the ground up. Enabling this merge of synthetic and systems biology will require greater predictive capability for modeling the behavior of cellular systems, and more comprehensive data sets for building and calibrating these models. The so-called "-omics" data sets can now be generated via high throughput techniques in the form of genomic, proteomic, transcriptomic, and metabolomic information on the engineered biological system. Of particular interest with respect to the engineering of microbes capable of producing biofuels and other chemicals economically and at scale are metabolomic datasets, and their insights into intracellular metabolic fluxes. Metabolic fluxes provide a rapid and easy to understand picture of how carbon and energy flow throughout the cell. Here, we present a detailed guide to performing metabolic flux analysis and modeling using the open source JBEI Quantitative Metabolic Modeling (jQMM) library. This library allows the user to transform metabolomics data in the form of isotope labeling data from a 13C labeling experiment into a determination of cellular fluxes that can be used to develop genetic engineering strategies for metabolic engineering.The jQMM library presents a complete toolbox for performing a range of different tasks of interest in metabolic engineering. Various different types of flux analysis and modeling can be performed such as flux balance analysis, 13C metabolic flux analysis, and two-scale 13C metabolic flux analysis (2S-13C MFA). 2S-13C MFA is a novel method that determines genome-scale fluxes without the need of every single carbon transition in the metabolic network. In addition to several other capabilities, the jQMM library can make model based predictions for how various genetic engineering strategies can be incorporated toward bioengineering goals: it can predict the effects of reaction knockouts on metabolism using both the MoMA and ROOM methodologies. In this chapter, we will illustrate the use of the jQMM library through a step-by-step demonstration of flux determination and knockout prediction in a complex eukaryotic model organism: Saccharomyces cerevisiae (S. cerevisiae). Included with this chapter is a digital Jupyter Notebook file that provides a computable appendix showing a self-contained example of jQMM usage, which can be changed to fit the user's specific needs. As an open source software project, users can modify and extend the code base to make improvements at will, allowing them to share their development work and contribute back to the jQMM modeling community.

Constraining Genome-Scale Models to Represent the Bow Tie Structure of Metabolism for 13C Metabolic Flux Analysis.

Determination of internal metabolic fluxes is crucial for fundamental and applied biology because they map how carbon and electrons flow through metabolism to enable cell function. 13 C Metabolic Flux Analysis ( 13 C MFA) and Two-Scale 13 C Metabolic Flux Analysis (2S- 13 C MFA) are two techniques used to determine such fluxes. Both operate on the simplifying approximation that metabolic flux from peripheral metabolism into central "core" carbon metabolism is minimal, and can be omitted when modeling isotopic labeling in core metabolism. The validity of this "two-scale" or "bow tie" approximation is supported both by the ability to accurately model experimental isotopic labeling data, and by experimentally verified metabolic engineering predictions using these methods. However, the boundaries of core metabolism that satisfy this approximation can vary across species, and across cell culture conditions. Here, we present a set of algorithms that (1) systematically calculate flux bounds for any specified "core" of a genome-scale model so as to satisfy the bow tie approximation and (2) automatically identify an updated set of core reactions that can satisfy this approximation more efficiently. First, we leverage linear programming to simultaneously identify the lowest fluxes from peripheral metabolism into core metabolism compatible with the observed growth rate and extracellular metabolite exchange fluxes. Second, we use Simulated Annealing to identify an updated set of core reactions that allow for a minimum of fluxes into core metabolism to satisfy these experimental constraints. Together, these methods accelerate and automate the identification of a biologically reasonable set of core reactions for use with 13 C MFA or 2S- 13 C MFA, as well as provide for a substantially lower set of flux bounds for fluxes into the core as compared with previous methods. We provide an open source Python implementation of these algorithms at https://github.com/JBEI/limitfluxtocore.

Two-Scale 13C Metabolic Flux Analysis for Metabolic Engineering.

Accelerating the Design-Build-Test-Learn (DBTL) cycle in synthetic biology is critical to achieving rapid and facile bioengineering of organisms for the production of, e.g., biofuels and other chemicals. The Learn phase involves using data obtained from the Test phase to inform the next Design phase. As part of the Learn phase, mathematical models of metabolic fluxes give a mechanistic level of comprehension to cellular metabolism, isolating the principle drivers of metabolic behavior from the peripheral ones, and directing future experimental designs and engineering methodologies. Furthermore, the measurement of intracellular metabolic fluxes is specifically noteworthy as providing a rapid and easy-to-understand picture of how carbon and energy flow throughout the cell. Here, we present a detailed guide to performing metabolic flux analysis in the Learn phase of the DBTL cycle, where we show how one can take the isotope labeling data from a 13C labeling experiment and immediately turn it into a determination of cellular fluxes that points in the direction of genetic engineering strategies that will advance the metabolic engineering process.For our modeling purposes we use the Joint BioEnergy Institute (JBEI) Quantitative Metabolic Modeling (jQMM) library, which provides an open-source, python-based framework for modeling internal metabolic fluxes and making actionable predictions on how to modify cellular metabolism for specific bioengineering goals. It presents a complete toolbox for performing different types of flux analysis such as Flux Balance Analysis, 13C Metabolic Flux Analysis, and it introduces the capability to use 13C labeling experimental data to constrain comprehensive genome-scale models through a technique called two-scale 13C Metabolic Flux Analysis (2S-13C MFA) [1]. In addition to several other capabilities, the jQMM is also able to predict the effects of knockouts using the MoMA and ROOM methodologies. The use of the jQMM library is illustrated through a step-by-step demonstration, which is also contained in a digital Jupyter Notebook format that enhances reproducibility and provides the capability to be adopted to the user's specific needs. As an open-source software project, users can modify and extend the code base and make improvements at will, providing a base for future modeling efforts.

The Experiment Data Depot: A Web-Based Software Tool for Biological Experimental Data Storage, Sharing, and Visualization.

Although recent advances in synthetic biology allow us to produce biological designs more efficiently than ever, our ability to predict the end result of these designs is still nascent. Predictive models require large amounts of high-quality data to be parametrized and tested, which are not generally available. Here, we present the Experiment Data Depot (EDD), an online tool designed as a repository of experimental data and metadata. EDD provides a convenient way to upload a variety of data types, visualize these data, and export them in a standardized fashion for use with predictive algorithms. In this paper, we describe EDD and showcase its utility for three different use cases: storage of characterized synthetic biology parts, leveraging proteomics data to improve biofuel yield, and the use of extracellular metabolite concentrations to predict intracellular metabolic fluxes.

The JBEI quantitative metabolic modeling library (jQMM): a python library for modeling microbial metabolism.

BACKGROUND: Modeling of microbial metabolism is a topic of growing importance in biotechnology. Mathematical modeling helps provide a mechanistic understanding for the studied process, separating the main drivers from the circumstantial ones, bounding the outcomes of experiments and guiding engineering approaches. Among different modeling schemes, the quantification of intracellular metabolic fluxes (i.e. the rate of each reaction in cellular metabolism) is of particular interest for metabolic engineering because it describes how carbon and energy flow throughout the cell. In addition to flux analysis, new methods for the effective use of the ever more readily available and abundant -omics data (i.e. transcriptomics, proteomics and metabolomics) are urgently needed. RESULTS: The jQMM library presented here provides an open-source, Python-based framework for modeling internal metabolic fluxes and leveraging other -omics data for the scientific study of cellular metabolism and bioengineering purposes. Firstly, it presents a complete toolbox for simultaneously performing two different types of flux analysis that are typically disjoint: Flux Balance Analysis and 13C Metabolic Flux Analysis. Moreover, it introduces the capability to use 13C labeling experimental data to constrain comprehensive genome-scale models through a technique called two-scale 13C Metabolic Flux Analysis (2S-13C MFA). In addition, the library includes a demonstration of a method that uses proteomics data to produce actionable insights to increase biofuel production. Finally, the use of the jQMM library is illustrated through the addition of several Jupyter notebook demonstration files that enhance reproducibility and provide the capability to be adapted to the user's specific needs. CONCLUSIONS: jQMM will facilitate the design and metabolic engineering of organisms for biofuels and other chemicals, as well as investigations of cellular metabolism and leveraging -omics data. As an open source software project, we hope it will attract additions from the community and grow with the rapidly changing field of metabolic engineering.

Cooperative Interactions between Different Classes of Disordered Proteins Play a Functional Role in the Nuclear Pore Complex of Baker's Yeast.

Nucleocytoplasmic transport is highly selective, efficient, and is regulated by a poorly understood mechanism involving hundreds of disordered FG nucleoporin proteins (FG nups) lining the inside wall of the nuclear pore complex (NPC). Previous research has concluded that FG nups in Baker's yeast (S. cerevisiae) are present in a bimodal distribution, with the "Forest Model" classifying FG nups as either di-block polymer like "trees" or single-block polymer like "shrubs". Using a combination of coarse-grained modeling and polymer brush modeling, the function of the di-block FG nups has previously been hypothesized in the Di-block Copolymer Brush Gate (DCBG) model to form a higher-order polymer brush architecture which can open and close to regulate transport across the NPC. In this manuscript we work to extend the original DCBG model by first performing coarse grained simulations of the single-block FG nups which confirm that they have a single block polymer structure rather than the di-block structure of tree nups. Our molecular simulations also demonstrate that these single-block FG nups are likely cohesive, compact, collapsed coil polymers, implying that these FG nups are generally localized to their grafting location within the NPC. We find that adding a layer of single-block FG nups to the DCBG model increases the range of cargo sizes which are able to translocate the pore through a cooperative effect involving single-block and di-block FG nups. This effect can explain the puzzling connection between single-block FG nup deletion mutants in S. cerevisiae and the resulting failure of certain large cargo transport through the NPC. Facilitation of large cargo transport via single-block and di-block FG nup cooperativity in the nuclear pore could provide a model mechanism for designing future biomimetic pores of greater applicability.

13C Metabolic Flux Analysis for Systematic Metabolic Engineering of S. cerevisiae for Overproduction of Fatty Acids.

Efficient redirection of microbial metabolism into the abundant production of desired bioproducts remains non-trivial. Here, we used flux-based modeling approaches to improve yields of fatty acids in Saccharomyces cerevisiae. We combined 13C labeling data with comprehensive genome-scale models to shed light onto microbial metabolism and improve metabolic engineering efforts. We concentrated on studying the balance of acetyl-CoA, a precursor metabolite for the biosynthesis of fatty acids. A genome-wide acetyl-CoA balance study showed ATP citrate lyase from Yarrowia lipolytica as a robust source of cytoplasmic acetyl-CoA and malate synthase as a desirable target for downregulation in terms of acetyl-CoA consumption. These genetic modifications were applied to S. cerevisiae WRY2, a strain that is capable of producing 460 mg/L of free fatty acids. With the addition of ATP citrate lyase and downregulation of malate synthase, the engineered strain produced 26% more free fatty acids. Further increases in free fatty acid production of 33% were obtained by knocking out the cytoplasmic glycerol-3-phosphate dehydrogenase, which flux analysis had shown was competing for carbon flux upstream with the carbon flux through the acetyl-CoA production pathway in the cytoplasm. In total, the genetic interventions applied in this work increased fatty acid production by ~70%.

Cytoskeletal Network Morphology Regulates Intracellular Transport Dynamics.

Intracellular transport is essential for maintaining proper cellular function in most eukaryotic cells, with perturbations in active transport resulting in several types of disease. Efficient delivery of critical cargos to specific locations is accomplished through a combination of passive diffusion and active transport by molecular motors that ballistically move along a network of cytoskeletal filaments. Although motor-based transport is known to be necessary to overcome cytoplasmic crowding and the limited range of diffusion within reasonable timescales, the topological features of the cytoskeletal network that regulate transport efficiency and robustness have not been established. Using a continuum diffusion model, we observed that the time required for cellular transport was minimized when the network was localized near the nucleus. In simulations that explicitly incorporated network spatial architectures, total filament mass was the primary driver of network transit times. However, filament traps that redirect cargo back to the nucleus caused large variations in network transport. Filament polarity was more important than filament orientation in reducing average transit times, and transport properties were optimized in networks with intermediate motor on and off rates. Our results provide important insights into the functional constraints on intracellular transport under which cells have evolved cytoskeletal structures, and have potential applications for enhancing reactions in biomimetic systems through rational transport network design.

Cooperative protofilament switching emerges from inter-motor interference in multiple-motor transport.

Within living cells, the transport of cargo is accomplished by groups of molecular motors. Such collective transport could utilize mechanisms which emerge from inter-motor interactions in ways that are yet to be fully understood. Here we combined experimental measurements of two-kinesin transport with a theoretical framework to investigate the functional ramifications of inter-motor interactions on individual motor function and collective cargo transport. In contrast to kinesin's low sidestepping frequency when present as a single motor, with exactly two kinesins per cargo, we observed substantial motion perpendicular to the microtubule. Our model captures a surface-associated mode of kinesin, which is only accessible via inter-motor interference in groups, in which kinesin diffuses along the microtubule surface and rapidly "hops" between protofilaments without dissociating from the microtubule. Critically, each kinesin transitions dynamically between the active stepping mode and this weak surface-associated mode enhancing local exploration of the microtubule surface, possibly enabling cellular cargos to overcome macromolecular crowding and to navigate obstacles along microtubule tracks without sacrificing overall travel distance.

Nuclear pore complex protein sequences determine overall copolymer brush structure and function.

The transport of cargo across the nuclear membrane is highly selective and accomplished by a poorly understood mechanism involving hundreds of nucleoporins lining the inside of the nuclear pore complex (NPC). Currently, there is no clear picture of the overall structure formed by this collection of proteins within the pore, primarily due to their disordered nature. We perform coarse-grained simulations of both individual nucleoporins and grafted rings of nups mimicking the in vivo geometry of the NPC and supplement this with polymer brush modeling. Our results indicate that different regions or blocks of an individual NPC protein can have distinctly different forms of disorder and that this property appears to be a conserved functional feature. Furthermore, this block structure at the individual protein level is critical to the formation of a unique higher-order polymer brush architecture that can exist in distinct morphologies depending on the effective interaction energy between the phenylalanine glycine (FG) domains of different nups. Because the interactions between FG domains may be modulated by certain forms of transport factors, our results indicate that transitions between brush morphologies could play an important role in regulating transport across the NPC, suggesting novel forms of gated transport across membrane pores with wide biomimetic applicability.

Physical motif clustering within intrinsically disordered nucleoporin sequences reveals universal functional features.

Bioinformatics of disordered proteins is especially challenging given high mutation rates for homologous proteins and that functionality may not be strongly related to sequence. Here we have performed a novel bioinformatic analysis, based on the spatial clustering of physically relevant features such as binding motifs and charges within disordered proteins, on thousands of Nuclear Pore Complex (NPC) FG motif containing proteins (FG nups). The biophysical mechanism by which FG nups regulate nucleocytoplasmic transport has remained elusive. Our analysis revealed a set of highly conserved spatial features in the sequence structure of individual FG nups, such as the separation, localization, and ordering of FG motifs and charged residues along the protein chain. These functionally conserved features provide insight into the particular biophysical mechanisms responsible for regulation of nucleocytoplasmic traffic in the NPC, strongly constraining current models. Additionally this method allows us to identify potentially functionally analogous disordered proteins across distantly related species.