Statistical mechanics in biology: how ubiquitous are long-range correlations?

doi:10.1016/0378-4371(94)90502-9

Physica A: Statistical Mechanics and its Applications

Volume 205, Issues 1–3, 1 April 1994, Pages 214-253

https://doi.org/10.1016/0378-4371(94)90502-9 Get rights and content

Abstract

The purpose of this opening talk is to describe examples of recent progress in applying statistical mechanics to biological systems. We first briefly review several biological systems, and then focus on the fractal features characterized by the long-range correlations found recently in DNA sequences containing non-coding material. We discuss the evidence supporting the finding that for sequences containing only coding regions, there are no long-range correlations. We also discuss the recent finding that the exponent α characterizing the long-range correlations increases with evolution, and we discuss two related models, the insertion model and the insertion-deletion model, that may account for the presence of long-range correlations. Finally, we summarize the analysis of long-term data on human heartbeats (up to 10⁴ heart beats) that supports the possibility that the successive increments in the cardiac beat-to-beat intervals of healthy subjects display scale-invariant, long-range “anti-correlations” (a tendency to beat faster is balanced by a tendency to beat slower later on). In contrast, for a group of subjects with severe heart disease, long-range correlations vanish. This finding suggests that the classical theory of homeostasis, according to which stable physiological processes seek to maintain “constancy,” should be extended to account for this type of dynamical, far from equilibrium, behavior.

References (66)

T. Vicsek et al.
Physica A
(1990)
S. Matsuura et al.
Physica A
(1992)
D. Stauffer et al.
From Newton to Mandelbrot: A Primer in Theoretical Physics
(1990)
E. Guyon et al.
T. VicsekJ. Feder
Fractals
(1988)
Proc. 1985 Cargèse NATO ASI
P. Bak et al.
C.-K. Peng et al.
Physica A
(1991)
C.-K. Peng et al.
Phys. Rev. A
(1991)
R.N. Mantegna
Physica A
(1991)
C.-K. Peng et al.
Nature
(1992)
A.S. Borovik, A.Yu. Grosberg and M.D. Frank-Kamenetskii, Nature, in...

H. Fujikawa et al.

J. Phys. Soc. Jpn.

(1989)

H. Fujikawa et al.

J. Phys. Soc. Jpn.

(1991)

Physica A

(1990)

B.B. Mandelbrot

(1992)

C.-K. Peng et al.

Phys. Rev. Lett.

(1993)

(1992)

I. Amato

Science

(1992)

Cited by (145)

