BZPEERReading, trust and collaboration for researchSource author record
S. K. Nechaev appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
I review few conceptual steps in analytic description of topological interactions, which constitute the basis of a new interdisciplinary branch in mathematical physics, "Statistical Topology", emerged at the edge of topology and statistical physics of fluctuating non-phantom rope-like objects. This new branch is called statistical (or probabilistic) topology.
We consider random non-directed networks subject to dynamics conserving vertex degrees and study analytically and numerically equilibrium three-vertex motif distributions in the presence of an external field, $h$, coupled to one of the motifs. For small $h$ the numerics is well described by the "chemical kinetics" for the concentrations of motifs based on the law of mass action. For larger $h$ a transition into some trapped motif state occurs in Erdős-Rényi networks. We explain the existence of the transition by employing the notion of the entropy of the motif distribution and describe it in terms of a phenomenological Landau-type theory with a non-zero cubic term. A localization transition should always occur if the entropy function is non-convex. We conjecture that this phenomenon is the origin of the motifs' pattern formation in real evolutionary networks.
We present a model of an AB-diblock random copolymer sequential self-packaging with local quenched interactions on a one-dimensional infinite sticky substrate. It is assumed that the A-A and B-B contacts are favorable, while A-B are not. The position of a newly added monomer is selected in view of the local contact energy minimization. The model demonstrates a self-organization behavior with the nontrivial dependence of the total energy, $E$ (the number of unfavorable contacts), on the number of chain monomers, $N$: $E\sim N^{3/4}$ for quenched random equally probable distribution of A- and B-monomers along the chain. The model is treated by mapping it onto the "lamplighter" random walk and the diffusion-controlled chemical reaction of $X+X\to 0$ type with the subdiffusive motion of reagents.
We consider the ensemble of $N\times N$ ($N\gg 1$) symmetric random matrices with the bimodal independent distribution of matrix elements: each element could be either "1" with the probability $p$, or "0" otherwise. We pay attention to the "diluted" sparse regime, taking $p=1/N +ε$, where $0<ε\ll 1/N$. In this limit the eigenvalue density, $ρ(λ)$, is essentially singular, consisting of a hierarchical ultrametric set of peaks. We provide two conjectures concerning the structure of $ρ(λ)$: (i) we propose an equation for the position of sequential (in heights) peaks, and (ii) we give an expression for the shape of an outbound enveloping curve. We point out some similarities of $ρ(λ)$ with the shapes constructed on the basis of the Dedekind modular $η$-function.
The relaxation of an elastic network, constructed by a contact map of a fractal (crumpled) polymer globule is investigated. We found that: i) the slowest mode of the network is separated from the rest of the spectrum by a wide gap, and ii) the network quickly relaxes to a low--dimensional (one-dimensional, in our demonstration) manifold spanned by slowest degrees of freedom with a large basin of attraction, and then slowly approaches the equilibrium not escaping this manifold. By these dynamic properties, the fractal globule elastic network is similar to real biological molecular machines, like myosin. We have demonstrated that unfolding of a fractal globule can be described as a cascade of equilibrium phase transitions in a hierarchical system. Unfolding manifests itself in a sequential loss of stability of hierarchical levels with the temperature change.
Connectivity networks have recently become widely used in biology due to increasing amounts of information on the physical and functional links between individual proteins. This connectivity data provides valuable material for expanding our knowledge far beyond the experimentally validated via mathematical analysis and theoretical predictions of new functional interactions. In this paper we demonstrate an application of several algorithms developed for the ranking of potential gene-expression regulators within the context of an associated network. We analyze how different types of connectivity between genes and proteins affect the topology of the integral C.elegans functional network and thereby validate algorithmic performance. We demonstrate the possible definition of co-expression gene clusters within a network context from their specific motif distribution signatures. We also show that the method based on the shortest path function (SPF) applied to gene interactions sub-network of the co-expression gene cluster, efficiently predicts novel regulatory transcription factors (TFs). Simultaneous application of other methods, including only interactions with neighborhood genes, allows rapid ranking of potential regulators that could be functionally linked with the group of co-expressed genes. Predicting functions of regulators for a cluster of ribosomal/mRNA metabolic genes we highlight a role of mRNA translation and decay in a longevity of organisms.
The statistical properties of intra-chromosome maps obtained by a genome-wide chromosome conformation capture method (Hi-C) are described in the framework of the hierarchical crumpling model of heteropolymer chain with quenched disorder in the primary sequence. We conjecture that the observed Hi-C maps are statistical averages over many different ways of hierarchical genome folding, and show that the existence of quenched primary structure coupled with hierarchical folding can induce the observed fine structure of intra-chromosome contact maps.
We study the fraction $f$ of nucleotides involved in the formation of a cactus--like secondary structure of random heteropolymer RNA--like molecules. In the low--temperature limit we study this fraction as a function of the number $c$ of different nucleotide species. We show, that with changing $c$, the secondary structures of random RNAs undergo a morphological transition: $f(c)\to 1$ for $c \le c_{\rm cr}$ as the chain length $n$ goes to infinity, signaling the formation of a virtually "perfect" gapless secondary structure; while $f(c)<1$ for $c>c_{\rm cr}$, what means that a non-perfect structure with gaps is formed. The strict upper and lower bounds $2 \le c_{\rm cr} \le 4$ are proven, and the numerical evidence for $c_{\rm cr}$ is presented. The relevance of the transition from the evolutional point of view is discussed.
We propose a new toy model of a heteropolymer chain capable of forming planar secondary structures typical for RNA molecules. In this model the sequential intervals between neighboring monomers along a chain are considered as quenched random variables. Using the optimization procedure for a special class of concave--type potentials, borrowed from optimal transport analysis, we derive the local difference equation for the ground state free energy of the chain with the planar (RNA--like) architecture of paired links. We consider various distribution functions of intervals between neighboring monomers (truncated Gaussian and scale--free) and demonstrate the existence of a topological crossover from sequential to essentially embedded (nested) configurations of paired links.
We consider the dynamical system described by the area--preserving standard mapping. It is known for this system that $P(t)$, the normalized number of recurrences staying in some given domain of the phase space at time $t$ (so-clled "survival probability") has the power--law asymptotics, $P(t)\sim t^{-ν}$. We present new semi--phenomenological arguments which enable us to map the dynamical system near the chaos border onto the effective "ultrametric diffusion" on the boundary of a tree--like space with hierarchically organized transition rates. In the frameworks of our approach we have estimated the exponent $ν$ as $ν=\ln 2/\ln (1+r_g)\approx 1.44$, where $r_g=(\sqrt{5}-1)/2$ is the critical rotation number.
The distribution of motifs in random hierarchical networks defined by nonsymmetric random block--hierarchical adjacency matrices, is constructed for the first time. According to the classification of U. Alon et al of network superfamilies by their motifs distributions, our artificial directed random hierarchical networks falls into the superfamily of natural networks to which the class of neuron networks belongs. This is the first example of ``handmade'' networks with the motifs distribution as in a special class of natural networks of essential biological importance.
A new statistical method of alignment of two heteropolymers which can form hierarchical cloverleaf-like secondary structures is proposed. This offers a new constructive algorithm for quantitative determination of binding free energy of two noncoding RNAs with arbitrary primary sequences. The alignment of ncRNAs differs from the complete alignment of two RNA sequences: in ncRNA case we align only the sequences of nucleotides which constitute pairs between two different RNAs, while the secondary structure of each RNA comes into play only by the combinatorial factors affecting the entropc contribution of each molecule to the total cost function. The proposed algorithm is based on two observations: i) the standard alignment problem is considered as a zero-temperature limit of a more general statistical problem of binding of two associating heteropolymer chains; ii) this last problem is generalized onto the sequences with hierarchical cloverleaf-like structures (i.e. of RNA-type). Taking zero-temperature limit at the very end we arrive at the desired "cost function" of the system with account for entropy of side cactus-like loops. Moreover, we have demonstrated in detail how our algorithm enables to solve the "structure recovery" problem. Namely, we can predict in zero-temperature limit the cloverleaf-like (i.e. secondary) structure of interacting ncRNAs by knowing only their primary sequences.
In this paper the question about statistical properties of block--hierarchical random matrices is raised for the first time in connection with structural characteristics of random hierarchical networks obtained by mipmapping procedure. In particular, we compute numerically the spectral density of large random adjacency matrices defined by a hierarchy of the Bernoulli distributions $\{q_1,q_2,...\}$ on matrix elements, where $q_γ$ depends on hierarchy level $γ$ as $q_γ=p^{-μγ}$ ($μ>0$). For the spectral density we clearly see the free--scale behavior. We show also that for the Gaussian distributions on matrix elements with zero mean and variances $σ_γ=p^{-νγ}$, the tail of the spectral density, $ρ_G(λ)$, behaves as $ρ_G(λ) \sim |λ|^{-(2-ν)/(1-ν)}$ for $|λ|\to\infty$ and $0<ν<1$, while for $ν\ge 1$ the power--law behavior is terminated. We also find that the vertex degree distribution of such hierarchical networks has a poly--scale fractal behavior extended to a very broad range of scales.
We consider a ${\rm A}_m{\rm B}_n$ diblock copolymer, whose links are capable of forming local reversible bonds with each other. We assume that the resulting structure of the bonds is RNA--like, i.e. topologically isomorphic to a tree. We show that, depending on the relative strengths of A--A, A--B and B--B contacts, such a polymer can be in one of two different states. Namely, if a self--association is preferable (i.e., A--A and B--B bonds are comparatively stronger than A--B contacts) then the polymer forms a typical randomly branched cloverleaf structure. On the contrary, if alternating association is preferable (i.e. A--B bonds are stronger than A--A and B--B contacts) then the polymer tends to form a generally linear necklace structure (with, probably, some rear side branches and loops, which do not influence the overall characteristics of the chain). The transition between cloverleaf and necklace states is studied in details and it is shown that it is a 2nd order phase transition.