Quantitative lid dynamics of MDM2 reveals differential ligand binding modes of the p53-binding cleft

The oncoprotein MDM2 regulates the activity and stability of the tumor suppressor p53 through protein-protein interaction involving their N-terminal domains. The N-terminal lid of MDM2 has been implicated in p53 regulation; however, due to its flexible nature, limited data are available concerning its role in ligand binding. The quantitative dynamics study using NMR reported here shows, for the first time, that the lid in apo-MDM2 slowly interconverts between a "closed" state that is associated with the p53-binding cleft and an "open" state that is highly flexible. Our results reveal that apo-MDM2 predominantly populates the closed state, whereas the p53-bound MDM2 exclusively populates the open state. Unlike p53 binding, the small molecule MDM2 antagonist nutlin-3 binds to the cleft essentially without perturbing the closed lid state. The lid dynamics thereby represents a signature for the experimental and virtual screening of therapeutic antagonists that target the p53-MDM2 interaction.     

Prospective virtual screening for novel p53–MDM2 inhibitors using ultrafast shape recognition

The p53 protein, known as the guardian of genome, is mutated or deleted in approximately 50 % of human tumors. In the rest of the cancers, p53 is expressed in its wild-type form, but its function is inhibited by direct binding with the murine double minute 2 (MDM2) protein. Therefore, inhibition of the p53–MDM2 interaction, leading to the activation of tumor suppressor p53 protein presents a fundamentally novel therapeutic strategy against several types of cancers. The present study utilized ultrafast shape recognition (USR), a virtual screening technique based on ligand–receptor 3D shape complementarity, to screen DrugBank database for novel p53–MDM2 inhibitors. Specifically, using 3D shape of one of the most potent crystal ligands of MDM2, MI-63, as the query molecule, six compounds were identified as potential p53–MDM2 inhibitors. These six USR hits were then subjected to molecular modeling investigations through flexible receptor docking followed by comparative binding energy analysis. These studies suggested a potential role of the USR-selected molecules as p53–MDM2 inhibitors. This was further supported by experimental tests showing that the treatment of human colon tumor cells with the top USR hit, telmisartan, led to a dose-dependent cell growth inhibition in a p53-dependent manner. It is noteworthy that telmisartan has a long history of safe human use as an approved anti-hypertension drug and thus may present an immediate clinical potential as a cancer therapeutic. Furthermore, it could also serve as a structurally-novel lead molecule for the development of more potent, small-molecule p53–MDM2 inhibitors against variety of cancers. Importantly, the present study demonstrates that the adopted USR-based virtual screening protocol is a useful tool for hit identification in the domain of small molecule p53–MDM2 inhibitors.     

Calculation of hot spots for protein–protein interaction in p53/PMI-MDM2/MDMX complexes

The recently developed MM/GBSA_IE method is applied to computing hot and warm spots in p53/PMI-MDM2/MDMX protein–protein interaction systems. Comparison of the calculated hot (>2 kcal/mol) and warm spots (>1 kcal/mol) in P53 and PMI proteins interacting with MDM2 and MDMX shows a good quantitative agreement with the available experimental data. Further, our calculation predicted hot spots in MDM2 and MDMX proteins in their interactions with P53 and PMI and they help elucidate the interaction mechanism underlying this important PPI system. In agreement with the experimental result, the present calculation shows that PMI has more hot and warm spots and binds stronger to MDM2/MDMX. The analysis of these hot and warm spots helps elucidate the fundamental difference in binding between P53 and PMI to the MDM2/MDMX systems. Specifically, for p53/PMI-MDM2 systems, p53 and PMI use essentially the same residues (L54, I61, Y67, Q72, V93, H96, and I99) of MDM2 for binding. However, PMI enhanced interactions with residues L54, Y67, and Q72 of MDM2. For the p53/PMI-MDMX system, p53 and PMI use similar residues (M53, I60, Y66, Q71, V92, and Y99) of MDMX for binding. However, PMI exploited three extra residues (M61, K93, and L98) of MDMX for enhanced binding. In addition, PMI enhanced interaction with four residues (M53, Y66, Q71, and Y99) of MDMX. These results gave quantitative explanation on why the binding affinities of PMI-MDM2/MDMX interactions are stronger than that of p53-MDM2/MDMX although their binding modes are similar. © 2018 Wiley Periodicals, Inc.     

Forces mediating protein–protein interactions: a computational study of p53 “approaching” MDM2

