Synergistic effects in the designs of neuraminidase ligands: Analysis from docking and molecular dynamics studies

Docking and molecular dynamics were used to study the nine ligands (see Scheme 1) at the neuraminidase (NA) active sites. Their binding modes are structurally and energetically different, with details given in the text. Compared with 1A (oseltamivir carboxylate), the changes of core template or/and functional groups in the other ligands cause the reductions of interaction energies and numbers of H-bonds with the NA proteins. Nonetheless, all these ligands occupy the proximity space at the NA active sites and share some commonness in their binding modes. The fragment approach was then used to analyze and understand the binding specificities of the nine ligands. The contributions of each core template and functional group were evaluated. It was found that the core templates rather than functional groups play a larger role during the binding processes; in addition, the binding qualities are determined by the synergistic effects of the core templates and functional groups. Among the nine ligands, 1A (oseltamivir carboxylate) has the largest synergistic energy and its functional groups fit perfectly with the NA active site, consistent with the largest interaction energy, numerous H-bonds with the NA active-site residues as well as experimentally lowest IC50 value. Owing to the poorer metabolizability than oseltamivir, large contribution of the benzene core template and fine synergistic effects of the functional groups, the 4-(N-acetylamino)-5-guanidino-3-(3-pentyloxy)benzoic acid should be an ideal lead compound for optimizing NA drugs.