Your new experience awaits. Try the new design now and help us make it even better

METHODS article

Front. Genet., 07 January 2026

Sec. Computational Genomics

Volume 16 - 2025 | https://doi.org/10.3389/fgene.2025.1613475

This article is part of the Research TopicFrom Graphs to Genes: Harnessing the Power of Machine Learning and BioengineeringView all articles

Resilience and evolutionary insights in PPI networks: comparative analysis of node resilience and centrality measures

  • School of Artificial Intelligence and Information Engineering, Yichun University, Yichun, China

Introduction: Protein-protein interaction (PPI) networks serve as the central framework for deciphering the modular structure of cellular functions and signal transduction mechanisms. While established network topological Measures (such as degree centrality, betweenness centrality, and closeness centrality) can statically characterize nodal connectivity density or pathway intermediation capacity, they fail to dynamically capture cascade following node failure.

Method: This study employs systems biology approaches to quantitatively analyze network resilience based on bacterial PPI network data obtained from the Stanford Network Analysis Platform (SNAP). First, a progressive node removal strategy was implemented to simulate cascading failure propagation and evaluate system-level resilience degradation dynamics. Subsequently, single-node knockout experiments were systematically conducted to quantify local topological disruption effects, with network fragmentation metrics (e.g., giant component size decay rate) being integrated to establish the Node Resilience (NR) index. To validate the biological relevance of NR, we developed a multidimensional analytical framework that performs cross-correlation analysis between NR and classical centrality measures [Degree centrality (DC), Betweenness centrality (BC), Closeness centrality (CC), Eigenvector centrality (EC)], enabling systematic revelation of consensus vital nodes identified by both approaches, and unique sensitive nodes detectable only through resilience-oriented perturbation analysis.

Results and Discussion: Our systematic node removal simulations revealed biphasic resilience degradation across bacterial PPI networks: progressive node failure induced gradual resilience decay whereas exceeding a critical threshold for each network triggered accelerated collapse. This phase transition aligns with evolutionary design principles - modular architectures buffer localized perturbations through functional redundancy, but inter-modular bridge depletion beyond criticality propagates cascading failures via weakly coupled connections. Notably, NR exhibited a strong negative correlation with BC, contrasting with weak associations for DC, CC, and EC. This dichotomy arises because BC quantifies cross-modular information brokerage - high-BC nodes act as structural keystones whose removal disconnects functional modules, drastically reducing global entropy. Conversely, for DC, CC, and EC primarily reflect local connectivity patterns with limited cascade propagation potential.

1 Introduction

Biological networks, as abstract representations of molecular interactions in living systems, encompass diverse subtypes including metabolic pathways, gene regulatory networks, and protein-protein interaction (PPI) networks. Research on general biological networks has laid the groundwork for understanding molecular interaction systems. As early as the 1990s, the analysis of metabolic pathways originated in the 1990s and established essential methods for network structure and function analysis (Fell and Cornish-Bowden, 1997). Subsequently, Barabási and Oltvai formalized the network biology theory, revealing that biological networks (including PPI networks) exhibit scale-free, hierarchical modular, and ultra-small-world topological principles that underpin their functional resilience and evolutionary dynamics (Barabasi and Oltvai, 2004).

Protein-protein interaction (PPI) networks serve as foundational frameworks in systems biology, encoding the intricate functional relationships that sustain cellular processes (Ciliberti et al., 2007). PPI networks enable systematic interrogation of protein complexes, signaling cascades, and functional modules, elucidating regulatory architectures governing cellular processes (Gavin et al., 2002; Von Mering et al., 2002; Ideker and Krogan, 2012). Their analysis has thus become indispensable for deciphering biological complexity in proteomics and systems biology (Wodak et al., 2013; Szklarczyk et al., 2019).