The protein MDM2 forms a complex with the tumor suppressing protein p53 and targets it for proteolysis in order to down-regulate p53 in normal cells. Inhibition of this interaction is of therapeutic importance. Molecular dynamics simulations of the association between p53 and MDM2 have revealed mutual modulation of the two surfaces. Analysis of the simulations of the two species approaching each other in solution shows how long range electrostatics steers these two proteins together. The net electrostatics is controlled largely by a few cationic residues that surround the MDM2 binding site. There is an overall separation in electrostatics of MDM2 and p53 that are mutually complementary and drive association. Upon close approach, there is significant energetic strain as the charges are occluded from water (desolvated). However, the complexation is driven by packing interactions that lead to highly favorable van der Waals interactions. Although the complementarity of the electrostatics of the two surfaces is essential for the two partners to form a complex, steric collisions of Y100 and short ranged van der Waals interactions of F19, W23, L26 of p53 determine the final steps of native complex formation. The electrostatics seem to be evolutionarily conserved, including variations in both partners.     

Long-range modulation of chain motions within the intrinsically disordered transactivation domain of tumor suppressor p53

The tumor suppressor p53 is a hub protein with a multitude of binding partners, many of which target its intrinsically disordered N-terminal domain, p53-TAD. Partners, such as the N-terminal domain of MDM2, induce formation of local structure and leave the remainder of the domain apparently disordered. We investigated segmental chain motions in p53-TAD using fluorescence quenching of an extrinsic label by tryptophan in combination with fluorescence correlation spectroscopy (PET-FCS). We studied the loop closure kinetics of four consecutive segments within p53-TAD and their response to protein binding and phosphorylation. The kinetics was multiexponential, showing that the conformational ensemble of the domain deviates from random coil, in agreement with previous findings from NMR spectroscopy. Phosphorylations or binding of MDM2 changed the pattern of intrachain kinetics. Unexpectedly, we found that upon binding and phosphorylation chain motions were altered not only within the targeted segments but also in remote regions. Long-range interactions can be induced in an intrinsically disordered domain by partner proteins that induce apparently only local structure or by post-translational modification.     

Protein–protein recognition: a computational mutagenesis study of the MDM2–P53 complex

Protein P53 is involved in more than 50% of the human cancers and the P53–MDM2 complex is a target for anticancer drug design. It is possible to engineer small P53 mimics that would be expected to disrupt the P53–MDM2 complex, and release P53 to initiate cell-cycle arrest or apoptosis. These small peptides should bind to the functional epitopes of the protein–protein interface, and prevent the interaction between P53 and MDM2. Here, we apply an improved computational alanine scanning mutagenesis method, which allows the determination of the hot spots present in both monomers, P53 and MDM2, of three protein complexes (the P53-binding domain of human MDM2, its analogue from Xenopus laevis, and the structure of human MDM2 in complex with an optimized P53 peptide). The importance of the hydrogen bonds formed by the protein backbone has been neglected due to the difficulty of measuring experimentally their contribution to the binding free energy. In this study we present a computational approach that allows the estimation of the contribution to the binding free energy of the C=O and N–H groups in the backbone of the P53 and MDM2 proteins. We have noticed that the hydrogen bond between the HE1 atom of the hot spot Trp23 and the O atom of the residue Leu54, as well as the NH-pi hydrogen bond between the Ile57 and Met58 were observed in the Molecular dynamics simulation, and their contribution to the binding free energy measured. This study not only shows the reliability of the computational mutagenesis method to detect hot spots but also demonstrates an excellent correlation between the quantitative calculated binding free energy contribution of the C=O and N–H backbone groups of the interfacial residues and the qualitative values expected for this kind of interaction. The study also increases our understanding of the P53–MDM2 interaction.     

A computational analysis of the binding model of MDM2 with inhibitors

Abstract  

It is a new and promising strategy for anticancer drug design to block the MDM2-p53 interaction using a non-peptide small-molecule inhibitor. We carry out molecular dynamics simulations to study the binding of a set of six non-peptide small-molecule inhibitors with the MDM2. The relative binding free energies calculated using molecular mechanics Poisson–Boltzmann surface area method produce a good correlation with experimentally determined results. The study shows that the van der Waals energies are the largest component of the binding free energy for each complex, which indicates that the affinities of these inhibitors for MDM2 are dominated by shape complementarity. The A-ligands and the B-ligands are the same except for the conformation of 2,2-dimethylbutane group. The quantum mechanics and the binding free energies calculation also show the B-ligands are the more possible conformation of ligands. Detailed binding free energies between inhibitors and individual protein residues are calculated to provide insights into the inhibitor-protein binding model through interpretation of the structural and energetic results from the simulations. The study shows that G1, G2 and G3 group mimic the Phe19, Trp23 and Leu26 residues in p53 and their interactions with MDM2, but the binding model of G4 group differs from the original design strategy to mimic Leu22 residue in p53.     

