Acta Phys. -Chim. Sin. ›› 2016, Vol. 32 ›› Issue (11): 2693-2708.doi: 10.3866/PKU.WHXB201608121

• ARTICLE • Previous Articles     Next Articles

QSAR, Molecular Docking and Molecular Dynamics of 3C-like Protease Inhibitors

Feng LIN1,Xin-Mei FU2,*(),Chao WANG1,Si-Yu JIANG1,Jing-Hui WANG1,Shu-Wei ZHANG1,Ling YANG3,Yan LI1,*()   

  1. 1 School of Chemical Engineering, Dalian University of Technology, Dalian 116024, Liaoning Province, P. R. China
    2 State Key Laboratory of Fine Chemicals, Dalian University of Technology, Dalian 116024, Liaoning Province, P. R. China
    3 Laboratory of Pharmaceutical Resource Discovery, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, Liaoning Province, P. R. China
  • Received:2016-05-31 Published:2016-11-08
  • Contact: Xin-Mei FU,Yan LI E-mail:fuxinmei@dlut.edu.cn;yanli@dlut.edu.cn
  • Supported by:
    the National Natural Science Foundation of China(11201049)

Abstract:

3C-like protease is an extremely important protease involved in the multiplicative process of coronaviruses, including the deadlyMiddle East respiratory syndrome coronavirus (MERS-CoV). 3C-like protease has become a hot research topic in the field of coronavirology. For the first time, a set of ligand-and receptorbased three-dimensional quantitative structure-activity relationships (3D-QSAR) models were carried out via comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) to explore the structure-activity correlation of 43 peptidomimetic inhibitors of the 3C-like protease of the bat coronavirus HKU4 (HKU4-CoV), which belongs to the same 2c lineage as MERS-CoV and shows high sequence similarity with MERS-CoV. Based on the ligand-based alignment, an optimal CoMSIAmodel (yielded by steric, electrostatic, H-bond donor and H-bond acceptor fields) was obtained with good predictive power of Q2=0.522, R2 ncv=0.996 and R2 pre=0.904 (Q2:cross-validated correlation coefficient, R2 ncv:non-cross-validated correlation coefficient, R2 pre:predicted correlation coefficient for the test set of compounds). Molecular docking and molecular dynamics simulations were performed according to this model to further determine the interaction mechanism between ligands and the receptor. The experimental results show:(1) based on the optimal CoMSIAmodel, the 3D contour maps vividly illustrate that the molecular biological activity is influenced by the steric, electrostatic, H-bond donor and H-bond acceptor interactions of molecular groups. (2) Based on the docking analysis, hydrophobicity, crystal water, His166 andGlu169 have important roles in the ligands and receptor binding process. (3) Molecular dynamics (MD) simulations were carried out for further verification of the reliability of the docking model, and provide two new key residues, Ser24 and Gln192, which have two strong hydrogen bonds with the ligands. Some new compounds were obtained based on the modeling that are potential peptidomimetic inhibitors of 3C-like protease. These results help establish the binding mechanism between 3C-like protease and peptidomimetic inhibitors, and provide a valuable reference for future anti-MERS-CoV drug design.

Key words: MERS-CoV, 3C-like protease, Peptidomimetic inhibitor, 3D-QSAR, Molecular docking, Molecular dynamics