External perturbations, sucn as genetic mutations, pathogenic invasions, or pharmacological agents, can destabilize PPI network topology and functionality (Albert et al., 2000; Rodrigues et al., 2021). Node removal (e.g., protein dysfunction) disrupts associated edges, isolating the node from its interactome. Progressive node loss exacerbates topological disintegration, compromising critical functional modules and system-wide stability (Jeong et al., 2000; Motter, 2004). Quantifying these dynamics is pivotal for understanding disease pathogenesis, therapeutic efficacy, and evolutionary adaptation (Goh et al., 2007).

Zitnik et al. (2019) conducted evolutionary analyses of 1,840 species and found that phylogenetically advanced organisms exhibit enhanced PPI network resilience, suggesting that natural selection favors robust topologies to buffer environmental stresses. However, the precise mechanisms governing resilience decay under sustained intensifying perturbations, particularly in microbial systems, remain poorly understood. Addressing this, we present a computational framework to model progressive node removal in bacterial PPI networks, integrating entropy-based resilience measures and iterative simulations to dissect fragility dynamics. Our study identifies biphasic resilience degradation across bacterial interactomes: gradual decay under low-intensity perturbations transitions abruptly to systemic collapse when node failure surpasses a critical threshold. This phase transition, characterized by steepened decay slopes, underscores the evolutionary trade-off between metabolic efficiency and redundancy allocation.

Furthermore, we introduce Node Resilience (NR), a novel metric quantifying individual proteins’ contributions to global network stability. Comparative analyses with traditional centrality measures reveal that NR inversely correlates with betweenness centrality (Joyce et al., 2010), highlighting bridge nodes as structural keystones whose removal triggers modular disarticulation.

To illustrate the significance of resilience-critical nodes, we conducted a case study on Listeria monocytogenes, demonstrating that these nodes frequently orchestrate virulence regulation and metabolic adaptation. These findings advance our understanding of PPI network evolvability, offering insights into therapeutic targeting and the evolutionary drivers of interactome robustness.

2 Methods

2.1 Network resilience

PPI networks represent functional associations between proteins within a biological system (Barabasi and Oltvai, 2004). Nodes denote individual proteins, while edges signify pairwise interactions. These networks are typically undirected and unweighted. Notably, the PPI networks in this study are derived from the largest connected component (LCC). Zitnik et al. formalized resilience as the network’s ability to sustain connectivity under iterative perturbations in 2019.

The network is progressively fragmented by iteratively isolating nodes, that is, removing all their links. Links represent protein-protein interactions, and their removal result in isolated proteins and smaller non-interacting components. The improved Shannon diversity index is employed to quantify the diversity of components in a network after nodes removal. The process involves removing nodes from the originally connected network G, which subsequently fragments into k disconnected components. The index is calculated as Equation 1:

HmshGf=1logNi=1kpilogpi(1)

where N is the total number of nodes in the original network G, k is the total number of components in the fragmented network Gf, pi = |ci|/N is the proportion of nodes in component ci relative to N. Therefore, pi represents the probability that a randomly selected node belongs to the connected component ci.

The modified Shannon diversity Hmsh was employed to quantify the component diversity of the protein interactome under increasing node failure levels f. At f = 0, the network is intact and Hmsh = 0. As f increases and more proteins are removed, the network disintegrates and Hmsh rises. When all proteins are isolated (f = 1), Hmsh peaks at 1.

The resilience metric integrates entropy across perturbation intensities. This integration quantifies the network’s ability to maintain stability under different disturbance levels, and its calculation is defined by Equation 2:

RG=101HmshGfdf(2)

where Hmsh(Gf) denotes entropy at failure rate f. A high R(G) value reflects delayed fragmentation during progressive node removal, indicating networks that retain a giant connected component despite cumulative perturbations.

Broad-scale evolutionary studies across diverse taxa have demonstrated that organisms positioned higher on the phylogenetic scale develop reinforced interactome robustness (Zitnik et al., 2019). Experimental validation from microbial evolution studies, including the Escherichia coli Long-Term Evolution Experiment, demonstrates that evolved strains exhibit PPI networks with significantly higher resilience than randomized counterparts, demonstrating selective maintenance of network stability under prolonged environmental pressures (Maddamsetti, 2021). These findings posit resilience as both a product of functional constraints and an evolutionary driver, enabling survival in fluctuating ecological niches.

