Combinatorial in Vitro and in Silico Approach To Describe Shear-Force Dependent Uptake of Nanoparticles in Microfluidic Vascular Models.

In the present work, we combine experimental and computational methods to define the critical shear stress as an alternative parameter for nanotoxicological and nanomedical evaluations using an in vitro microfluidic vascular model. We demonstrate that our complementary in vitro and in silico approach is well suited to assess the fluid flow velocity above which clathrin-mediated (active) nanoparticle uptake per cell decreases drastically although higher numbers of nanoparticles per cell are introduced. Results of our study revealed a critical shear stress of 1.8 dyn/cm2, where maximum active cellular nanoparticle uptake took place, followed by a 70% decrease in uptake of 249 nm nanoparticles at 10 dyn/cm2, respectively. The observed nonlinear relationship between flow velocity and nanoparticle uptake strongly suggests that fluid mechanical forces also need to be considered in order to predict potential in vivo distribution, bioaccumulation, and clearance of nanomaterials and novel nanodrugs.

Designing Optimized Production Hosts by Metabolic Modeling.

Many of the complex and expensive production steps in the chemical industry are readily available in living cells. In order to overcome the metabolic limits of these cells, the optimal genetic intervention strategies can be computed by the use of metabolic modeling. Elementary flux mode analysis (EFMA) is an ideal tool for this task, as it does not require defining a cellular objective function. We present two EFMA-based methods to optimize production hosts: (1) the standard approach that can only be used for small and medium scale metabolic networks and (2) the advanced dual system approach that can be utilized to directly compute intervention strategies in a genome-scale metabolic model.

From elementary flux modes to elementary flux vectors: Metabolic pathway analysis with arbitrary linear flux constraints.

Elementary flux modes (EFMs) emerged as a formal concept to describe metabolic pathways and have become an established tool for constraint-based modeling and metabolic network analysis. EFMs are characteristic (support-minimal) vectors of the flux cone that contains all feasible steady-state flux vectors of a given metabolic network. EFMs account for (homogeneous) linear constraints arising from reaction irreversibilities and the assumption of steady state; however, other (inhomogeneous) linear constraints, such as minimal and maximal reaction rates frequently used by other constraint-based techniques (such as flux balance analysis [FBA]), cannot be directly integrated. These additional constraints further restrict the space of feasible flux vectors and turn the flux cone into a general flux polyhedron in which the concept of EFMs is not directly applicable anymore. For this reason, there has been a conceptual gap between EFM-based (pathway) analysis methods and linear optimization (FBA) techniques, as they operate on different geometric objects. One approach to overcome these limitations was proposed ten years ago and is based on the concept of elementary flux vectors (EFVs). Only recently has the community started to recognize the potential of EFVs for metabolic network analysis. In fact, EFVs exactly represent the conceptual development required to generalize the idea of EFMs from flux cones to flux polyhedra. This work aims to present a concise theoretical and practical introduction to EFVs that is accessible to a broad audience. We highlight the close relationship between EFMs and EFVs and demonstrate that almost all applications of EFMs (in flux cones) are possible for EFVs (in flux polyhedra) as well. In fact, certain properties can only be studied with EFVs. Thus, we conclude that EFVs provide a powerful and unifying framework for constraint-based modeling of metabolic networks.

Optimal knockout strategies in genome-scale metabolic networks using particle swarm optimization.

BACKGROUND: Knockout strategies, particularly the concept of constrained minimal cut sets (cMCSs), are an important part of the arsenal of tools used in manipulating metabolic networks. Given a specific design, cMCSs can be calculated even in genome-scale networks. We would however like to find not only the optimal intervention strategy for a given design but the best possible design too. Our solution (PSOMCS) is to use particle swarm optimization (PSO) along with the direct calculation of cMCSs from the stoichiometric matrix to obtain optimal designs satisfying multiple objectives. RESULTS: To illustrate the working of PSOMCS, we apply it to a toy network. Next we show its superiority by comparing its performance against other comparable methods on a medium sized E. coli core metabolic network. PSOMCS not only finds solutions comparable to previously published results but also it is orders of magnitude faster. Finally, we use PSOMCS to predict knockouts satisfying multiple objectives in a genome-scale metabolic model of E. coli and compare it with OptKnock and RobustKnock. CONCLUSIONS: PSOMCS finds competitive knockout strategies and designs compared to other current methods and is in some cases significantly faster. It can be used in identifying knockouts which will force optimal desired behaviors in large and genome scale metabolic networks. It will be even more useful as larger metabolic models of industrially relevant organisms become available.

