Spectral analysis and the dynamic response of complex networks

Spectral properties of hyperbolic nanonetworks with tunable aggregation of simplexes

Cooperative self-assembly is a ubiquitous phenomenon found in natural systems which is used for designing nanostructured materials with new functional features. Its origin and mechanisms, leading to improved functionality of the assembly, have attracted much attention from researchers in many branches of science and engineering. These complex structures often come with hyperbolic geometry; however, the relation between the hyperbolicity and their spectral and dynamical properties remains unclear. Using the model of aggregation of simplexes introduced by Šuvakov et al. [Sci. Rep. 8, 1987 (2018)2045-232210.1038/s41598-018-20398-x], here we study topological and spectral properties of a large class of self-assembled structures or nanonetworks consisting of monodisperse building blocks (cliques of size n=3,4,5,6) which self-assemble via sharing the geometrical shapes of a lower order. The size of the shared substructure is tuned by varying the chemical affinity ν such that for significant positive ν sharing the largest face is the most probable, while for ν<0, attaching via a single node dominates. Our results reveal that, while the parameter of hyperbolicity remains δ_{max}=1 across the assemblies, their structure and spectral dimension d_{s} vary with the size of cliques n and the affinity when ν≥0. In this range, we find that d_{s}>4 can be reached for n≥5 and sufficiently large ν. For the aggregates of triangles and tetrahedra, the spectral dimension remains in the range d_{s}∈[2,4), as well as for the higher cliques at vanishing affinity. On the other end, for ν<0, we find d_{s}≂1.57 independently on n. Moreover, the spectral distribution of the normalized Laplacian eigenvalues has a characteristic shape with peaks and a pronounced minimum, representing the hierarchical architecture of the simplicial complexes. These findings show how the structures compatible with complex dynamical properties can be assembled by controlling the higher-order connectivity among the building blocks.

Inferring pathogen-host interactions between Leptospira interrogans and Homo sapiens using network theory

Leptospirosis is the most emerging zoonotic disease of epidemic potential caused by pathogenic species of Leptospira. The bacterium invades the host system and causes the disease by interacting with the host proteins. Analyzing these pathogen-host protein interactions (PHPIs) may provide deeper insight into the disease pathogenesis. For this analysis, inter-species as well as intra-species protein interactions networks of Leptospira interrogans and human were constructed and investigated. The topological analyses of these networks showed lesser connectivity in inter-species network than intra-species, indicating the perturbed nature of the inter-species network. Hence, it can be one of the reasons behind the disease development. A total of 35 out of 586 PHPIs were identified as key interactions based on their sub-cellular localization. Two outer membrane proteins (GpsA and MetXA) and two periplasmic proteins (Flab and GlyA) participating in PHPIs were found conserved in all pathogenic, intermediate and saprophytic spp. of Leptospira. Furthermore, the bacterial membrane proteins involved in PHPIs were found playing major roles in disruption of the immune systems and metabolic processes within host and thereby causing infectious disease. Thus, the present results signify that the membrane proteins participating in such interactions hold potential to serve as effective immunotherapeutic candidates for vaccine development.

Network spectra for drug-target identification in complex diseases: new guns against old foes

The fundamental understanding of altered complex molecular interactions in a diseased condition is the key to its cure. The overall functioning of these molecules is kind of jugglers play in the cell orchestra and to anticipate these relationships among the molecules is one of the greatest challenges in modern biology and medicine. Network science turned out to be providing a successful and simple platform to understand complex interactions among healthy and diseased tissues. Furthermore, much information about the structure and dynamics of a network is concealed in the eigenvalues of its adjacency matrix. In this review, we illustrate rapid advancements in the field of network science in combination with spectral graph theory that enables us to uncover the complexities of various diseases. Interpretations laid by network science approach have solicited insights into molecular relationships and have reported novel drug targets and biomarkers in various complex diseases.

Interaction paths promote module integration and network-level robustness of spliceosome to cascading effects

