Molecular dynamics simulations studies and free energy analysis on inhibitors of MDM2–p53 interaction

•The MD simulations and MM-GBSA analyze the differences of three potent inhibitors binding to MDM2.•The residue located at 22′ is also important in inhibitor–MDM2 binding except F19′, W23′ and L26′.•Conformation variability of residue Y100 influences the binding conformation of inhibitor–MDM2 complex.

The oncoprotein MDM2 (murine double minute 2) negatively regulates the activity and stability of tumor suppressor p53. Inactivation of the MDM2–p53 interaction by potent inhibitors offers new possibilities for anticancer therapy. Molecular dynamics (MD) simulations were performed on three inhibitors–MDM2 complexes to investigate the stability and structural transitions. Simulations show that the backbone of MDM2 maintains stable during the whole time. However, slightly structural changes of inhibitors and MDM2 are observed. Furthermore, the molecular mechanics generalized Born surface area (MM-GBSA) approach was introduced to analyze the interactions between inhibitors and MDM2. The results show that binding of inhibitor pDIQ to MDM2 is significantly stronger than that of pMI and pDI to MDM2. The side chains of residues have more contribution than backbone of residues in energy decomposition. The structure–affinity analyses show that L54, I61, M62, Y67, Q72, H73 and V93 produce important interaction energy with inhibitors. The residue W/Y22′ is also very important to the interaction between the inhibitors and MDM2. The three-dimensional structures at different times indicate that the mobility of Y100 influences on the binding of inhibitors to MDM2, and its change has important role in conformations of inhibitors and MDM2.

Graphical abstractFigure optionsDownload full-size imageDownload high-quality image (354 K)Download as PowerPoint slide