2.2 Node resilience

In PPI networks, spatially heterogeneous nodes exhibit marked functional heterogeneity, while the systematic identification of topological hubs under dynamic environmental perturbations, remains a central unresolved challenge. Although traditional centrality measures [e.g., degree centrality, closeness centrality (Freeman, 1978), eigenvector centrality (Bonacich, 1987), and betweenness centrality (Freeman, 1977)] can effectively predict essential proteins, they exhibit systematic limitations in quantifying the cascade disruption effects induced by node removal on network connectivity.

To quantify the contribution of individual nodes to network resilience, we introduce the concept of Node Resilience (NR). This metric evaluates the criticality of a node by measuring the change in network resilience after its removal. The calculation process is as follows: Consider a PPI network G, represented as G=(V,E), where V is the set of N nodes (vi∈V, i = 1,2,…,N), and E is the set of edges representing interactions between nodes.

For each node vi in G, its removal results in a new network G'=(V|{vi},E′), where E′ is the set of edges after removing all interactions involving vi. The resilience of the residual network G′ is then calculated, yielding the resilience value for node vi, denoted as NRi.

The complete set of node resilience values {NR1,NR2,…,NRn } is generated by iteratively removing each node vi from original network G and quantifying the resilience R of the resulting residual network G|{vi}, formally expressed as Equation 3:

NR1,...,NRn=RG|viviV(3)

3 Results

3.1 Progressive degradation of network resilience under node removal

Systematic interrogation of how PPI networks degrade under escalating perturbations (a critical determinant of cellular robustness) remains limited. Here, we dissect resilience dynamics in bacterial interactomes through iterative node failure simulations, probing both evolutionary design rules and catastrophic fragility thresholds.

We developed a computational framework to model progressive network degradation via stochastic node removal, wherein node removal rates (f) increased incrementally from 0% to 100% in 1% steps, simulating evolutionary stressors such as pathogenic gene silencing. Network resilience was quantified using the entropy-integrated, which captures a network’s capacity to buffer perturbations. To ensure statistical reliability, each failure rate underwent 100 independent iterative simulations.

We analyzed six bacterial protein-protein interaction (PPI) networks as described in Table 1, which were retrieved from the SNAP database (Figure 1A). The six bacteria cover multiple key evolutionary branches of the bacterial domain, with distinct biological functions and survival strategies, enabling comprehensive verification of network resilience trend across different functional bacterial types.

Table 1
www.frontiersin.org

Table 1. Network and computational characteristics of bacterial PPI networks.

Figure 1
Scatter plots show resilience versus node failure ratio for six bacterial species with breakpoints indicated: Rhodobacter sphaeroides (51.62%), Bacillus subtilis (63.76%), Treponema pallidum (59.53%), Burkholderia pseudomallei (63.32%), Synechococcus (49.98%), and Stenotrophomonas maltophilia (53.73%). A boxplot summarizes breakpoint positions with a mean of 52.26% and a 95% confidence interval.

Figure 1. Resilience degradation dynamics in bacterial PPI networks under progressive node failure rate. (A) The six scatter plots on the left illustrate the resilience variation trends of six distinct bacterial species during the progressive node removal process (with 1% removal per step). In each plot, the red fitted line intuitively depicts the change pattern of resilience with the Node Failure Ratio. A breakpoint, marked by a green dashed line, exists in all plots. (B) The box plot on the right summarizes the distribution range of resilience breakpoints across all bacteria. The results show that the breakpoints are concentrated in the Node Failure Ratio interval of 40%–65%.

The resilience trajectories of six representative bacterial PPI networks are visualized through progressive node removal simulations. The results reveal the network resilience decreases relatively gently to the left of the breakpoint, while the decline rate accelerates significantly to the right, reflecting the “resilience threshold” characteristic of the network structure. Expanded analysis of bacterial PPI networks with node counts exceeding 100 from the SNAP Database demonstrated measurable criticality in resilience dynamics (Figure 1B). The metabolic networks of these bacteria can maintain resilience through structural redundancy before approximately 40%–65% of nodes are removed, and then lose stability rapidly beyond this threshold.

