PyCoMo: a python package for community metabolic model creation and analysis.

SUMMARY: PyCoMo is a python package for quick and easy generation of genome-scale compartmentalized community metabolic models that are compliant with current openCOBRA file formats. The resulting models can be used to predict (i) the maximum growth rate at a given abundance profile, (ii) the feasible community compositions at a given growth rate, and (iii) all exchange metabolites and cross-feeding interactions in a community metabolic model independent of the abundance profile; we demonstrate PyCoMo's capability by analysing methane production in a previously published simplified biogas community metabolic model. AVAILABILITY AND IMPLEMENTATION: PyCoMo is freely available under an MIT licence at http://github.com/univieCUBE/PyCoMo, the Python Package Index, and Zenodo.

Logistic PCA explains differences between genome-scale metabolic models in terms of metabolic pathways.

Genome-scale metabolic models (GSMMs) offer a holistic view of biochemical reaction networks, enabling in-depth analyses of metabolism across species and tissues in multiple conditions. However, comparing GSMMs Against each other poses challenges as current dimensionality reduction algorithms or clustering methods lack mechanistic interpretability, and often rely on subjective assumptions. Here, we propose a new approach utilizing logisitic principal component analysis (LPCA) that efficiently clusters GSMMs while singling out mechanistic differences in terms of reactions and pathways that drive the categorization. We applied LPCA to multiple diverse datasets, including GSMMs of 222 Escherichia-strains, 343 budding yeasts (Saccharomycotina), 80 human tissues, and 2943 Firmicutes strains. Our findings demonstrate LPCA's effectiveness in preserving microbial phylogenetic relationships and discerning human tissue-specific metabolic profiles, exhibiting comparable performance to traditional methods like t-distributed stochastic neighborhood embedding (t-SNE) and Jaccard coefficients. Moreover, the subsystems and associated reactions identified by LPCA align with existing knowledge, underscoring its reliability in dissecting GSMMs and uncovering the underlying drivers of separation.

Tailored Mass Spectral Data Exploration Using the SpecXplore Interactive Dashboard.

Untargeted metabolomics promises comprehensive characterization of small molecules in biological samples. However, the field is hampered by low annotation rates and abstract spectral data. Despite recent advances in computational metabolomics, manual annotations and manual confirmation of in-silico annotations remain important in the field. Here, exploratory data analysis methods for mass spectral data provide overviews, prioritization, and structural hypothesis starting points to researchers facing large quantities of spectral data. In this research, we propose a fluid means of dealing with mass spectral data using specXplore, an interactive Python dashboard providing interactive and complementary visualizations facilitating mass spectral similarity matrix exploration. Specifically, specXplore provides a two-dimensional t-distributed stochastic neighbor embedding embedding as a jumping board for local connectivity exploration using complementary interactive visualizations in the form of partial network drawings, similarity heatmaps, and fragmentation overview maps. SpecXplore makes use of state-of-the-art ms2deepscore pairwise spectral similarities as a quantitative backbone while allowing fast changes of threshold and connectivity limitation settings, providing flexibility in adjusting settings to suit the localized node environment being explored. We believe that specXplore can become an integral part of mass spectral data exploration efforts and assist users in the generation of structural hypotheses for compounds of interest.

ecmtool: fast and memory-efficient enumeration of elementary conversion modes.

MOTIVATION: Characterizing all steady-state flux distributions in metabolic models remains limited to small models due to the explosion of possibilities. Often it is sufficient to look only at all possible overall conversions a cell can catalyze ignoring the details of intracellular metabolism. Such a characterization is achieved by elementary conversion modes (ECMs), which can be conveniently computed with ecmtool. However, currently, ecmtool is memory intensive, and it cannot be aided appreciably by parallelization. RESULTS: We integrate mplrs-a scalable parallel vertex enumeration method-into ecmtool. This speeds up computation, drastically reduces memory requirements and enables ecmtool's use in standard and high-performance computing environments. We show the new capabilities by enumerating all feasible ECMs of the near-complete metabolic model of the minimal cell JCVI-syn3.0. Despite the cell's minimal character, the model gives rise to 4.2×109 ECMs and still contains several redundant sub-networks. AVAILABILITY AND IMPLEMENTATION: ecmtool is available at https://github.com/SystemsBioinformatics/ecmtool. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Homologue series detection and management in LC-MS data with homologueDiscoverer.

