Investigation on the Binding Properties of N1 Neuraminidase of H5N1 Influenza Virus in Complex with Fluorinated Sialic Acid Analog Compounds—a Study by Molecular Docking and Molecular Dynamics Simulations

Worldwide, only two types of antiviral inhibitors (M2 ion channel protein inhibitor and Neuraminidase inhibitors) are approved to treat the influenza viral infection. But the mutation of amino acid sequence in the viral membrane proteins creates the viral resistance to existing antiviral drugs or inhibitors. So the corresponding antiviral drugs have to be reformulated to match these antigenic variations. Fluorination on the carbon–based molecule significantly enriches its biological properties. Hence this study is motivated to design the fluorinated sialic acid (SIA) analog inhibitors for the neuraminidase of H5N1 influenza A virus by substituting fluorine atom at different hydroxyls (O2, O4, O7, O8, and O9) of sialic acid. 100 ns molecular dynamics simulations are carried out for each protein–ligand complex system. NAMD pair interaction energy and MM–PBSA binding free energy calculations predict two possible binding modes for N1–SIA_F2, N1–SIA_F4, and N1–SIA_F7 complexes and single binding mode for N1–SIA_F8 and N1–SIA_F9 complexes. RMSD, RMSF, and hydrogen bonding analyses are used to understand the conformational flexibility and structural stability of each complex system. It has been concluded that the fluorinated sialic acid drug candidates SIA_F2 and SIA_F7 have better inhibiting potency against the N1 neuraminidase of H5N1 influenza virus.