The functionality of distinct types of protein networks depends on the patterns of protein-protein interactions. A problem to solve is understanding the fragility of protein networks to predict system malfunctioning due to mutations and other errors. Spectral graph theory provides tools to understand the structural and dynamical properties of a system based on the mathematical properties of matrices associated with the networks. We combined two of such tools to explore the fragility to cascading effects of the network describing protein interactions within a key macromolecular complex, the spliceosome. Using S. cerevisiae as a model system we show that the spliceosome network has more indirect paths connecting proteins than random networks. Such multiplicity of paths may promote routes to cascading effects to propagate across the network. However, the modular network structure concentrates paths within modules, thus constraining the propagation of such cascading effects, as indicated by analytical results from the spectral graph theory and by numerical simulations of a minimal mathematical model parameterized with the spliceosome network. We hypothesize that the concentration of paths within modules favors robustness of the spliceosome against failure, but may lead to a higher vulnerability of functional subunits, which may affect the temporal assembly of the spliceosome. Our results illustrate the utility of spectral graph theory for identifying fragile spots in biological systems and predicting their implications.

Spectral properties of complex networks

This review presents an account of the major works done on spectra of adjacency matrices drawn on networks and the basic understanding attained so far. We have divided the review under three sections: (a) extremal eigenvalues, (b) bulk part of the spectrum, and (c) degenerate eigenvalues, based on the intrinsic properties of eigenvalues and the phenomena they capture. We have reviewed the works done for spectra of various popular model networks, such as the Erdős-Rényi random networks, scale-free networks, 1-d lattice, small-world networks, and various different real-world networks. Additionally, potential applications of spectral properties for natural processes have been reviewed.

Understanding cancer complexome using networks, spectral graph theory and multilayer framework

Cancer complexome comprises a heterogeneous and multifactorial milieu that varies in cytology, physiology, signaling mechanisms and response to therapy. The combined framework of network theory and spectral graph theory along with the multilayer analysis provides a comprehensive approach to analyze the proteomic data of seven different cancers, namely, breast, oral, ovarian, cervical, lung, colon and prostate. Our analysis demonstrates that the protein-protein interaction networks of the normal and the cancerous tissues associated with the seven cancers have overall similar structural and spectral properties. However, few of these properties implicate unsystematic changes from the normal to the disease networks depicting difference in the interactions and highlighting changes in the complexity of different cancers. Importantly, analysis of common proteins of all the cancer networks reveals few proteins namely the sensors, which not only occupy significant position in all the layers but also have direct involvement in causing cancer. The prediction and analysis of miRNAs targeting these sensor proteins hint towards the possible role of these proteins in tumorigenesis. This novel approach helps in understanding cancer at the fundamental level and provides a clue to develop promising and nascent concept of single drug therapy for multiple diseases as well as personalized medicine.

The network organization of protein interactions in the spliceosome is reproduced by the simple rules of food-web models

The network structure of biological systems provides information on the underlying processes shaping their organization and dynamics. Here we examined the structure of the network depicting protein interactions within the spliceosome, the macromolecular complex responsible for splicing in eukaryotic cells. We show the interactions of less connected spliceosome proteins are nested subsets of the connections of the highly connected proteins. At the same time, the network has a modular structure with groups of proteins sharing similar interaction patterns. We then investigated the role of affinity and specificity in shaping the spliceosome network by adapting a probabilistic model originally designed to reproduce food webs. This food-web model was as successful in reproducing the structure of protein interactions as it is in reproducing interactions among species. The good performance of the model suggests affinity and specificity, partially determined by protein size and the timing of association to the complex, may be determining network structure. Moreover, because network models allow building ensembles of realistic networks while encompassing uncertainty they can be useful to examine the dynamics and vulnerability of intracelullar processes. Unraveling the mechanisms organizing the spliceosome interactions is important to characterize the role of individual proteins on splicing catalysis and regulation.

Ground-state phase-space structures of two-dimensional ±J spin glasses: A network approach