SUMMARY: Untargeted metabolomics data analysis is highly labour intensive and can be severely frustrated by both experimental noise and redundant features. Homologous polymer series is a particular case of features that can either represent large numbers of noise features or alternatively represent features of interest with large peak redundancy. Here, we present homologueDiscoverer, an R package that allows for the targeted and untargeted detection of homologue series as well as their evaluation and management using interactive plots and simple local database functionalities. AVAILABILITY AND IMPLEMENTATION: homologueDiscoverer is freely available at GitHub https://github.com/kevinmildau/homologueDiscoverer. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Elementary vectors and autocatalytic sets for resource allocation in next-generation models of cellular growth.

Traditional (genome-scale) metabolic models of cellular growth involve an approximate biomass "reaction", which specifies biomass composition in terms of precursor metabolites (such as amino acids and nucleotides). On the one hand, biomass composition is often not known exactly and may vary drastically between conditions and strains. On the other hand, the predictions of computational models crucially depend on biomass. Also elementary flux modes (EFMs), which generate the flux cone, depend on the biomass reaction. To better understand cellular phenotypes across growth conditions, we introduce and analyze new classes of elementary vectors for comprehensive (next-generation) metabolic models, involving explicit synthesis reactions for all macromolecules. Elementary growth modes (EGMs) are given by stoichiometry and generate the growth cone. Unlike EFMs, they are not support-minimal, in general, but cannot be decomposed "without cancellations". In models with additional (capacity) constraints, elementary growth vectors (EGVs) generate a growth polyhedron and depend also on growth rate. However, EGMs/EGVs do not depend on the biomass composition. In fact, they cover all possible biomass compositions and can be seen as unbiased versions of elementary flux modes/vectors (EFMs/EFVs) used in traditional models. To relate the new concepts to other branches of theory, we consider autocatalytic sets of reactions. Further, we illustrate our results in a small model of a self-fabricating cell, involving glucose and ammonium uptake, amino acid and lipid synthesis, and the expression of all enzymes and the ribosome itself. In particular, we study the variation of biomass composition as a function of growth rate. In agreement with experimental data, low nitrogen uptake correlates with high carbon (lipid) storage.

Sulfate limitation increases specific plasmid DNA yield and productivity in E. coli fed-batch processes.

Plasmid DNA (pDNA) is a key biotechnological product whose importance became apparent in the last years due to its role as a raw material in the messenger ribonucleic acid (mRNA) vaccine manufacturing process. In pharmaceutical production processes, cells need to grow in the defined medium in order to guarantee the highest standards of quality and repeatability. However, often these requirements result in low product titer, productivity, and yield. In this study, we used constraint-based metabolic modeling to optimize the average volumetric productivity of pDNA production in a fed-batch process. We identified a set of 13 nutrients in the growth medium that are essential for cell growth but not for pDNA replication. When these nutrients are depleted in the medium, cell growth is stalled and pDNA production is increased, raising the specific and volumetric yield and productivity. To exploit this effect we designed a three-stage process (1. batch, 2. fed-batch with cell growth, 3. fed-batch without cell growth). The transition between stage 2 and 3 is induced by sulfate starvation. Its onset can be easily controlled via the initial concentration of sulfate in the medium. We validated the decoupling behavior of sulfate and assessed pDNA quality attributes (supercoiled pDNA content) in E. coli with lab-scale bioreactor cultivations. The results showed an increase in supercoiled pDNA to biomass yield by 33% and an increase of supercoiled pDNA volumetric productivity by 13 % upon limitation of sulfate. In conclusion, even for routinely manufactured biotechnological products such as pDNA, simple changes in the growth medium can significantly improve the yield and quality.

