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.     

Small molecule inhibitors of the MDM2-p53 interaction discovered by ensemble-based receptor models

Five nonpeptide, small-molecule inhibitors of the human MDM2-p53 interaction are presented, and each inhibitor represents a new scaffold. The most potent compound exhibited a Ki of 110 +/- 30 nM. These compounds were identified using our multiple protein structure (MPS) method which incorporates protein flexibility into a receptor-based pharmacophore model that identifies appropriate hotspots of binding. Docking the inhibitors with an induced-fit docking protocol suggested that the inhibitors mimicked the three critical binding residues of p53 (Phe19, Trp23, and Leu26). Docking also predicted a new orientation of the scaffolds that more fully fills the binding cleft, enabling the inhibitors to take advantage of additional hydrogen-bonding possibilities not explored by other small molecule inhibitors. One inhibitor in particular was proposed to probe the hydrophobic core of the protein by taking advantage of the flexibility of the binding cleft floor. These results show that the MPS technique is a promising advance for structure-based drug discovery and that the method can truly explore broad chemical space efficiently in the quest to discover potent, small-molecule inhibitors of protein-protein interactions. Our MPS technique is one of very few ensemble-based techniques to be proven through experimental verification of the discovery of new inhibitors.     

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.     

Structure-based design of potent non-peptide MDM2 inhibitors

A successful structure-based design of a class of non-peptide small-molecule MDM2 inhibitors targeting the p53-MDM2 protein-protein interaction is reported. The most potent compound 1d binds to MDM2 protein with a Ki value of 86 nM and is 18 times more potent than a natural p53 peptide (residues 16-27). Compound 1d is potent in inhibition of cell growth in LNCaP prostate cancer cells with wild-type p53 and shows only a weak activity in PC-3 prostate cancer cells with a deleted p53. Importantly, 1d has a minimal toxicity to normal prostate epithelial cells. Our studies provide a convincing example that structure-based strategy can be employed to design highly potent, non-peptide, cell-permeable, small-molecule inhibitors to target protein-protein interaction, which remains a very challenging area in chemical biology and drug design.     

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.     

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.     

Bioorthogonal Probes for the Study of MDM2‐p53 Inhibitors in Cells and Development of High‐Content Screening Assays for Drug Discovery

To study the behavior of MDM2‐p53 inhibitors in a disease‐relevant cellular model, we have developed and validated a set of bioorthogonal probes that can be fluorescently labeled in cells and used in high‐content screening assays. By using automated image analysis with single‐cell resolution, we could visualize the intracellular target binding of compounds by co‐localization and quantify target upregulation upon MDM2‐p53 inhibition in an osteosarcoma model. Additionally, we developed a high‐throughput assay to quantify target occupancy of non‐tagged MDM2‐p53 inhibitors by competition and to identify novel chemical matter. This approach could be expanded to other targets for lead discovery applications.     

Crystallographic analysis of an 8-mer p53 peptide analogue complexed with MDM2

The most potent inhibitor of the p53-MDM2 interaction reported to date is an 8-mer p53 peptide analogue (Novartis peptide), which contains 6-chlorotryptophane (Cl-Trp) and phosphonomethylphenylalanine (Pmp) as key residues for the enhanced activity. We report here a crystal structure of the co-complex between MDM2 and the Novartis peptide solved at 1.8 A resolution. The structural basis for the role of the two aromatic residues are delineated by comparing the present structure with crystal structures of the MDM2 co-complex bound to other inhibitors including the wt-p53 peptide itself.     

A computational analysis of binding modes and conformation changes of MDM2 induced by p53 and inhibitor bindings

Molecular dynamics (MD) simulations followed by principal component analysis were performed to study the conformational change of MDM2 induced by p53 and two inhibitor (P4 and MI63a) bindings. The results show that the hydrophobic cleft of MDM2 is very flexible and adaptive to different structural binding partners. The cleft tends to become wider and more stable as MDM2 binds to the three binding partners, while unbound MDM2 shows a narrower and pretty flexible cleft, which agrees with recent experimental data and theoretical studies. It was also found that the binding of P4 and p53 stabilizes the motion of the loop L2 linking the helix α2 and β strand (β3), but the presence of MI63a makes the motion of L2 disordered. In addition, the binding free energies of the three partners to MDM2 were calculated using molecular mechanics generalized Born surface area to explain the binding modes of these three partners to MDM2. This study will be helpful not only for better understanding the functional, concerted motion of MDM2, but also for the rational design of potent anticancer drugs targeting the p53–MDM2 interaction.     

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.     

Multicomponent Peptide Stapling as a Diversity‐Driven Tool for the Development of Inhibitors of Protein–Protein Interactions

