Glycan binding and specificity of viral influenza neuraminidases by classical molecular dynamics and replica exchange

Jiraphorn Phanich1, Siraphob Threeracheep2, Nawee Kungwan3,4

  • 1a Faculty of Science, Department of Chemistry, Center of Excellence in Computational Chemistry (CECC) , Chulalongkorn University , Bangkok , Thailand.

Two important glycoproteins on the influenza virus membrane, hemagglutinin (HA) and neuraminidase (NA), are relevant to virus replication. As previously reported, HA has a substrate specificity towards SIA-2,3-GAL-1,4-NAG (3SL) and SIA-2,6-GAL-1,4-NAG (6SL) glycans, while NA can cleave both types of linkages. However, the substrate binding into NA and its preference are not well understood. In this work, the glycan binding and specificity of human and avian NAs were evaluated by classical molecular dynamics (MD) simulations, whilst the conformational diversity of 3SL avian and 6SL human glycans in an unbound state was investigated by replica exchange MD simulations. The results indicated that the 3SL avian receptor fits well in the binding cavity of all NAs and does not require a conformational change for such binding compared to the flexible shape of the 6SL human receptor. From the QM/MM-GBSA binding free energy and decomposition free energy data, 6SL showed a much stronger binding towards human NAs (H1N1, H2N2 and H3N2) than to avian NAs (H5N1 and H7N9). This suggests that influenza NAs have a substrate specificity corresponding to their HA, indicating the functional balance between the two important glycoproteins. Both linkages show distinct glycan topologies when complexed with NAs, while the flexibility of torsion angles between GAL and NAG in 6SL results in the various shapes of glycan and different binding patterns. Lower conformational diversities of both glycans when bound to NA compared to the unbound state were found, and were required in order to be accommodated within the NA cavity. Communicated by Ramaswamy H. Sarma.

Related Concept Videos

Molecular Models02:00

Molecular Models

Physical models representing molecular architectures of chemical compounds play essential roles in understanding chemistry. The use of molecular models makes it easier to visualize the structures and shapes of atoms and molecules.
43.8K
Dynamic Equilibrium02:20

Dynamic Equilibrium

A reversible chemical reaction represents a chemical process that proceeds in both forward (left to right) and reverse (right to left) directions. When the rates of the forward and reverse reactions are equal, the concentrations of the reactant and product species remain constant over time and the system is at equilibrium. A special double arrow is used to emphasize the reversible nature of the reaction. The relative concentrations of reactants and products in equilibrium systems vary greatly;...
63.0K
Molecular Orbital Theory I02:35

Molecular Orbital Theory I

Overview of Molecular Orbital Theory
47.7K
Molecular Orbital Theory II03:51

Molecular Orbital Theory II

Molecular Orbital Energy Diagrams
27.6K
Predicting Molecular Geometry02:27

Predicting Molecular Geometry

VSEPR Theory for Determination of Electron Pair Geometries
46.0K
Molecular and Ionic Solids02:54

Molecular and Ionic Solids

Crystalline solids are divided into four types: molecular, ionic, metallic, and covalent network based on the type of constituent units and their interparticle interactions.
Molecular Solids
Molecular crystalline solids, such as ice, sucrose (table sugar), and iodine, are solids that are composed of neutral molecules as their constituent units. These molecules are held together by weak intermolecular forces such as London dispersion forces, dipole-dipole interactions, or hydrogen bonds, which...
20.2K