Stability of Signaling Pathways during Aging-A Boolean Network Approach

Leveraging quantum computing for dynamic analyses of logical networks in systems biology

The dynamics of cellular mechanisms can be investigated through the analysis of networks. One of the simplest but most popular modeling strategies involves logic-based models. However, these models still face exponential growth in simulation complexity compared with a linear increase in nodes. We transfer this modeling approach to quantum computing and use the upcoming technique in the field to simulate the resulting networks. Leveraging logic modeling in quantum computing has many benefits, including complexity reduction and quantum algorithms for systems biology tasks. To showcase the applicability of our approach to systems biology tasks, we implemented a model of mammalian cortical development. Here, we applied a quantum algorithm to estimate the tendency of the model to reach particular stable conditions and further revert dynamics. Results from two actual quantum processing units and a noisy simulator are presented, and current technical challenges are discussed.

Review and assessment of Boolean approaches for inference of gene regulatory networks

Boolean descriptions of gene regulatory networks can provide an insight into interactions between genes. Boolean networks hold predictive power, are easy to understand, and can be used to simulate the observed networks in different scenarios. We review fundamental and state-of-the-art methods for inference of Boolean networks. We introduce a methodology for a straightforward evaluation of Boolean inference approaches based on the generation of evaluation datasets, application of selected inference methods, and evaluation of performance measures to guide the selection of the best method for a given inference problem. We demonstrate this procedure on inference methods REVEAL (REVerse Engineering ALgorithm), Best-Fit Extension, MIBNI (Mutual Information-based Boolean Network Inference), GABNI (Genetic Algorithm-based Boolean Network Inference) and ATEN (AND/OR Tree ENsemble algorithm), which infers Boolean descriptions of gene regulatory networks from discretised time series data. Boolean inference approaches tend to perform better in terms of dynamic accuracy, and slightly worse in terms of structural correctness. We believe that the proposed methodology and provided guidelines will help researchers to develop Boolean inference approaches with a good predictive capability while maintaining structural correctness and biological relevance.

Discrete Logic Modeling of Cell Signaling Pathways

Cell signaling pathways often crosstalk generating complex biological behaviors observed in different cellular contexts. Frequently, laboratory experiments focus on a few putative regulators, alone unable to predict the molecular mechanisms behind the observed phenotypes. Here, systems biology complements these approaches by giving a holistic picture to complex signaling crosstalk. In particular, Boolean network models are a meaningful tool to study large network behaviors and can cope with incomplete kinetic information. By introducing a model describing pathways involved in hematopoietic stem cell maintenance, we present a general approach on how to model cell signaling pathways with Boolean network models.

Systems approaches to investigate the role of NF-κB signaling in aging

The nuclear factor-κB (NF-κB) signaling pathway is one of the most well-studied pathways related to inflammation, and its involvement in aging has attracted considerable attention. As aging is a complex phenomenon and is the result of a multi-step process, the involvement of the NF-κB pathway in aging remains unclear. To elucidate the role of NF-κB in the regulation of aging, different systems biology approaches have been employed. A multi-omics data-driven approach can be used to interpret and clarify unknown mechanisms but cannot generate mechanistic regulatory structures alone. In contrast, combining this approach with a mathematical modeling approach can identify the mechanistics of the phenomena of interest. The development of single-cell technologies has also helped clarify the heterogeneity of the NF-κB response and underlying mechanisms. Here, we review advances in the understanding of the regulation of aging by NF-κB by focusing on omics approaches, single-cell analysis, and mathematical modeling of the NF-κB network.

Reconstructing Boolean network ensembles from single-cell data for unraveling dynamics in the aging of human hematopoietic stem cells

Regulatory dependencies in molecular networks are the basis of dynamic behaviors affecting the phenotypical landscape. With the advance of high throughput technologies, the detail of omics data has arrived at the single-cell level. Nevertheless, new strategies are required to reconstruct regulatory networks based on populations of single-cell data. Here, we present a new approach to generate populations of gene regulatory networks from single-cell RNA-sequencing (scRNA-seq) data. Our approach exploits the heterogeneity of single-cell populations to generate pseudo-timepoints. This allows for the first time to uncouple network reconstruction from a direct dependency on time series measurements. The generated time series are then fed to a combined reconstruction algorithm. The latter allows a fast and efficient reconstruction of ensembles of gene regulatory networks. Since this approach does not require knowledge on time-related trajectories, it allows us to model heterogeneous processes such as aging. Applying the approach to the aging-associated NF-κB signaling pathway-based scRNA-seq data of human hematopoietic stem cells (HSCs), we were able to reconstruct eight ensembles, and evaluate their dynamic behavior. Moreover, we propose a strategy to evaluate the resulting attractor patterns. Interaction graph-based features and dynamic investigations of our model ensembles provide a new perspective on the heterogeneity and mechanisms related to human HSCs aging.

Network Simulations Reveal Molecular Signatures of Vulnerability to Age-Dependent Stress and Tau Accumulation

