The team’s approach rests on a protein side chain-based network method, a technique in which each amino acid residue is represented as a node, and edges are drawn between residues whose side chains come into close non-covalent contact. Such contacts include the hydrogen bonds, hydrophobic interactions, and electrostatic attractions that hold proteins together without forming chemical bonds. Rather than analyzing a single static structure, the researchers averaged their network construction over many snapshots drawn from molecular dynamics simulations produced by D. E. Shaw Research, an unusually long and detailed set of simulations that capture the complex in motion at atomic resolution. Averaging across these frames filters out fleeting, transient contacts and reveals the persistent connectivity patterns that define the complex’s functional architecture.

Once the network is built, the analysis proceeds at multiple hierarchical levels. At the most local scale, the researchers examined cliques—fully connected subgroups of residues in which every member interacts with every other member. Cliques at the viral–human interface act as tightly knit hubs of interaction, and their abundance and composition reveal how densely the two proteins interlock. At the global scale, the team analyzed clusters of connected residues and performed percolation analysis, a concept borrowed from statistical physics that asks how a network holds together as connections are progressively removed. In a percolation framework, the strongest substructures are identified by stripping away weaker links and observing which clusters survive.

The central finding of the study is unambiguous: the SARS-CoV-2 spike protein forms a more dominant and robust connection with the human ACE2 receptor than SARS-CoV-1 does. At the global level, this superiority shows up clearly in the percolation results. When the researchers restricted the networks to only the strongest interfacial bonds, the interface cluster of the SARS-CoV-2–ACE2 complex persisted, remaining a single connected unit, whereas the corresponding cluster in the SARS-CoV-1–ACE2 complex fragmented. In other words, the stronger non-covalent interactions between the newer virus and its receptor form a self-reinforcing web that cannot be dismantled by discarding weak contacts—a hallmark of structural robustness.

The atomic-scale analysis tells a complementary story. In the SARS-CoV-2 complex, network cliques at the interface form in a manner that reflects denser and better-organized inter-protein communication. Cliques are demanding structures: because every residue in a clique must contact every other, their emergence signals a level of mutual reinforcement that looser collections of contacts lack. The fact that these clique structures at the SARS-CoV-2 interface consistently encode stronger connectivity provides a mechanistic, residue-level explanation for the virus’s enhanced receptor engagement, one that emerges naturally from the geometry of interactions rather than from any single pinned atom or bond.

Beyond confirming what virologists had suspected from binding-affinity measurements, the network framework delivers something more practical: an objective method to map how atomic-level non-covalent interactions give rise to global connectivity between any two amino acids in the complex. Bridging these scales has long been a challenge in molecular biophysics. A single hydrogen bond at an interface may seem inconsequential in isolation, yet through the network it can contribute to pathways of communication that span the entire protein assembly. The side chain-network approach makes this propagation explicit, allowing researchers to trace exactly which local interactions feed into which large-scale structural couplings, and to rank residues by their contribution to overall interface integrity.

Using this ranking ability, the team pinpointed key functional residues in SARS-CoV-2 that play important roles in establishing the higher connectivity with ACE2. These residues, identified through their network centrality and their participation in robust interfacial clusters and cliques, constitute the structural backbone of the viral grip on human cells. The authors also examined specific individual snapshots from the molecular dynamics trajectories rather than only time-averaged networks, and this analysis highlighted prominent variations in network topology across different conformations. The spike–ACE2 complex does not sit in a single rigid state; it explores a diverse conformational landscape, and the connectivity architecture shifts meaningfully from one moment to another. Capturing these fluctuations matters because a receptor-binding interface is a dynamic entity, and transient configurations may expose vulnerabilities or reveal alternative interaction modes.

Perhaps the most striking connection to real-world virology emerges when the network map is overlaid with the evolutionary record of the pandemic. The researchers found that most mutations occurring in the SARS-CoV-2 spike protein in variants of concern and variants of interest—including the currently active JN.1 variant and its subvariants—have been observed at the interface with the ACE2 receptor. This is precisely the region that the network analysis identifies as the seat of the virus’s enhanced connectivity. The convergence is telling: as the virus evolves to evade immunity generated by prior infection or vaccination, natural selection keeps returning to the same interfacial neighborhood that network science flags as functionally dominant. Mutations there can simultaneously alter antibody recognition and tune receptor binding, which explains why this region is both the engine of immune escape and the object of intense surveillance by public health laboratories worldwide.