We illustrate a complex-network approach to study the phase spaces of spin glasses. By mapping the whole ground-state phase spaces of two-dimensional Edwards-Anderson bimodal (±J) spin glasses exactly into networks for analysis, we discovered various phase-space properties. The Gaussian connectivity distribution of the phase-space networks demonstrates that both the number of free spins and the visiting frequency of all microstates follow the Gaussian distribution. The spectra of phase-space networks are Gaussian, which is proven to be exact when the system is infinitely large. The phase-space networks exhibit community structures. By coarse graining to the community level, we constructed a network representing the entropy landscape of the ground state and discovered its scale-free property. The phase-space networks exhibit fractal structures, as a result of the rugged entropy landscape. Moreover, we show that the connectivity distribution, community structures, and fractal structures change drastically at the ferromagnetic-to-glass phase transition. These quantitative measurements of the ground states provide new insight into the study of spin glasses. The phase-space networks of spin glasses share a number of common features with those of lattice gases and geometrically frustrated spin systems and form a new class of complex networks with unique topology.

Origin and implications of zero degeneracy in networks spectra

The spectra of many real world networks exhibit properties which are different from those of random networks generated using various models. One such property is the existence of a very high degeneracy at the zero eigenvalue. In this work, we provide all the possible reasons behind the occurrence of the zero degeneracy in the network spectra, namely, the complete and partial duplications, as well as their implications. The power-law degree sequence and the preferential attachment are the properties which enhances the occurrence of such duplications and hence leading to the zero degeneracy. A comparison of the zero degeneracy in protein-protein interaction networks of six different species and in their corresponding model networks indicates importance of the degree sequences and the power-law exponent for the occurrence of zero degeneracy.

Assortative and disassortative mixing investigated using the spectra of graphs

We investigate the impact of degree-degree correlations on the spectra of networks. Even though density distributions exhibit drastic changes depending on the (dis)assortative mixing and the network architecture, the short-range correlations in eigenvalues exhibit universal random matrix theory predictions. The long-range correlations turn out to be a measure of randomness in (dis)assortative networks. The analysis further provides insight into the origin of high degeneracy at the zero eigenvalue displayed by a majority of biological networks.

Randomness and preserved patterns in cancer network

Breast cancer has been reported to account for the maximum cases among all female cancers till date. In order to gain a deeper insight into the complexities of the disease, we analyze the breast cancer network and its normal counterpart at the proteomic level. While the short range correlations in the eigenvalues exhibiting universality provide an evidence towards the importance of random connections in the underlying networks, the long range correlations along with the localization properties reveal insightful structural patterns involving functionally important proteins. The analysis provides a benchmark for designing drugs which can target a subgraph instead of individual proteins.

Uncovering randomness and success in society

An understanding of how individuals shape and impact the evolution of society is vastly limited due to the unavailability of large-scale reliable datasets that can simultaneously capture information regarding individual movements and social interactions. We believe that the popular Indian film industry, “Bollywood”, can provide a social network apt for such a study. Bollywood provides massive amounts of real, unbiased data that spans more than 100 years, and hence this network has been used as a model for the present paper. The nodes which maintain a moderate degree or widely cooperate with the other nodes of the network tend to be more fit (measured as the success of the node in the industry) in comparison to the other nodes. The analysis carried forth in the current work, using a conjoined framework of complex network theory and random matrix theory, aims to quantify the elements that determine the fitness of an individual node and the factors that contribute to the robustness of a network. The authors of this paper believe that the method of study used in the current paper can be extended to study various other industries and organizations.

Universality in network dynamics

Despite significant advances in characterizing the structural properties of complex networks, a mathematical framework that uncovers the universal properties of the interplay between the topology and the dynamics of complex systems continues to elude us. Here we develop a self-consistent theory of dynamical perturbations in complex systems, allowing us to systematically separate the contribution of the network topology and dynamics. The formalism covers a broad range of steady-state dynamical processes and offers testable predictions regarding the system's response to perturbations and the development of correlations. It predicts several distinct universality classes whose characteristics can be derived directly from the continuum equation governing the system's dynamics and which are validated on several canonical network-based dynamical systems, from biochemical dynamics to epidemic spreading. Finally, we collect experimental data pertaining to social and biological systems, demonstrating that we can accurately uncover their universality class even in the absence of an appropriate continuum theory that governs the system's dynamics.

Algorithms for systematic identification of small subgraphs

