# Structural analysis of antigenic variation and adaptive evolution of the H5N1 neuraminidase gene

## TL;DR

This paper studies how the H5N1 influenza virus's neuraminidase enzyme evolves and functions, identifying key residues like K207 that control its dynamic behavior and suggesting new antiviral strategies.

## Contribution

The study identifies K207 as a master regulator of neuraminidase dynamics and proposes a novel antiviral strategy targeting this dynamic vulnerability.

## Key findings

- Structural analysis of nine NA subtypes reveals four phylogenetic subgroups and 40 conserved functional residues.
- K207 substitutions (K207W, K207H, K207I) cause distinct structural and dynamic outcomes, re-wiring allosteric communication.
- Phylogenetic analysis shows host-specific selection pressures, with gallus isolates showing hyper-reassortment and human isolates showing sporadic spillover.

## Abstract

The concern regarding H5N1 outbreak, particularly the accelerated mutagenesis of its core genomic elements, underscores the persistent threat of influenza to global health. Neuraminidase (NA), a pivotal sialidase integral to virion egress and propagation, comprises nine distinct isoforms, exhibiting unique evolutionary trajectories and structural adaptations. Despite extensive characterization of hemagglutinin subtypes, the functional divergence of the nine NA subtypes remains inadequately understood. To address this gap, we conducted a structural analysis of NA subtypes, employing structural superimposition and motif-guided sequence alignment to delineate subtype-specific residues. Hierarchical clustering stratified the nine NA subtypes into four distinct subgroups: NA2 (subgroup I), NA1 and NA4 (subgroup II), NA9/NA7/NA6/NA3 (subgroup III), and NA8/NA5 (subgroup 4). We identified 40 highly conserved and functionally significant amino acid loci, likely modulating enzymatic activity and substrate specificity across subtypes. To investigate the structural basis of adaptation in H5N1, we generated NA1 mutants by swapping family specific position (FSP) residues and analyzed their dynamics using Molecular Dynamics (MD) simulations, complemented by a deep phylogenetic analysis across six host reservoirs. MD simulation parameters reveal a dynamic paradox: the Wild-Type (WT) NA1 maintains a conserved global compactness Rg, which masks a complex, bi-modal switching mechanism essential for its catalytic function, validated by multi-basin free energy landscape (FEL) topography. We identify Lysine-207 (K207) as the master determinant of this switching mechanism and the enzyme’s dynamic fate. Substitutions at this conserved nexus produced diametrically opposite outcomes: K207W imposed structural rigidification (abolishing the switch), K207H achieved dynamic preservation, and K207I drove expanded disorder and collapse. Furthermore, dynamic correlation analysis shows that these single-point substitutions function as molecular switches that significantly re-wire the enzyme’s allosteric communication networks, extending far beyond the active site. To assess the role of NA1 in host tropism and adaptive evolution, we conducted a phylogenetic analysis of NA1 genes from H5N1 isolates across multiple host reservoirs; H. sapiens, G. gallus, Anser anser domesticus, M. gallopavo, B. taurus, and C. olor. Notably, we observed opposing selection pressures and diversification patterns: G. gallus isolates showed signatures of positive selection consistent with hyper-reassortment, while human isolates displayed highly diverse, sporadic spillover events. We conclude that the evolutionary contribution of NA1 to H5N1 host adaptation is not encoded in static structure, but certain residues such as K207 defines a pivotal mechanism for regulating the enzyme’s function through dynamic states. Our MD data thus proposes a novel strategy for next-generation antivirals by targeting this dynamic vulnerability—the Nexus for Dynamic Ablation—to permanently entrain the enzyme in a non-functional conformation.

The persistent threat of H5N1 influenza demands a deeper understanding of its critical enzyme, Neuraminidase (NA). We first conducted a structural analysis of all nine NA subtypes, stratifying them into four distinct phylogenetic subgroups and identifying 40 highly conserved functional residues. Focusing on the H5N1 NA1, Molecular Dynamics simulations revealed a dynamic paradox: while the enzyme maintains conserved global compactness, its catalytic function is regulated by a complex, bi-modal switching mechanism. We identified Lysine-207 (K207) as the master determinant of this dynamic fate. Substitutions at K207 produced diametrically opposite outcomes—from structural rigidification to expanded disorder—and significantly re-wired the enzyme’s allosteric communication networks far beyond the active site. Phylogenetic analysis of NA1 across host reservoirs further showed that NA1’s adaptive evolution is defined by these dynamic states, not static structure. Our data proposes a novel strategy for next-generation antivirals by targeting this dynamic vulnerability—the Nexus for Dynamic Ablation—to permanently lock the enzyme into a non-functional state.