Simulating molecular mechanisms of the MDM2-mediated regulatory interactions: a conformational selection model of the MDM2 lid dynamics

Deciphering the Multi-state Conformational Equilibrium of HDM2 in the Regulation of p53 Binding: Perspectives from Molecular Dynamics Simulation and NMR Analysis

HDM2 negatively regulates the activity of the tumor suppressor p53. Previous NMR studies have shown that apo-HDM2 interconverts between an "open" state in which the N-terminal "lid" is disordered and a "closed" state in which the lid covers the p53-binding site in the core region. Molecular dynamics (MD) simulation studies have been performed to elucidate the conformational dynamics of HDM2, but the direct relevance of the experimental and computational analyses is unclear. In addition, how the phosphorylation of S17 in the lid contributes to the inhibition of p53 binding remains controversial. Here, we used both NMR and MD simulations to investigate the conformational dynamics of apo-HDM2. The NMR analysis revealed that apo-HDM2 exists in a fast-exchanging equilibrium within two closed states, closed 1 and closed 2, in addition to a previously demonstrated slow-exchanging "open-closed" equilibrium. MD simulations visualized two characteristic closed states, where the spatial orientation of the key residues corresponds well to the chemical shift changes of the NMR spectra. Furthermore, the phosphorylation of S17 induced an equilibrium shift toward closed 1, thereby suppressing the binding of p53 to HDM2. This study reveals a multi-state equilibrium of apo-HDM2 and provides new insights into the regulation mechanism of HDM2-p53 interactions.

Role of the N-terminal lid in regulating the interaction of phosphorylated MDMX with p53

Murine double minute 4 protein (MDMX) is crucial for the regulation of the tumor suppressor protein p53. Phosphorylation of the N-terminal domain of MDMX is thought to affect its binding with the transactivation domain of p53, thus playing a role in p53 regulation. In this study, the effects of MDMX phosphorylation on the binding of p53 were investigated using molecular dynamics simulations. It is shown that in addition to the previously proposed mechanism in which phosphorylated Y99 of MDMX inhibits p53 binding through steric clash with P27 of p53, the N-terminal lid of MDMX also appears to play an important role in regulating the phosphorylation-dependent interactions between MDMX and p53. In the proposed mechanism, phosphorylated Y99 aids in pulling the lid into the p53-binding pocket, thus inhibiting the binding between MDMX and p53. Rebinding of p53 appears to be facilitated by the subsequent phosphorylation of Y55, which draws the lid away from the binding pocket by electrostatic attraction of the lid's positively charged N-terminus. The ability to target these mechanisms for the proper regulation of p53 could have important implications for understanding cancer biology and for drug development.

Protein-peptide molecular docking with large-scale conformational changes: the p53-MDM2 interaction

Protein-peptide interactions are often associated with large-scale conformational changes that are difficult to study either by classical molecular modeling or by experiment. Recently, we have developed the CABS-dock method for flexible protein-peptide docking that enables large-scale rearrangements of the protein chain. In this study, we use CABS-dock to investigate the binding of the p53-MDM2 complex, an element of the cell cycle regulation system crucial for anti-cancer drug design. Experimental data suggest that p53-MDM2 binding is affected by significant rearrangements of a lid region - the N-terminal highly flexible MDM2 fragment; however, the details are not clear. The large size of the highly flexible MDM2 fragments makes p53-MDM2 intractable for exhaustive binding dynamics studies using atomistic models. We performed extensive dynamics simulations using the CABS-dock method, including large-scale structural rearrangements of MDM2 flexible regions. Without a priori knowledge of the p53 peptide structure or its binding site, we obtained near-native models of the p53-MDM2 complex. The simulation results match well the experimental data and provide new insights into the possible role of the lid fragment in p53 binding. The presented case study demonstrates that CABS-dock methodology opens up new opportunities for protein-peptide docking with large-scale changes of the protein receptor structure.

Medicinal Chemistry Strategies to Disrupt the p53-MDM2/MDMX Interaction

