The control of complex networks is one of the most challenging problems in the fields of biology and engineering. In this study, the authors explored the controllability and control energy of several signalling networks, which consisted of many interconnected pathways, including networks with a bow-tie architecture. On the basis of the theory of structure controllability, they revealed that biological mechanisms, such as cross-pathway interactions, compartmentalisation and so on make the networks easier to fully control. Furthermore, using numerical simulations for two realistic examples, they demonstrated that the control energy of normal networks with crosstalk is lower than in networks without crosstalk. These results indicate that the biological networks are optimally designed to achieve their normal functions from the viewpoint of the control theory. The authors’ work provides a comprehensive understanding of the impact of network structures and properties on controllability.

Parameter estimation for high dimension complex dynamic system is a hot topic. However, the current statistical model and inference approach is known as a large p small n problem. How to reduce the dimension of the dynamic model and improve the accuracy of estimation is more important. To address this question, the authors take some known parameters and structure of system as priori knowledge and incorporate it into dynamic model. At the same time, they decompose the whole dynamic model into subset network modules, based on different modules, and then they apply different estimation approaches. This technique is called Rao-Blackwellised particle filters decomposition methods. To evaluate the performance of this method, the authors apply it to synthetic data generated from repressilator model and experimental data of the JAK-STAT pathway, but this method can be easily extended to large-scale cases.

Xiamenmycins, a series of prenylated benzopyran compounds with anti-fibrotic bioactivities, were isolated from a mangrove-derived Streptomyces xiamenensis. To fulfil the requirements of pharmaceutical investigations, a high production of xiamenmycin is needed. In this study,, the authors present a kinetic metabolic model to evaluate fluxes in an engineered Streptomyces lividans with xiamenmycin-oriented genetic modification based on generic enzymatic rate equations and stability constraints. Lyapunov function was used for a viability optimisation. From their kinetic model, the flux distributions for the engineered S. lividans fed on glucose and glycerol as carbon sources were calculated. They found that if the bacterium can utilise glucose simultaneously with glycerol, xiamenmycin production can be enhanced by 40% theoretically, while maintaining the same growth rate. Glycerol may increase the flux for phosphoenolpyruvate synthesis without interfering citric acid cycle. They therefore believe this study demonstrates a possible new direction for bioengineering of S. lividans.

A structure-based statistical potential is developed for transcription factor binding site (TFBS) prediction. Besides the direct contact between amino acids from TFs and DNA bases, the authors also considered the influence of the neighbouring base. This three-body potential showed better discriminate powers than the two-body potential. They validate the performance of the potential in TFBS identification, binding energy prediction and binding mutation prediction.

In this study, the authors studied the protein structure prediction problem by the two-dimensional hydrophobic–polar model on triangular lattice. Particularly the non-compact conformation was modelled to fold the amino acid sequence into a relatively larger triangular lattice, which is more biologically realistic and significant than the compact conformation. Then protein structure prediction problem was abstracted to match amino acids to lattice points. Mathematically, the problem was formulated as an integer programming and they transformed the biological problem into an optimisation problem. To solve this problem, classical particle swarm optimisation algorithm was extended by the single point adjustment strategy. Compared with square lattice, conformations on triangular lattice are more flexible in several benchmark examples. They further compared the authors’ algorithm with hybrid of hill climbing and genetic algorithm. The results showed that their method was more effective in finding solution with lower energy and less running time.

The study of biology and medicine in a noise environment is an evolving direction in biological data analysis. Among these studies, analysis of electrocardiogram (ECG) signals in a noise environment is a challenging direction in personalized medicine. Due to its periodic characteristic, ECG signal can be roughly regarded as sparse biomedical signals. This study proposes a two-stage recovery algorithm for sparse biomedical signals in time domain. In the first stage, the concentration subspaces are found in advance. Then by exploiting these subspaces, the mixing matrix is estimated accurately. In the second stage, based on the number of active sources at each time point, the time points are divided into different layers. Next, by constructing some transformation matrices, these time points form a row echelon-like system. After that, the sources at each layer can be solved out explicitly by corresponding matrix operations. It is noting that all these operations are conducted under a weak sparse condition that the number of active sources is less than the number of observations. Experimental results show that the proposed method has a better performance for sparse ECG signal recovery problem.

The big omics data are challenging translational bioinformatics in an unprecedented way for its complexities and volumes. How to employ big omics data to achieve a rivalling-clinical, reproducible disease diagnosis from a systems approach is an urgent problem to be solved in translational bioinformatics and machine learning. In this study, the authors propose a novel transcriptome marker diagnosis to tackle this problem using big RNA-seq data by viewing whole transcriptome as a profile marker systematically. The systems diagnosis not only avoids the reproducibility issue of the existing gene-/network-marker-based diagnostic methods, but also achieves rivalling-clinical diagnostic results by extracting true signals from big RNA-seq data. Their method demonstrates a better fit for personalised diagnostics by attaining exceptional diagnostic performance via using systems information than its competitive methods and prepares itself as a good candidate for clinical usage. To the best of their knowledge, it is the first study on this topic and will inspire the more investigations in big omics data diagnostics.