This study reveals that PPI networks exhibit biphasic dynamics: they buffer low-intensity perturbations, until damage exceeds a critical threshold, where cascading failures drive abrupt robustness collapse (Cohen et al., 2000; Motter, 2004). This phenomenon reflects the robustness-efficiency trade-off strategy evolved in biological systems. The modular architecture achieves effective buffering against stochastic noise under metabolic cost constraints through spatial compartmentalization that restricts local disturbance propagation and redundant functional backups (sparing critical components). When incurred damage surpasses the critical threshold, weakly coupled inter-modular connections may transform into propagation channels for cascade failures, ultimately triggering systemic collapse (Callaway et al., 2000; Gao et al., 2016).

3.2 Node resilience versus centrality metrics

We calculated the six bacterial species mentioned above and found that most nodes in the network have a limited impact on overall resilience when removed. However, the removal of a small subset of key nodes, causes significant declines in network resilience and severe fragmentation of the network structure (Figure 2). From Figure 2, we can distinctly observe that across all six bacterial PPI networks, the node resilience (NR) values exhibit a remarkable pattern: the vast majority of nodes cluster showing only marginal fluctuations. In contrast, a small subset of nodes deviates significantly from this cluster, displaying notably higher or lower resilience values. This distribution pattern implies that the network’s overall resilience is robust to the removal of most nodes—their elimination barely perturbs the network’s ability to maintain structural and functional integrity. Conversely, the removal of those few outlying nodes, which possess exceptionally high or low node resilience. This observation aligns with the modular and hierarchical nature of PPI networks, where most nodes participate in redundant or auxiliary interactions, while a select few mediate critical pathways or structural connections that underpin the network’s resilience.

Figure 2
Scatter plots of resilience versus node number for six bacteria: a) Rhodobacter sphaeroides, b) Bacillus subtilis, c) Treponema pallidum, d) Burkholderia pseudomallei, e) Synechococcus, f) Stenotrophomonas maltophilia. Each plot shows mean, 95% data range, and outlier nodes highlighted.

Figure 2. Node resilience of the six bacterial PPI networks. In each subgraph, the horizontal axis represents the order of nodes’ appearance in the PPI network, while the vertical axis denotes node resilience (quantifying the network’s resilience after the removal of individual nodes). As shown in the figure, most nodes have resilience values falling within the 95% confidence interval, with only a small fraction lying outside this range-these outlier nodes are highlighted in red. This indicates that the resilience values of most nodes show little variation, while the red-highlighted subset exhibits significant deviations.

Comprehensive correlation analyses demonstrated that NR captures distinct aspects of topological vulnerability compared to traditional centrality paradigms. A pronounced negative correlation emerged between NR and BC, indicating that nodes functioning as inter-modular bridges—those with high BC values—exhibit disproportionately low resilience. In contrast, NR displayed substantially weaker correlations with DC, CC, and EC.

To assess the generalizability of these findings, we expanded our analysis to all 84 bacterial species from the SNAP database. The violin plot in Figure 3E summarizes the distribution of correlation coefficients across this extended dataset. BC maintained a robust negative correlation with NR, while DC, CC, and EC showed consistently weak anti-correlations.

Figure 3
Scatter plots labeled A to D show relationships between Node Resilience and various centrality measures, each plot displaying a trend line and data points with different colors indicating specific data sets. Statistical information, including correlation coefficients and p-values, is provided in legends. A violin plot labeled E compares distributions of different resilience measures across four categories represented in various colors.

Figure 3. Comparison of NR with Four Centrality measures. (A–D) Scatter plots representing the relationship between node resilience and various centrality measures, including DC,CC, BC, and EC. (E) Violin plot displaying the distribution of correlation coefficients between node resilience and the centrality measures across all 84 bacterial species from the SNAP database.

