A modified Ising model of Barabási–Albert network with gene-type spins

Abstract

The central question of systems biology is to understand how individual components of a biological system such as genes or proteins cooperate in emerging phenotypes resulting in the evolution of diseases. As living cells are open systems in quasi-steady state type equilibrium in continuous exchange with their environment, computational techniques that have been successfully applied in statistical thermodynamics to describe phase transitions may provide new insights to the emerging behavior of biological systems. Here we systematically evaluate the translation of computational techniques from solid-state physics to network models that closely resemble biological networks and develop specific translational rules to tackle problems unique to living systems. We focus on logic models exhibiting only two states in each network node. Motivated by the apparent asymmetry between biological states where an entity exhibits boolean states i.e. is active or inactive, we present an adaptation of symmetric Ising model towards an asymmetric one fitting to living systems here referred to as the modified Ising model with gene-type spins. We analyze phase transitions by Monte Carlo simulations and propose a mean-field solution of a modified Ising model of a network type that closely resembles a real-world network, the Barabási–Albert model of scale-free networks. We show that asymmetric Ising models show similarities to symmetric Ising models with the external field and undergoes a discontinuous phase transition of the first-order and exhibits hysteresis. The simulation setup presented herein can be directly used for any biological network connectivity dataset and is also applicable for other networks that exhibit similar states of activity. The method proposed here is a general statistical method to deal with non-linear large scale models arising in the context of biological systems and is scalable to any network size.

Critical transitions in degree mixed networks: A discovery of forbidden tipping regions in networked spin systems

“…What is more, many networked systems cannot realistically be assumed to have such simple grid like network structures. Components in these systems may instead be more irregularly connected forming inhomogeneous network structures with varying degrees [13].…”

A Review of Mathematical and Computational Methods in Cancer Dynamics

“…Further, many complex cancer processes including phosphoproteomics (on/off switches of protein conformations) and epigenetic-chromatin states (e.g., acetylation/methylation dynamics of histones) can be defined as binary states. As such, Ising models are simple yet powerful tools to study cell fate transitions from such complex gene and protein network dynamics ( 54 , 127 ). However, even the optimization of a two-dimensional Ising model is NP-hard.…”

SSI: A Statistical Sensitivity Index for Chemical Reaction Networks in cancer

“…Following an assignment of 419 initial conditions, the solution of the system of ODEs provides the simulated time course of the concentrations, and the (possibly allowed) asymptotic stationary state, to be identified with an equilibrium state. Numerical simulations may be used to predict changes in the model response resulting from varied external signals or internal conditions (Zi, 2011), to put into evidence feedback effects (Santra, 2018; Sommariva and others , 2021 b ), and to determine the impact of the introduction of targeted drugs (Santra, 2018; Krishnan and others , 2020; Sommariva and others , 2021 b ).…”