Microfluidic Migration and Wound Healing Assay Based on Mechanically Induced Injuries of Defined and Highly Reproducible Areas.

All cell migration and wound healing assays are based on the inherent ability of adherent cells to move into adjacent cell-free areas, thus providing information on cell culture viability, cellular mechanisms and multicellular movements. Despite their widespread use for toxicological screening, biomedical research and pharmaceutical studies, to date no satisfactory technological solutions are available for the automated, miniaturized and integrated induction of defined wound areas. To bridge this technological gap, we have developed a lab-on-a-chip capable of mechanically inducing circular cell-free areas within confluent cell layers. The microdevices were fabricated using off-stoichiometric thiol-ene-epoxy (OSTEMER) polymer resulting in hard-polymer devices that are robust, cost-effective and disposable. We show that the pneumatically controlled membrane deflection/compression method not only generates highly reproducible (RSD 4%) injuries but also allows for repeated wounding in microfluidic environments. Performance analysis demonstrated that applied surface coating remains intact even after multiple wounding, while cell debris is simultaneously removed using laminar flow conditions. Furthermore, only a few injured cells were found along the edge of the circular cell-free areas, thus allowing reliable and reproducible cell migration of a wide range of surface sensitive anchorage dependent cell types. Practical application is demonstrated by investigating healing progression and endothelial cell migration in the absence and presence of an inflammatory cytokine (TNF-α) and a well-known cell proliferation inhibitor (mitomycin-C).

Trend analysis of performance parameters of pre-packed columns for protein chromatography over a time span of ten years.