The resilience-connectivity duality in PPI networks operates through distinct mechanistic pathways. BC quantifies a node’s role as a bridge in shortest paths—a topological property intrinsically linked to network modular architecture. These bridge nodes facilitate communication between network modules, and their removal induces catastrophic modular disarticulation, effectively fragmenting the network into isolated components. High-BC nodes act as structural keystones whose removal disrupts inter-module communication, effectively decoupling functional units like metabolic pathways or signal transduction cascades (Guimera and Nunes Amaral, 2005; Newman, 2006). This explains their disproportionate impact on network resilience, as measured by NR. Conversely, peripheral nodes with low BC values primarily contribute to intra-modular connectivity. Their removal induces localized rather than whole effects, resulting in minimal NR deviations. While DC, CC, and EC effectively map nodes’ local connectivity influence, proximity to network cores, or leadership within densely connected clusters (Rayfield et al., 2011), they fail to capture the cross-modular integration capacity epitomized by BC. This analytical gap positions NR as a complementary metric that specifically highlights nodes orchestrating global network coherence—those whose removal triggers cascading failure through modular decoupling rather than localized connectivity loss.

3.3 A comparative study of NR and BC in the Listeria monocytogenes

Our topological analysis of the PPI network of Listeria monocytogenes focused on its largest connected compenent (187 nodes, 341 edges) to elucidate core architectural features. Two complementary topological ranking strategies were implemented: NR identified Key Component Nodes (KCNs) — the 10% of nodes whose removal caused maximal resilience loss—while BC pinpointed Key Bridge Nodes (KBNs) — the 10% most frequent mediators of shortest paths. The Venn diagram (Figure 4) illustrates that there is an 89.4% overlap (16 out of 19 nodes) between KCNs and KBNs, highlighting their dual critical roles in structural robustness and information routing. The details of the Key Component Nodes are provided in Table 2.

Figure 4
Venn diagram with two overlapping circles. The left circle represents Key Component Nodes (KCNs) with three unique elements. The right circle represents Key Bridge Nodes (KBNs) also with three unique elements. The overlapping area shows sixteen elements common to both sets.

Figure 4. Intersection of KCNs and KBNs datasets in the Listeria PPI network: the blue area represents KCNs, and the red area represents KBNs. The overlapping area highlights the dual core nodes.

Table 2
www.frontiersin.org

Table 2. Key Nodes information.

The core nodes shared by KCNs and KBNs (e.g., pykA (lmo1570),pdhB (lmo1053), gcvT (lmo1348), thyA (lmo1874), ndk (lmo1929), among others) are extensively involved in critical pathways. Specifically, PykA maintains cell membrane integrity by enhancing fatty acid synthesis (Liu et al., 2017); pdhB conserves resources by downregulating energy metabolism (Stasiewicz et al., 2011); gcvT participates in amino acid metabolic reprogramming (Mujahid et al., 2013); thyA ensures genetic stability by supporting nucleotide synthesis (Kimman et al., 2008); and ndk acts as a pleiotropic hub, simultaneously coordinating nucleotide metabolic balance and virulence expression (Yu et al., 2017). Notably, folD is exclusively present in KCNS and the bifunctional enzyme encoded by the gene is a central component of folate-mediated one-carbon metabolism and directly supports rapid bacterial growth and proliferation (Zhang et al., 2022). Together, these genes form a core network enabling bacterial stress response and enhanced pathogenicity.

4 Discussion

This study introduced a new network science approach to evaluate the resilience of PPI networks under simulated disruptions. By constructing a model with randomly removed nodes, we investigated the ability of these networks to maintain topological stability in the face of faults or attacks. Our findings revealed a biphasic resilience degradation pattern: in the pre-critical phase, the network demonstrated an evolutionary adaptation to buffer random failures; however, once node failure reaches the critical threshold, the network exhibits accelerated collapse and inherent vulnerability under continued perturbations. The identification of this critical failure threshold aligns with theoretical models positing that biological networks evolve to maximize robustness against common perturbations while tolerating rare, catastrophic events. This dual-phase behavior underscores the delicate balance between resilience and fragility in biological systems, reflecting their evolutionary optimization for both stability and adaptability.