The growth inhibitory activity of p53 tumor suppressor is tightly regulated by interaction with two negative regulatory proteins, murine double minute 2 (MDM2) and X (MDMX), which are overexpressed in about half of all human tumors. The elucidation of crystallographic structures of MDM2/MDMX complexes with p53 has been pivotal for the identification of several classes of inhibitors of the p53-MDM2/MDMX interaction. The present review provides in silico strategies and screening approaches used in drug discovery as well as an overview of the most relevant classes of small-molecule inhibitors of the p53-MDM2/MDMX interaction, their progress in pipeline, and highlights particularities of each class of inhibitors. Most of the progress made with high-throughput screening has led to the development of inhibitors belonging to the cis-imidazoline, piperidinone, and spiro-oxindole series. However, novel potent and selective classes of inhibitors of the p53-MDM2 interaction with promising antitumor activity are emerging. Even with the discovery of the 3D structure of complex p53-MDMX, only two small molecules were reported as selective p53-MDMX antagonists, WK298 and SJ-172550. Dual inhibition of the p53-MDM2/MDMX interaction has shown to be an alternative approach since it results in full activation of the p53-dependent pathway. The knowledge of structural requirements crucial to the development of small-molecule inhibitors of the p53-MDMs interactions has enabled the identification of novel antitumor agents with improved in vivo efficacy.

Impact of Ser17 Phosphorylation on the Conformational Dynamics of the Oncoprotein MDM2

MDM2 is an important oncoprotein that downregulates the activity of the tumor suppressor protein p53 via binding of its N-terminal domain to the p53 transactivation domain. The first 24 residues of the MDM2 N-terminal domain form an intrinsically disordered "lid" region that interconverts on a millisecond time scale between "open" and "closed" states in unliganded MDM2. While the former conformational state is expected to facilitate p53 binding, the latter competes in a pseudo-substrate manner with p53 for its binding site. Phosphorylation of serine 17 in the MDM2 lid region is thought to modulate the equilibrium between "open" and "closed" lid states, but contradictory findings on the favored lid conformational state upon phosphorylation have been reported. Here, the nature of the conformational states of MDM2 pSer17 and Ser17Asp variants was addressed by means of enhanced sampling molecular dynamics simulations. Detailed analyses of the computed lid conformational ensembles indicate that both lid variants stabilize a "closed" state, with respect to wild type. Nevertheless, the nature of the closed-state conformational ensembles differs significantly between the pSer17 and Ser17Asp variants. Thus, care should be applied in the interpretation of biochemical experiments that use phosphomimetic variants to model the effects of phosphorylation on the structure and dynamics of this disordered protein region.

Probing Origin of Binding Difference of inhibitors to MDM2 and MDMX by Polarizable Molecular Dynamics Simulation and QM/MM-GBSA Calculation

