Unscented Kalman filter with parameter identifiability analysis for the estimation of multiple parameters in kinetic models

In systems biology, experimentally measured parameters are not always available, necessitating the use of computationally based parameter estimation. In order to rely on estimated parameters, it is critical to first determine which parameters can be estimated for a given model and measurement set. This is done with parameter identifiability analysis. A kinetic model of the sucrose accumulation in the sugar cane culm tissue developed by Rohwer et al. was taken as a test case model. What differentiates this approach is the integration of an orthogonal-based local identifiability method into the unscented Kalman filter (UKF), rather than using the more common observability-based method which has inherent limitations. It also introduces a variable step size based on the system uncertainty of the UKF during the sensitivity calculation. This method identified 10 out of 12 parameters as identifiable. These ten parameters were estimated using the UKF, which was run 97 times. Throughout the repetitions the UKF proved to be more consistent than the estimation algorithms used for comparison

13C-assisted metabolic flux analysis to investigate heterotrophic and mixotrophic metabolism in Cupriavidus necator H16

Introduction. Cupriavidus necator H16 is a gram-negative bacterium, capable of lithoautotrophic growth by utilizing hydrogen as an energy source and fixing carbon dioxide (CO2) through Calvin-Benson-Bassham (CBB) cycle. The potential to utilize synthesis gas (Syngas) and the prospects of rerouting carbon from polyhydroxybutyrate synthesis to value-added compounds makes C. necator an excellent chassis for industrial application. Objectives. In the context of lack of sufficient quantitative information of the metabolic pathways and to advance in rational metabolic engineering for optimized product synthesis in C. necator H16, we carried out a metabolic flux analysis based on steady-state 13C-labelling. Methods. In this study, steady-state carbon labelling experiments, using either D-[1-13C]fructose or [1,2-13C]glycerol, were undertaken to investigate the carbon flux through the central carbon metabolism in C. necator H16 under heterotrophic and mixotrophic growth conditions, respectively. Results. We found that the CBB cycle is active even under heterotrophic condition, and growth is indeed mixotrophic. While Entner-Doudoroff (ED) pathway is shown to be the major route for sugar degradation, tricarboxylic acid (TCA) cycle is highly active in mixotrophic condition. Enhanced flux is observed in reductive pentose phosphate pathway (redPPP) under the mixotrophic condition to supplement the precursor requirement for CBB cycle. The flux distribution was compared to the mRNA abundance of genes encoding enzymes involved in key enzymatic reactions of the central carbon metabolism. Conclusion. This study leads the way to establishing 13C-based quantitative fluxomics for rational pathway engineering in C. necator H16

<it>iMS2Flux</it> – a high–throughput processing tool for stable isotope labeled mass spectrometric data used for metabolic flux analysis

<p>Abstract</p> <p>Background</p> <p>Metabolic flux analysis has become an established method in systems biology and functional genomics. The most common approach for determining intracellular metabolic fluxes is to utilize mass spectrometry in combination with stable isotope labeling experiments. However, before the mass spectrometric data can be used it has to be corrected for biases caused by naturally occurring stable isotopes, by the analytical technique(s) employed, or by the biological sample itself. Finally the MS data and the labeling information it contains have to be assembled into a data format usable by flux analysis software (of which several dedicated packages exist). Currently the processing of mass spectrometric data is time-consuming and error-prone requiring peak by peak cut-and-paste analysis and manual curation. In order to facilitate high-throughput metabolic flux analysis, the automation of multiple steps in the analytical workflow is necessary.</p> <p>Results</p> <p>Here we describe <it>iMS2Flux,</it> software developed to automate, standardize and connect the data flow between mass spectrometric measurements and flux analysis programs. This tool streamlines the transfer of data from extraction via correction tools to ¹³C-Flux software by processing MS data from stable isotope labeling experiments. It allows the correction of large and heterogeneous MS datasets for the presence of naturally occurring stable isotopes, initial biomass and several mass spectrometry effects. Before and after data correction, several checks can be performed to ensure accurate data. The corrected data may be returned in a variety of formats including those used by metabolic flux analysis software such as <it>13CFLUX</it>, <it>OpenFLUX</it> and <it>13CFLUX2</it>.</p> <p>Conclusion</p> <p><it>iMS2Flux</it> is a versatile, easy to use tool for the automated processing of mass spectrometric data containing isotope labeling information. It represents the core framework for a standardized workflow and data processing. Due to its flexibility it facilitates the inclusion of different experimental datasets and thus can contribute to the expansion of flux analysis applications.</p

Separating the wheat from the chaff - a strategy to utilize plant genetic resources from <em>ex situ</em> genebanks.

The need for higher yielding and better-adapted crop plants for feeding the world&#39;s rapidly growing population has raised the question of how to systematically utilize large genebank collections with their wide range of largely untouched genetic diversity. Phenotypic data that has been recorded for decades during various rounds of seed multiplication provides a rich source of information. Their usefulness has remained limited though, due to various biases induced by conservation management over time or changing environmental conditions. Here, we present a powerful procedure that permits an unbiased trait-based selection of plant samples based on such phenotypic data. Applying this technique to the wheat collection of one of the largest genebanks worldwide, we identified groups of plant samples displaying contrasting phenotypes for selected traits. As a proof of concept for our discovery pipeline, we resequenced the entire major but conserved flowering time locus Ppd-D1 in just a few such selected wheat samples - and nearly doubled the number of hitherto known alleles

An improved constraint filtering technique for inferring hidden states and parameters of a biological model

Motivation: In systems biology, kinetic models represent the biological system using a set of ordinary differential equations (ODEs). The correct values of the parameters within these ODEs are critical for a reliable study of the dynamic behaviour of such systems. Typically, it is only possible to experimentally measure a fraction of these parameter values. The rest must be indirectly determined from measurements of other quantities. In this article, we propose a novel statistical inference technique to computationally estimate these unknown parameter values. By characterizing the ODEs with non-linear state-space equations, this inference technique models the unknown parameters as hidden states, which can then be estimated from noisy measurement data.Results: Here we extended the square-root unscented Kalman filter SR-UKF proposed by Merwe and Wan to include constraints with the state estimation process. We developed the constrained square-root unscented Kalman filter (CSUKF) to estimate parameters of non-linear state-space models. This probabilistic inference technique was successfully used to estimate parameters of a glycolysis model in yeast and a gene regulatory network. We showed that our method is numerically stable and can reliably estimate parameters within a biologically meaningful parameter space from noisy observations. When compared with the two common non-linear extensions of Kalman filter in addition to four widely used global optimization algorithms, CSUKF is shown to be both accurate and computationally efficient. With CSUKF, statistical analysis is straightforward, as it directly provides the uncertainty on the estimation result