We proposed a novel metric node resilience (NR), which calculates the resilience value of the network after the removal of each node. This metric quantifies the impact of individual nodes on overall network stability and provides insights into the structural and functional dynamics of PPI networks. The results demonstrate that most nodes in the network have a limited impact on overall resilience when removed. However, the removal of a small subset of key nodes, causes significant declines in network resilience and severe fragmentation of the network structure. These critical nodes play essential roles in connecting different components of the network and maintaining its global connectivity. Their removal disrupts network cohesion, underscoring their importance for preserving biological system functions.

NR is a valuable complement to traditional centrality measures, with distinct correlations tied to network connectivity and modular isolation: a strong negative correlation with BC, as high-BC inter-modular bridges drive systemic fragmentation when removed; a weak negative correlation with DC, since high-DC local hubs only affect intra-modular connectivity and cause minimal isolation; a weak negative correlation with CC, as proximity to cores does not equate to maintaining inter-module links; and a negligible negative correlation with EC, as EC overlooks cross-modular integration—NR’s key focus for resisting modular isolation. While measures like betweenness centrality identify key bridging nodes, node resilience quantifies how the removal of specific nodes affects the entire network’s structural integrity. This dual perspective enriches our understanding of the roles of individual nodes in PPI networks and their contributions to maintaining stability and connectivity.

Despite the promising insights offered by this approach, there are limitations to this study. First, the dataset used is restricted to bacterial organisms, and future research should validate these findings across diverse species and experimental conditions. Second, as network scale increases, computational time grows significantly, which is a scalability limitation of our current framework. Third, for large-scale PPI networks, focusing solely on the impact of individual node failure on network stability has obvious limitations. Our subsequent research will employ protein grouping to systematically categorize nodes, followed by exploring the dynamic changes in network resilience with functional groups as the research unit, which not only addresses the inadequacy of single-node analysis but also provides a more biologically meaningful perspective on how functional modules collectively influence network stability while mitigating the computational burden associated with large-scale networks.

In summary, the use of node resilience as a new metric for assessing the contribution of individual nodes to network stability represents a significant advancement in understanding and analyzing complex biological systems. This approach complements existing network analysis techniques and opens new avenues for studying PPI networks, with potential applications in systems biology, disease modeling, and therapeutic target identification.

Data availability statement

Publicly available datasets were analyzed in this study. This data can be found here: http://snap.stanford.edu/tree-of-life/.

Author contributions

JZ: Conceptualization, Methodology, Software, Formal analysis, Visualization, Writing – original draft, Writing – review and editing.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This work was supported by the National Natural Science Foundation of China (Grant No. 62161050), Science and Technology Research Project of Jiangxi Provincial Department of Education (Grant No. GJJ201605), and Jiangxi Provincial Natural Science Foundation (Grant No. 20252BAC240210).

Conflict of interest

The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Albert, R., Jeong, H., and Barabási, A. L. (2000). Error and attack tolerance of complex networks. Nature 406 (6794), 378–382. doi:10.1038/35019019

PubMed Abstract | CrossRef Full Text | Google Scholar

Barabasi, A.-L., and Oltvai, Z. N. (2004). Network biology: understanding the cell's functional organization. Nat. Reviews. Genet. 5 (2), 101–113. doi:10.1038/nrg1272

PubMed Abstract | CrossRef Full Text | Google Scholar

Bonacich, P. (1987). Power and centrality: a family of measures. Am. Journal Sociology 92 (5), 1170–1182. doi:10.1086/228631

CrossRef Full Text | Google Scholar

Callaway, D. S., Newman, M. E., Strogatz, S. H., and Watts, D. J. (2000). Network robustness and fragility: percolation on random graphs. Phys. Review Letters 85 (25), 5468–5471. doi:10.1103/PhysRevLett.85.5468

PubMed Abstract | CrossRef Full Text | Google Scholar