Targeting the conformational transitions of MDM2 and MDMX: insights into dissimilarities and similarities of p53 recognition

MDM2 and MDMX are oncogenic homologue proteins that regulate the activity and stability of p53, a tumor suppressor protein involved in more than 50% of human cancers. While the large body of experiments so far accumulated has validated MDM2 as a therapeutically important target for the development of anticancer drugs, it is only recently that MDMX has also become an attractive target for the treatment of tumor cells expressing wild type p53. The availability of structural information of the N-terminal domain of MDM2 in complex with p53-derived peptides and inhibitors, and the very recent disclosure of the crystal structure of the N-terminal domain of MDMX bound to a p53 peptide, offer an unprecedented opportunity to provide insight into the molecular basis of p53 recognition and the identification of discriminating features affecting the binding of the tumor suppressor protein at MDM2 and MDMX. By using coarse graining simulations, in this study we report the exploration of the conformational transitions featured in the pathway leading from the apo-MDM2 and apo-MDMX states to the p53-bound MDM2 and p53-bound MDMX states, respectively. The results have enabled us to identify a pool of diverse conformational states of the oncogenic proteins that affect the binding of p53 and the presence of conserved and non-conserved interactions along the conformational transition pathway that may be exploited in the design of selective and dual modulators of MDM2 and MDMX activity.     

Interrogation of MDM2 phosphorylation in p53 activation using native chemical ligation: the functional role of Ser17 phosphorylation in MDM2 reexamined

The E3 ubiquitin ligase MDM2 functions as a crucial negative regulator of the p53 tumor suppressor protein by antagonizing p53 transactivation activity and targeting p53 for degradation. Cellular stress activates p53 by alleviating MDM2-mediated functional inhibition, even though the molecular mechanisms of stress-induced p53 activation still remain poorly understood. Two opposing models have been proposed to describe the functional and structural role in p53 activation of Ser17 phosphorylation in the N-terminal "lid" (residues 1-24) of MDM2. Using the native chemical ligation technique, we synthesized the p53-binding domain (1-109)MDM2 and its Ser17-phosphorylated analogue (1-109)MDM2 pS17 as well as (1-109)MDM2 S17D and (25-109)MDM2, and comparatively characterized their interactions with a panel of p53-derived peptide ligands using surface plasmon resonance, fluorescence polarization, and NMR and CD spectroscopic techniques. We found that the lid is partially structured in apo-MDM2 and occludes p53 peptide binding in a ligand size-dependent manner. Binding of (1-109)MDM2 by the (15-29)p53 peptide fully displaces the lid and renders it completely disordered in the peptide-protein complex. Importantly, neither Ser17 phosphorylation nor the phospho-mimetic mutation S17D has any functional impact on p53 peptide binding to MDM2. Although Ser17 phosphorylation or its mutation to Asp contributes marginally to the stability of the lid conformation in apo-MDM2, neither modification stabilizes apo-MDM2 globally or the displaced lid locally. Our findings demonstrate that Ser17 phosphorylation is functionally neutral with respect to p53 binding, suggesting that MDM2 phosphorylation at a single site is unlikely to play a dominant role in stress-induced p53 activation.     

Systematic mutational analysis of an ubiquitin ligase (MDM2)-binding peptide: computational studies