Pre-packed small scale chromatography columns are increasingly used for process development, for determination of design space in bioprocess development, and for post-licence process verifications. The packing quality of 30,000 pre-packed columns delivered to customers over a period 10 years has been analyzed by advanced statistical tools. First, the data were extracted and checked for inconsistencies, and then were tabulated and made ready for statistical processing using the programming language Perl (https://www.perl.org/) and the statistical computing environment R (https://www.r-project.org/). Reduced HETP and asymmetry were plotted over time to obtain a trend of packing quality over 10 years. The obtained data were used as a visualized coefficient of variation analysis (VCVA), a process that has often been applied in other industries such as semiconductor manufacturing. A typical fluctuation of reduced HETP was seen. A Tsunami effect in manufacturing, the effect of propagation of manufacturing deviations leading to out-of-specification products, was not observed with these pre-packed columns. Principal component analysis (PCA) showed that all packing materials cluster. Our data analysis showed that the current commercially available chromatography media used for biopharmaceutical manufacturing can be reproducibly and uniformly packed in polymer-based chromatography columns, which are designed for ready-to-use purposes. Although the number of packed columns has quadrupled over one decade the packing quality has remained stable.

Which sets of elementary flux modes form thermodynamically feasible flux distributions?

Elementary flux modes (EFMs) are non-decomposable steady-state fluxes through metabolic networks. Every possible flux through a network can be described as a superposition of EFMs. The definition of EFMs is based on the stoichiometry of the network, and it has been shown previously that not all EFMs are thermodynamically feasible. These infeasible EFMs cannot contribute to a biologically meaningful flux distribution. In this work, we show that a set of thermodynamically feasible EFMs need not be thermodynamically consistent. We use first principles of thermodynamics to define the feasibility of a flux distribution and present a method to compute the largest thermodynamically consistent sets (LTCSs) of EFMs. An LTCS contains the maximum number of EFMs that can be combined to form a thermodynamically feasible flux distribution. As a case study we analyze all LTCSs found in Escherichia coli when grown on glucose and show that only one LTCS shows the required phenotypical properties. Using our method, we find that in our E. coli model < 10% of all EFMs are thermodynamically relevant.

Exact quantification of cellular robustness in genome-scale metabolic networks.

MOTIVATION: Robustness, the ability of biological networks to uphold their functionality in spite of perturbations, is a key characteristic of all living systems. Although several theoretical approaches have been developed to formalize robustness, it still eludes an exact quantification. Here, we present a rigorous and quantitative approach for the structural robustness of metabolic networks by measuring their ability to tolerate random reaction (or gene) knockouts. RESULTS: In analogy to reliability theory, based on an explicit consideration of all possible knockout sets, we exactly quantify the probability of failure for a given network function (e.g. growth). This measure can be computed if the network's minimal cut sets (MSCs) are known. We show that even in genome-scale metabolic networks the probability of (network) failure can be reliably estimated from MSCs with lowest cardinalities. We demonstrate the applicability of our theory by analyzing the structural robustness of multiple Enterobacteriaceae and Blattibacteriaceae and show a dramatically low structural robustness for the latter. We find that structural robustness develops from the ability to proliferate in multiple growth environments consistent with experimentally found knowledge. CONCLUSION: The probability of (network) failure provides thus a reliable and easily computable measure of structural robustness and redundancy in (genome-scale) metabolic networks. AVAILABILITY AND IMPLEMENTATION: Source code is available under the GNU General Public License at https://github.com/mpgerstl/networkRobustnessToolbox CONTACT: juergen.zanghellini@boku.ac.at SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

ICT: isotope correction toolbox.

SUMMARY: Isotope tracer experiments are an invaluable technique to analyze and study the metabolism of biological systems. However, isotope labeling experiments are often affected by naturally abundant isotopes especially in cases where mass spectrometric methods make use of derivatization. The correction of these additive interferences--in particular for complex isotopic systems--is numerically challenging and still an emerging field of research. When positional information is generated via collision-induced dissociation, even more complex calculations for isotopic interference correction are necessary. So far, no freely available tools can handle tandem mass spectrometry data. We present isotope correction toolbox, a program that corrects tandem mass isotopomer data from tandem mass spectrometry experiments. Isotope correction toolbox is written in the multi-platform programming language Perl and, therefore, can be used on all commonly available computer platforms. AVAILABILITY AND IMPLEMENTATION: Source code and documentation can be freely obtained under the Artistic License or the GNU General Public License from: https://github.com/jungreuc/isotope_correction_toolbox/ CONTACT: {christian.jungreuthmayer@boku.ac.at,juergen.zanghellini@boku.ac.at} SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The 3D pore structure and fluid dynamics simulation of macroporous monoliths: High permeability due to alternating channel width.

Polymethacrylate-based monoliths have excellent flow properties. Flow in the wide channel interconnected with narrow channels is theoretically assumed to account for favorable permeability. Monoliths were cut into 898 slices in 50nm distances and visualized by serial block face scanning electron microscopy (SBEM). A 3D structure was reconstructed and used for the calculation of flow profiles within the monolith and for calculation of pressure drop and permeability by computational fluid dynamics (CFD). The calculated and measured permeabilities showed good agreement. Small channels clearly flowed into wide and wide into small channels in a repetitive manner which supported the hypothesis describing the favorable flow properties of these materials. This alternating property is also reflected in the streamline velocity which fluctuated. These findings were corroborated by artificial monoliths which were composed of regular (interconnected) cells where narrow cells followed wide cells. In the real monolith and the artificial monoliths with interconnected flow channels similar velocity fluctuations could be observed. A two phase flow simulation showed a lateral velocity component, which may contribute to the transport of molecules to the monolith wall. Our study showed that the interconnection of small and wide pores is responsible for the excellent pressure flow properties. This study is also a guide for further design of continuous porous materials to achieve good flow properties.

Designing minimal microbial strains of desired functionality using a genetic algorithm.

BACKGROUND: The rational, in silico prediction of gene-knockouts to turn organisms into efficient cell factories is an essential and computationally challenging task in metabolic engineering. Elementary flux mode analysis in combination with constraint minimal cut sets is a particularly powerful method to identify optimal engineering targets, which will force an organism into the desired metabolic state. Given an engineering objective, it is theoretically possible, although computationally impractical, to find the best minimal intervention strategies. RESULTS: We developed a genetic algorithm (GA-MCS) to quickly find many (near) optimal intervention strategies while overcoming the above mentioned computational burden. We tested our algorithm on Escherichia coli metabolic networks of three different sizes to find intervention strategies satisfying three different engineering objectives. CONCLUSIONS: We show that GA-MCS finds all practically relevant targets for any (non)-linear engineering objective. Our algorithm also found solutions comparable to previously published results. We show that for large networks optimal solutions are found within a fraction of the time used for a complete enumeration.

Avoiding the Enumeration of Infeasible Elementary Flux Modes by Including Transcriptional Regulatory Rules in the Enumeration Process Saves Computational Costs.

Despite the significant progress made in recent years, the computation of the complete set of elementary flux modes of large or even genome-scale metabolic networks is still impossible. We introduce a novel approach to speed up the calculation of elementary flux modes by including transcriptional regulatory information into the analysis of metabolic networks. Taking into account gene regulation dramatically reduces the solution space and allows the presented algorithm to constantly eliminate biologically infeasible modes at an early stage of the computation procedure. Thereby, computational costs, such as runtime, memory usage, and disk space, are extremely reduced. Moreover, we show that the application of transcriptional rules identifies non-trivial system-wide effects on metabolism. Using the presented algorithm pushes the size of metabolic networks that can be studied by elementary flux modes to new and much higher limits without the loss of predictive quality. This makes unbiased, system-wide predictions in large scale metabolic networks possible without resorting to any optimization principle.

Predicting genetic engineering targets with Elementary Flux Mode Analysis: a review of four current methods.

Elementary flux modes (EFMs) are a well-established tool in metabolic modeling. EFMs are minimal, feasible, steady state pathways through a metabolic network. They are used in various approaches to predict targets for genetic interventions in order to increase production of a molecule of interest via a host cell. Here we give an introduction to the concept of EFMs, present an overview of four methods which use EFMs in order to predict engineering targets and lastly use a toy model and a small-scale metabolic model to demonstrate and compare the capabilities of these methods.

Metabolomics integrated elementary flux mode analysis in large metabolic networks.

Elementary flux modes (EFMs) are non-decomposable steady-state pathways in metabolic networks. They characterize phenotypes, quantify robustness or identify engineering targets. An EFM analysis (EFMA) is currently restricted to medium-scale models, as the number of EFMs explodes with the network's size. However, many topologically feasible EFMs are biologically irrelevant. We present thermodynamic EFMA (tEFMA), which calculates only the small(er) subset of thermodynamically feasible EFMs. We integrate network embedded thermodynamics into EFMA and show that we can use the metabolome to identify and remove thermodynamically infeasible EFMs during an EFMA without losing biologically relevant EFMs. Calculating only the thermodynamically feasible EFMs strongly reduces memory consumption and program runtime, allowing the analysis of larger networks. We apply tEFMA to study the central carbon metabolism of E. coli and find that up to 80% of its EFMs are thermodynamically infeasible. Moreover, we identify glutamate dehydrogenase as a bottleneck, when E. coli is grown on glucose and explain its inactivity as a consequence of network embedded thermodynamics. We implemented tEFMA as a Java package which is available for download at https://github.com/mpgerstl/tEFMA.

tEFMA: computing thermodynamically feasible elementary flux modes in metabolic networks.

UNLABELLED: : Elementary flux modes (EFMs) are important structural tools for the analysis of metabolic networks. It is known that many topologically feasible EFMs are biologically irrelevant. Therefore, tools are needed to find the relevant ones. We present thermodynamic tEFM analysis (tEFMA) which uses the cellular metabolome to avoid the enumeration of thermodynamically infeasible EFMs. Specifically, given a metabolic network and a not necessarily complete metabolome, tEFMA efficiently returns the full set of thermodynamically feasible EFMs consistent with the metabolome. Compared with standard approaches, tEFMA strongly reduces the memory consumption and the overall runtime. Thus tEFMA provides a new way to analyze unbiasedly hitherto inaccessible large-scale metabolic networks. AVAILABILITY AND IMPLEMENTATION: https://github.com/mpgerstl/tEFMA CONTACT: : christian.jungreuthmayer@boku.ac.at or juergen.zanghellini@boku.ac.at SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Design of optimally constructed metabolic networks of minimal functionality.

BACKGROUND: Metabolic engineering aims to design microorganisms that will generate a product of interest at high yield. Thus, a variety of in silico modeling strategies has been applied successfully, including the concepts of elementary flux modes (EFMs) and constrained minimal cut sets (cMCSs). The EFMs (minimal, steady state pathways through the system) can be calculated given a metabolic model. cMCSs are sets of reaction deletions in such a network that will allow desired pathways to survive and disable undesired ones (e.g., those with low product secretion or low growth rates). Grouping the modes into desired and undesired categories had to be done manually until now. RESULTS: Although the optimal solution for a given set of pathways will always be found with the currently available tools, manual selection may lead to a sub-optimal solution with respect to a metabolic engineering target. A small change in the selection of modes can reduce the number of necessary deletions while only slightly reducing production. Based on our recently introduced formulation of cut set calculations using binary linear programming, we suggest an algorithm that does not require manual selection of the desired pathways. CONCLUSIONS: We demonstrated the principle of our algorithm with the help of a small toy network and applied it to a model of E. coli using different design objectives. Furthermore we validated our method by reproducing previously obtained results without requiring manual grouping of modes.

Comparison and improvement of algorithms for computing minimal cut sets.

BACKGROUND: Constrained minimal cut sets (cMCSs) have recently been introduced as a framework to enumerate minimal genetic intervention strategies for targeted optimization of metabolic networks. Two different algorithmic schemes (adapted Berge algorithm and binary integer programming) have been proposed to compute cMCSs from elementary modes. However, in their original formulation both algorithms are not fully comparable. RESULTS: Here we show that by a small extension to the integer program both methods become equivalent. Furthermore, based on well-known preprocessing procedures for integer programming we present efficient preprocessing steps which can be used for both algorithms. We then benchmark the numerical performance of the algorithms in several realistic medium-scale metabolic models. The benchmark calculations reveal (i) that these preprocessing steps can lead to an enormous speed-up under both algorithms, and (ii) that the adapted Berge algorithm outperforms the binary integer approach. CONCLUSIONS: Generally, both of our new implementations are by at least one order of magnitude faster than other currently available implementations.

Fast computation of minimal cut sets in metabolic networks with a Berge algorithm that utilizes binary bit pattern trees.

Minimal cut sets are a valuable tool for analyzing metabolic networks and for identifying optimal gene intervention strategies by eliminating unwanted metabolic functions and keeping desired functionality. Minimal cut sets rely on the concept of elementary flux modes, which are sets of indivisible metabolic pathways under steady-state condition. However, the computation of minimal cut sets is nontrivial, as even medium-sized metabolic networks with just 100 reactions easily have several hundred million elementary flux modes. We developed a minimal cut set tool that implements the well-known Berge algorithm and utilizes a novel approach to significantly reduce the program run time by using binary bit pattern trees. By using the introduced tree approach, the size of metabolic models that can be analyzed and optimized by minimal cut sets is pushed to new and considerably higher limits.

Elementary flux modes in a nutshell: properties, calculation and applications.

Elementary flux mode (EFM) analysis allows the unbiased decomposition of a metabolic network into minimal functional units, making it a powerful tool for metabolic engineering. While the use of EFM analysis (EFMA) is still limited by the size of the models it can handle, EFMA has been successfully applied to solve real-world metabolic engineering problems. Here we provide a user-oriented introduction to EFMA, provide examples of recent applications, analyze current research strategies to overcome the computational restrictions and give an overview over current approaches, which aim to identify and calculate only biologically relevant EFMs.

regEfmtool: speeding up elementary flux mode calculation using transcriptional regulatory rules in the form of three-state logic.

Despite the considerable progress made in recent years, the computation of the complete set of elementary flux modes of large or even genome-scale metabolic networks is still impossible. We present regEfmtool which is an extension to efmtool that utilizes transcriptional regulatory networks for the computation of elementary flux modes. The implemented extension significantly decreases the computational costs for the calculation of elementary flux modes, such as runtime, memory usage and disk space by omitting biologically infeasible solutions. Hence, using the presented regEfmtool pushes the size of metabolic networks that can be studied by elementary flux modes to new limits.