This convergence between evolutionary observation and network topology gives the study practical significance for therapeutic design. If a defined set of residues sustains the robust connectivity between the viral spike and ACE2, then molecules that mimic the relevant receptor surface could, in principle, compete with ACE2 for spike binding and neutralize the virus. The authors note that their analyses highlight the importance of interface interactions and provide a rationale for designing receptor-like peptides and proteins to combat immunity-escaping variants. Because such decoys target the physical grip on ACE2 rather than the antibody-recognized epitopes that mutate most rapidly under immune pressure, they may retain effectiveness even as the virus drifts genetically—a property that conventional antibody drugs increasingly struggle to maintain against JN.1-descended lineages.

More broadly, the work demonstrates the value of complex-systems thinking in molecular biology. A virus–receptor complex is not merely a collection of atoms or a static ribbon diagram; it is a network whose resilience, communication pathways, and failure points determine biological function. By providing an objective, quantitative bridge from atomic non-covalent interactions to global connectivity between any two amino acids, the framework developed by Subramanyan and colleagues can in principle be applied to any protein complex, well beyond coronaviruses. For SARS-CoV-2 specifically, it delivers a durable conceptual map of why this virus binds human cells so effectively, why evolution keeps redrawing the same interfacial territory, and where next-generation antiviral designs might most productively aim as the virus continues to diversify.

Subject of Research: Network analysis of the SARS-CoV-2 spike–human ACE2 receptor complex connectivity

Article Title: From atomic to global connectivity in the structure of the SARS-CoV2-human ACE2 receptor complex

Article References: Subramanyan, V., Anto, A., Bhattacharyya, M., Vishveshwara, S., & Vishveshwara, S. (2026). From atomic to global connectivity in the structure of the SARS-CoV2-human ACE2 receptor complex. PLOS Complex Systems, 3(4), e0000096. https://doi.org/10.1371/journal.pcsy.0000096

Image Credits: AI Generated

DOI: 10.1371/journal.pcsy.0000096

Keywords: SARS-CoV-2, ACE2 receptor, spike protein, SARS-CoV-1, molecular dynamics, protein structure networks, percolation analysis, network cliques, variants of concern, JN.1 variant, immune escape, peptide design

The scale of the underlying dataset is what sets this analysis apart. The investigators assembled a cohort of every individual registered with the Lombardy Regional Health Service as of January 2020, aggregating monthly data by sex, age, country of origin, municipality or ZIP code area, and pre-existing comorbidities. The region’s eight local public health authorities, which administer surveillance and prevention across distinct catchment areas, contributed records that allowed the researchers to track three separate outcomes: confirmed COVID-19 incidence, hospital admissions, and mortality. Crucially, the team also had access to SARS-CoV-2 lineage data from the regional genomic repository, which enabled them to divide the study period into four sub-periods defined by the predominant viral variants circulating at the time. This variant-based segmentation meant the researchers could ask not only how the epidemic changed over time, but how the relationship between risk factors and outcomes itself evolved as the virus mutated.

The demographic profile of the population under study helps explain why Lombardy was hit so hard. Of the more than ten million residents included in the cohort, 45 percent were aged 50 or older, and 22 percent carried a history of cardiac disease. These figures describe a region with a substantial reservoir of older adults and cardiovascular vulnerability, precisely the characteristics epidemiologists associate with severe COVID-19 outcomes. The study’s multivariable regression models, run separately for each of the four variant-defined sub-periods, confirmed that advancing age remained the dominant risk factor for both infection and death throughout the entire pandemic. Neurologic and renal comorbidities also emerged as consistent predictors of mortality, alongside cardiac conditions, underscoring that the burden of severe disease fell disproportionately on people with chronic organ dysfunction.

One of the most striking findings concerns the geography of the early epidemic. Contrary to the intuitive expectation that a respiratory virus would ignite first in the largest urban centers, the outbreak in Lombardy was initially concentrated in a handful of Local Health Districts, none of which contained large cities. This observation aligns with the now well-understood story of the virus’s silent spread through smaller provincial towns in the Po Valley before the first cases were even recognized. The geographic heterogeneity persisted throughout the study period: incidence, hospitalization rates, and mortality varied substantially across districts, and the researchers found that geographical location itself remained an independent predictor of outcomes even after accounting for age, sex, comorbidities, and other individual-level characteristics.

That persistent geographic signal carries an important interpretive message. Because the models adjusted for measurable individual risk factors, the residual effect of location hints at the influence of local epidemic dynamics, such as the timing and intensity of transmission waves in each community, and of health system factors, including differences in testing capacity, care pathways, and hospital load between districts. The authors suggest that areas impacted early in the pandemic experienced distinct trajectories compared with those hit later, and that these local dynamics shaped outcomes in ways that individual-level data alone cannot capture. In practical terms, the same person’s risk of severe outcomes depended not only on their age and health status but also on where and when they lived through the pandemic, a finding with clear implications for how regional health systems should allocate resources during future epidemic emergencies.