To understand the importance of amino acids that comprise the peptide PMI (p53-MDM2/MDMX inhibitor), a p53-mimicking peptide with high affinity for the ubiquitin ligase MDM2, computational alanine scanning has been carried out using various protocols. This approach is very useful for identifying regions of a peptide that can be mutated to yield peptides that bind to their targets with higher affinities. Computational alanine scanning is a very useful technique that involves mutating each amino acid of the peptide in its complex with its target (MDM2 in the current study) to alanine, running short simulations on the mutated complex and computing the difference in interaction energies between the mutant peptides and the target protein (MDM2 in the current study) relative to the interaction energy of the original (wild-type) peptide and the target protein (MDM2 in the current study). We find that running multiple short simulations yield values of computed binding affinities (enthalpies) that are similar to those obtained from a long simulation and are well correlated with the trends in the data available from experiments that used Surface Plasmon Resonance to obtain dissociation constants. The p53-mimicking peptides contain three amino acids (F19, W23 and L26) that are major determinants of the interactions between the peptides and MDM2 and form an essential motif. We find in the current study that the trends amongst the contributions to experimental binding affinities of the hydrophobic residues F19, W23 and L26 are the best reproduced in all the computational protocols examined here. This study suggests that running such short simulations may provide a rapid method to redesign peptides to obtain high-affinity variants against a target protein. We further observe that modelling an extended conformation at the C-terminus of the helical PMI peptides, in accord with the conformation of the p53-peptide complexed to MDM2, reproduces the trends seen amongst the experimental affinities of the peptides that carry the alanine mutations at their C-termini. This suggests that some of the mutant peptides possibly interconvert between helical and extended states and can bind to MDM2 in either conformation. This novel feature, not obvious from the crystallographic data, if factored into modelling protocols, may yield novel high-affinity peptides. Our findings suggest that such protocols may enable rapid investigations of at least certain types of amino acid mutations, notably from large to small amino acids.     

Analysis of chemical shift changes reveals the binding modes of isoindolinone inhibitors of the MDM2-p53 interaction

In this study we present a method for defining the binding modes of a set of structurally related isoindolinone inhibitors of the MDM2-p53 interaction. This approach derives the location and orientation of isoindolinone binding, based on an analysis of the patterns of magnitude and direction of chemical shift perturbations for a series of inhibitors of the MDM2-p53 interaction. The MDM2-p53 complex is an attractive target for therapeutic intervention in cancer cells with intact tumor suppressor p53, as it offers the possibility of releasing p53 by blocking the MDM2-p53 binding site with a small molecule antagonist to promote apoptosis. Isoindolinones are a novel class of MDM2-antagonists of moderate affinity, which still require the development of more potent candidates for clinical applications. As the applicability of conventional structural methods to this system is limited by a number of fundamental factors, the exploitation of the information contained in chemical shift perturbations has offered a useful route to obtaining structural information to guide the development of more potent compounds. For a set of 12 structurally related isoindolinones, the data suggests 4 different orientations of binding, caused by subtle changes in the chemical structure of the inhibitors.     

Direct observations of conformational distributions of intrinsically disordered p53 peptides using UV Raman and explicit solvent simulations

We report the first experimental measurements of Ramachandran Ψ-angle distributions for intrinsically disordered peptides: the N-terminal peptide fragment of tumor suppressor p53 and its P27S mutant form. To provide atomically detailed views of the conformational distributions, we performed classical, explicit-solvent molecular dynamics simulations on the microsecond time scale. Upon binding its partner protein, MDM2, wild-type p53 peptide adopts an α-helical conformation. Mutation of Pro27 to serine results in the highest affinity yet observed for MDM2-binding of the p53 peptide. Both UV resonance Raman spectroscopy (UVRR) and simulations reveal that the P27S mutation decreases the extent of PPII helical content and increases the probability for conformations that are similar to the α-helical MDM2-bound conformation. In addition, UVRR measurements were performed on peptides that were isotopically labeled at the Leu26 residue preceding the Pro27 in order to determine the conformational distributions of Leu26 in the wild-type and mutant peptides. The UVRR and simulation results are in quantitative agreement in terms of the change in the population of non-PPII conformations involving Leu26 upon mutation of Pro27 to serine. Finally, our simulations reveal that the MDM2-bound conformation of the peptide is significantly populated in both the wild-type and mutant isolated peptide ensembles in their unbound states, suggesting that MDM2 binding of the p53 peptides may involve conformational selection.     

Use of a retroinverso p53 peptide as an inhibitor of MDM2

An N-terminal helical region of the tumor suppressor p53 binds in a hydrophobic cleft of the oncoprotein MDM2. A retroinverso isomer of the natural N-terminal helical peptide was found to interact with MDM2 using the same hydrophobic residues, Phe, Trp, and Leu. We propose that the retroinverso d-peptide adopts a right-handed helical conformation to achieve functional mimicry of the p53 peptide.     

Unveiling the “Three‐Finger Pharmacophore” Required for p53–MDM2 Inhibition by Saturation‐Transfer Difference (STD) NMR Initial Growth‐Rates Approach

