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.     

Turning a scorpion toxin into an antitumor miniprotein

The oncoproteins MDM2 and MDMX negatively regulate the activity and stability of the tumor suppressor protein p53 and are important molecular targets for anticancer therapy. Grafting four residues of p53 critical for MDM2/MDMX binding to the N-terminal alpha-helix of BmBKTx1, a scorpion toxin isolated from the venom of the Asian scorpion Buthus martensi Karsch, converts the miniature protein into an effective inhibitor of p53 interactions with MDM2 and MDMX. Additional mutations enable the 27-residue miniprotein inhibitor to traverse the cell membrane and selectively kill tumor cells in a p53 dependent manner.     

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.     

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

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.     

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.     

Native PAGE to study the interaction between the oncosuppressor p53 and its protein ligands

In the present study, we investigated a new approach for studying the interaction between p53 and MDM2/X (where MDM is murine double minute protein). The method is based on the different mobility between the interacting domains of the oncosuppressor p53 and its protein ligands MDM2/X on polyacrylamide gels under native conditions. While the two proteins MDM2/X alone were able to enter the gel, the formation of a binary complex between p53 and MDM2/X prevented the gel entry. The novel technique is reliable for determining the different affinity elicited by MDM2 or MDMX toward p53, and can be useful for analyzing the dissociation power exerted by other molecules on the p53–MDM2/X complex.     

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.     

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.     

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.     

Nondenaturing polyacrylamide gel electrophoresis to study the dissociation of the p53·MDM2/X complex by potentially anticancer compounds

A new analytical method to study the dissociation of the complexes between the oncosuppressor p53 and its negative modulators murine double‐minute protein 2 (MDM2) or MDMX, is proposed. This technique is reliable to determine the dissociative power exerted by small molecules on the complex taking advantage of the appearance of migrating MDM2 or MDMX in a native polyacrylamide gel, when inhibitors are added to the complex mixture. Therefore, we propose this new approach to easily screen library of compounds, with potential pharmacological anticancer activity.     

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.     

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.     

p53突变体Y220C小分子稳定剂的虚拟筛选

为了获得p53突变体的稳定剂,依次利用利宾斯基五原则,通过2次分子对接和全原子分子动力学(MD)模拟从Drug Bank 4.0数据库中筛选获得了潜在的稳定剂他克林.利用MD模拟进一步验证他克林和目标蛋白质之间的亲和作用.结果表明,他克林能够紧密结合到Y220C突变所形成的疏水空腔之中;他克林和目标蛋白质之间的主要作用力为疏水和静电相互作用,其中疏水相互作用占主导地位.此外,他克林分别与目标蛋白质的残基Leu145,Val147和Asp228形成3个氢键.基于MD模拟轨迹分析了他克林与p53CY220C的结合过程.由硫黄素T荧光光谱进一步证明他克林能够提高p53C-Y220C突变体的稳定性.     

The Mechanism of p53 Rescue by SUSP4

p53 is an important tumor‐suppressor protein deactivation of which by mdm2 results in cancers. A SUMO‐specific protease 4 (SUSP4) was shown to rescue p53 from mdm2‐mediated deactivation, but the mechanism is unknown. The discovery by NMR spectroscopy of a “p53 rescue motif” in SUSP4 that disrupts p53‐mdm2 binding is presented. This 29‐residue motif is pre‐populated with two transient helices connected by a hydrophobic linker. The helix at the C‐terminus binds to the well‐known p53‐binding pocket in mdm2 whereas the N‐terminal helix serves as an affinity enhancer. The hydrophobic linker binds to a previously unidentified hydrophobic crevice in mdm2. Overall, SUSP4 appears to use two synergizing modules, the p53 rescue motif described here and a globular‐structured SUMO‐binding catalytic domain, to stabilize p53. A p53 rescue motif peptide exhibits an anti‐tumor activity in cancer cell lines expressing wild‐type p53. A pre‐structures motif in the intrinsically disordered proteins is thus important for target recognition.     

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.     

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.     

TEAD–YAP Interaction Inhibitors and MDM2 Binders from DNA-Encoded Indole-Focused Ugi Peptidomimetics

DNA-encoded combinatorial synthesis provides efficient and dense coverage of chemical space around privileged molecular structures. The indole side chain of tryptophan plays a prominent role in key, or “hot spot”, regions of protein–protein interactions. A DNA-encoded combinatorial peptoid library was designed based on the Ugi four-component reaction by employing tryptophan-mimetic indole side chains to probe the surface of target proteins. Several peptoids were synthesized on a chemically stable hexathymidine adapter oligonucleotide “hexT”, encoded by DNA sequences, and substituted by azide-alkyne cycloaddition to yield a library of 8112 molecules. Selection experiments for the tumor-relevant proteins MDM2 and TEAD4 yielded MDM2 binders and a novel class of TEAD-YAP interaction inhibitors that perturbed the expression of a gene under the control of these Hippo pathway effectors.     

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.     

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.