Ciliberti, S., Martin, O. C., and Wagner, A. (2007). Robustness can evolve gradually in complex regulatory gene networks with varying topology. Plos Comput. Biol. 3 (2), 164–173. doi:10.1371/journal.pcbi.0030015

PubMed Abstract | CrossRef Full Text | Google Scholar

Cohen, E., Ben, A., and Havlin, (2000). Resilience of the internet to random breakdowns. Phys. Review Letters 85 (21), 4626–4628. doi:10.1103/PhysRevLett.85.4626

CrossRef Full Text | Google Scholar

Fell, D., and Cornish-Bowden, A. (1997). Understanding the control of metabolism. London: Portland press.

Google Scholar

Freeman, L. C. (1977). A set of measures of centrality based on betweenness. Sociometry 40, 35–41. doi:10.2307/3033543

CrossRef Full Text | Google Scholar

Freeman, L. C. (1978). Centrality in social networks conceptual clarification. Soc. Networks 1 (3), 215–239. doi:10.1016/0378-8733(78)90021-7

CrossRef Full Text | Google Scholar

Gao, J., Barzel, B., and Barabási, A. L. (2016). Universal resilience patterns in complex networks (vol 530, pg 307, 2016). Nature 536 (7615), 238. doi:10.1038/nature18019

PubMed Abstract | CrossRef Full Text | Google Scholar

Gavin, A. C., Bosche, M., Krause, R., Grandi, P., Marzioch, M., Bauer, A., et al. (2002). Functional organization of the yeast proteome by systematic analysis of protein complexes. Nat. 415, 141–147. doi:10.1038/415141a

PubMed Abstract | CrossRef Full Text | Google Scholar

Goh, K.-I., Cusick, M. E., Valle, D., Childs, B., Vidal, M., and Barabási, A.-L. (2007). The human disease network. Proc. Natl. Acad. Sci. 104(21): 8685–8690. doi:10.1073/pnas.0701361104

PubMed Abstract | CrossRef Full Text | Google Scholar

Guimera, R., and Nunes Amaral, L. A. (2005). Functional cartography of complex metabolic networks. Nature 433 (7028), 895–900. doi:10.1038/nature03288

PubMed Abstract | CrossRef Full Text | Google Scholar

Ideker, T., and Krogan, N. J. (2012). Differential network biology. Mol. Syst. Biol. 8 (1), 565. doi:10.1038/msb.2011.99

PubMed Abstract | CrossRef Full Text | Google Scholar

Jeong, H., Tombor, B., Albert, R., Oltvai, Z. N., and Barabasi, A. L. (2000). The large-scale organization of metabolic networks. Nature 407 (6804), 651–654. doi:10.1038/35036627

PubMed Abstract | CrossRef Full Text | Google Scholar

Joyce, K. E., Laurienti, P. J., Burdette, J. H., and Hayasaka, S. (2010). A new measure of centrality for brain networks. Plos One 5 (8), e12200. doi:10.1371/journal.pone.0012200

PubMed Abstract | CrossRef Full Text | Google Scholar

Kimman, T. G., Smit, E., and Klein, M. R. (2008). Evidence-based biosafety: a review of the principles and effectiveness of microbiological containment measures. Clin. Microbiol. Rev. 21(3): 403–425. doi:10.1128/cmr.00014-08

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, X., Basu, U., Miller, P., and McMullen, L. M. (2017). Differential gene expression and filamentation of Listeria monocytogenes 08-5923 exposed to sodium lactate and sodium diacetate. Food Microbiology 63, 153–158. doi:10.1016/j.fm.2016.11.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Maddamsetti, R. (2021). Selection maintains protein interactome resilience in the long-term evolution experiment with Escherichia coli. Genome Biol. Evol. 13 (6), evab074. doi:10.1093/gbe/evab074

PubMed Abstract | CrossRef Full Text | Google Scholar

Motter, A. E. (2004). Cascade control and defense in complex networks. Phys. Review Letters 93 (9), 098701. doi:10.1103/PhysRevLett.93.098701

PubMed Abstract | CrossRef Full Text | Google Scholar

