Understanding of known drug-target interactions in the catalytic pocket of neuraminidase subtype N1

To provide detailed information and insight into the drug-target interaction, structure, solvation, and dynamic and thermodynamic properties, the three known-neuraminidase inhibitors-oseltamivir (OTV), zanamivir (ZNV), and peramivir (PRV)-embedded in the catalytic site of neuraminidase (NA) subtype N1 were studied using molecular dynamics simulations. In terms of ligand conformation, there were major differences in the structures of the guanidinium and the bulky groups. The atoms of the guanidinium group of PRV were observed to form many more hydrogen bonds with the surrounded residues and were much less solvated by water molecules, in comparison with the other two inhibitors. Consequently, D151 lying on the 150-loop (residues 147-152) of group-1 neuraminidase (N1, N4, N5, and N8) was considerably shifted to form direct hydrogen bonds with the --OH group of the PRV, which was located rather far from the 150-loop. For the bulky group, direct hydrogen bonds were detected only between the hydrophilic side chain of ZNV and residues R224, E276, and E277 of N1 with rather weak binding, 20-70% occupation. This is not the case for OTV and PRV, in which flexibility and steric effects due to the hydrophobic side chain lead to the rearrangement of the surrounded residues, that is, the negatively charged side chain of E276 was shifted and rotated to form hydrogen bonds with the positively charged moiety of R224. Taking into account all the ligand-enzyme interaction data, the gas phase MM interaction energy of -282.2 kcal/mol as well as the binding free energy (DeltaG(binding)) of -227.4 kcal/mol for the PRV-N1 are significantly lower than those of the other inhibitors. The ordering of DeltaG(binding) of PRV < ZNV < OTV agrees well with the ordering of experimental IC(50) value.     

Molecular prediction of oseltamivir efficiency against probable influenza A (H1N1-2009) mutants: molecular modeling approach

To predict the susceptibility of the probable 2009 influenza A (H1N1-2009) mutant strains to oseltamivir, MD/LIE approach was applied to oseltamivir complexed with the most frequent drug-resistant strains of neuraminidase subtypes N1 and N2: two mutations on the framework residues (N294S and H274Y) and the two others on the direct-binding residues (E119V and R292K) of oseltamivir. Relative to those of the wild type (WT), loss of drug–target interaction energies, especially in terms of electrostatic contributions and hydrogen bonds were dominantly established in the E119V and R292K mutated systems. The inhibitory potencies of oseltamivir towards the WT and mutants were predicted according to the ordering of binding-free energies: WT (−12.3 kcal mol⁻¹) > N294S (−10.4 kcal mol⁻¹) > H274Y (−9.8 kcal mol⁻¹) > E119 V (−9.3 kcal mol⁻¹) > R292K (−7.7 kcal mol⁻¹), suggesting that the H1N1-2009 influenza with R292K substitution, perhaps, conferred a high level of oseltamivir resistance, while the other mutants revealed moderate resistance levels. This result calls for an urgent need to develop new potent anti-influenza agents against the next pandemic of potentially higher oseltamivir-resistant H1N1-2009 influenza.