Elucidation of Conformational Dynamics of MDM2 and Alterations Induced Upon Inhibitor Binding Using Elastic Network Simulations and Molecular Docking

Abstract

Elastic network model simulations were performed to investigate the conformational changes of MDM2 protein induced by its native substrate p53 and two small molecule inhibitor (NVP-CGM097 and HDM 201) bindings. Residues Phe 19, Trp 23, Leu 26 were observed to reside in the minima of slowest modes of p53, pointing to the accepted three finger binding model. Pro 27 displays the most significant hinge present in p53 and comes out to be another functionally important residue. Three distinct conserved regions are identified in MDM2. Regions I (residues 50–77) and III (residues 90–105) correspond to the binding interface of MDM2 including [Formula: see text] helix-2 ([Formula: see text]), Loop-2 (L2) and [Formula: see text] helix-4 ([Formula: see text]) domains which are stabilized during complex formation. Region II (residues 77–90) is a highly flexible region in both unbound and bound forms exhibiting high amplitude collective motion. MDM2 exhibits a scattered profile in the fastest modes of motion, while binding of p53 and inhibitors puts restraints on MDM2 pointing to induced-fit mechanism. Flexible docking using AutoDock Vina is performed to account for the flexible nature of the receptor and to elucidate the essential interactions in the binding cleft. [Formula: see text] domain controls the size of the cleft by keeping the cleft narrow in unbound MDM2; and open in the bound states. Inhibitors studied in this work (NVP-CGM097 and HDM201), which are recently undergoing clinical trials, succeed in mimicking p53 behavior which would shed light on the rational design of novel anticancer drugs.