Infections caused by the Zika virus (ZIKV) have detrimental effects on human health, in particular on infants. As no potent drug or vaccine is available to date to contain this viral disease, it is necessary to design inhibitors that can target the NS2B-NS3 protease of the ZIKV, which is mainly responsible for the proliferation of the virus inside the host cells . Here, molecular dynamics (MD) simulation and molecular mechanics energies combined with the generalized Born and surface area continuum solvation model (MM/GBSA) are used to understand the binding modes and stabilities of R, KR, KKR, WKR, WKKR, YKKR, and FKKR peptide inhibitors bound to the NS3-NS2B protease. The results are compared with the corresponding results obtained for covalent (compound ) and non-covalent (compound ) peptidomimetic inhibitors . It is revealed that peptide inhibitors can bind strongly with the ZIKV protease with the ΔG ranging from -12 kcal/mol to -73 kcal/mol. Among these peptides, YKKR is found to make the most stable complex with the protease and fully occupy the electrostatically active substrate binding site. Hence, it would inhibit the protease activities of ZIKV strongly. The residue-wise decomposition of ΔG indicates that Asp75, Asp129, Tyr130, Ser135, Gly151, Asn152, Glys153, and Tyr161 of NS3 and Ser81, Asp83, and Phe84 of NS2B play a prominent role in the inhibitor binding. Therefore, any future design of inhibitors should be aimed to target these residues.

Author