Molecular modeling studies on dengue and West Nile Virus NS2B/NS3 protease inhibitor interaction

Abstract

Dengue Virus (DV) and West Nile Virus (WNV) are causes of severe public-health problems in many countries. Neither effective human vaccines nor drugs against both DV and WNV are available. NS2B (an essential cofactor) together with NS3 protease, which is a vital enzyme for the viruses in the replication cycle, is currently focused as the promising enzyme target for drug development. However, the active enzyme-inhibitor complex of DV is not yet experimentally available. Therefore, 3D structure of the DV's enzyme was built using the related WNV NS2B/NS3 structure as a template (2FP7.pdb) [1]. The DV model and 2FP7.pdb (WNV-X) were subjected to molecular dynamics (MD) simulations. MD simulations reveal the important of NS2B for interacting with both NS3 protease and P2 subsite of inhibitor. The final MD snapshot of DV model and WNV-X were employed for GRID field calculations. Results show bigger N+ and OH fields at the S1 pocket of DV model than those of WNV-X implying the importance of the electrostatic or H-bond interaction at the S1 pocket than that of WNV. On the other hand, the Dry field at the S1 pocket of WNV-X is larger than that of DV model indicating the relevance of hydrophobic interaction at this pocket for WNV. In addition, MD simulations of the NS2B/NS3 protease of WNV and DV complexed with Bz-Nle-Lys-Arg-Ala-H (most selective inhibitor, Ki = 193.00 μM for DV and IC50 = 4.6 μM for WNV) [2,3] were also carried out. Results reveal that AlaP1 of the inhibitor is not able to make favourable electrostatic interactions at the S1/S2 pockets of DV. The simulations showed also that both S1-P1 and S2-P2 interactions are important for inhibiting the DV protease whereas in case of the WNV protease the interaction at the S1-pocket is less important.