Inhibitors of the p53‐MDM2 protein–protein interaction are emerging as a new and validated approach to treating cancer. Herein, we describe the synthesis and inhibitory evaluation of a series of isoquinolin‐1‐one analogues, and highlight the utility of an initial growth‐rates saturation‐transfer difference (STD) NMR approach supported by protein–ligand docking to investigate p53‐MDM2 inhibition. The approach is illustrated by the study of compound 1 , providing key insights into the binding mode of this kind of MDM2 ligands and, more importantly, readily unveiling the previously proposed three‐finger pharmacophore requirement for p53‐MDM2 inhibition.     

An integrated in silico screening strategy for identifying promising disruptors of p53-MDM2 interaction

The p53 protein, also called guardian of the genome, plays a critical role in the cell cycle regulation and apoptosis. This protein is frequently inactivated in several types of human cancer by abnormally high levels of its negative regulator, mouse double minute 2 (MDM2). As a result, restoration of p53 function by inhibiting p53-MDM2 protein–protein interaction has been pursued as a compelling strategy for cancer therapy. To date, a limited number of small-molecules have been reported as effective p53−MDM2 inhibitors. X-ray structures of MDM2 in complex with some ligands are available in Protein Data Bank and herein, these data have been exploited to efficiently identify new p53-MDM2 interaction antagonists through a hierarchical virtual screening strategy. For this purpose, the first step was aimed at compiling a focused library of 686,630 structurally suitable compounds, from PubChem database, similar to two known effective inhibitors, Nutlin-3a and DP222669. These compounds were subjected to the subsequent structure-based approaches (quantum polarized ligand docking and molecular dynamics simulation) to select potential compounds with highest binding affinity for MDM2 protein. Additionally, ligand binding energy, ADMET properties and PAINS analysis were also considered as filtering criteria for selecting the most promising drug-like molecules. On the basis of these analyses, three top-ranked hit molecules, CID_118439641, CID_60452010 and CID_3106907, were found to have acceptable pharmacokinetics properties along with superior in silico inhibitory ability towards the p53-MDM2 interaction compared to known inhibitors. Molecular docking and molecular dynamics results well confirmed the interactions of the final selected compounds with critical residues within p53 binding site on the MDM2 hydrophobic clefts with satisfactory thermodynamics stability. Consequently, the new final scaffolds identified by the presented computational approach could offer a set of guidelines for designing promising anti-cancer agents targeting p53-MDM2 interaction.     

Dual Targeting of MDM4 and FTH1 by MMRi71 for Induced Protein Degradation and p53-Independent Apoptosis in Leukemia Cells

MDM2 and MDM4 are cancer drug targets validated in multiple models for p53-based cancer therapies. The RING domains of MDM2 and non-p53-binder MDM2 splice isoforms form RING domain heterodimer polyubiquitin E3 ligases with MDM4, which regulate p53 stability in vivo and promote tumorigenesis independent of p53. Despite the importance of the MDM2 RING domain in p53 regulation and cancer development, small molecule inhibitors targeting the E3 ligase activity of MDM2-MDM4 are poorly explored. Here, we describe the synthesis and characterization of quinolinol derivatives for the identification of analogs that are capable of targeting the MDM2-MDM4 heterodimer E3 ligase and inducing apoptosis in cells. The structure-activity-relationship (SAR) study identified structural moieties critical for the inhibitory effects toward MDM2-MDM4 E3 ligase, the targeted degradation of MDM4 and FTH1 in cells, and anti-proliferation activity. Lead optimization led to the development of compound MMRi71 with improved activity. In addition to accumulating p53 proteins in wt-p53 bearing cancer cells as expected of any MDM2 inhibitors, MMRi71 effectively kills p53-null leukemia cells, an activity that conventional MDM2-p53 disrupting inhibitors lack. This study provides a prototype structure for developing MDM4/FTH1 dual-targeting inhibitors as potential cancer therapeutics.     

Proteomimetic libraries: design, synthesis, and evaluation of p53-MDM2 interaction inhibitors

The p53-MDM2 interaction regulates p53-mediated cellular responses to DNA damage, and MDM2 is overexpressed in 7% of all cancers. Structure-based computational design was applied to this system to design libraries centered on a scaffold that projects side chain functionalities with distance and angular relationships equivalent to those seen in the MDM2 interacting motif of p53. A library of 173 such compounds was synthesized using solution phase parallel chemistry. The in vitro competitive ability of the compounds to block p53 peptide binding to MDM2 was determined using a fluorescence polarization competition assay. The most active compound bound with K(d) = 12 microM, and its binding was characterized by (15)N-(1)H HSQC NMR.