The ability to analyze large biological networks proves to be a computationally expensive task, but the information one can gain is worth the cost and effort. In cancer research for example, one is able to derive knowledge about putative drug targets by revealing the strengths and weaknesses inherent in a protein-protein interaction (PPI) network. Further, network analyses can be used to optimize high-throughput genetic and proteomic experiments. In addition, the study of biological networks is now an active part of molecular biology. In this chapter, we review techniques for studying biological networks in general but with a focus on PPI networks, including an example of a bacterial PPI network. After a brief introduction, we concentrate on methods based on the analysis of subnetworks, namely, graph motifs and graphlets.

Spectral properties of directed random networks with modular structure

We study spectra of directed networks with inhibitory and excitatory couplings. We investigate in particular eigenvector localization properties of various model networks for different values of correlation among their entries. Spectra of random networks with completely uncorrelated entries show a circular distribution with delocalized eigenvectors, whereas networks with correlated entries have localized eigenvectors. In order to understand the origin of localization we track the spectra as a function of connection probability and directionality. As connections are made directed, eigenstates start occurring in complex-conjugate pairs and the eigenvalue distribution combined with the localization measure shows a rich pattern. Moreover, for a very well distinguished community structure, the whole spectrum is localized except few eigenstates at the boundary of the circular distribution. As the network deviates from the community structure there is a sudden change in the localization property for a very small value of deformation from the perfect community structure. We search for this effect for the whole range of correlation strengths and for different community configurations. Furthermore, we investigate spectral properties of a metabolic network of zebrafish and compare them with those of the model networks.

Random matrix analysis of localization properties of gene coexpression network

We analyze gene coexpression network under the random matrix theory framework. The nearest-neighbor spacing distribution of the adjacency matrix of this network follows Gaussian orthogonal statistics of random matrix theory (RMT). Spectral rigidity test follows random matrix prediction for a certain range and deviates afterwards. Eigenvector analysis of the network using inverse participation ratio suggests that the statistics of bulk of the eigenvalues of network is consistent with those of the real symmetric random matrix, whereas few eigenvalues are localized. Based on these IPR calculations, we can divide eigenvalues in three sets: (a) The nondegenerate part that follows RMT. (b) The nondegenerate part, at both ends and at intermediate eigenvalues, which deviates from RMT and expected to contain information about important nodes in the network. (c) The degenerate part with zero eigenvalue, which fluctuates around RMT-predicted value. We identify nodes corresponding to the dominant modes of the corresponding eigenvectors and analyze their structural properties.

Phase-space networks of the six-vertex model under different boundary conditions

The six-vertex model is mapped to three-dimensional sphere stacks and different boundary conditions corresponding to different containers. The shape of the container provides a qualitative visualization of the boundary effect. Based on the sphere-stacking picture, we map the phase spaces of the six-vertex models to discrete networks. A node in the network represents a state of the system, and an edge between two nodes represents a zero-energy spin flip, which corresponds to adding or removing a sphere. The network analysis shows that the phase spaces of systems with different boundary conditions share some common features. We derived a few formulas for the number and the sizes of the disconnected phase-space subnetworks under the periodic boundary conditions. The sphere stacking provides new challenges in combinatorics and may cast light on some two-dimensional models.

Phase-space networks of geometrically frustrated systems

We illustrate a network approach to the phase-space study by using two geometrical frustration models: antiferromagnet on triangular lattice and square ice. Their highly degenerated ground states are mapped as discrete networks such that the quantitative network analysis can be applied to phase-space studies. The resulting phase spaces share some comon features and establish a class of complex networks with unique Gaussian spectral densities. Although phase-space networks are heterogeneously connected, the systems are still ergodic due to the random Poisson processes. This network approach can be generalized to phase spaces of some other complex systems.

Spectral analysis of deformed random networks

