Reconstruction of arbitrary biochemical reaction networks: A compressive sensing approach

Reconstruction of biochemical reaction networks (BRN) and genetic regulatory networks (GRN) in particular is a central topic in systems biology which raises crucial theoretical challenges in system identification. Nonlinear Ordinary Differential Equations (ODEs) that involve polynomial and rational functions are typically used to model biochemical reaction networks. Such nonlinear models make the problem of determining the connectivity of biochemical networks from time-series experimental data quite difficult. In this paper, we present a network reconstruction algorithm that can deal with ODE model descriptions containing polynomial and rational functions. Rather than identifying the parameters of linear or nonlinear ODEs characterised by pre-defined equation structures, our methodology allows us to determine the nonlinear ODEs structure together with their associated parameters. To solve the network reconstruction problem, we cast it as a compressive sensing (CS) problem and use sparse Bayesian learning (SBL) algorithms as a computationally efficient and robust way to obtain its solution. © 2012 IEEE.

Dynamical structure function identifiability conditions enabling signal structure reconstruction

Networks of controlled dynamical systems exhibit a variety of interconnection patterns that could be interpreted as the structure of the system. One such interpretation of system structure is a system's signal structure, characterized as the open-loop causal dependencies among manifest variables and represented by its dynamical structure function. Although this notion of structure is among the weakest available, previous work has shown that if no a priori structural information is known about the system, not even the Boolean structure of the dynamical structure function is identifiable. Consequently, one method previously suggested for obtaining the necessary a priori structural information is to leverage knowledge about target specificity of the controlled inputs. This work extends these results to demonstrate precisely the a priori structural information that is both necessary and sufficient to reconstruct the network from input-output data. This extension is important because it significantly broadens the applicability of the identifiability conditions, enabling the design of network reconstruction experiments that were previously impossible due to practical constraints on the types of actuation mechanisms available to the engineer or scientist. The work is motivated by the proteomics problem of reconstructing the Per-Arnt-Sim Kinase pathway used in the metabolism of sugars. © 2012 IEEE.

Quantitative convergence analysis of multi-agent systems

We introduce a characterization of contraction for bounded convex sets. For discrete-time multi-agent systems we provide an explicit upperbound on the rate of convergence to a consensus under the assumptions of contractiveness and (weak) connectedness (across an interval.) Convergence is shown to be exponential when either the system or the function characterizing the contraction is linear. Copyright © 2007 IFAC.

Implementation of shading effect for reconstruction of smooth layerbased 3D holographic images

A holographic rendering algorithm using a layer-based structure with angular tiling supports view-dependent shading and accommodation cues. This approach also has the advantages of rapid computation speed and visual reduction of layer gap artefacts compared to other approaches. Holograms rendered with this algorithm are displayed using an SLM to demonstrate view-dependent shading and occlusion. © 2013 SPIE-IS&T.

Expressive visual text-to-speech using active appearance models

This paper presents a complete system for expressive visual text-to-speech (VTTS), which is capable of producing expressive output, in the form of a 'talking head', given an input text and a set of continuous expression weights. The face is modeled using an active appearance model (AAM), and several extensions are proposed which make it more applicable to the task of VTTS. The model allows for normalization with respect to both pose and blink state which significantly reduces artifacts in the resulting synthesized sequences. We demonstrate quantitative improvements in terms of reconstruction error over a million frames, as well as in large-scale user studies, comparing the output of different systems. © 2013 IEEE.

Free energy surface reconstruction from umbrella samples using Gaussian process regression

We demonstrate how the Gaussian process regression approach can be used to efficiently reconstruct free energy surfaces from umbrella sampling simulations. By making a prior assumption of smoothness and taking account of the sampling noise in a consistent fashion, we achieve a significant improvement in accuracy over the state of the art in two or more dimensions or, equivalently, a significant cost reduction to obtain the free energy surface within a prescribed tolerance in both regimes of spatially sparse data and short sampling trajectories. Stemming from its Bayesian interpretation the method provides meaningful error bars without significant additional computation. A software implementation is made available on www.libatoms.org.

Free energy surface reconstruction from umbrella samples using Gaussian process regression II: Multiple collective variables

We demonstrate how a prior assumption of smoothness can be used to enhance the reconstruction of free energy profiles from multiple umbrella sampling simulations using the Bayesian Gaussian process regression approach. The method we derive allows the concurrent use of histograms and free energy gradients and can easily be extended to include further data. In Part I we review the necessary theory and test the method for one collective variable. We demonstrate improved performance with respect to the weighted histogram analysis method and obtain meaningful error bars without any significant additional computation. In Part II we consider the case of multiple collective variables and compare to a reconstruction using least squares fitting of radial basis functions. We find substantial improvements in the regimes of spatially sparse data or short sampling trajectories. A software implementation is made available on www.libatoms.org.

A high throughput Nile red method for quantitative measurement of neutral lipids in microalgae

Isolation of high neutral lipid-containing microalgae is key to the commercial success of microalgae-based biofuel production. The Nile red fluorescence method has been successfully applied to the determination of lipids in certain microalgae, but has been unsuccessful in many others, particularly those with thick, rigid cell walls that prevent the penetration of the fluorescence dye. The conventional "one sample at a time" method was also time-consuming. In this study, the solvent dimethyl sulfoxide (DMSO) was introduced to microalgal samples as the stain carrier at an elevated temperature. The cellular neutral lipids were determined and quantified using a 96-well plate on a fluorescence spectrophotometer with an excitation wavelength of 530 nm and an emission wavelength of 575 run. An optimized procedure yielded a high correlation coefficient (R-2 = 0.998) with the lipid standard triolein and repeated measurements of replicates. Application of the improved method to several green algal strains gave very reproducible results with relative standard errors of 8.5%, 3.9% and 8.6%, 4.5% for repeatability and reproducibility at two concentration levels (2.0 mu g/mL and 20 mu g/mL), respectively. Moreover, the detection and quantification limits of the improved Nile red staining method were 0.8 mu g/mL and 2.0 mu g/mL for the neutral lipid standard triolein, respectively. The modified method and a conventional gravimetric determination method provided similar results on replicate samples. The 96-well plate-based Nile red method can be used as a high throughput technique for rapid screening of a broader spectrum of naturally-occurring and genetically-modified algal strains and mutants for high neutral lipid/oil production. (C) 2009 Published by Elsevier B.V.

Network reconstruction using knock-out and over-expression data

This paper outlines necessary and sufficient conditions for network reconstruction of linear, time-invariant systems using data from either knock-out or over-expression experiments. These structural system perturbations, which are common in biological experiments, can be formulated as unknown system inputs, allowing the network topology and dynamics to be found. We assume that only partial state measurements are available and propose an algorithm that can reconstruct the network at the level of the measured states using either time-series or steady-state data. A simulated example illustrates how the algorithm successfully reconstructs a network from data. © 2013 EUCA.