Steady state approach to model gene regulatory networks—Simulation of microarray experiments

Biomaterialomics: Data science-driven pathways to develop fourth-generation biomaterials

Conventional approaches to developing biomaterials and implants require intuitive tailoring of manufacturing protocols and biocompatibility assessment. This leads to longer development cycles, and high costs. To meet existing and unmet clinical needs, it is critical to accelerate the production of implantable biomaterials, implants and biomedical devices. Building on the Materials Genome Initiative, we define the concept 'biomaterialomics' as the integration of multi-omics data and high-dimensional analysis with artificial intelligence (AI) tools throughout the entire pipeline of biomaterials development. The Data Science-driven approach is envisioned to bring together on a single platform, the computational tools, databases, experimental methods, machine learning, and advanced manufacturing (e.g., 3D printing) to develop the fourth-generation biomaterials and implants, whose clinical performance will be predicted using 'digital twins'. While analysing the key elements of the concept of 'biomaterialomics', significant emphasis has been put forward to effectively utilize high-throughput biocompatibility data together with multiscale physics-based models, E-platform/online databases of clinical studies, data science approaches, including metadata management, AI/ Machine Learning (ML) algorithms and uncertainty predictions. Such integrated formulation will allow one to adopt cross-disciplinary approaches to establish processing-structure-property (PSP) linkages. A few published studies from the lead author's research group serve as representative examples to illustrate the formulation and relevance of the 'Biomaterialomics' approaches for three emerging research themes, i.e. patient-specific implants, additive manufacturing, and bioelectronic medicine. The increased adaptability of AI/ML tools in biomaterials science along with the training of the next generation researchers in data science are strongly recommended. STATEMENT OF SIGNIFICANCE: This leading opinion review paper emphasizes the need to integrate the concepts and algorithms of the data science with biomaterials science. Also, this paper emphasizes the need to establish a mathematically rigorous cross-disciplinary framework that will allow a systematic quantitative exploration and curation of critical biomaterials knowledge needed to drive objectively the innovation efforts within a suitable uncertainty quantification framework, as embodied in 'biomaterialomics' concept, which integrates multi-omics data and high-dimensional analysis with artificial intelligence (AI) tools, like machine learning. The formulation of this approach has been demonstrated for patient-specific implants, additive manufacturing, and bioelectronic medicine.

Computational modelling of the inflammatory response in trauma, sepsis and wound healing: implications for modelling resilience

Resilience refers to the ability to recover from illness or adversity. At the cell, tissue, organ and whole-organism levels, the response to perturbations such as infections and injury involves the acute inflammatory response, which in turn is connected to and controlled by changes in physiology across all organ systems. When coordinated properly, inflammation can lead to the clearance of infection and healing of damaged tissues. However, when either overly or insufficiently robust, inflammation can drive further cell stress, tissue damage, organ dysfunction and death through a feed-forward process of inflammation → damage → inflammation. To address this complexity, we have obtained extensive datasets regarding the dynamics of inflammation in cells, animals and patients, and created data-driven and mechanistic computational simulations of inflammation and its recursive effects on tissue, organ and whole-organism (patho)physiology. Through this approach, we have discerned key regulatory mechanisms, recapitulated in silico key features of clinical trials for acute inflammation and captured diverse, patient-specific outcomes. These insights may allow for the determination of individual-specific tolerances to illness and adversity, thereby defining the role of inflammation in resilience.

Steady state analysis of the genetic regulatory network incorporating underlying molecular mechanisms for anaerobic metabolism in Escherichia coli

A Gene Regulatory Network (GRN) represents complex connections between genes in a cell which interact with each other through their RNA and protein expression products, thereby determining the expression levels of mRNA and proteins required for functioning of the cell. Microarray experiments yield the log fold change in mRNA abundance and quantify the expression levels for a GRN at the genome level. While Boolean or Bayesian modeling along with expression and location data are useful in analyzing microarray data, they lack underlying mechanistic details present in GRNs. Our objective is to understand the role of molecular mechanisms in quantifying a GRN. To that effect, we analyze under steady state, the complete GRN for the central metabolic pathway during anaerobiosis in Escherichia coli. We simulate the microarray experiments using a steady state gene expression simulator (SSGES) that models molecular mechanistic details such as dimerization, multiple-site binding, auto-regulation and feedback. Given a GRN, the SSGES provided the log fold change in mRNA expression values as the output, which can be compared to data from microarray experiments. We predict the log fold changes for mutants obtained by knocking out crucial transcriptional regulators such as FNR (F), ArcA (A), IHFA-B (I) and DpiA (D) and observe a high degree of correlation with previously reported experimental data. We also predict the microarray expression values for hitherto unknown combinations of deletion mutants. We hierarchically cluster the predicted log fold change values for these mutants and postulate that E. coli has evolved from a predominantly lactate secreting (FAID mutant) into a mixed acid secreting phenotype as seen in the wild type (WT) during anaerobiosis. Upon simulating a model without incorporating the mechanistic details, not only the correlation with the experimental data reduced considerably, but also the clustering of expression data indicated WT to be closer to the quadruple mutant FAID. This clearly demonstrates the significance of incorporating mechanistic data while quantifying the expression profile of a GRN which can help in predicting the effect of a gene mutant and understanding the evolution of transcriptional control.