Mujahid, S., Orsi, R., Vangay, P., Boor, K., and Wiedmann, M. (2013). Refinement of the Listeria monocytogenes σB regulon through quantitative proteomic analysis. Microbiology 159 (Pt_6), 1109–1119. doi:10.1099/mic.0.066001-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Newman, M. E. J. (2006). Modularity and community structure in networks. Proc. Natl. Acad. Sci. U. S. A., 103(23): 8577–8582. doi:10.1073/pnas.0601602103

PubMed Abstract | CrossRef Full Text | Google Scholar

Rayfield, B., Fortin, M. J., and Fall, A. (2011). Connectivity for conservation: a framework to classify network measures. Ecology 92 (4), 847–858. doi:10.1890/09-2190.1

PubMed Abstract | CrossRef Full Text | Google Scholar

Rodrigues, C. H. M., Pires, D. E. V., and Ascher, D. B. (2021). mmCSM-PPI: predicting the effects of multiple point mutations on protein-protein interactions. Nucleic Acids Res. 49 (W1), W417–W424. doi:10.1093/nar/gkab273

PubMed Abstract | CrossRef Full Text | Google Scholar

Stasiewicz, M. J., Wiedmann, M., and Bergholz, T. M. (2011). The transcriptional response of Listeria monocytogenes during adaptation to growth on lactate and diacetate includes synergistic changes that increase fermentative acetoin production. Appl. Environmental Microbiology 77 (15), 5294–5306. doi:10.1128/AEM.02976-10

PubMed Abstract | CrossRef Full Text | Google Scholar

Szklarczyk, D., Gable, A. L., Lyon, D., Junge, A., Wyder, S., Huerta-Cepas, J., et al. (2019). STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 47 (D1), D607–D613. doi:10.1093/nar/gky1131

PubMed Abstract | CrossRef Full Text | Google Scholar

von Mering, C., Krause, R., Snel, B., Cornell, M., Oliver, S. G., Fields, S., et al. (2002). Comparative assessment of large-scale data sets of protein-protein interactions. Nature 417 (6887), 399–403. doi:10.1038/nature750

PubMed Abstract | CrossRef Full Text | Google Scholar

Wodak, S. J., Vlasblom, J., Turinsky, A. L., and Pu, S. Y. (2013). Protein-protein interaction networks: the puzzling riches. Curr. Opin. Struct. Biol. 23 (6), 941–953. doi:10.1016/j.sbi.2013.08.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, H., Rao, X., and Zhang, K. (2017). Nucleoside diphosphate kinase (ndk): a pleiotropic effector manipulating bacterial virulence and adaptive responses. Microbiol. Res. 205, 125–134. doi:10.1016/j.micres.2017.09.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Anaya-Sanchez, A., and Portnoy, D. A. (2022). Para-aminobenzoic acid biosynthesis is required for Listeria monocytogenes growth and pathogenesis. Infect. Immun. 90 (11), e00207–e00222. doi:10.1128/iai.00207-22

PubMed Abstract | CrossRef Full Text | Google Scholar

Zitnik, M., Sosic, R., Feldman, M. W., and Leskovec, J. (2019). Evolution of resilience in protein interactomes across the tree of life. Proc. Natl. Acad. Sci. U. S. A. 116 (10), 4426–4433. doi:10.1073/pnas.1818013116

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: biological networks, network structure, network resilience, node resilience, centrality measures

Citation: Zhang J (2026) Resilience and evolutionary insights in PPI networks: comparative analysis of node resilience and centrality measures. Front. Genet. 16:1613475. doi: 10.3389/fgene.2025.1613475

Received: 17 April 2025; Accepted: 10 December 2025;
Published: 07 January 2026.

Edited by:

Oksana Sorokina, University of Edinburgh, United Kingdom

Reviewed by:

Yuriy L. Orlov, I. M. Sechenov First Moscow State Medical University, Russia
Colin David Mclean, University of Edinburgh, United Kingdom

Copyright © 2026 Zhang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Jiarui Zhang, amlhemhhbmcyMDIwQDEyNi5jb20=

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.