Probabilistic quotient's work and pharmacokinetics' contribution: countering size effect in metabolic time series measurements.

Metabolomic time course analyses of biofluids are highly relevant for clinical diagnostics. However, many sampling methods suffer from unknown sample sizes, commonly known as size effects. This prevents absolute quantification of biomarkers. Recently, several mathematical post acquisition normalization methods have been developed to overcome these problems either by exploiting already known pharmacokinetic information or by statistical means. Here we present an improved normalization method, MIX, that combines the advantages of both approaches. It couples two normalization terms, one based on a pharmacokinetic model (PKM) and the other representing a popular statistical approach, probabilistic quotient normalization (PQN), in a single model. To test the performance of MIX, we generated synthetic data closely resembling real finger sweat metabolome measurements. We show that MIX normalization successfully tackles key weaknesses of the individual strategies: it (i) reduces the risk of overfitting with PKM, and (ii), contrary to PQN, it allows to compute sample volumes. Finally, we validate MIX by using real finger sweat as well as blood plasma metabolome data and demonstrate that MIX allows to better and more robustly correct for size effects. In conclusion, the MIX method improves the reliability and robustness of quantitative biomarker detection in finger sweat and other biofluids, paving the way for biomarker discovery and hypothesis generation from metabolomic time course data.

Inclusion of maintenance energy improves the intracellular flux predictions of CHO.

Chinese hamster ovary (CHO) cells are the leading platform for the production of biopharmaceuticals with human-like glycosylation. The standard practice for cell line generation relies on trial and error approaches such as adaptive evolution and high-throughput screening, which typically take several months. Metabolic modeling could aid in designing better producer cell lines and thus shorten development times. The genome-scale metabolic model (GSMM) of CHO can accurately predict growth rates. However, in order to predict rational engineering strategies it also needs to accurately predict intracellular fluxes. In this work we evaluated the agreement between the fluxes predicted by parsimonious flux balance analysis (pFBA) using the CHO GSMM and a wide range of 13C metabolic flux data from literature. While glycolytic fluxes were predicted relatively well, the fluxes of tricarboxylic acid (TCA) cycle were vastly underestimated due to too low energy demand. Inclusion of computationally estimated maintenance energy significantly improved the overall accuracy of intracellular flux predictions. Maintenance energy was therefore determined experimentally by running continuous cultures at different growth rates and evaluating their respective energy consumption. The experimentally and computationally determined maintenance energy were in good agreement. Additionally, we compared alternative objective functions (minimization of uptake rates of seven nonessential metabolites) to the biomass objective. While the predictions of the uptake rates were quite inaccurate for most objectives, the predictions of the intracellular fluxes were comparable to the biomass objective function.

EFMlrs: a Python package for elementary flux mode enumeration via lexicographic reverse search.

BACKGROUND: Elementary flux mode (EFM) analysis is a well-established, yet computationally challenging approach to characterize metabolic networks. Standard algorithms require huge amounts of memory and lack scalability which limits their application to single servers and consequently limits a comprehensive analysis to medium-scale networks. Recently, Avis et al. developed mplrs-a parallel version of the lexicographic reverse search (lrs) algorithm, which, in principle, enables an EFM analysis on high-performance computing environments (Avis and Jordan. mplrs: a scalable parallel vertex/facet enumeration code. arXiv:1511.06487 , 2017). Here we test its applicability for EFM enumeration. RESULTS: We developed EFMlrs, a Python package that gives users access to the enumeration capabilities of mplrs. EFMlrs uses COBRApy to process metabolic models from sbml files, performs loss-free compressions of the stoichiometric matrix, and generates suitable inputs for mplrs as well as efmtool, providing support not only for our proposed new method for EFM enumeration but also for already established tools. By leveraging COBRApy, EFMlrs also allows the application of additional reaction boundaries and seamlessly integrates into existing workflows. CONCLUSION: We show that due to mplrs's properties, the algorithm is perfectly suited for high-performance computing (HPC) and thus offers new possibilities for the unbiased analysis of substantially larger metabolic models via EFM analyses. EFMlrs is an open-source program that comes together with a designated workflow and can be easily installed via pip.