We study spectral behavior of sparsely connected random networks under the random matrix framework. Subnetworks without any connection among them form a network having perfect community structure. As connections among the subnetworks are introduced, the spacing distribution shows a transition from the Poisson statistics to the Gaussian orthogonal ensemble statistics of random matrix theory. The eigenvalue density distribution shows a transition to the Wigner's semicircular behavior for a completely deformed network. The range for which spectral rigidity, measured by the Dyson-Mehta Delta3 statistics, follows the Gaussian orthogonal ensemble statistics depends upon the deformation of the network from the perfect community structure. The spacing distribution is particularly useful to track very slight deformations of the network from a perfect community structure, whereas the density distribution and the Delta3 statistics remain identical to the undeformed network. On the other hand the Delta3 statistics is useful for the larger deformation strengths. Finally, we analyze the spectrum of a protein-protein interaction network for Helicobacter, and compare the spectral behavior with those of the model networks.

Approximating the largest eigenvalue of network adjacency matrices

The largest eigenvalue of the adjacency matrix of a network plays an important role in several network processes (e.g., synchronization of oscillators, percolation on directed networks, and linear stability of equilibria of network coupled systems). In this paper we develop approximations to the largest eigenvalue of adjacency matrices and discuss the relationships between these approximations. Numerical experiments on simulated networks are used to test our results.

Random matrix analysis of complex networks

We study complex networks under random matrix theory (RMT) framework. Using nearest-neighbor and next-nearest-neighbor spacing distributions we analyze the eigenvalues of the adjacency matrix of various model networks, namely, random, scale-free, and small-world networks. These distributions follow the Gaussian orthogonal ensemble statistic of RMT. To probe long-range correlations in the eigenvalues we study spectral rigidity via the Delta_{3} statistic of RMT as well. It follows RMT prediction of linear behavior in semilogarithmic scale with the slope being approximately 1pi;{2} . Random and scale-free networks follow RMT prediction for very large scale. A small-world network follows it for sufficiently large scale, but much less than the random and scale-free networks.

Universality in complex networks: random matrix analysis

We apply random matrix theory to complex networks. We show that nearest neighbor spacing distribution of the eigenvalues of the adjacency matrices of various model networks, namely scale-free, small-world, and random networks follow universal Gaussian orthogonal ensemble statistics of random matrix theory. Second, we show an analogy between the onset of small-world behavior, quantified by the structural properties of networks, and the transition from Poisson to Gaussian orthogonal ensemble statistics, quantified by Brody parameter characterizing a spectral property. We also present our analysis for a protein-protein interaction network in budding yeast.

Eigenfunction statistics of complex systems: a common mathematical formulation

We derive a common mathematical formulation for the eigenfunction statistics of Hermitian operators, represented by a multiparametric probability density. The system information in the formulation enters through two parameters only, namely, system size and the complexity parameter, a function of all system parameters including size. The behavior is contrary to the eigenvalue statistics which is sensitive to the complexity parameter only and shows a single parametric scaling. The existence of a mathematical formulation of both eigenfunctions and eigenvalues common to a wide range of complex systems indicates the possibility of a similar formulation for many physical properties. This also suggests the possibility to classify them in various universality classes defined by the complexity parameter.

Functional integration of ecological networks through pathway proliferation

Large-scale structural patterns commonly occur in network models of complex systems including a skewed node degree distribution and small-world topology. These patterns suggest common organizational constraints and similar functional consequences. Here, we investigate a structural pattern termed pathway proliferation. Previous research enumerating pathways that link species determined that as pathway length increases, the number of pathways tends to increase without bound. We hypothesize that this pathway proliferation influences the flow of energy, matter, and information in ecosystems. In this paper, we clarify the pathway proliferation concept, introduce a measure of the node-node proliferation rate, describe factors influencing the rate, and characterize it in 17 large empirical food-webs. During this investigation, we uncovered a modular organization within these systems. Over half of the food-webs were composed of one or more subgroups that were strongly connected internally, but weakly connected to the rest of the system. Further, these modules had distinct proliferation rates. We conclude that pathway proliferation in ecological networks reveals subgroups of species that will be functionally integrated through cyclic indirect effects.

Synchronizabilities of networks: a new index

The random matrix theory is used to bridge the network structures and the dynamical processes defined on them. We propose a possible dynamical mechanism for the enhancement effect of network structures on synchronization processes, based upon which a dynamic-based index of the synchronizability is introduced in the present paper.