A Modified Ising Model of Barab\'asi-Albert Network with Gene-type Spins

TL;DR

This paper introduces a modified asymmetric Ising model tailored for biological networks, particularly scale-free Barabási-Albert networks, to analyze phase transitions and emergent behaviors in gene regulation systems.

Contribution

It develops a novel asymmetric Ising model for biological systems, adapting statistical physics techniques to network models that resemble real biological networks.

Findings

Asymmetric Ising models exhibit first-order phase transitions and hysteresis.

The model's simulation setup is applicable to real biological network data.

The approach is scalable to large biological networks.

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 emerging behavior of biological systems. Here we will 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. Hence we will 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 mean-field solution of modified Ising model of a network type that closely resembles real-world network, the Barab\'{a}si-Albert model of scale-free networks. We show that asymmetric Ising models show similarities to symmetric Ising models with external field and undergoes a discontinuous phase transition of the first-order and exhibits hysteresis. The simulation setup presented here can be directly used for any biological network connectivity dataset and is also applicable for other networks that exhibit similar states of activity. This 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.