Alzheimer’s disease (AD) is the leading cause of dementia and one of the most common causes of death worldwide. As an age-dependent multifactorial disease, the causative triggers of AD are rooted in spontaneous declines in cellular function and metabolic capacity with increases in protein stressors such as the tau protein. This multitude of age-related processes that cause neurons to change from healthy states to ones vulnerable to the damage seen in AD are difficult to simultaneously investigate and even more difficult to quantify. Here we aimed to diminish these gaps in our understanding of neuronal vulnerability in AD development by using simulation methods to theoretically quantify an array of cellular stress responses and signaling molecules. This temporally-descriptive molecular signature was produced using a novel multimethod simulation approach pioneered by our laboratory for biological research; this methodology combines hierarchical agent-based processes and continuous equation-based modeling in the same interface, all while maintaining intrinsic distributions that emulate natural biological stochasticity. The molecular signature was validated for a normal organismal aging trajectory using experimental longitudinal data from Caenorhabditis elegans and rodent studies. In addition, we have further predicted this aging molecular signature for cells impacted by the pathogenic tau protein, giving rise to distinct stress response conditions needed for cytoprotective aging. Interestingly, our simulation experiments showed that oxidative stress signaling (via daf-16 and skn-1 activities) does not substantially protect cells from all the early stressors of aging, but that it is essential in preventing a late-life degenerative cellular phenotype. Together, our simulation experiments aid in elucidating neurodegenerative triggers in the onset of AD for different genetic conditions. The long-term goal of this work is to provide more detailed diagnostic and prognostic tools for AD development and progression, and to provide more comprehensive preventative measures for this disease.

Concepts in Boolean network modeling: What do they all mean?

Boolean network models are one of the simplest models to study complex dynamic behavior in biological systems. They can be applied to unravel the mechanisms regulating the properties of the system or to identify promising intervention targets. Since its introduction by Stuart Kauffman in 1969 for describing gene regulatory networks, various biologically based networks and tools for their analysis were developed. Here, we summarize and explain the concepts for Boolean network modeling. We also present application examples and guidelines to work with and analyze Boolean network models.

A Robustness Analysis of Dynamic Boolean Models of Cellular Circuits

With ever growing amounts of omics data, the next challenge in biological research is the interpretation of these data to gain mechanistic insights about cellular function. Dynamic models of cellular circuits that capture the activity levels of proteins and other molecules over time offer great expressive power by allowing the simulation of the effects of specific internal or external perturbations on the workings of the cell. However, the study of such models is at its infancy and no large-scale analysis of the robustness of real models to changing conditions has been conducted to date. Here we provide a computational framework to study the robustness of such models using a combination of stochastic simulations and integer linear programming techniques. We apply our framework to a large collection of cellular circuits and benchmark the results against randomized models. We find that the steady states of real circuits tend to be more robust in multiple aspects compared with their randomized counterparts.

Identification via Numerical Computation of Transcriptional Determinants of a Cell Phenotype Decision Making

Complex cellular processes, such as phenotype decision making, are exceedingly difficult to analyze experimentally, due to the multiple-layer regulation of gene expression and the intercellular variability referred to as biological noise. Moreover, the heterogeneous experimental approaches used to investigate distinct macromolecular species, and their intrinsic differential time-scale dynamics, add further intricacy to the general picture of the physiological phenomenon. In this respect, a computational representation of the cellular functions of interest can be used to extract relevant information, being able to highlight meaningful active markers within the plethora of actors forming an active molecular network. The multiscale power of such an approach can also provide meaningful descriptions for both population and single-cell level events. To validate this paradigm a Boolean and a Markov model were combined to identify, in an objective and user-independent manner, a signature of genes recapitulating epithelial to mesenchymal transition in-vitro. The predictions of the model are in agreement with experimental data and revealed how the expression of specific molecular markers is related to distinct cell behaviors. The presented method strengthens the evidence of a role for computational representation of active molecular networks to gain insight into cellular physiology and as a general approach for integrating in-silico/in-vitro study of complex cell population dynamics to identify their most relevant drivers.

Automatic Screening for Perturbations in Boolean Networks

A common approach to address biological questions in systems biology is to simulate regulatory mechanisms using dynamic models. Among others, Boolean networks can be used to model the dynamics of regulatory processes in biology. Boolean network models allow simulating the qualitative behavior of the modeled processes. A central objective in the simulation of Boolean networks is the computation of their long-term behavior—so-called attractors. These attractors are of special interest as they can often be linked to biologically relevant behaviors. Changing internal and external conditions can influence the long-term behavior of the Boolean network model. Perturbation of a Boolean network by stripping a component of the system or simulating a surplus of another element can lead to different attractors. Apparently, the number of possible perturbations and combinations of perturbations increases exponentially with the size of the network. Manually screening a set of possible components for combinations that have a desired effect on the long-term behavior can be very time consuming if not impossible. We developed a method to automatically screen for perturbations that lead to a user-specified change in the network's functioning. This method is implemented in the visual simulation framework ViSiBool utilizing satisfiability (SAT) solvers for fast exhaustive attractor search.