Integration of steady-state and temporal gene expression data for the inference of gene regulatory networks

We develop a new regression algorithm, cMIKANA, for inference of gene regulatory networks from combinations of steady-state and time-series gene expression data. Using simulated gene expression datasets to assess the accuracy of reconstructing gene regulatory networks, we show that steady-state and time-series data sets can successfully be combined to identify gene regulatory interactions using the new algorithm. Inferring gene networks from combined data sets was found to be advantageous when using noisy measurements collected with either lower sampling rates or a limited number of experimental replicates. We illustrate our method by applying it to a microarray gene expression dataset from human umbilical vein endothelial cells (HUVECs) which combines time series data from treatment with growth factor TNF and steady state data from siRNA knockdown treatments. Our results suggest that the combination of steady-state and time-series datasets may provide better prediction of RNA-to-RNA interactions, and may also reveal biological features that cannot be identified from dynamic or steady state information alone. Finally, we consider the experimental design of genomics experiments for gene regulatory network inference and show that network inference can be improved by incorporating steady-state measurements with time-series data.

Inference of biological pathway from gene expression profiles by time delay boolean networks

One great challenge of genomic research is to efficiently and accurately identify complex gene regulatory networks. The development of high-throughput technologies provides numerous experimental data such as DNA sequences, protein sequence, and RNA expression profiles makes it possible to study interactions and regulations among genes or other substance in an organism. However, it is crucial to make inference of genetic regulatory networks from gene expression profiles and protein interaction data for systems biology. This study will develop a new approach to reconstruct time delay boolean networks as a tool for exploring biological pathways. In the inference strategy, we will compare all pairs of input genes in those basic relationships by their corresponding p-scores for every output gene. Then, we will combine those consistent relationships to reveal the most probable relationship and reconstruct the genetic network. Specifically, we will prove that O(log n) state transition pairs are sufficient and necessary to reconstruct the time delay boolean network of n nodes with high accuracy if the number of input genes to each gene is bounded. We also have implemented this method on simulated and empirical yeast gene expression data sets. The test results show that this proposed method is extensible for realistic networks.

A glance at DNA microarray technology and applications

INTRODUCTION

Because of huge impacts of "OMICS" technologies in life sciences, many researchers aim to implement such high throughput approach to address cellular and/or molecular functions in response to any influential intervention in genomics, proteomics, or metabolomics levels. However, in many cases, use of such technologies often encounters some cybernetic difficulties in terms of knowledge extraction from a bunch of data using related softwares. In fact, there is little guidance upon data mining for novices. The main goal of this article is to provide a brief review on different steps of microarray data handling and mining for novices and at last to introduce different PC and/or web-based softwares that can be used in preprocessing and/or data mining of microarray data.

METHODS

To pursue such aim, recently published papers and microarray softwares were reviewed.

RESULTS

It was found that defining the true place of the genes in cell networks is the main phase in our understanding of programming and functioning of living cells. This can be obtained with global/selected gene expression profiling.

CONCLUSION

Studying the regulation patterns of genes in groups, using clustering and classification methods helps us understand different pathways in the cell, their functions, regulations and the way one component in the system affects the other one. These networks can act as starting points for data mining and hypothesis generation, helping us reverse engineer.

Short time-series microarray analysis: methods and challenges

The detection and analysis of steady-state gene expression has become routine. Time-series microarrays are of growing interest to systems biologists for deciphering the dynamic nature and complex regulation of biosystems. Most temporal microarray data only contain a limited number of time points, giving rise to short-time-series data, which imposes challenges for traditional methods of extracting meaningful information. To obtain useful information from the wealth of short-time series data requires addressing the problems that arise due to limited sampling. Current efforts have shown promise in improving the analysis of short time-series microarray data, although challenges remain. This commentary addresses recent advances in methods for short-time series analysis including simplification-based approaches and the integration of multi-source information. Nevertheless, further studies and development of computational methods are needed to provide practical solutions to fully exploit the potential of this data.