Search

Background: Recent advances in omics technologies have raised great opportunities to study large-scale regulatory networks inside the cell. In addition, single-cell experiments have measured the gene and protein activities in a large number of cells under the same experimental conditions. However, a significant challenge in computational biology and bioinformatics is how to derive quantitative information from the single-cell observations and how to develop sophisticated mathematical models to describe the dynamic properties of regulatory networks using the derived quantitative information., Methods: This work designs an integrated approach to reverse-engineer gene networks for regulating early blood development based on singel-cell experimental observations. The wanderlust algorithm is initially used to develop the pseudo-trajectory for the activities of a number of genes. Since the gene expression data in the developed pseudo-trajectory show large fluctuations, we then use Gaussian process regression methods to smooth the gene express data in order to obtain pseudo-trajectories with much less fluctuations. The proposed integrated framework consists of both bioinformatics algorithms to reconstruct the regulatory network and mathematical models using differential equations to describe the dynamics of gene expression., Results: The developed approach is applied to study the network regulating early blood cell development. A graphic model is constructed for a regulatory network with forty genes and a dynamic model using differential equations is developed for a network of nine genes. Numerical results suggests that the proposed model is able to match experimental data very well. We also examine the networks with more regulatory relations and numerical results show that more regulations may exist. We test the possibility of auto-regulation but numerical simulations do not support the positive auto-regulation. In addition, robustness is used as an importantly additional criterion to select candidate networks., Conclusion: The research results in this work shows that the developed approach is an efficient and effective method to reverse-engineer gene networks using single-cell experimental observations.

Background: Influenza A virus (IAV) infection-induced inflammatory regulatory networks (IRNs) are extremely complex and dynamic. Specific biological experiments for investigating the interactions between individual inflammatory factors cannot provide a detailed and insightful multidimensional view of IRNs. Recently, data from high-throughput technologies have permitted system-level analyses. The construction of large and cell-specific IRNs from high-throughput data is essential to understanding the pathogenesis of IAV infection., Results: In this study, we proposed a computational method, which combines nonlinear ordinary differential equation (ODE)-based optimization with mutual information, to construct a cell-specific optimized IRN during IAV infection by integrating gene expression data with a prior knowledge of network topology. Moreover, we used the average relative error and sensitivity analysis to evaluate the effectiveness of our proposed approach. Furthermore, from the optimized IRN, we confirmed 45 interactions between proteins in biological experiments and identified 37 new regulatory interactions and 8 false positive interactions, including the following interactions: IL1β regulates TLR3, TLR3 regulates IFN-β and TNF regulates IL6. Most of these regulatory interactions are statistically significant by Z-statistic. The functional annotations of the optimized IRN demonstrated clearly that the defense response, immune response, response to wounding and regulation of cytokine production are the pivotal processes of IAV-induced inflammatory response. The pathway analysis results from the Kyoto Encyclopaedia of Genes and Genomes (KEGG) showed that 8 pathways are enriched significantly. The 5 pathways were validated by experiments, and 3 other pathways, including the intestinal immune network for IgA production, the cytosolic DNA-sensing pathway and the allograft rejection pathway, are the predicted novel pathways involved in the inflammatory response., Conclusions: Integration of knowledge-driven and data-driven methods allows us to construct an effective IRN during IAV infection. Based on the constructed network, we have identified new interactions among inflammatory factors and biological pathways. These findings provide new insight into our understanding of the molecular mechanisms in the inflammatory network in response to IAV infection. Further characterization and experimental validation of the interaction mechanisms identified from this study may lead to a novel therapeutic strategy for the control of infections and inflammatory responses.

Background: The biochemical oscillator that controls periodic events during the Xenopus embryonic cell cycle is centered on the activity of CDKs, and the cell cycle is driven by a protein circuit that is centered on the cyclin-dependent protein kinase CDK1 and the anaphase-promoting complex (APC). Many studies have been conducted to confirm that the interactions in the cell cycle can produce oscillations and predict behaviors such as synchronization, but much less is known about how the various elaborations and collective behavior of the basic oscillators can affect the robustness of the system. Therefore, in this study, we investigate and model a multi-cell system of the Xenopus embryonic cell cycle oscillators that are coupled through a common complex protein, and then analyze their synchronization ability under four different external stimuli, including a constant input signal, a square-wave periodic signal, a sinusoidal signal and a noise signal., Results: Through bifurcation analysis and numerical simulations, we obtain synchronization intervals of the sensitive parameters in the individual oscillator and the coupling parameters in the coupled oscillators. Then, we analyze the effects of these parameters on the synchronization period and amplitude, and find interesting phenomena, e.g., there are two synchronization intervals with activation coefficient in the Hill function of the activated CDK1 that activates the Plk1, and different synchronization intervals have distinct influences on the synchronization period and amplitude. To quantify the speediness and robustness of the synchronization, we use two quantities, the synchronization time and the robustness index, to evaluate the synchronization ability. More interestingly, we find that the coupled system has an optimal signal strength that maximizes the synchronization index under different external stimuli. Simulation results also show that the ability and robustness of the synchronization for the square-wave periodic signal of cyclin synthesis is strongest in comparison to the other three different signals., Conclusions: These results suggest that the reaction process in which the activated cyclin-CDK1 activates the Plk1 has a very important influence on the synchronization ability of the coupled system, and the square-wave periodic signal of cyclin synthesis is more conducive to the synchronization and robustness of the coupled cell-cycle oscillators. Our study provides insight into the internal mechanisms of the cell cycle system and helps to generate hypotheses for further research.