Binding abilities of current inhibitors to MDMX are weaker than to MDM2. Polarizable molecular dynamics simulations (MD) followed by Quantum mechanics/molecular mechanics generalized Born surface area (QM//MM-GBSA) calculations were performed to investigate the binding difference of inhibitors to MDM2 and MDMX. The predicted binding free energies not only agree well with the experimental results, but also show that the decrease in van der Walls interactions of inhibitors with MDMX relative to MDM2 is a main factor of weaker bindings of inhibitors to MDMX. The analyses of dihedral angles based on MD trajectories suggest that the closed conformation formed by the residues M53 and Y99 in MDMX leads to a potential steric clash with inhibitors and prevents inhibitors from arriving in the deep of MDMX binding cleft, which reduces the van der Waals contacts of inhibitors with M53, V92, P95 and L98. The calculated results using the residue-based free energy decomposition method further prove that the interaction strength of inhibitors with M53, V92, P95 and L98 from MDMX are obviously reduced compared to MDM2. We expect that this study can provide significant theoretical guidance for designs of potent dual inhibitors to block the p53-MDM2/MDMX interactions.

Elucidation of Ligand-Dependent Modulation of Disorder-Order Transitions in the Oncoprotein MDM2

Numerous biomolecular interactions involve unstructured protein regions, but how to exploit such interactions to enhance the affinity of a lead molecule in the context of rational drug design remains uncertain. Here clarification was sought for cases where interactions of different ligands with the same disordered protein region yield qualitatively different results. Specifically, conformational ensembles for the disordered lid region of the N-terminal domain of the oncoprotein MDM2 in the presence of different ligands were computed by means of a novel combination of accelerated molecular dynamics, umbrella sampling, and variational free energy profile methodologies. The resulting conformational ensembles for MDM2, free and bound to p53 TAD (17-29) peptide identify lid states compatible with previous NMR measurements. Remarkably, the MDM2 lid region is shown to adopt distinct conformational states in the presence of different small-molecule ligands. Detailed analyses of small-molecule bound ensembles reveal that the ca. 25-fold affinity improvement of the piperidinone family of inhibitors for MDM2 constructs that include the full lid correlates with interactions between ligand hydrophobic groups and the C-terminal lid region that is already partially ordered in apo MDM2. By contrast, Nutlin or benzodiazepinedione inhibitors, that bind with similar affinity to full lid and lid-truncated MDM2 constructs, interact additionally through their solubilizing groups with N-terminal lid residues that are more disordered in apo MDM2.

Life as we know it depends on interactions between proteins. There is substantial evidence that many interactions between proteins involve very flexible protein regions. These disordered regions may undergo disorder/order transitions upon forming an interaction with another protein. Many successful approaches to medicinal chemistry are based on mimicking the interactions of biological molecules with man-made small molecules. However how drug-like small-molecules may modulate protein disorder is currently poorly understood, largely because it is difficult to measure in details this type of interaction with experimental methods. Here we have used computer simulations to resolve with great details the process by which different small-molecules modulate the flexibility of a disordered region of the protein MDM2. This protein is overexpressed in many cancers and small-molecules that recognize MDM2 have been developed over the last decade as possible novel anti-cancer agents. We show that the flexible MDM2 “lid” region adopts different conformational states in the presence of different small-molecules. Our results suggest why some classes of small-molecules form favorable interactions with the lid region, whereas others do not. These findings may prove crucial to develop new and more effective MDM2 inhibitors, and more generally to help drug designers target disordered proteins regions with small-molecules.

Binding modes of three inhibitors 8CA, F8A and I4A to A-FABP studied based on molecular dynamics simulation

Adipocyte fatty-acid binding protein (A-FABP) is an important target of drug designs treating some diseases related to lipid-mediated biology. Molecular dynamics (MD) simulations coupled with solvated interaction energy method (SIE) were carried out to study the binding modes of three inhibitors 8CA, F8A and I4A to A-FABP. The rank of our predicted binding affinities is in accordance with experimental data. The results show that the substitution in the position 5 of N-benzyl and the seven-membered ring of N-benzyl-indole carboxylic acids strengthen the I4A binding, while the substitution in the position 2 of N-benzyl weakens the F8A binding. Computational alanine scanning and dynamics analyses were performed and the results suggest that the polar interactions of the positively charged residue R126 with the three inhibitors provide a significant contribution to inhibitor bindings. This polar interaction induces the disappearance of the correlated motion of the C terminus of A-FABP relative to the N terminus and favors the stability of the binding complex. This study is helpful for the rational design of potent inhibitors within the fields of metabolic disease, inflammation and atherosclerosis.

Inhibition of nutlin-resistant HDM2 mutants by stapled peptides

Pharmacological modulation of p53 activity is an attractive therapeutic strategy in cancers with wild-type p53. Presently in clinical trials, the small molecule Nutlin-3A competitively binds to HDM2, a key negative regulator of p53 and blocks its activity. We have described resistance mutations in HDM2 that selectively reduce affinity for Nutlin but not p53. In the present communication, we show that stapled peptides targeting the same region of HDM2 as Nutlin are refractory to these mutations, and display reduced discrimination between the wild-type and mutant HDM2s with regards to functional abrogation of interaction with p53. The larger interaction footprint afforded by stapled peptides suggests that this class of ligands may prove comparatively more resilient to acquired resistance in a clinical setting.

Effective screening strategy using ensembled pharmacophore models combined with cascade docking: application to p53-MDM2 interaction inhibitors

Protein-protein interactions (PPIs) play a crucial role in cellular function and form the backbone of almost all biochemical processes. In recent years, protein-protein interaction inhibitors (PPIIs) have represented a treasure trove of potential new drug targets. Unfortunately, there are few successful drugs of PPIIs on the market. Structure-based pharmacophore (SBP) combined with docking has been demonstrated as a useful Virtual Screening (VS) strategy in drug development projects. However, the combination of target complexity and poor binding affinity prediction has thwarted the application of this strategy in the discovery of PPIIs. Here we report an effective VS strategy on p53-MDM2 PPI. First, we built a SBP model based on p53-MDM2 complex cocrystal structures. The model was then simplified by using a Receptor-Ligand complex-based pharmacophore model considering the critical binding features between MDM2 and its small molecular inhibitors. Cascade docking was subsequently applied to improve the hit rate. Based on this strategy, we performed VS on NCI and SPECS databases and successfully discovered 6 novel compounds from 15 hits with the best, compound 1 (NSC 5359), K(i) = 180 ± 50 nM. These compounds can serve as lead compounds for further optimization.