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.     

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.     

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.     

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.     

Single-molecule analysis of interaction between p53TAD and MDM2 using aerolysin nanopores

Protein–protein interactions (PPIs) are regarded as important, but undruggable targets. Intrinsically disordered p53 transactivation domain (p53TAD) mediates PPI with mouse double minute 2 (MDM2), which is an attractive anticancer target for therapeutic intervention. Here, using aerolysin nanopores, we probed the p53TAD peptide/MDM2 interaction and its modulation by small-molecule PPI inhibitors or p53TAD phosphorylation. Although the p53TAD peptide showed short-lived (<100 ms) translocation, the protein complex induced the characteristic extraordinarily long-lived (0.1 s ∼ tens of min) current blockage, indicating that the MDM2 recruitment by p53TAD peptide almost fully occludes the pore. Simultaneously, the protein complex formation substantially reduced the event frequency of short-lived peptide translocation. Notably, the addition of small-molecule PPI inhibitors, Nutlin-3 and AMG232, or Thr18 phosphorylation of p53TAD peptide, were able to diminish the extraordinarily long-lived events and restore the short-lived translocation of the peptide rescued from the complex. Taken together, our results elucidate a novel mechanism of single-molecule sensing for analyzing PPIs and their inhibitors using aerolysin nanopores. This novel methodology may contribute to remarkable improvements in drug discovery targeted against undruggable PPIs.

Using aerolysin nanopores, we probed protein–protein interaction (PPI) between p53TAD and MDM2 and its modulation by small-molecule PPI inhibitors and p53TAD phosphorylation.     

Design,synthesis and structure-activity relationship of 4,5-dihydropyrrolo[3,4-c]pyrazol-6(1H)-ones as potent p53-MDM2 inhibitors

In the past decade, the p53-MDM2 protein-protein interaction by small molecules has been confirmed as a successful strategy for cancer therapy. In our previous work, pyrrolo[3,4-c]pyrazol-6(1H)-ones were found to be potent p53-MDM2 inhibitors. Further optimization and structure-activity relationship studies were described in the present work. The result revealed that benzyl group on position N1 of imidazole and bromine on C4-phenyl of pyrrolidone showed higher inhibitory activities. In vitro antiproliferative assay demonstrated the potent p53-MDM2 inhibitor 5c with 4-fold selectivity for U2 OS and Saos-2 cells. These data indicated that 4,5-dihydropyrrolo[3,4-c]pyrazol-6(1H)-one moiety is a valuable scaffold for further development of p53-MDM2 inhibitors.     

MDM2-Associated Clusterization-Triggered Emission and Apoptosis Induction Effectuated by a Theranostic Spiropolymer

Heteroatom-containing spiropolymers were constructed in a facile manner by a catalyst-free multicomponent spiropolymerization route. P1a2b as the most potent of these spiropolymers, demonstrates cluster-triggered emission resulting from strong interactions with the MDM2 protein. By preventing the anti-apoptotic p53/MDM2 interaction, P1a2b triggers apoptosis in cancerous cells, while demonstrating a good biocompatibility and non-toxicity in non-cancerous cells. The combined results from solution and cell-based cluster-triggered emission studies, docking, protein expression experiments and cytotoxicity data strongly support the MDM2-binding hypothesis and indicate a potential application as a fluorescent cancer marker as well as therapeutic for this spiropolymer.     

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.     

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.     

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.     

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.     

MDM2‐Associated Clusterization‐Triggered Emission and Apoptosis Induction Effectuated by a Theranostic Spiropolymer

Heteroatom‐containing spiropolymers were constructed in a facile manner by a catalyst‐free multicomponent spiropolymerization route. P1a2b as the most potent of these spiropolymers, demonstrates cluster‐triggered emission resulting from strong interactions with the MDM2 protein. By preventing the anti‐apoptotic p53/MDM2 interaction, P1a2b triggers apoptosis in cancerous cells, while demonstrating a good biocompatibility and non‐toxicity in non‐cancerous cells. The combined results from solution and cell‐based cluster‐triggered emission studies, docking, protein expression experiments and cytotoxicity data strongly support the MDM2‐binding hypothesis and indicate a potential application as a fluorescent cancer marker as well as therapeutic for this spiropolymer.     

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.     

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...     

Discovery of Monoacylglycerol Lipase (MAGL) Inhibitors Based on a Pharmacophore-Guided Virtual Screening Study

Monoacylglycerol lipase (MAGL) is an important enzyme of the endocannabinoid system that catalyzes the degradation of the major endocannabinoid 2-arachidonoylglycerol (2-AG). MAGL is associated with pathological conditions such as pain, inflammation and neurodegenerative diseases like Parkinson’s and Alzheimer’s disease. Furthermore, elevated levels of MAGL have been found in aggressive breast, ovarian and melanoma cancer cells. Due to its different potential therapeutic implications, MAGL is considered as a promising target for drug design and the discovery of novel small-molecule MAGL inhibitors is of great interest in the medicinal chemistry field. In this context, we developed a pharmacophore-based virtual screening protocol combined with molecular docking and molecular dynamics simulations, which showed a final hit rate of 50% validating the reliability of the in silico workflow and led to the identification of two promising and structurally different reversible MAGL inhibitors, VS1 and VS2. These ligands represent a valuable starting point for structure-based hit-optimization studies aimed at identifying new potent MAGL inhibitors.     

Monitoring the effects of antagonists on protein-protein interactions with NMR spectroscopy

We describe an NMR method that directly monitors the influence of ligands on protein-protein interactions. For a two-protein interaction complex, the size of one component should be small enough (less than ca. 15 kDa) to provide a good quality (15)N((13)C) HSQC spectrum after (15)N((13)C) labeling. The size of the second unlabeled component should be large enough so that the molecular weight of the preformed complex is larger than ca. 40 kDa. When the smaller protein binds to a larger one, broadening of NMR resonances results in the disappearance of most of its cross-peaks in the HSQC spectrum. Addition of an antagonist that can dissociate the complex would restore the HSQC spectrum of the smaller component. The method directly shows whether an antagonist releases proteins in their wild-type folded states or whether it induces their denaturation, partial unfolding, or precipitation. We illustrate the method by studying lead compounds that have recently been reported to block the MDM2-p53 interaction. Activation of p53 in tumor cells by inhibiting its interaction with MDM2 offers new strategy for cancer 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.