Evaluation of stochastic effects on biomolecular networks using the generalised Nyquist stability criterion

Abstract—Stochastic differential equations are now commonly used to model biomolecular networks in systems biology, and much recent research has been devoted to the development of methods to analyse their stability properties. Stability analysis of such systems may be performed using the Laplace transform, which requires the calculation of the exponential matrix involving time symbolically. However, the calculation of the symbolic exponential matrix is not feasible for problems of even moderate size, as the required computation time increases exponentially with the matrix order. To address this issue, we present a novel method for approximating the Laplace transform which does not require the exponential matrix to be calculated explicitly. The calculation time associated with the proposed method does not increase exponentially with the size of the system, and the approximation error is shown to be of the same order as existing methods. Using this approximation method, we show how a straightforward application of the generalized Nyquist stability criterion provides necessary and sufficient conditions for the stability of stochastic biomolecular networks. The usefulness and computational efficiency of the proposed method is illustrated through its application to the problem of analysing a model for limit-cycle oscillations in cAMP during aggregation of Dictyostelium cells

Domain-mediated interactions for protein subfamily identification

Within a protein family, proteins with the same domain often exhibit different cellular functions, despite the shared evolutionary history and molecular function of the domain. We hypothesized that domain-mediated interactions (DMIs) may categorize a protein family into subfamilies because the diversified functions of a single domain often depend on interacting partners of domains. Here we systematically identified DMI subfamilies, in which proteins share domains with DMI partners, as well as with various functional and physical interaction networks in individual species. In humans, DMI subfamily members are associated with similar diseases, including cancers, and are frequently co-associated with the same diseases. DMI information relates to the functional and evolutionary subdivisions of human kinases. In yeast, DMI subfamilies contain proteins with similar phenotypic outcomes from specific chemical treatments. Therefore, the systematic investigation here provides insights into the diverse functions of subfamilies derived from a protein family with a link-centric approach and suggests a useful resource for annotating the functions and phenotypic outcomes of proteins.11Ysciescopu

Robustness analysis of magnetic torquer controlled spacecraft attitude dynamics

This paper describes a systematic approach to the robustness analysis of linear periodically time-varying (LPTV) systems. The method uses the technique known as Lifting to transform the original time-varying uncertain system into linear fractional transformation (LFT) form. The stability and performance robustness of the system to structured parametric uncertainty can then be analysed non-conservatively using the structured singular value μ. The method is applied to analyse the stability robustness of an attitude control law for a spacecraft controlled by magnetic torquer bars, whose linearised dynamics can naturally be written in linear periodically time-varying form. The proposed method allows maximum allowable levels of uncertainty, as well as worst-case uncertainty combinations to be computed. The destabilising effect of these uncertain parameter combinations is verified in time-domain simulations

Validation of a model of regulation in the tryptophan operon against multiple experiment data using global optimisation

This paper is concerned with validating a mathematical model of regulation in the tryptophan operon using global optimization. Although a number of models for this biochemical network are proposed, in many cases only qualitative agreement between the model output and experimental data was demonstrated, since very little information is currently available to guide the selection of parameter values for the models. This paper presents a model validating method using both multiple experimental data and global optimization

Renormalization analysis of intermittency in two coupled maps

The critical behavior for intermittency is studied in two coupled one-dimensional (1D) maps. We find two fixed maps of an approximate renormalization operator in the space of coupled maps. Each fixed map has a common relavant eigenvaule associated with the scaling of the control parameter of the uncoupled one-dimensional map. However, the relevant ``coupling eigenvalue'' associated with coupling perturbation varies depending on the fixed maps. These renormalization results are also confirmed for a linearly-coupled case.Comment: 11 pages, RevTeX, 2 eps figure