Flux tope analysis: studying the coordination of reaction directions in metabolic networks.

Motivation: Elementary flux mode (EFM) analysis allows an unbiased description of metabolic networks in terms of minimal pathways (involving a minimal set of reactions). To date, the enumeration of EFMs is impracticable in genome-scale metabolic models. In a complementary approach, we introduce the concept of a flux tope (FT), involving a maximal set of reactions (with fixed directions), which allows one to study the coordination of reaction directions in metabolic networks and opens a new way for EFM enumeration. Results: A FT is a (nontrivial) subset of the flux cone specified by fixing the directions of all reversible reactions. In a consistent metabolic network (without unused reactions), every FT contains a 'maximal pathway', carrying flux in all reactions. This decomposition of the flux cone into FTs allows the enumeration of EFMs (of individual FTs) without increasing the problem dimension by reaction splitting. To develop a mathematical framework for FT analysis, we build on the concepts of sign vectors and hyperplane arrangements. Thereby, we observe that FT analysis can be applied also to flux optimization problems involving additional (inhomogeneous) linear constraints. For the enumeration of FTs, we adapt the reverse search algorithm and provide an efficient implementation. We demonstrate that (biomass-optimal) FTs can be enumerated in genome-scale metabolic models of B.cuenoti and E.coli, and we use FTs to enumerate EFMs in models of M.genitalium and B.cuenoti. Availability and implementation: The source code is freely available at https://github.com/mpgerstl/FTA. Supplementary information: Supplementary data are available at Bioinformatics online.

What CHO is made of: Variations in the biomass composition of Chinese hamster ovary cell lines.

BACKGROUND: Cell line-specific, genome-scale metabolic models enable rigorous and systematic in silico investigation of cellular metabolism. Such models have recently become available for Chinese hamster ovary (CHO) cells. However, a key ingredient, namely an experimentally validated biomass function that summarizes the cellular composition, was so far missing. Here, we close this gap by providing extensive experimental data on the biomass composition of 13 parental and producer CHO cell lines under various conditions. RESULTS: We report total protein, lipid, DNA, RNA and carbohydrate content, cell dry mass, and detailed protein and lipid composition. Furthermore, we present meticulous data on exchange rates between cells and environment and provide detailed experimental protocols on how to determine all of the above. The biomass composition is converted into cell line- and condition-specific biomass functions for use in cell line-specific, genome-scale metabolic models of CHO. Finally, flux balance analysis (FBA) is used to demonstrate consistency between in silico predictions and experimental analysis. CONCLUSIONS: Our study reveals a strong variability of the total protein content and cell dry mass across cell lines. However, the relative amino acid composition is independent of the cell line and condition and thus needs not be explicitly measured for each new cell line. In contrast, the lipid composition is strongly influenced by the growth media and thus will have to be determined in each case. These cell line-specific variations in biomass composition have a small impact on growth rate predictions with FBA, as inaccuracies in the predictions are rather dominated by inaccuracies in the exchange rate spectra. Cell-specific biomass variations only become important if the experimental errors in the exchange rate spectra drop below twenty percent.

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.

A mathematical framework for yield (vs. rate) optimization in constraint-based modeling and applications in metabolic engineering.