The temporal trends documented in the study reveal a pattern that public health officials observed qualitatively but had rarely been able to quantify so systematically. As overall SARS-CoV-2 transmission increased in the population, driven by successive waves of more transmissible variants and eventually by widespread immunity from infection and vaccination, hospital admissions and mortality declined over time. The decline was especially pronounced in areas that had been impacted early, such as the districts at the epicenter of the first wave, where the combination of prior infection, vaccination coverage, and adapted clinical management appears to have progressively decoupled case numbers from severe outcomes. The researchers also tracked test positivity rates and vaccination coverage across the local public health authorities and districts, providing a parallel account of how diagnostic intensity and immunization campaigns varied geographically and interacted with the epidemic’s course.

Methodologically, the study exemplifies the power of routinely collected administrative health data when linked at population scale. Rather than relying on sampled surveys or voluntary reporting, the team exploited the complete registry of Regional Health Service beneficiaries, pseudonymized according to the ISO 25237:2017 standard by each health protection agency before analysis. Only fully anonymized, aggregated data were made available to the study group, a design that satisfied Italian data protection regulations and the EU General Data Protection Regulation while preserving the statistical power of a true total-population cohort. The multivariable regression framework, applied within each variant-defined sub-period, allowed the investigators to estimate how the hazard of infection and death associated with each risk factor changed as the virus and the population’s immunity evolved, a level of temporal resolution that few previous studies have achieved.

The findings also carry lessons that extend well beyond Lombardy. The confirmation that neurologic and renal comorbidities predict mortality, in addition to the more widely recognized cardiac and pulmonary conditions, refines the clinical risk stratification that guided vaccination prioritization and protective measures during the pandemic. The demonstration that geographic location independently shapes outcomes even in a region with a single, relatively uniform health system suggests that epidemic management cannot rely on one-size-fits-all policies; sub-regional surveillance and locally tailored interventions matter. And the observation that severity declined as transmission rose, reflecting the transition from a naive population facing a lethal novel virus to an immunized population facing endemic circulation, provides a quantitative template for anticipating how future respiratory pathogens might behave as they make the same transition.

Lombardy’s experience was in many respects the harshest test any European health system faced in 2020, when hospitals in Bergamo, Brescia, and Cremona became global symbols of the pandemic’s toll. This new analysis converts that traumatic experience into structured evidence, mapping exactly where and when the virus struck, who was most vulnerable, and how the determinants of severe disease shifted across the ancestral, Alpha, Delta, and Omicron eras. The work was funded by Fondazione Cariplo as part of a broader initiative on post-COVID syndrome research, and it forms part of the PASCNET network’s effort to build lasting analytical capacity from the pandemic’s data legacy. As the authors conclude, the findings contribute essential insights into the dynamics of COVID-19 in one of Europe’s earliest and most severely affected regions, evidence that will inform pandemic preparedness, health system planning, and epidemiological modeling for years to come. For a region that served as the pandemic’s European ground zero, the three-year map now drawn represents both a record of what was endured and a guide for what may come.

Subject of Research: Population-based epidemiological mapping of COVID-19 incidence, hospitalization, and mortality in Lombardy, Italy, from 2020 to 2022

Article Title: Mapping the COVID-19 epidemic in Lombardy, Italy, years 2020–2022: a population-based study of temporal trends, geographic distribution, and risk factors for incidence, hospitalization, and mortality

Article References: Villa, S., Magnoni, P., Mazzali, C., Zucchi, A., Maifredi, G., Cavalieri d’Oro, L., Gambino, M. L., Fanetti, A. C., Perotti, P. G., Villa, M., Valsecchi, M. G., Vigani, D., Lucifora, C., Pregliasco, F. E., Cereda, D., Leoni, O., Russo, A. G., the PASCNET study group, Franzoni, F., … Zuccaro, V. (2026). Mapping the COVID-19 epidemic in Lombardy, Italy, years 2020–2022: a population-based study of temporal trends, geographic distribution, and risk factors for incidence, hospitalization, and mortality. BMC Public Health. https://doi.org/10.1186/s12889-026-29682-2

Image Credits: AI Generated

DOI: 10.1186/s12889-026-29682-2

Keywords: COVID-19, Lombardy, SARS-CoV-2, epidemiology, spatial analysis, risk factors, hospitalization, mortality, variants of concern, vaccination coverage, public health, population-based study