Nonextensive realizations in interacting ion channels: Implications for mechano-electrical transducer mechanisms
2023, BioSystems
We propose a theoretical model to investigate the thermodynamics of single and coupled two-state ion channels, associated with mechanoelectrical transduction ( $M E T$ ) and hair cell biophysics. The modeling was based on the Tsallis nonextensive statistical mechanics. The choice for a nonextensive framework in modeling ion channels is encouraged on the fact that we take into account the presence of interactions or long-range correlations in the dynamics of single and coupled ion channels. However, the basic assumptions that support Boltzmann–Gibbs statistics, traditionally used to model ion channel dynamics, state that the system is formed by independent or weakly interacting elements. Despite being well studied in many biological systems, the literature has not yet addressed the study of both entropy and mutual information related to isolated or physically interacting pairs of $M E T$ channels. Inspired by hair cell biophysics, we show how the presence of nonextensivity, or subadditivity and superadditivity modulates the nonextensive entropy and mutual information as functions of stereocilia displacements. We also observe that the magnitude of the interaction between the two channels, given by a nonextensive parameter, influences the amplitude of the nonextensive joint entropy and mutual information as functions of the hair cell displacements. Finally, we show how nonextensivity regulates the current versus displacement curve for a single and a pair of interacting two-state channels. The present findings shed light on the thermodynamic process involved in the molecular mechanisms of the auditory system.
The q-exponentials do not maximize the Rényi entropy
2021, Physica A: Statistical Mechanics and its Applications
It is generally assumed that the Rényi entropy is maximized by the $q$ -exponentials and is hence useful to construct a generalized statistical mechanics. However, to the best of our knowledge, this assumption has never been explicitly checked. In this work, we consider the Rényi entropy with the linear and escort mean value constraints and check whether it is indeed maximized by $q$ -exponentials. We show, both theoretically and numerically, that the Rényi entropy yields erroneous inferences concerning the optimum distributions of the $q$ -exponential form and moreover exhibits high estimation errors in the regime of long range correlations. Finally, we note that the Shannon entropy successfully detects the power law distributions when the logarithmic mean value constraint is used.
Spatial constrains and information content of sub-genomic regions of the human genome
2021, iScience
Complexity metrics and machine learning (ML) models have been utilized to analyze the lengths of segmental genomic entities of DNA sequences (exonic, intronic, intergenic, repeat, unique) with the purpose to ask questions regarding the segmental organization of the human genome within the size distribution of these sequences. For this we developed an integrated methodology that is based upon the reconstructed phase space theorem, the non-extensive statistical theory of Tsallis, ML techniques, and a technical index, integrating the generated information, which we introduce and named complexity factor (COFA). Our analysis revealed that the size distribution of the genomic regions within chromosomes are not random but follow patterns with characteristic features that have been seen through its complexity character, and it is part of the dynamics of the whole genome. Finally, this picture of dynamics in DNA is recognized using ML tools for clustering, classification, and prediction with high accuracy.
An alternative description of power law correlations in DNA sequences
2020, Physica A: Statistical Mechanics and its Applications
We analyze the coding sequence for the Homo Sapiens via a model which naturally embraces power law correlations (PLC) among the bases in DNA sequences of living organisms. This model is based on a principle of universal optimization, which is the core of all statistical arguments, being associated with the power law distribution function of the length of DNA, measured in base pairs (bp). This distribution provides a PLC parameter introduced through a nonadditive framework in which such parameter measures the PLC in the DNA sequence. The results show that the Short-Range-Correlations (SRC), always present in coding DNA sequences, are appropriately captured through the power law distribution, adequately describing the cumulative length distribution of DNA bases, in contrast with the case of the traditional exponential statistical model. We use an Empirical cumulative distribution function and the database of proteins compiled by the Ensembl Project to show that the power law distribution provides the best description of the data. A Bayesian analysis of the data further confirms this result.
Walsh code based numerical mapping method for the identification of protein coding regions in eukaryotes
2020, Biomedical Signal Processing and Control
The protein coding regions play a significant role for gene applications in genomic signal processing. Unlike prokaryotes, the coding regions in eukaryotes are arranged in a random manner. Owing to unequal lengths and low volume density of coding regions, the identification of coding regions makes cumbersome. In this work, a new numerical mapping method based on Walsh codes is proposed to detect the coding regions in eukaryotes. The Walsh code for each nucleotide is obtained using the statistical features of a DNA sequence. The proposed method uses static type of mapping to convert a string of DNA nucleotides into a numerical sequence. The numerical sequence is given as an input to the digital signal processing based spectrum identification tool to detect the existence of quasi-periodic components within the coding region. The advantage of our method is that it is simple to design and easy to represent. The performance of the proposed method has been tested on four benchmark databases and a random set of sequences collected from the National Center for Biological Information (NCBI) database. Furthermore, it has been compared with other state-of-the-art spectrum based numerical mapping methods for statistical features such as sensitivity, specificity and accuracy. The proposed method is efficient as it attains 94 % accuracy, 85 % sensitivity and 96 % specificity when tested on the benchmark C. Elegans gene sequence.
Assessing information content and interactive relationships of subgenomic DNA sequences of the MHC using complexity theory approaches based on the non-extensive statistical mechanics
2018, Physica A: Statistical Mechanics and its Applications
Citation Excerpt :
Many researchers have tried to understand a series of fundamental problems such as the identification of motifs in the noncoding and coding regions or the role of intergenic regions in the function of DNA by means of methods, which generally fall into the category of the complexity theory concepts [12–14]. Furthermore, many scientists [12,15–21] have introduced and studied the DNA random walk process as a basic physical process-model for the detection and understanding of the observed long-range correlations of nucleotides in DNA sequences. According to Voss [15] there is no significant difference between long-range correlation properties of coding and noncoding sequences.
This study combines two independent domains of science, the high throughput DNA sequencing capabilities of Genomics and complexity theory from Physics, to assess the information encoded by the different genomic segments of exonic, intronic and intergenic regions of the Major Histocompatibility Complex (MHC) and identify possible interactive relationships. The dynamic and non-extensive statistical characteristics of two well characterized MHC sequences from the homozygous cell lines, PGF and COX, in addition to two other genomic regions of comparable size, used as controls, have been studied using the reconstructed phase space theorem and the non-extensive statistical theory of Tsallis. The results reveal similar non-linear dynamical behavior as far as complexity and self-organization features. In particular, the low-dimensional deterministic nonlinear chaotic and non-extensive statistical character of the DNA sequences was verified with strong multifractal characteristics and long-range correlations. The nonlinear indices repeatedly verified that MHC sequences, whether exonic, intronic or intergenic include varying levels of information and reveal an interaction of the genes with intergenic regions, whereby the lower the number of genes in a region, the less the complexity and information content of the intergenic region. Finally we showed the significance of the intergenic region in the production of the DNA dynamics. The findings reveal interesting content information in all three genomic elements and interactive relationships of the genes with the intergenic regions. The results most likely are relevant to the whole genome and not only to the MHC. These findings are consistent with the ENCODE project, which has now established that the non-coding regions of the genome remain to be of relevance, as they are functionally important and play a significant role in the regulation of expression of genes and coordination of the many biological processes of the cell.

View all citing articles on Scopus

View full text