BACKGROUND: The optimization of metabolic rates (as linear objective functions) represents the methodical core of flux-balance analysis techniques which have become a standard tool for the study of genome-scale metabolic models. Besides (growth and synthesis) rates, metabolic yields are key parameters for the characterization of biochemical transformation processes, especially in the context of biotechnological applications. However, yields are ratios of rates, and hence the optimization of yields (as nonlinear objective functions) under arbitrary linear constraints is not possible with current flux-balance analysis techniques. Despite the fundamental importance of yields in constraint-based modeling, a comprehensive mathematical framework for yield optimization is still missing. RESULTS: We present a mathematical theory that allows one to systematically compute and analyze yield-optimal solutions of metabolic models under arbitrary linear constraints. In particular, we formulate yield optimization as a linear-fractional program. For practical computations, we transform the linear-fractional yield optimization problem to a (higher-dimensional) linear problem. Its solutions determine the solutions of the original problem and can be used to predict yield-optimal flux distributions in genome-scale metabolic models. For the theoretical analysis, we consider the linear-fractional problem directly. Most importantly, we show that the yield-optimal solution set (like the rate-optimal solution set) is determined by (yield-optimal) elementary flux vectors of the underlying metabolic model. However, yield- and rate-optimal solutions may differ from each other, and hence optimal (biomass or product) yields are not necessarily obtained at solutions with optimal (growth or synthesis) rates. Moreover, we discuss phase planes/production envelopes and yield spaces, in particular, we prove that yield spaces are convex and provide algorithms for their computation. We illustrate our findings by a small example and demonstrate their relevance for metabolic engineering with realistic models of E. coli. CONCLUSIONS: We develop a comprehensive mathematical framework for yield optimization in metabolic models. Our theory is particularly useful for the study and rational modification of cell factories designed under given yield and/or rate requirements.

Metabolic engineering of Pichia pastoris.

Besides its use for efficient production of recombinant proteins the methylotrophic yeast Pichia pastoris (syn. Komagataella spp.) has been increasingly employed as a platform to produce metabolites of varying origin. We summarize here the impressive methodological developments of the last years to model and analyze the metabolism of P. pastoris, and to engineer its genome and metabolic pathways. Efficient methods to insert, modify or delete genes via homologous recombination and CRISPR/Cas9, supported by modular cloning techniques, have been reported. An outstanding early example of metabolic engineering in P. pastoris was the humanization of protein glycosylation. More recently the cell metabolism was engineered also to enhance the productivity of heterologous proteins. The last few years have seen an increased number of metabolic pathway design and engineering in P. pastoris, mainly towards the production of complex (secondary) metabolites. In this review, we discuss the potential role of P. pastoris as a platform for metabolic engineering, its strengths, and major requirements for future developments of chassis strains based on synthetic biology principles.

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.

Comprehensive assessment of measurement uncertainty in 13C-based metabolic flux experiments.

In the field of metabolic engineering 13C-based metabolic flux analysis experiments have proven successful in indicating points of action. As every step of this approach is affected by an inherent error, the aim of the present work is the comprehensive evaluation of factors contributing to the uncertainty of nonnaturally distributed C-isotopologue abundances as well as to the absolute flux value calculation. For this purpose, a previously published data set, analyzed in the course of a 13C labeling experiment studying glycolysis and the pentose phosphate pathway in a yeast cell factory, was used. Here, for isotopologue pattern analysis of these highly polar metabolites that occur in multiple isomeric forms, a gas chromatographic separation approach with preceding derivatization was used. This rendered a natural isotope interference correction step essential. Uncertainty estimation of the resulting C-isotopologue distribution was performed according to the EURACHEM guidelines with Monte Carlo simulation. It revealed a significant increase for low-abundance isotopologue fractions after application of the necessary correction step. For absolute flux value estimation, isotopologue fractions of various sugar phosphates, together with the assessed uncertainties, were used in a metabolic model describing the upper part of the central carbon metabolism. The findings pinpointed the influence of small isotopologue fractions as sources of error and highlight the need for improved model curation. Graphical abstract ᅟ.

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).

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.