Stapled peptides are chemical entities in‐between biologics and small molecules, which have proven to be the solution to high affinity protein–protein interaction antagonism, while keeping control over pharmacological performance such as stability and membrane penetration. We demonstrate that the multicomponent reaction‐based stapling is an effective strategy for the development of α‐helical peptides with highly potent dual antagonistic action of MDM2 and MDMX binding p53. Such a potent inhibitory activity of p53‐MDM2/X interactions was assessed by fluorescence polarization, microscale thermophoresis, and 2D NMR, while several cocrystal structures with MDM2 were obtained. This MCR stapling protocol proved efficient and versatile in terms of diversity generation at the staple, as evidenced by the incorporation of both exo‐ and endo‐cyclic hydrophobic moieties at the side chain cross‐linkers. The interaction of the Ugi‐staple fragments with the target protein was demonstrated by crystallography.     

A direct arylation-oxidation route to 3-arylisoindolinone inhibitors of MDM2-p53 interaction

A route to the MDM2-p53 inhibitor isoindolinone pharmacophore from a pre-formed phthalimide is detailed. The route involves treatment of 3-hydroxy-2-(n-propyl)isoindolinone with a substituted benzene in the presence of triflic acid. The resulting 3-aryl-2-(n-propyl)isoindolinones are then oxidized to the corresponding 3-hydroxy-3-aryl-2-(n-propyl) isoindolinones by treatment with 2,2’-bipyridinium chlorochromate. The benzylic oxidation represents a rather rare oxochromium (VI)-mediated reaction in which a selective C-H to C-OH transformation occurs.     

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.     

Multicomponent Peptide Stapling as a Diversity-Driven Tool for the Development of Inhibitors of Protein–Protein Interactions

Stapled peptides are chemical entities in-between biologics and small molecules, which have proven to be the solution to high affinity protein–protein interaction antagonism, while keeping control over pharmacological performance such as stability and membrane penetration. We demonstrate that the multicomponent reaction-based stapling is an effective strategy for the development of α-helical peptides with highly potent dual antagonistic action of MDM2 and MDMX binding p53. Such a potent inhibitory activity of p53-MDM2/X interactions was assessed by fluorescence polarization, microscale thermophoresis, and 2D NMR, while several cocrystal structures with MDM2 were obtained. This MCR stapling protocol proved efficient and versatile in terms of diversity generation at the staple, as evidenced by the incorporation of both exo- and endo-cyclic hydrophobic moieties at the side chain cross-linkers. The interaction of the Ugi-staple fragments with the target protein was demonstrated by crystallography.     

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.     

Discovery of spirooxindole–ferrocene hybrids as novel MDM2 inhibitors

A series of spirooxindole-ferrocene hybrids bearing five or four contiguous chiral centers were designed and synthesized via organocatalysis. In vitro protein binding and cellular proliferation assays suggested that compound 5 d was the most potent mouse double minute 2 homolog(MDM2) inhibitor. In addition,mechanistic studies indicated that compound 5 d suppressed MDM2-mediated p53 degradation, induced apoptosis and promoted oxidative damage. Molecular docking studies have suggested that 5 d binds to MDM2 by mimicking the Trp23 and Leu26 residues of p53. This work can provide a basis for the development of novel multifunctional MDM2 inhibitors. The further exploration of more derivatives from this library and additional investigation of organocatalysis application in the development of new molecules may generate new potential lead compounds for cancer-targeted therapy.     

Co-operative intra-protein structural response due to protein–protein complexation revealed through thermodynamic quantification: study of MDM2-p53 binding

The p53 protein activation protects the organism from propagation of cells with damaged DNA having oncogenic mutations. In normal cells, activity of p53 is controlled by interaction with MDM2. The well understood p53-MDM2 interaction facilitates design of ligands that could potentially disrupt or prevent the complexation owing to its emergence as an important objective for cancer therapy. However, thermodynamic quantification of the p53-peptide induced structural changes of the MDM2-protein remains an area to be explored. This study attempts to understand the conformational free energy and entropy costs due to this complex formation from the histograms of dihedral angles generated from molecular dynamics simulations. Residue-specific quantification illustrates that, hydrophobic residues of the protein contribute maximum to the conformational thermodynamic changes. Thermodynamic quantification of structural changes of the protein unfold the fact that, p53 binding provides a source of inter-element cooperativity among the protein secondary structural elements, where the highest affected structural elements (α2 and α4) found at the binding site of the protein affects faraway structural elements (β1 and Loop1) of the protein. The communication perhaps involves water mediated hydrogen bonded network formation. Further, we infer that in inhibitory F19A mutation of P53, though Phe19 is important in the recognition process, it has less prominent contribution in the stability of the complex. Collectively, this study provides vivid microscopic understanding of the interaction within the protein complex along with exploring mutation sites, which will contribute further to engineer the protein function and binding affinity.     

Bicyclic stapled peptides based on p53 as dual inhibitors for the interactions of p53 with MDM2 and MDMX

In recent years, the strategy of inhibiting the interactions of p53 with murine double minute 2(MDM2)and murine double minute X(MDMX) has been proved to be a promising approach for tumor therapy.However, the poor proteolytical stability and low intracellular delivery efficiency of peptide inhibitors limit their clinical application. Here, we designed and synthesized the bicyclic stapled peptides based on p53 by combining all-hydrocarbon stapling and lactam stapling strategies. We demonstrated th...