Molecular Dynamic Simulation Insights into the Normal State and Restoration of p53 Function
2012

Ting Fu email@example.com
Dalian Institute of Chemical Physics
Laboratory of Molecular Modeling and Design
State key Laboratory of Molecular Reaction Dynamics
Chinese Academy of Sciences
457 Zhongshan Road116023DalianChina

Department of Bioscience and Biotechnology
Dalian University of Technology
116024DalianChina

Graduate University of Chinese Academy of Sciences
100049BeijingChina

Hanyi Min
Department of Ophthalmology
Peking Union Medical College Hospital, Chinese Academy of Medical Sciences & Peking Union Medical College
100730BeijingChina

Yong Xu email@example.com
Institute of Chemical Biology
Guangzhou Institutes of Biomedicine and Health
Chinese Academy of Sciences
190 Kaiyuan Avenue510530GuangzhouChina

Jianzhong Chen firstname.lastname@example.com
Department of Mathematics and Physics
Shandong Jiaotong University
250031JinanChina

Guohui Li
Dalian Institute of Chemical Physics
Laboratory of Molecular Modeling and Design
State key Laboratory of Molecular Reaction Dynamics
Chinese Academy of Sciences
457 Zhongshan Road116023DalianChina

Molecular Dynamic Simulation Insights into the Normal State and Restoration of p53 Function

Int. J. Mol. Sci
13201210.3390/ijms13089709Received: 11 May 2012; in revised form: 6 July 2012 / Accepted: 11 July 2012 /International Journal of Molecular Sciences Review † These authors contributed equally to this work. * Author to whom correspondence should be addressed; E-Mail: email@example.com;p53MDM2MDMXmolecular dynamic simulationconformational changeprotein-protein interactionprotein-ligand interaction
As a tumor suppressor protein, p53 plays a crucial role in the cell cycle and in cancer prevention. Almost 50 percent of all human malignant tumors are closely related to a deletion or mutation in p53. The activity of p53 is inhibited by over-active celluar antagonists, especially by the over-expression of the negative regulators MDM2 and MDMX. Protein-protein interactions, or post-translational modifications of the C-terminal negative regulatory domain of p53, also regulate its tumor suppressor activity. Restoration of p53 function through peptide and small molecular inhibitors has become a promising strategy for novel anti-cancer drug design and development. Molecular dynamics simulations have been extensively applied to investigate the conformation changes of p53OPEN ACCESSInt. J. Mol. Sci. 2012, 13    9710 induced by protein-protein interactions and protein-ligand interactions, including peptide and small molecular inhibitors. This review focuses on the latest MD simulation research, to provide an overview of the current understanding of interactions between p53 and its partners at an atomic level.

Introduction

p53 is a tumor suppressor protein, encoded by the P53 gene initially reported in 1979 [1]. p53 is commonly referred to as "the guardian of the genome", as it plays a vitally important role in multi-cellular organisms, where it regulates the cell cycle (promoting cell apoptosis, maintaining genomic stability and inhibiting tumor angiogenesis) and functions as a tumor suppressor involved in slowing or monitoring the cell division and preventing genome mutation in normal circumstances [2][3][4][5]. When cells are damaged, p53 determines the degree of DNA damage and either promotes cellular repair when the degree of damage is not severe, and cell apoptosis when the damage is too severe to be repaired. Damaged cells with a P53 deficiency do not have this functionality and will continue to begin uncontrolled division even in adverse conditions. Under normal circumstances, p53 is strictly controlled at minimal levels through a fast continuous degradation mechanism with a half-life of only about 20 min. In abnormal conditions, however, p53 with extended half-life is stabilized and activated in response to stressful stimuli such as DNA damage, ionizing radiation, UV-radiation, heat shock, nitric oxide, hypoxia and oncogene activation [6][7][8].

p53 is a homo-tetramer, with each flexible monomer containing 393 residues that can be divided into several regions, as follows [9][10][11]: The N-terminal transcriptional activation domain (TAD) includes an acidic activation domain 1 (AD1, residues: 1-42) and domain 2 (AD2, residues: 43-62). This region is a highly charged area including a large number of acidic amino acid residues and phosphorylation sites, and is involved in the regulation of apoptotic activity. The adjacent proline-rich domain at residues 63-97, is essential for the apoptotic activity of p53. The core DNA-binding domain (DBD, residues: 102-292), contains one zinc atom and several electropositive arginine amino acids which interact with DNA. This region can bind to sequence-specific DNA, and is often unstable and contains most of the inactivating single missense mutations occurring in human tumors [12]. The nuclear localization signaling domain at residues 316-325, is involved in the intracellular localization of p53. The oligomerization domain (OD, residues: 326-356) is responsible for tetramerization which is essential for the activity of p53 in vivo. The intrinsically disordered C-terminal regulatory domain (CTD), a flexible region is involved in down-regulation of the central DNA binding domain at residues 363-393. The three-dimensional structures of the entire DNA-binding domain [13][14][15] (PDB ID: 3KMD, Figure 1B) [15] and oligomerization domain (PDB ID: 1PES, Figure 1C) [16] have been solved while the structures of the N-terminal transcriptional activation domain and C-terminal regulatory domain have not been determined ( Figure 1). Bell et al. [17] argued that full-length wild-type p53 protein contains large unstructured regions in its N-and C-terminal domains, is destabilized and easy to unfold and loses its biological activities in the absence of modifications or stabilizing partners. The three-dimensional structures of p53 TAD fragment bound to MDM2 (PDB ID: 1YCR, Figure 1A) [18] and p53 CTD fragment bound to S100 calcium-binding protein B (PDB 1DT7, Figure 1D) [19] are shown in Figure 1. All the figures were created with Pymol [20].  [18] is shown in cartoon, p53 TAD fragment (residues [17][18][19][20][21][22][23][24][25][26][27][28][29] is shown in magenta and the three most important residues are shown in stick, MDM2 (residues 25-109) is shown in green; (B) The tetramer of the DBD of p53 (PDB 3KMD) [15] is shown in cartoon and the four monomers (residues 92-291) are colored in green, cyan, magenta and yellow, respectively; Zn 2+ is shown in sphere and dirtyviolet, and the DNA is shown in stick; (C) The tetramer of oligomerization domain (OD) of p53 (PDB 1PES) [16] is shown in cartoon and the four monomers (residues 325-355) are colored in green, cyan, magenta and yellow, respectively; (D) The complex of p53 C-terminal regulatory domain (CTD) fragment bound to S100 calcium-binding protein B (PDB 1DT7) [19] is shown in cartoon, p53 CTD fragment (residues 377-387) is shown in magenta and yellow, S100B (residues 1-91) is shown in green and cyan and the two Ca 2+ are shown in sphere and are colored in, consistent with the S100B protein for the two subunits, respectively. Figures were created with Pymol (http://pymol.org) [20].

It is clear that the stability and transcriptional activity of p53 are regulated through a complex cascade of post-translational modifications, such as phosphorylation (the 17 known phosphorylation sites in human p53 are Ser6, Ser9, Ser15, Thr18, Ser20, Ser33, Ser37, Ser46, Thr55, Thr81, Ser149, Ser150, Ser155, Ser315, Ser376, Ser378 and Ser392), and acetylation of critical lysines (AcLys382), methylation (MeLys382) and ubiquitination [21][22][23][24]. Furthermore, the destabilized structure may allow the physiological interaction of p53 with numerous protein partners and regulation of its turnover [14]. Many biological, structural, mutagenesis and computational studies showed that the pro-apoptotic activity of p53 is complicated, and affected by protein-protein interactions [25,26]. For example, the TAD fragment of p53 involving residues 12-26, has high probability of forming a short α-helix that is capable of interacting with protein partners, such as the transformed mouse 3T3 cell double minute 2 (MDM2, or HDM2 for the human congener, PDB ID: 1YCR, Figure 1A) [18] and MDM2-related protein (MDMX, also named MDM4) [27]. As a negative regulator, MDM2/X can induce inactivation of over-expressed p53 in a normal cell. In addition to the key regulators MDM2 and MDMX which interact with the target p53 through TAD, some other partners have been found in recent years. Bcl-XL, one member of the Bcl-2 family proteins, is identified as a binding target of p53 via TAD and results in transcription-independent apoptotic activity [28][29][30]. Azurin, a copper-containing protein with electron transfer activity, has been reported to bind p53 via either the TAD or the DBD domains of p53 [31][32][33]. The single-stranded DNA-binding protein, replication protein A (RPA) (PDB ID: 2G3B) [34] and the RNA polymerase II transcription factor B subunit 1 are also found to interact with p53 TAD (PDB ID: 2GS0) [35]. The DBD of p53 is mainly responsible for sequence-specific DNA binding (PDB ID: 3KMD, Figure 1B) [15] and some protein-protein interactions. The large T-antigen of Simian Virus 40 binds to DBD and induces the dramatic conformational changes at the DBD of p53 (PDB ID: 2H1L) [36]. Moreover, the extreme CTD not only binds to DNA and RNA sequences , but also is critical for regulation of p53 function [37] and is capable of adopting multiple folded conformations upon binding to different partners such as S100 calcium-binding protein B (S100B) (PDB ID: 1DT7, Figure 1D) [19], sirtuin protein (Sir2) (PDB ID: 1MA3) [38], cAMP response element-binding (CREB) binding protein (CBP) (PDB ID: 1JSP) [39], the histone methyltransferase SET9 (also known as SET7/9) (PDB ID: 1XQH) [40] and the cyclin A/cyclin-dependent protein kinase 2 complex (PDB ID: 1H26) [41].

Mutation of TP53 is the most common genetic change in human cancers, which mainly occurs in the highly conserved residues, particularly Arg175, Gly245, Arg248, Arg249, Arg273 and Arg282 hot spots. More than 50 percent of all the human malignant tumors, for instance liver cancer, breast cancer, bladder cancer, gastric cancer, colon cancer, prostate cancer, soft tissue sarcoma, ovarian cancer, brain tumor, lymphatic tumors, esophageal cancer, lung cancer, osteogenic sarcoma, are the result of a mutation or deletion of the TP53 gene [42][43][44]. Due to a conformational change, the mutational TP53 gene loses regulatory control of cell growth, apoptosis and DNA repair. For the other half of human cancers, p53 retains its wild-type form but the activity is greatly reduced due to protein-protein interactions with key regulatory partners, such as MDM2 and MDMX. Therefore, restoring endogenous p53 functions by the disruption of p53-MDM2/X interaction using non-peptide small molecule inhibitors [45,46] or peptide [47][48][49] currently plays one of the most promising strategies for anti-cancer drug design and development. Several related reviews [46,50-58] have reported, a large number of potent and selective inhibitors that have been identified and can be divided into two categories ( Figure 2). One type is the non-peptide small molecule inhibitors, which include but are not limited to, Spiro-oxindoles, piperazine-4-phenyl derivatives, chlorofusin, Cis-imidazolines (nutlins, Figure 2A) [59], imidazole-indoles (K23, Figure 2B) [45], norbornanes, benzodiazepinediones (DIZ, Figure 2C) [60], chromenotriazolopyrimidines, terphenyls, isoindolinones, sulfonamides, chalcones and boronic chalcones. All the above mentioned small molecule inhibitors are used to fill the Trp23 sub-pocket with either a p-halogenated phenyl group or a combination of 6-chloro(ox)indole and the Phe19 and Leu26 sub-pockets with multiple hydrophobic interactions. Another type includes peptide inhibitors such as p53-like peptides (PMI, Figure 2D) [49], miniproteins and peptidominetic inhibitors, which mimic the identical hot-spot triad of Phe19, Trp23 and Leu26 of p53. To avoid proteolytic degradation, a series of stabilized alpha-helix of p53 (SAH-p53) cyclic peptides are generated by adopting the "peptide stapling" strategy ( Figure 2E) [61]. In addition, S100B has been confirmed to interact with the CTD of p53 protein and inhibit p53 aggregation and down-regulate p53 protein levels and activity [62,63]. Blocking the S100B-p53 interaction and protecting p53 from S100B-dependent down-regulation by ligands is another way to restore the p53 tumor suppressor function [64]. Several other proteins bind to or modify the p53 CTD, for instance, Sir2 specifically binds and deacetylates p53 protein, and activate the tumor suppressor [38]. , the parts imitating the three most important residues of p53 (F19′, W23′ and L26′) are presented as stick; (E) MDM2 with SAH-p53-8 stapled-peptide (from PDB 3V3B) [61], the all-hydrocarbon staple is shown in stick and magenta. Figures were created with Pymol (http://pymol.org) [20].

Computational methods, including molecular docking, molecular dynamics (MD) simulation and binding free energy prediction are usually applied to supplement the experimental studies in order to understand the mechanism at an atomic level or to improve the efficiency of molecular design or discovery. Over the past decade, there has been a dramatic increase in the number of theoretical studies of p53, along with the experimental data and structural information. Previous reviews [50][51][52][53][54][55][56][57] primarily have been concerned with the interactions of p53-MDM2/X and MDM2/X-inhibitors based on the three-dimensional structures available. Many computational studies were carried out to understand p53's intrinsic flexibility, the molecular mechanism of interactions between p53 and partners, and the conformation changes of p53 induced by protein-protein interactions and protein-ligand interactions, including peptide and small molecular inhibitors. More than 40 articles discussing the dynamic simulations of p53 have been published in the latest three years. The dynamics, flexibility and ligand-induced conformational changes of p53 have not been completely summarized. Additionally, the dynamics characteristic of other domains of p53 has not been reviewed. In this review, the latest MD simulations discussing the protein-protein interactions and the flexibility of native p53 in TAD, DBD, OD and CTD domains, and the inhibitors (small molecule and peptide inhibitors) are discussed.

Advance in Molecular Simulations

Molecular Docking, MD Simulation and Binding Free Energy Prediction

Molecular docking, a useful computational method for rational drug design, can predict the preferred binding orientation and conformations of small molecule, and the strength of association or binding affinity of the complex through scoring functions. Docking software includes AutoDock, Dock, 3D-Dock, Affinity, LigandFit, FRED, Surflex, HEX, Glide, GOLD, ICM, and MVD as examples. Molecular Dynamics (MD) is a computer simulation method that relies mainly on Newton mechanics to simulate the physical interactions and movement of atoms and molecular systems. To obtain the dynamic characteristics and the understanding of interaction mechanism at an atomistic scale, MD simulation packages like AMBER, CHARMM, GROMACS and NAMD are applied widely in the modeling of biomolecules. Binding free energy prediction has been regarded as a powerful and valuable tool to explore the binding mechanisms and binding affinity of inhibitors to proteins. The molecular mechanics Possion-Boltzmann/generalized Born surface area (MM-PB/GBSA) method has been used successfully in explaining protein-protein and protein-inhibitor interactions [65]. Several computational strategies have been adopted to get a deep insight into the interactions between p53 and its partners.

Analysis of p53-Protein Interaction

Protein-Protein Interaction Located in TAD of p53

As the principal negative protein regulators of p53, MDM2 binds a short amphipathic α-helix to the TAD domain of p53 (residues 17-29) (Figure 3) [18], and inhibits the transcriptional activity of p53. MDM2 acts as an E3 ubiquitin ligase and is able to target p53 for proteasomal degradation. Since MDM2 is the target gene of p53, the p53-MDM2 system forms a negative feedback loop, which precisely regulates and controls the protein level and activity of p53 [66][67][68][69]. MDMX, highly homologous to MDM2 but lacking ubiquitin ligase activity, negatively regulates the transactivation function of p53 via interaction with the p53 TAD domain [70,71]. Studies have shown that the over-expression or hyper-activation of MDM2 and MDMX directly leads to the negative consequences of loss of p53 normal function during the development of nearly half of human cancers [72]. Expressional levels of MDM2/X and the interactions between p53-MDM2/X are important for p53 functionality. The above evidence illustrates that p53 TAD is one of the most pivotal regions in p53 function.

Small lipophilic molecules or p53 peptide analogues, which disrupt or prevent the p53-MDM2/X interaction and restore the p53 tumor suppression, have become an important strategy for anti-cancer drug design and discovery [50,51,73,74]. The stability of p53 TAD fragment and conformational changes due to protein-protein interactions, in particular between p53 and MDM2/X, and protein-ligand interactions, including peptide and small molecular inhibitors have been most extensively studied [58, [75][76][77][78][79]. The three-dimensional structure of p53-MDM2 (PDB ID: 1YCR) ( Figures 1A and 3) [18] reveals that the peptide segment (Glu17-Asn29) of p53 TAD as an amphipathic α-helix binds to the deep hydrophobic cleft, which is formed by the N-terminal domain of MDM2 (Glu25-Val109) containing of two pairs of α-helices and an anti-parallel β-sheets. The four key residues Phe19, Leu22, Trp23 and Leu26 of p53, and at least residues Leu54/Met53, Leu57/Leu56, Gly58/Gly57, Ile61/Ile60, Met62/Met61, Tyr67/Tyr66, Gln72/Gln71, His73/His72, Val75/Val74, Phe91/Phe90, Val93/Val92, His96/Pro95, Ile99/Leu98 and Tyr100/Tyr99 of MDM2/X play the dominant roles for the p53-MDM2/X interactions. In particular, three key residues Phe19, Trp23, and Leu26 of p53 insert deeply into the hydrophobic cleft of MDM2. Similarly, the crystal structure of the N-terminal domain of human MDMX (Gln23-Ile111) bound to a 15-residue p53 TAD fragment (Ser15-Asn29) was solved at 1.90 Å (PDB ID: 3DAB) [28] (Figure 3). Although the sequence is highly conserved, and the structure and binding mode are very analogous, the MDMX binding pocket is differently shaped and smaller than that of MDM2, due to the side chains of Met53 and Tyr99 (corresponds to Leu46 and Tyr100 in MDM2) protruding into the hydrophobic cleft of MDMX. Analysis of the dynamics reveals that the most flexible region of MDM2 is the hydrophobic cleft, which becomes more wider and stable when bound to p53 while narrower and highly flexible in the absence of p53 [80]. Moreover, based on available structural data of the MDM2/X protein in complex with ligand, including peptide and small molecular inhibitors [45,47,48,59,60], it is possible to understand the detailed molecular interactions between p53 and MDM2/X.

It is well known that p53 TAD in the absence of partners in solution lacks persistent structural order, which is similar to a classical random coil except for a small helix, and contains numerous phosphorylation sites. This intrinsically disordered structural characteristic enables it to interact with multiple protein partners. To explore the native conformational properties, Espinoza-Fonseca and co-workers [76] carried out all-atom MD simulations of the extended p53 TAD fragment (Met1-Ala39) in solution and found that the leucine-rich hydrophobic clusters play an important role in the formation and stabilization of the p53 TAD fragment (Met1-Ala39) globule-like folded structure. Terakawa et al. [81] tested diverse structural ensembles of p53 TAD by a systematic multiscale computational method. They used all-atom replica-exchange MD simulations of a coarse-grained model to generate the structural ensemble, and their results showed that the p53 TAD fragment was preferable and formed a kinked structure in the AD1 region and extended conformation in a proline-rich region. Xiong et al. [82] experimentally measured the distributions of Ramachandran Ψ angles for this intrinsically disordered p53 TAD and P27S mutation peptide in unbound state using UV resonance Raman spectroscopy, and they also performed explicit solvent simulations to provide detailed views of these conformational distributions. Both experiment and simulations revealed that the number of PPII-like helical conformation content was decreased and the non-PPII content was increased upon P27S mutation. It turns out that binding of MDM2 induces a helical pattern in the fragment Phe19-Leu25 of p53 TAD. Espinoza-Fonseca et al. [83] confirmed that the fragment Phe19-Leu22 displayed a stable, transient helical pattern at sub-microsecond periods. Another MD simulation was carried out by Mavinahalli et al. [84] using an implicit solvent method to explore the folding ability of the TAD domains of p53, p63 and p73. The binding region folds into a stable helix when p53 TAD peptide is embedded into MDM2. Residues Phe19, Leu22, Trp23, Leu25 and Leu26 constitute a hydrophobic cluster, which controls the kinetics of helix formation. An ionic interaction between Asp20 and Lys24 along with the other three hydrogen bonds between Leu22 and Glu17 (42% occupancy), Phe19-Trp23 (29% occupancy) and Ser20-Lys24 (21% occupancy) are responsible for the stability of the helix. It is known that the helical conformation of the p53 TAD is necessary for optimal binding of p53 to MDM2. It was once believed that [85] the loss of a key hbond between Thr18 and Asp21 would lead to the destabilized helical structure and reduced affinity for MDM2. Mavinahalli et al. [84] shows that phosphorylation of Thr18 results in perturbation of this helicity, whereas phosphorylation of Ser20 actually enhances the helical propensity of p53, and the electrostatic interaction of the phosphate and Lys24 play a vital role in the thermodynamic stability of the helix. MD simulations [84], in accord with the earlier findings [86], suggest that upon phosphorylation of either Thr18 or Ser20, its helical structure is still retained but the affinities for MDM2 are reduced. The computational interpretation [86] is that though the important hbond between Thr18 and Asp21 is broken, there are other alternate patterns of hbonds formed to maintain the helicity; a charge-charge repulsion between the negatively charged phosphate group on Thr18 and the anionic patches in the region of the MDM2 surface, result in the weakening p53-MDM2 interactions. This idea has been confirmed by experiment study [87].

From 1999 to the present, p53-MDM2 interactions have been extensively studied using MD simulations. In recent years researchers are paying more attentions to the binding free energy calculations by using the MM-PB/GBSA method, and estimations are in excellent agreement with the experimental results. This has been well reviewed by Lauria et al. [79] and will not be discussed in the present paper.

Although MDM2 and MDMX proteins are highly homologous, they display different dynamic behavior. More recently, Dastidar et al. [88] emphasized the role of the dynamics of Tyr100, which was recognized as a putative gatekeeper in p53-MDM2 interactions by MD simulations. It appears that the transition of the Tyr100 side chain from "closed" to "open" enlarges the binding site, thus enhancing the binding of MDM2. It is well known that electrostatic complementarity plays a crucial role in the complex formation. The p53 TAD fragment (Glu17-Asn29) has a strong anionic potential with a net -2e charge, in contrast to the largely cationic potential in MDM2 with net charge of +5e while MDMX with a net charge of +1e. The binding affinity of p53-MDMX is much lower than that of p53-MDM2. Then, the same author [89] examined in full atomic detail of the role of electrostatic interaction on the formation of the p53-MDM2 complex by varying the distance of separation between MDM2 and p53. This study reveals that electrostatic interaction controls p53-MDM2 complex formation at long range, while van der Waals interactions of Phe19, Trp23, Leu26 of p53 determine the complex formation at short range.

In order to investigate the conformational events of free and p53-bound states of MDM2 and MDMX, MD simulations were performed by Carotti et al. [90]. The analysis of the trajectories shows that the diverse conformations of Tyr99/Tyr100 (MDMX/MDM2) hold crucial roles in affecting the p53-MDM2 and p53-MDMX recognition, and different interactions in the apo and p53-bound states might affect the shape of the p53 binding cleft. The van der Waals interaction was found to be favorable in the p53-MDM2 simulation while the electrostatic interactions mainly dominate in the p53-MDMX simulation. However, the gain in electrostatic energy could not offset the loss of solvation energy, van der Waals interaction between both proteins MDM2/X and Phe19 and Trp23 of p53 was identified as the dominate event in the peptide recognition.

A systematic analysis of single mutations in the 12-residue fragment of the p53 TAD (Q16ETFSDLWKLLP27) and an equivalent phage optimized peptide (12/1, M1PRFMDYWEGLN12) were undertaken by Brown et al. [91]. The 10 complexes between MDM2 and the wild type p53 TAD and mutants P27S, P27A, P27T and P27N, and wild type 12/1 and N12S, N12A, N12T, N12P, were studied by using MD simulations to elucidate the differences in the mechanism and thermodynamics of their interactions. The analysis revealed that single residue substitutions dramatically increase the affinity of the MDM2 binding peptides. In addition, the modification of Ser17 in the N-terminal "lid" of MDM2 is expected to regulate the p53-MDM2 binding, such as S17D can promote the binding of p53 to MDM2. A similar study was performed by Dastidar et al. [92], where they carried out detailed MD simulations on the p53 binding domain of wild type (WT) MDM2 (residues 1-119), S17D mutant and the phosphorylated Ser17 (pSer17) with the lid in "open" and "closed" conformations to evaluate the effects of Ser17 modification on p53 binding.

Recently, the recognition process between TAD fragment of p53and Bcl-X L , has been studied by Bharatham et al. [93], finding that the residues Glu96, Phe97, Arg100, Phe105, Leu108, Val126, Leu130, Ala142 of Bcl-X L play important roles in the interaction with p53 TAD, in which the three conserved hydrophobic residues (Phe19, Trp23 and Leu26) were found to be anchored into three hydrophobic pockets of Bcl-X L .

It has been observed that the copper-containing protein azurin (AZ), with electron transfer activity in Pseudomonas aeruginosa, binds p53 protein via either the TAD or the DBD and enhances the structural stability of p53 to activate its tumor suppressor function [31][32][33]. To explore the molecular interaction mechanism of the p53-AZ complex, Taranta et al. [94] modeled the complex structure of TAD peptide of p53 (residues 17-56) and the bacterial protein AZ, performing MD simulations followed by binding free energy estimation, and validating the predicted complex structure by mutagenesis studies. The short-range hydrophobic interactions and a lot of hydrogen bonds were proposed to be the main driving force for the stabilization of the complex.

These results shed further light on different features of recognitions between p53 TAD and protein partners and maybe helpful for rational design and development of new potent and selective inhibitors which can restore the function of p53.

Protein-Protein Interactions Located in DBD and OD of p53

The three-dimensional structures of the entire DNA-binding domain [15] and oligomerization domain [16] have been determined ( Figure 1B,C). Several studies found that the metal cation located in the DBD is important for aggregation and sequence-specific DNA binding [95,96]. DBD is often unstable and contains most of the inactivating single missense mutations occurring in human tumors [12].

The protein Azurin (AZ) has been demonstrated to enhance p53 stability via binding to not only TAD but also the p53 DBD. Computational methods including molecular docking, MD simulations and free binding energy estimation by MM-PBSA, were carried out to understand the interaction between AZ and p53 DBD [97]. Recent research by Santini et al. [98] focused on the interaction between p28 (the peptide fragment of Azurin, residues 50 to 77) and the p53 DBD peptide using computational docking coupled with MD simulations and binding free energy estimations. The p28 fragment should retain the cellular penetration ability and anti-proliferative activity of the whole protein. Both the folded and the unfolded structures of the p28 peptide have been found to form a stable complex with the p53 DBD by favorable low binding free energy, high shape complementarity, predominant polar interfaces and hydrogen bond interactions.

Tetramerization of the p53 protein is essential for its function, such as for DNA binding, protein-protein interactions, post-translational modifications, and p53 degradation [99]. The small oligomerization domain (OD) comprising one β-strand, one tight turn and one α-helix ( Figure 1C), is responsible for the formation of the stable native tetramer, which is generally formed by a dimer (consisting of an anti-parallel β-sheet and an anti-parallel helix-helix interface) of dimers [16,100,101].

Dimerization of the p53 OD and unfolding of the tetramerization domain have been studied by running MD simulations [102,103], revealing that formation of the folding nucleus precedes and appears to be necessary for the formation of hydrogen bonds in the β-sheet as well as other inter-monomer contacts. The tetramerization is stabilized by hydrophobic interactions while disruption of the native tetrameric hydrophobic core results in the beginning of unfolding.

Moreover, single mutation R337H can destabilize the tetramer and inactivate its function [104,105]. A strong salt bridge formed by Arg337 and Asp352 from another monomer at four equivalent sites, is very important to the overall tetramer stability [100]. Interestingly, a fluid salt-bridging cluster of Arg337, Arg333, Glu349 and Asp352 was observed during the MD simulations of wild type p53 and the R337H mutant at several different pH and salt conditions [106], suggesting that these four residues contribute significantly to the tetramer stability. The mutation R337H still interacts with Asp352, but loses much of its interaction with Glu349. When His becomes deprotonated at high pH, the interaction between this deprotonated residue and Asp352 is also lost, which was considered the reason for strong destabilization of mutation R337H observed at high pH.

Additionally, published observations [107][108][109][110][111][112] suggest that S100B can binds to the tetramerization domain of p53 when it is exposed in lower oligomerization states, and S100B not only inhibits p53 tetramer formation but also promotes disassembly of the p53 oligomers. These findings imply that reduced the levels of S100B in tumors could restore p53 functions and contribute to cancer prevention. However, up till now its molecule mechanism of the interactions between S100B and the p53 OD has not been clearly understood and will require further research.

Protein-Protein Interactions Located in CTD of p53

The C-terminal region is subject to extensive post-translational modifications including acetylation, methylation, phosphorylation and ubiquitination that regulate function and cellular protein levels.

Unlike other proteins that bind to or modify the CTD of p53, proteins that activate the tumor suppressor, inactivation was observed for S100B, a member of the best characterized proteins of the S100 protein family involved in cell cycle progression and cell differentiation [62,63,[107][108][109][110][111][112]. Researchers have reported that the tumor suppressor protein p53 interacts in a calcium-dependent manner with dimeric S100B [107], and consequently the PKC-dependent phosphorylation of p53 in its CTD is suppressed by protein-protein interactions [62,63,111,112]. Interestingly, S100B not only binds to the p53 tetramerization domain and controls its oligomerization state, but also binds to the p53 CTD, which is unstructured in its native form, as revealed by nuclear magnetic resonance (NMR) spectroscopy. It forms an α-helix secondary structure that spans residues 377-387 when binding to the Ca 2+ -bound S100B(betabeta) partner (PDB ID: 1DT7) ( Figure 1D) [19]. Two Ca 2+ at the loop regions (residues 20-25, and 61-66) are bound to each monomer of the homo-dimeric S100B in solution and trigger the binding of S100B with p53. It is believed that the p53 CTD recognizes protein S100B through a "fly-casting"-like process [113]. The S100B sterically hinders sites of phosphorylation (at positions Ser376 and Thr377) and acetylation (at position Lys382) on p53 which are crucial for the activation of transcription [19,114].

The hydrophobic surface of S100B is mainly formed by the helix 2 (residues 29-39), the hinge region (residues 40-49), the helix 3 (residues 50-60), and the helix 4 (residues 70-88) ( Figure 1D). The comparative structural studies reveal that S100B undergoes a conformational change upon binding Ca 2+ and generates the hydrophobic surface through changing the orientation of helix 3 and helix 4 and placing helix 3 nearly perpendicular to helix 4 and elongating the helix 4 by one turn. This conformational change results in the formation of several favorable hydrophobic and electrostatic interactions and promotes its interaction with the target p53 protein as well as disassembly of the p53 oligomerization [19,62,107,111,115]. Since p53 protein functions as a tumor suppressor only in its tetrameric form, the binding of S100B disrupts the tetramerization equilibrium of p53 and thus closely relates to the genesis and progression of malignant tumors [116][117][118]. Therefore, it is necessary to develop high-efficient and low-toxic ligands to block the p53-S100B interaction in order to protect p53 from S100B-dependent down-regulation and restore the tumor suppressor function. Several theoretical studies were performed to understand the p53-S100B interaction [119][120][121]. In this section, latest molecular dynamic simulation studies are discussed.

Recently, the binding process between the p53 CTD and five partners including the S100B was investigated using MD simulation by Allen et al. [120]. Each complex in addition to the p53 CTD fragment alone and the S100B alone were simulated in triplicate for 100 ns with GROMACS. When bound to the S100B, the α-helix structure of the p53 CTD comprising residues 380-387 underwent little change from the starting structure in the first and third replicates during each 100 ns simulation while the α helix was found to form a turn in the second replicate. When unbound to the S100B, the secondary structure of p53 CTD changes from α-helix to other forms, most notably β-sheets and turns. The RMSD of the S100B and p53 CTD fragment fluctuates less in bound state, especial for the p53 CTD compared with its free state in solution. The hinge region (residues 40-49), helix 3 and helix 4, shaping the p53 CTD binding site, all fluctuate less but the hinge region located in calcium binding site (residues 20-25) fluctuated more for the p53 CTD fragment-bound state. The observation was in agreement with the structural research [19] that S100B underwent a large conformational change in helix 3 upon binding Ca 2+ and this change was advantageous to the opening of the large hydrophobic core, including S100B residues Met79, Val80, Leu44, and Val56, which were necessary for the binding of p53 CTD. The configurational principal components analysis (PCA) which describes large collective motions, also confirmed a stabilization effect of S100B on the p53 CTD fragment upon binding to the S100B. The electropositive p53 CTD residue Arg379 forms stable hydrogen bonds with His42 and Glu45, so does Lys382 with Glu86 and Phe88. Other residues His368, Ser378, and Lys381 also form many favorable electrostatic contacts with S100B residues Glu45, Glu46, Lys48, Glu49, Glu51, Lys55, and Glu89. In addition, the van der Waals contacts formed by Arg379, Lys382, Leu383, Met384 and Phe385 of the p53 CTD fragment and S100B hydrophobic core defined by residues Leu44, Val52, Val56, Met79, Val80, and Phe87 stabilize the protein-protein complex. The author Allen [120] offered a hypothesis about the specific α-helix conformation adopted by the p53 CTD peptide when bound to S100B, which enables alignment of the hydrophobic residues (Leu383, Met384, and Phe385) toward the hydrophobic binding site. The Ser376 and Thr377 of the p53 CTD were believed to be buried against the surface of S100B, and the MD simulations confirmed that Thr377 tends to stay buried and forms hydrogen bonds with residues Glu45, Glu46, and Glu49, but Ser376 becomes solvent exposed.

To investigate the structural effects of a novel single point mutation (G389E), MD simulations of the free wild type (WT) and mutant (G389E) peptides corresponding to the extreme CTD, in solution and in complex with S100B were carried out in 2010 by Pirolli and co-workers [121]. In these studies, the helical stability of the CTD fragment of p53 was analyzed. The p53 α-helix G389E mutant was found to be unfolded rapidly in the absence of S100B in solution and WT peptide gradually loses its helical conformation which confirmed that the free p53 CTD exists in a disordered structure. The G389E mutation has negligible effect on the phosphorylation (residues Ser376, Thr377) and acetylation (residue Lys382) sites when it is free in the solution. The p53 CTD peptide is very flexible when bound to S100B, especially in the G389E complex. The mutant residue increases its α-helix content, while S100B is stable in both simulations. The α-helix of WT p53 was found to be gradually lost and undergo a large conformational change and form a closed peptide at the end of the simulation. Snapshots from simulations showed that the negatively charged residue Asp391 and Asp393 alternatively interact with positively charged residues Lys370, Lys372 and Lys373. The persistent direct salt bridge between Asp391 and Lys373 and one hydrogen bond between the Glu388 and Thr377 contributed to this conformational change. In the presence of this novel G389E mutation, the binding interface is significantly altered, and Lys382 and Lys370 form a direct salt bridge with Thr82 and Glu89 located in helix 4 of S100B, respectively, thus the normal function of p53 is prevented.

It has been known the acetylation of Lysine of p53 is necessary to destabilize the p53-MDM2 interaction and facilitate p53-mediated stress response [122]. The activity and stability of p53 is also regulated by the lysine acetyltransferase and transcriptional coactivator cAMP response element-binding (CREB)-binding protein (CBP)/p300 [123][124][125]. Modulation of acetylation-mediated interactions by small inhibitors may be a new approach for anti-cancer drug design and research [125][126][127][128]. NMR structural analysis reveals that bromo-domain of CBP (residues 1081-1196) binds specifically to the acetylated lysine 382 in the CTD (residues 367-386) [39]. The p53-CBP interaction was investigated by MD simulations and the energetics of the binding complex was also estimated by Free Energy Perturbation (FEP) [120,129].

Allen et al. [120] observed that a moderate change of bend conformation occurred along the middle of the chain of p53 CTD fragment in the complex during MD simulation. Drastic changes in the N-and C-terminal ends cause them to come together as β-strands to form a stable and persisted anti-parallel β-sheet. These occurred after 20 ns simulation for the free CTD fragment in solution. On the other side, the opening and closing of the hydrophobic binding site between the ZA loop (residues 1115-1138) and the BC loop (residues 1170-1190) in the bromo-domain represented 23.8% when bound to p53 CTD and up to 58.4% of the total motion for CBP protein in the absence of the p53 CTD according to the PCs analysis. In the initial structure, only six residues (381-386) located in the C-terminal of p53 CTD are in close contact with CBP, while the rest of the fragment quickly folds over and the polar residues His368, Ser371, Lys373, and Gln375 of p53 CTD interact with Asp1124 and Asp1127 of the CBP during the MD simulation. The residue AcLys382 consistently remains to be inserted into the hydrophobic cavity which is formed by Val1115, Ile1122, Tyr1125, Tyr1167, Val1174, and Phe1177 of CBP. The FEP simulations of the acetylated and non-acetylated peptides in complex with CBP bromo-domain showed that the relative contribution of the acetyl group to p53 binding is 4.8 ± 0.5 kcal/mol [129], which is a major portion of the binding free energy and in excellent agreement with experimental results spanning a range of 4.3-4.9 kcal/mol [130]. Phe385 and Lys386 were found to form contacts with Arg1112. Thus hydrophobic contacts are the dominant factor in the binding of CBP bromo-domain with p53 CTD. In addition, the salt bridges between Arg379 and Asp1124, Lys386 and Glu1105, and the H-bonding between AcLys382 and Tyr1167, AcLys382 and Asn1168, Leu383 and Ser1172 also contribute to the interaction [129]. All in all, it is clear that the AcLys382 plays an important role in the interaction between p53 and CBP.

Sirtuin-2 (SIRT2), a NAD + -dependent histone deacetylases, is crucial in cell cycle regulation [131][132][133][134], specifically binds and deacetylates p53 protein and functions in transcriptional silencing processes [135,136]. The X-ray co-crystal structure shows that the Lys382 acetylated p53 CTD peptide forms a β-strand from residues 379 to 387 (PDB ID: 1MA3), where the acetyl-lysine side chain inserts into a conserved, largely hydrophobic tunnel when it binds to a Sirtuin protein (Sir2) [38]. This β-strand of the acetylated peptide binds in the cleft and forms a stable β-sheet along with two other β-strands of Sir2. The binding process between the p53 CTD and Sir2 was investigated by MD simulation [120]. The results show that acetylated p53 CTD fragment exhibits very little secondary structure, with only a turn formed around residue Lys 381 in two of the three replicate simulations when bound to Sir2. However, the unbound acetylated p53 CTD fragment adopts many different ordered conformations, including a discontinuous turn and most notably a very short sheet-turn-sheet structural motif. The RMSD analysis reveals that the acetylated p53 CTD fragment becomes stable when bound to the Sir2. In the absence of p53 CTD, the p53 CTD-Sir2 interface (residues 45-60 and residues 150-175 in the Sir2) and most notably hydrophilic and acidic residues near the hydrophobic binding site of Sir2 fluctuate more. The backbone atoms of His380, AcLys382, Leu383 and Phe385 of the p53 CTD form persistent contacts with Gly166, Glu167, Leu169, Val195 and Tyr197 of Sir2 in all three replicate simulations. Van der Waals interactions between the AcLys382 of p53 CTD fragment and His118, Val63, Phe165, Leu169, and Val196 of Sir2 were found to be the major contributors to the interaction energy. Among all residues, the AcLys382 retains its position and plays the largest role in the acetylated p53 CTD-Sir2 interaction. Lysine acetylation of p53 is also required for its activation as a transcription factor.

Several studies have established that lysine methylation, like protein histone lysine methyltransferases KMT5 (SET9) methylates p53 at specific C-terminal lysines, and is a novel mechanism for modulating of p53 function [23,40,137,138]. In order to understand the underlying mechanisms, MD simulations were preformed to investigate the details of this interaction [120,139].

The crystal structure of the SET9 in complex with a mono-methylated (MeLys372) p53 CTD fragment (residues 369-374) was solved at 1.75 Å and used as the initial structure for MD simulations. It is different from the clear secondary structure when bound to S100B [19] or Sir2 [38], the C-terminal residues of p53 CTD were not solved and only 6-mer fragment left, it lacks any ordered secondary structure. According to Allen's study [120], only a bend centered at Ser 371 was detected when the p53 CTD in complex with SET9, whereas a sporadic turn and bend (residues 370-374) in the absence of SET9. The increased fluctuation was observed in the absence of the p53 CTD for the two important hydrophilic and acidic residues Asp256 and Thr266 in p53 CTD binding. PCA suggested a single vector corresponding to the opening and closing of the CTD binding site, described 72.0% of the overall motion in Set9 when it is free, but dropped to be only 38.3% when bound to p53 CTD. The MeLys372 inserts into a hydrophobic channel partly formed by the residues Val255 and Leu267. The residues Asp256, Thr266, and Ser268 of Set9 form stable hydrogen bonds with the p53 CTD fragment.

It is generally accepted that water plays an important role in the methylation process. Another study [139] was reported recently, finding the channel that surrounds p53-MeLys372 was not suitable for water movement, but a new channel composed of Gly292, Ala295, Tyr305 and Tyr335 was large enough for water moving in and out of the active site. The results indicated that SET9 is able to dimethylate p53-MeLys372.

Phosphorylation also modulates p53 function. The crystal structure of a complex containing a 9-mer p53 CTD fragment (residues 378-386) and cyclin A was solved at 2.24 Å [41]. To investigate the interactions between the two proteins, MD simulations were preformed [120]. The results pointed that the structure of the p53 CTD fragment became stabile when bound to cyclin A. The cyclin A forms a small hydrophobic binding pocket, in which Leu253 interacts with the hydrophobic residues of the p53 CTD. The residues Glu220, Glu224 and Asp283 of cyclin A contact with Ser378, Arg379, His380 and Lys381 of the p53 CTD comprising an important electrostatic contribution, Ile213, Leu214, and Leu253 of cyclin A interact with Leu383 and Phe385 of the p53 CTD creating the van der Waals contribution.

In summary, site specific post-translational modifications of acetylation, methylation and phosphorylation in CTD play an important role in the function of regulation of p53.

Small Molecule and Peptide Inhibitors Computational Modeling and Design

Recently, inhibition of the interaction between MDM2/X and p53 using non-peptide small-molecule or peptide inhibitor has received increasing attention and has been widely studied experimentally and computationally. For the well-defined structure of MDM2/X, inhibitors prefer to mimic p53 rather than MDM2/X, mainly mimic the side chain of three hydrophobic residues Phe19, Trp23 and Leu26, and the Leu22 which is also subsequently found to play a vital role in the p53-MDM2/X interaction, in the hope of gaining higher binding affinities than p53 TAD peptide [140]. The structures of p53-MDM2 and p53-MDMX (Figure 3) showed that the interactions of Phe19 and Trp23 of p53 are highly similar, but Leu26 is different. Besides, the three residues Leu54, His96 and Ile99 in MDM2 are different from the equivalent residues Met53, Pro95 and Leu98 in MDMX. These differences are possibly attributed to the variant binding of p53 peptide to MDM2 and MDMX. Hu et al. [141] carried out molecular docking and MD simulations to explore the dynamic behavior of MDM2 with a set of six non-peptide small-molecule inhibitors, extracted from the compounds designed by Ding et al. [142], and estimated the relative binding affinities by MM-PBSA. These inhibitors have a similar structure including 4 groups. For the inhibitor 8, G1 represents for 6-chloro-2-oxo-indole group, G2 for 3-chloro-2-fluoro-phenyl group, G3 for 2,2-dimethylbutane group and G4 for 2-morpholin-4yl-ethylamine group. The study showed that the van der Waals energies (−30 to −44 kcal/mol) are the dominant driving force for each complex, which indicated that the shape complementarity is the major contribution to the affinities. The MD simulations also showed that G1, G2 and G3 group in averaged structure mimic the Phe19, Trp23 and Leu26 residues in p53, respectively, with the hydrogen bond well kept. However, the G4 group does not act as the residue Leu22 of p53, interacting with the residues Phe55, Gln59 and Met62 of MDM2 by C-H…π and C-H…O interactions.

The inhibitory mechanisms of four small molecules on the MDM2-p53 interaction were investigated in our lab by 10 ns MD simulations and MM-GBSA analysis [143]. The binding free energies were −15.59, −15.63, −14.21 and −12.49 kcal/mol for IMZ, DIZ, YIN and K23 to MDM2, respectively. It can be seen that the estimated rank was in excellent agreement with the experimental rank. The van der Waals energy (spanning −37 to −45 kcal/mol) was the dominant factor for the binding of the four inhibitors among the binding free energy components. The R1, R2 and R3 (R4) groups of these inhibitors insert into three hydrophobic sub-pockets of MDM2, producing CH-π and π-π interactions, which are favorable for the stable binding. In this study, the six residues Leu54, Gly58, Ile61, Met62, Val93 and His96 in MDM2 were found to play important roles in hydrophobic interactions with the four inhibitors.

The inhibitory activity of 1,4-benzodiazepine-2,5-diones (BDPs) and N-Acylpolyamine (NAPA) derivatives to HDM2 was investigated by Wang et al. [144] through 3D-QSAR and MD-based docking. For BDPs, the key residues are Gln24, Lys51, Gly16, Ser17, Leu54, Gln72, Gly58, Leu57, Val93, Phe55 and His73. The small electronegative groups at R1 and R3 regions mainly involved H-bonds formation, while large groups at R2 and R4 regions are favorable for BDP's inhibitory activity. On the other hand, for the NAPA-protein complex, the key residues are Glu25, Thr26, Tyr100, Tyr104, His96, Leu54, Met50 and Phe55. Both huge and electropositive groups in B ring and A ring, and small groups at region P are to the benefit of NAPA derivatives inhibitory activity.

The p53-mimic peptide PMI (sequence TSFAEYWNLLSP, Figure 2D) has higher affinity than the p53 (the corresponding sequence ETFSDLWKLLPE) to the binding of MDM2. In order to explore the role of each residue in this peptide, Liu et al. [145] preformed multiple short simulations in explicit water environment of each alanine mutation. The binding free energies from the short MD simulations were in good agreement with that from 50 ns MD simulations. The F3A, Y6A, W7A, L10A mutants of PMI incur large loss of both van der Waals (8, 9, 12 and 6 kT for F3A, Y6A, W7A and L10A, respectively) and electrostatic interaction energies (12, 11 and 5 kT for F3A, Y6A and W7A, respectively) than other mutants.

The MD simulations followed by MM-GBSA method were performed to explore the binding mechanisms of peptide and small molecule inhibitors to MDMX (PMI-MDMX and WW8-MDMX complexes) by Cheng et al. [146] (Figure 2B,D). The non-peptide WW8 induces high flexibility, in particular, the α-helix in the C-terminal of MDMX is more obvious, and a stronger anti-correlate motion than peptide PMI for MDMX in most regions except the loop L2. Compared to the initial structure, the side chain conformation of residue Tyr99 was changed to generate a steric clash in both simulations. The binding free energies estimated by MM-GBSA are −15.49 and −9.26 kcal/mol, and van der Waals energies are −62.71 and −45.90 kcal/mol for PMI-MDMX and WW8-MDMX, respectively. The van der Waals energies are the major favorable contributors to the inhibitor-MDMX bindings. For its extra interactions with Lys50, Tyr66, Gln71 and His72, the peptide inhibitor PMI generate much stronger electrostatic and van der Waals interactions than non-peptide inhibitor WW8 with MDMX. Six common residues Met53, Ile60, Met61, Tyr66, Val92 and Leu98 are located in the hot spots of the interface between the inhibitors and MDMX determined by computational alanine scanning. These results lead to a conclusion that the hydrophobic CH-CH, CH-π and π-π interaction play important roles in the inhibitor-MDMX binding, and the peptide inhibitors stabilize the inhibitor-MDMX complex.

Though the two proteins MDM2 and MDMX share structural homology, their p53 binding sites are not exact same, and some inhibitors (like nutlin) show different affinity to MDM2 and MDMX. Consequently, to identify MDM2/MDMX dual-inhibitors opens a new era of anti-cancer drug development. Barakat et al. [147] conducted such a study by ensemble-based virtual screening to identify effective dual-inhibitors which could simultaneously disrupt the p53-MDM2/X interactions. Their results confirmed the three MDM2-inhibitors, Nutlin3, MI-219 and TDP, have low affinities to MDMX.

The first non-peptide small-molecule used to block p53-MDM2/X interactions was nutlin (Figure 2A). The different binding modes of p53 peptide and nutlin to MDM2 and MDMX were investigated by Joseph et al. [148] using MD simulations. Residues Glu69, Lys70 and Val93 in MDM2 and Gln69, Glu70, Asp79 and Asp94 in MDMX displayed the highest fluctuations. The p53 peptide appears to stable the two proteins so that the p53-bound states display less flexibility than their apo and nutlin-bound states. A longer helix was formed for p53 peptide when bound to MDM2 than to MDMX, for a salt bridge between Glu17 of p53 with Arg65 in MDM2 was formed. The hydrogen bond between Trp23 (p53) and Leu54 (MDM2) and the intra-molecular hydrogen bond between Glu17 (MDM2) and Ser20 (MDM2), in particular, the salt bridge between Glu17 and Arg65 plays an important role in maintaining the stability of complex. The hydrogen bond between Trp23 (p53) and Leu54 (MDMX) and the intra-molecular hydrogen bond between Phe19 (MDMX) and Gln17 (MDMX), are retained with the occupancy of 48% and 66% respectively during the simulations. However, the corresponding salt bridge Glu17-Arg65 in MDM2 could not form in MDMX for this Arg65 was replaced by Gln64. The binding free energies were estimated to be −6.4 and −3.7 kcal/mol for p53-MDM2 and p53-MDMX, respectively. Since the van der Waals interactions were similar, the driving force for the binding of p53-MDM2 was easily identified due to stronger electrostatic interactions than p53-MDMX. When bound to nutlin, Leu54 in MDM2 is replaced by the larger Met53 in MDMX, the binding pocket of MDMX is changed to be smaller than that in MDM2, resulting in reduced electrostatic and van der Waals interactions for MDMX-nutlin. The larger entropic penalty for p53 upon the binding leads to a lower affinity of p53 than nutlin for MDM2. The binding free energies were estimated to be −17.9 and −0.3 kcal/mol for MDM2 and MDMX in the nutlin-bound state, respectively. MDM2 binds both p53 and nutlin while MDMX interacts with nutlin with extremely low affinity.

Lu et al. [149] studied the pyrrolopyrimidine-based α-helix mimetic dual inhibitors of MDM2 and MDMX, including an active compound 3a and one inactive compound NC-1, using molecular docking and MD simulations followed by MM-PB/GBSA calculations. The R1, R2, and R3 groups of the pyrrolopyrimidine scaffold compound 3a could mimic the spatial orientation of the side chains of Phe19, Trp23, and Leu26 in p53 and their interactions with MDM2/X, respectively. However, the R1 group of NC-1 is a methyl sulfonyl moiety that could not mimic the side chain of Phe19 of p53. The compound 3a has similar binding affinity for both protein MDM2 and MDMX than that of NC-1, which departs from the active site of MDMX during MD simulations, with the binding free energies of The residues Leu54/Met53, Ile61/Ile60, Gln72/Gln71, His73/His72, and Val93/Val92 for MDM2 and MDMX, which locate in the Phe19 sub-pocket of p53 have made outstanding contributions to the binding in the hydrophobic interface cleft. This result also proved that the critical residue Phe19 is indispensable for the interaction between p53 and MDM2/X. The interactions between those two inhibitors and MDM2/X are dominated by the large van der Waals energy (−47.69, −34.35 and −48.87 kcal/mol), which indicates that the shape complementarity is important for designing high affinity inhibitors.

Another study to probe the difference in binding mechanisms of peptide (pDI6W, sequence LTFEHWWAQLTS) and non-peptide inhibitor (K23, with four aromatic groups) to MDM2/X was also investigated by our lab [150] (Figure 2C,E). The dynamic behavior analysis suggested that the inhibitor-MDM2 is much more stable than inhibitor-MDMX complex. Tyr99 in MDMX undergoes a larger conformational change for K23-bound than pDI6W-bound. The predicted binding free energies are −21.93, −19.74, −16.81 and −14.89 kcal/mol for pDI6W-MDM2, pDI6W-MDMX, K23-MDM2 and K23-MDMX complex, respectively, which agree well with the experimental rank. The van der Waals energies were found to dominate the bindings of the inhibitors in the hydrophobic cleft of MDM2/X. The peptide inhibitor pDI6W generates much stronger electrostatic and van der Waals interactions with the two proteins than the non-peptide inhibitor K23. The binding mode analysis for inhibitor-MDM2 complex revealed that four residues Tyr67, Met62, Ile61 and Gln72 form the first hydrophobic sub-pocket (Phe19 sub-pocket); Leu54, Gly58 and Val93 form the second hydrophobic sub-pocket (Trp23 sub-pocket); Leu54, Val93, His96 and Ile99 form the third hydrophobic sub-pocket (Leu26 sub-pocket). For the inhibitor-MDMX complex, residues Ile60, Met53, Tyr66 and Gln71 enclose the first hydrophobic sub-pocket; Met53, His54, Gly57 and Val92 shape the second hydrophobic sub-pocket; Met53, Val92, Pro95 and Leu98 shape the third hydrophobic sub-pocket. According to the computational alanine scanning analysis, it is clear that the π-π, CH-π and CH-CH interactions produced by eight residues (Lys51, Leu54, Ile61, Met62, Tyr67, Val93, His96 and Ile99) in the pDI6W-MDM2 complex, four residues (Leu54, Ile61, Val93 and Ile99) in the K23-MDM2 complex, six residues (Lys50, His54, Ile60, Met61, Tyr66 and Val92) in the pDI6W-MDMX complex and four residues (Met53, Ile60, Val92 and Leu98) in the K23-MDMX complex, play a key roles in the binding of the inhibitors in the hydrophobic cleft of MDM2/X.

In some sense peptide inhibitors have an advantage over non-peptide small molecules in targeting protein-protein interaction interfaces since they can be finely tuned to be much more specific. However, peptide inhibitors are usually unstable and subject to proteolytic degradation. In order to address this problem, the "peptide stapling" strategy using an all-hydrocarbon cross-link (staple) was applied to generate a series of stabilized alpha-helix of p53 (SAH-p53) cyclic peptides [151], in which the key interacting residues were kept and the amphiphilic α-helical structure was reinforced. Such SAH-p53 peptides exhibit significantly high affinity for MDM2/X in contrast to the corresponding native p53 peptide, and improve the cell permeability and increased helical stability compared to the peptidomimetics, therefore, they are becoming a new trend of anti-cancer drug development. The structure of the MDM2 bound to a stapled p53 peptide was determined recently [61] (Figure 2E). There have been some stapled peptide inhibitors, with increased helicity and cellular uptake without apparent cytotoxicity, are also reported [152][153][154].

To explore the structural and dynamic stability of the hydrocarbon cross-links stapled peptides, Guo et al. [155] studied 10 stapled α-helical peptides of the 16-residue of the p53 TAD fragment over a range of temperatures in solution by all-atom MD simulations. The overall propensity to form a-helical structure of stapled peptides was found to exhibit substantial variations, depending both on the relative position of the stapling and the specific residue substitutions. In order to understand their effects on the inhibition of MDM2/X, Joseph et al. [156] reported novel interactions of the staples with the target proteins by performing on all-atom MD simulations on the complexes of 9 stapled peptides to MDM2 and 2 stapled peptides to MDMX. The study confirmed that good binders extend the helicity in solution and increase the affinity when bound to MDM2 for the extensive interaction between hydrocarbon staples and the MDM2 surface. In addition, the contributions of stapled peptides are equivalent to Phe19 and Trp23 of p53 peptide. The stapling site, charge quantity and distribution, and overall hydrophobicity are important for the inhibitory activity. These studies provide new insights into the anti-cancer drug design of stapled peptides and the development of potent inhibitors preventing p53 from interacting with target proteins.

Like MDM2, S100B is activated at the transcriptional level by active p53 and then subsequently negatively regulates p53 function as a tumor suppressor via feedback loop control [112,157,158]. As one of the most effective ideas to develop anti-cancer drugs, there has been some researches focusing on the identification, characterization, design and in vitro screening of small molecule inhibitors of the interaction between calcium-dependent S100B and p53 tumor suppressor [64,[159][160][161]. In addition, nuclear magnetic resonance and X-ray crystallography studies have been reported about the structures of three small molecules (SBi132, SBi1279, and SBi523) and these three molecules were found to bind in distinct locations and orientations within the hydrophobic binding cleft of calcium-bound S100B [162]. The strategy of combining computer aided drug design method with NMR spectroscopy and X-ray crystallography is being used to identify inhibitors of the p53-S100B interaction. Several computational studies about the ligands with the potential to block the binding of p53-S100B have been published [163,164].

Whitlow JL and co-workers [163] studied the interaction between p53 and S100B (PDB ID: 1DT7) using MD simulations followed by the MM-PBSA method. Initially, organic molecules were identified as potential inhibitors to block the p53-S100B interaction using the LUDI program implemented in Cerius2 package and ZINC database, and then 14 ZINC compounds were selected for the S100B-ligand binding affinity calculation by MM-PBSA. The RMSD of the p53-S100B complex from MD demonstrated that the C-terminal peptide of p53 is very flexible whereas the S100B structure stays very stable. The same results also appeared in the subsequent MD studies [120,121], indicating that the α-helical structure of C-terminal peptide of p53 gradually loses its secondary structure during simulation. The analysis of the interaction between S100B and p53 showed that the van der Waals interaction is very weak and the electrostatic interaction dominates the binding of S100B with p53 CTD. This indicates that the residues Glu45 and Glu46 of S100B have the strongest interaction with p53. For the S100B-ligand systems, the 14 ZINC compounds were found to interact with S100B at residues Glu45, Glu46, Leu44, Val52, Lys55, Val56, Thr59, Phe76, Met79, Val80, Ala83 and Phe87, displaying the similar binding modes and interactions as p53 with S100B, and therefore it was expected that they have the potential to prevent S100B from binding to p53 [163].

Another paper exploring ligand blockage of the interaction with p53 [164] was published in 2009. Three compounds with the highest binding affinity to the S100B and the highest potential to block the protein interaction of S100B and p53 were selected from the NCI small molecular library using the Glide docking program. Based on the docked structures, it could be seen that the hydrophobic aromatic rings of the three selected compounds directly contact with S100B in the hydrophobic binding site mainly formed by the helix 2, 3, and 4. Besides, the hydrogen bonds between compounds and residues Glu45, Thr59, Gln71, and Lys111 of S100B were formed. MD simulations showed that the ligand induces the separation of p53 from S100B by about 2Å and the decreases of the binding affinity of S100B and p53 by ~8.5-14.6 kcal/mol compared to the original p53-S100B complex. Based on the binding free energy calculation done by MM-PBSA, the van der Waals energy could be the major driving force for the binding of the complexes of p53-S100B with or without the three ligands and the complexes of S100B-ligand. It could be perceived that the three compounds with hydrophobic aromatic rings propelled p53 away from S100B, and they are able to bind to free S100B and form stable S100B-ligand complexes with the binding affinities no weaker than p53. Therefore the availability of p53 is able to restore its biological function as a tumor suppressor. Design of compounds by effectively blocking p53-S100B interactions has been regarded as a new promising approach to the development of anti-cancer drugs.

Outlook

It is well known that p53 exerts its tumor suppressor function through a series of interactions of independently folded and intrinsically disordered functional domains. The full-length three-dimensional structure of this protein has not been determined to date. The unstructured TAD and CTD are especially regarded as the important modulatory domains for core interactions with specific sites for DNA binding.

Mutation, post-translational modifications and physiological interactions of vital residues in p53 with numerous negative regulator protein partners have a certain influence on the biological function of p53. At the same time, these negative regulator protein partners also interact with each other, such as MDM2-MDMX interaction [165,166] and MDM2-S100B interaction [167]. In addition, several proteins can bind p53 with more than one domain. All of these factors combine to create a complex regulation mechanism and network for p53 function.

Inhibitor design primarily mimics key interactions between p53 and protein partners, in order to achieve the purpose of releasing p53. Restoring endogenous p53 functions by the disruption of p53-protein interactions using peptide or non-peptide small molecule inhibitors holds a lot of promising strategies for anti-cancer drug design and development. Generally speaking, peptide inhibitors have an advantage over non-peptide small molecule inhibitors in targeting protein-protein interaction interfaces. Nevertheless, peptide inhibitors are unstructured in solution and highly susceptible to proteolytic degradation in vivo and therefore unsuitable for potential drugs. The "peptide stapling" strategy, using an all-hydrocarbon staple to enhance the stability and improve the cell permeability, may open a new era in the design and development of anti-cancer drug. Furthermore, traditional Chinese medicine and natural products attract a lot of attention [168][169][170][171][172][173] for their inhibitory activity of tumor growth and ability to induce apoptosis through restoring the p53 function to some extent. This research provides some trains of thought on the discovery and development of novel anti-cancer drugs.

With a better understanding of the interaction mechanism between p53 and its protein partner, it is possible to achieve efficient drug design. Today's MD simulations mainly adopt the classical molecular mechanics force fields, which rely on the representation of atomic electrostatic properties using fixed-point charge models. However, these crude models neither fully describe the real molecular electrostatic distribution, nor consider the influence of electrostatic polarization. The estimated binding free energies from MD simulations are often artificially high compared to measured values. Previous MD simulation studies suggested that electrostatic interactions and electrostatic component of solvation free energy provide a dominant role in p53-protein interactions. From the molecular simulation point of view, the polarizable force field is one of the hot topics considering the accurate electrostatic and polarization effect, which may be more suitable for this interaction study, thereby contributing to the rational development and design of new anti-cancer drugs.

Figure 1 .
1Structures of p53 protein. (A) The complex of p53 transcriptional activation domain (TAD) fragment bound to MDM2 (PDB 1YCR)

Figure 2 .
2Structures of the complexes of MDM2 with the non-peptide inhibitor and peptide inhibitors. MDM2 are shown in cartoon and cyan, inhibitors are shown in stick and green; (A) MDM2 with Nutlin (from PDB 1RV1) [59]; (B) MDM2 with K23 (from PDB 3LBK) [45]; (C) MDM2 with DIZ (from PDB 1T4E) [60]; (D) MDM2 with PMI (from PDB 3LNZ) [49]

Figure 3 .
3The superimposition of the structures of p53-MDM2 and p53-MDMX complex. The complex of p53-MDM2 (PDB 1YCR)[18] is shown in cartoon, p53 TAD fragment (residues 17-29) is shown in yellow and the three most important residues Phe19, Trp23 and Leu26 are shown in stick, MDM2 (residues 25-109) is shown in white and the four residues Leu54, His96, Ile99 and Tyr100 in MDM2 are shown in green and stick; The complex of p53-MDMX (PDB 3DAB)[27] is shown in cartoon, p53 TAD fragment (residues 17-29) is shown in orange and the three most important residues Phe19, Trp23 and Leu26 are shown in stick, MDMX (residues 23-110) is shown in grey and the equivalent residues Met53, Pro95, Leu98 and Tyr99 in MDMX are shown in light magenta and stick; the spatial orientation of the side chains of Phe19 and Trp23 of p53 are highly similar, but that of Leu26 is different in p53-MDM2 and p53-MDMX complex. (Insert) Magnified view of important residues of p53 and MDM2/X. Figures were created with Pymol (http://pymol.org)[20].

−15.16(PB)/ −14.74(GB) −3.01(PB)/ −2.23(GB) and −13.97(PB)/ −14.24(GB) kcal/mol for MDM2-3a, MDM2-NC-1 and MDMX-3a respectively.

42 .
42Greenblatt, M.S.; Bennett, W.P.; Hollstein, M.; Harris, C.C. Mutations in the p53 tumor suppressor gene: Clues to cancer etiology and molecular pathogenesis. Cancer Res. 1994, 54, 4855-4878. 43. Pfeifer, G.P.; Besaratinia, A. Mutational spectra of human cancer. Hum. Genet. 2009, 125, 493-506. 44. Goh, A.M.; Coffill, C.R.; Lane, D.P. The role of mutant p53 in human cancer. J. Pathol. 2011, 223, 116-126. 45. Popowicz, G.M.; Czarna, A.; Wolf, S.; Wang, K.; Wang, W.; Dömling, A.; Holak, T.A. Structures of low molecular weight inhibitors bound to MDMX and MDM2 reveal new approaches for p53-MDMX/MDM2 antagonist drug discovery. Cell Cycle 2010, 9, 1104-1111. 46. Dickens, M.P.; Fitzgerald, R.; Fischer, P.M. Small-molecule inhibitors of MDM2 as new anticancer therapeutics. Semin. Cancer Biol. 2010, 20, 10-18. 47. Pazgier, M.; Liu, M.; Zou, G.; Yuan, W.; Li, C.; Li, C.; Li, J.; Monbo, J.; Zella, D.; Tarasov, S.G.; Lu, W. Structural basis for high-affinity peptide inhibition of p53 interactions with MDM2 and MDMX. Proc. Natl. Acad. Sci. USA 2009, 106, 4665-4670. 48. Phan, J.; Li, Z.; Kasprzak, A.; Li, B.; Sebti, S.; Guida, W.; Schönbrunn, E.; Chen, J. Structure-based design of high affinity peptides inhibiting the interaction of p53 with MDM2 and MDMX. J. Biol. Chem. 2010, 285, 2174-2183. 49. Li, C.; Pazgier, M.; Li, C.; Yuan, W.; Liu, M.; Wei, G.; Lu, W.Y.; Lu, W. Systematic mutational analysis of peptide inhibition of the p53-MDM2/MDMX interactions. J. Mol. Biol. 2010, 398, 200-213. 50. Dudkina, A.S.; Lindsley, C.W. Small molecule protein-protein inhibitors for the p53-MDM2 interaction. Curr. Top. Med. Chem. 2007, 7, 952-960. 51. Zhan, C.; Lu, W. Peptide activators of the p53 tumor suppressor. Curr. Pharm. Des. 2011, 17, 603-609. 52. Hu, C.Q.; Hu, Y.Z. Small molecule inhibitors of the p53-MDM2. Curr. Med. Chem. 2008, 15, 1720-1730. 53. Patel, S.; Player, M.R. Small-molecule inhibitors of the p53-HDM2 interaction for the treatment of cancer. Expert Opin. Investig. Drugs 2008, 17, 1865-1882. 54. Riedinger, C.; McDonnell, J.M. Inhibitors of MDM2 and MDMX: A structural perspective. Future Med. Chem. 2009, 1, 1075-1094. 55. Millard, M.; Pathania, D.; Grande, F.; Xu, S.; Neamati, N. Small-molecule inhibitors of p53-MDM2 interaction: The 2006-2010 update. Curr. Pharm. Des. 2011, 17, 536-559. 56. Vu, B.T.; Vassilev, L. Small-molecule inhibitors of the p53-MDM2 interaction. Curr. Top. Microbiol. Immunol. 2011, 348, 151-172. 57. Popowicz, G.M.; Dömling, A.; Holak, T.A. The structure-based design of Mdm2/Mdmx-p53 inhibitors gets serious. Angew. Chem. Int. Ed. Engl. 2011, 50, 2680-2688. 58. Joerger, A.C.; Fersht, A.R. The tumor suppressor p53: From structures to drug discovery. Cold Spring Harb. Perspect. Biol. 2010, 2, doi:10.1101/cshperspect.a000919. 59. Vassilev, L.T.; Vu, B.T.; Graves, B.; Carvajal, D.; Podlaski, F.; Filipovic, Z.; Kong, N.; Kammlott, U.; Lukacs, C.; Klein, C.; Fotouhi, N.; Liu, E.A. In vivo activation of the p53 pathway by small-molecule antagonists of MDM2. Science 2004, 303, 844-848.
AcknowledgmentsThis work is supported by the National Natural Science Foundation of China (31070641), the National Basic Research Program of China (2012CB721002), and "Hundreds Talents Program" of the Chinese Academy of Sciences.
Detection of a transformation-related antigen in chemically induced sarcomas and other transformed cells of the mouse. A B Deleo, G Jay, E Appella, G C Dubois, L W Law, L J Old, Proc. Natl. Acad. Sci. Natl. Acad. SciUSA76DeLeo, A.B.; Jay, G.; Appella, E.; Dubois, G.C.; Law, L.W.; Old, L.J. Detection of a transformation-related antigen in chemically induced sarcomas and other transformed cells of the mouse. Proc. Natl. Acad. Sci. USA 1979, 76, 2420-2424.

Two faces of p53: Aging and tumor suppression. F Rodier, J Campisi, D Bhaumik, Nucleic Acids Res. 35Rodier, F.; Campisi, J.; Bhaumik. D. Two faces of p53: Aging and tumor suppression. Nucleic Acids Res. 2007, 35, 7475-7484.

Inhibition of tumor angiogenesis by p53: A new role for the guardian of the genome. J G Teodoro, S K Evans, M R Green, J. Mol. Med. (Berl. 85Teodoro, J.G.; Evans, S.K.; Green, M.R. Inhibition of tumor angiogenesis by p53: A new role for the guardian of the genome. J. Mol. Med. (Berl.) 2007, 85, 1175-1186.

p53 in health and disease. K H Vousden, D P Lane, Nat. Rev. Mol. Cell Biol. 8Vousden, K.H.; Lane, D.P. p53 in health and disease. Nat. Rev. Mol. Cell Biol. 2007, 8, 275-283.

The genetics of the p53 pathway, apoptosis and cancer therapy. A Vazquez, E E Bond, A J Levine, G L Bond, Nat. Rev. Drug Discov. 7Vazquez, A.; Bond, E.E.; Levine, A.J.; Bond, G.L. The genetics of the p53 pathway, apoptosis and cancer therapy. Nat. Rev. Drug Discov. 2008, 7, 979-987.

Participation of p53 protein in the cellular response to DNA damage. M B Kastan, O Onyekwere, D Sidransky, B Vogelstein, R W Craig, Cancer Res. 51Kastan, M.B.; Onyekwere, O.; Sidransky, D.; Vogelstein, B.; Craig, R.W. Participation of p53 protein in the cellular response to DNA damage. Cancer Res. 1991, 51, 6304-6311.

UV-radiation, apoptosis and skin. N Pustisek, M Situm, Coll. Antropol. 35Pustisek, N.; Situm, M. UV-radiation, apoptosis and skin. Coll. Antropol. 2011, 35, 339-341.

Activation of p53 following ionizing radiation, but not other stressors, is dependent on the proline-rich domain (PRD). H G Campbell, R Mehta, A A Neumann, C Rubio, M Baird, T L Slatter, A W Braithwaite, 10.1038/onc.2012.102Oncogene. Campbell, H.G.; Mehta, R.; Neumann, A.A.; Rubio, C.; Baird, M.; Slatter, T.L.; Braithwaite, A.W. Activation of p53 following ionizing radiation, but not other stressors, is dependent on the proline-rich domain (PRD). Oncogene 2012, doi:10.1038/onc.2012.102.

The amino acid sequence of murine p53 determined from a c-DNA clone. D Pennica, D V Goeddel, J S Hayflick, N C Reich, C W Anderson, A J Levine, Virology. 134Pennica, D.; Goeddel, D.V.; Hayflick, J.S.; Reich, N.C.; Anderson, C.W.; Levine, A.J. The amino acid sequence of murine p53 determined from a c-DNA clone. Virology 1984, 134, 477-482.

The DNA-binding domain of p53 contains the four conserved regions and the major mutation hot spots. N P Pavletich, K A Chambers, C O Pabo, Genes Dev. 7Pavletich, N.P.; Chambers, K.A.; Pabo, C.O. The DNA-binding domain of p53 contains the four conserved regions and the major mutation hot spots. Genes Dev. 1993, 7, 2556-2564.

Identification of a novel p53 functional domain that is necessary for efficient growth suppression. K K Walker, A J Levine, Proc. Natl. Acad. Sci. Natl. Acad. SciUSA93Walker, K.K.; Levine, A.J. Identification of a novel p53 functional domain that is necessary for efficient growth suppression. Proc. Natl. Acad. Sci. USA 1996, 93, 15335-15340.

Structure-function-rescue: The diverse nature of common p53 cancer mutants. A C Joerger, A R Fersht, Oncogene. 26Joerger, A.C.; Fersht, A.R. Structure-function-rescue: The diverse nature of common p53 cancer mutants. Oncogene 2007, 26, 2226-2242.

Crystal structure of a p53 tumor suppressor-DNA complex: Understanding tumorigenic mutations. Y Cho, S Gorina, P D Jeffrey, N P Pavletich, Science. 265Cho, Y.; Gorina, S.; Jeffrey, P.D.; Pavletich, N.P. Crystal structure of a p53 tumor suppressor-DNA complex: Understanding tumorigenic mutations. Science 1994, 265, 346-355.

. M Kitayner, H Rozenberg, R Rohs, O Suad, D Rabinovich, B Honig, Shakked, Z. Diversity in DNA recognition by p53 revealed by crystal structures with Hoogsteen base pairsKitayner, M.; Rozenberg, H.; Rohs, R.; Suad, O.; Rabinovich, D.; Honig, B.; Shakked, Z. Diversity in DNA recognition by p53 revealed by crystal structures with Hoogsteen base pairs.

. Nat. Struct. Mol. Biol. 17Nat. Struct. Mol. Biol. 2010, 17, 423-429.

Crystal structure of the p53 core domain bound to a full consensus site as a self-assembled tetramer. Y Chen, R Dey, L Chen, Structure. 18Chen, Y.; Dey, R.; Chen, L. Crystal structure of the p53 core domain bound to a full consensus site as a self-assembled tetramer. Structure 2010, 18, 246-256.

Solution structure of the tetrameric minimum transforming domain of p53. W Lee, T S Harvey, Y Yin, P Yau, D Litchfield, C H Arrowsmith, Nat. Struct. Biol. 1Lee, W.; Harvey, T.S.; Yin, Y.; Yau, P.; Litchfield, D.; Arrowsmith, C.H. Solution structure of the tetrameric minimum transforming domain of p53. Nat. Struct. Biol. 1994, 1, 877-890.

p53 contains large unstructured regions in its native state. S Bell, C Klein, L Müller, S Hansen, J Buchner, J. Mol. Biol. 322Bell, S.; Klein, C.; Müller, L.; Hansen, S.; Buchner, J. p53 contains large unstructured regions in its native state. J. Mol. Biol. 2002, 322, 917-927.

Structure of the MDM2 oncoprotein bound to the p53 tumor suppressor transactivation domain. P H Kussie, S Gorina, V Marechal, B Elenbaas, J Moreau, A J Levine, N P Pavletich, Science. 274Kussie, P.H.; Gorina, S.; Marechal, V.; Elenbaas, B.; Moreau, J.; Levine, A.J.; Pavletich, N.P. Structure of the MDM2 oncoprotein bound to the p53 tumor suppressor transactivation domain. Science 1996, 274, 948-953.

Structure of the negative regulatory domain of p53 bound to S100B(betabeta). R R Rustandi, D M Baldisseri, D J Weber, Nat. Struct. Biol. 7Rustandi, R.R.; Baldisseri, D.M.; Weber, D.J. Structure of the negative regulatory domain of p53 bound to S100B(betabeta). Nat. Struct. Biol. 2000, 7, 570-574.

The PyMOL Molecular Graphics System, Version 1.2r3pre. DeLano Scientific. San Carlos, CA, USAThe PyMOL Molecular Graphics System, Version 1.2r3pre. DeLano Scientific: San Carlos, CA, USA. Available online: http://pymol.org (accessed on 1 January 2012).

DNA damage activates p53 through a phosphorylation-acetylation cascade. K Sakaguchi, J E Herrera, S Saito, T Miki, M Bustin, A Vassilev, C W Anderson, E Appella, Genes Dev. 12Sakaguchi, K.; Herrera, J.E.; Saito, S.; Miki, T.; Bustin, M.; Vassilev, A.; Anderson, C.W.; Appella, E. DNA damage activates p53 through a phosphorylation-acetylation cascade. Genes Dev. 1998, 12, 2831-2841.

Influence of tetramerisation on site-specific post-translational modifications of p53: Comparison of human and murine p53 tumor suppressor protein. L J Warnock, A Knox, T R Mee, S A Raines, J Milner, Cancer Biol. Ther. 7Warnock, L.J.; Knox, A.; Mee, T.R.; Raines, S.A.; Milner, J. Influence of tetramerisation on site-specific post-translational modifications of p53: Comparison of human and murine p53 tumor suppressor protein. Cancer Biol. Ther. 2008, 7, 1481-1489.

Protein methylation: A new mechanism of p53 tumor suppressor regulation. A Scoumanne, X Chen, Histol. Histopathol. 23Scoumanne, A.; Chen, X. Protein methylation: A new mechanism of p53 tumor suppressor regulation. Histol. Histopathol. 2008, 23, 1143-1149.

Post-translational modifications of p53 tumor suppressor: Determinants of its functional targets. N Taira, K Yoshida, Histol. Histopathol. 27Taira, N.; Yoshida, K. Post-translational modifications of p53 tumor suppressor: Determinants of its functional targets. Histol. Histopathol. 2012, 27, 437-443.

Targeting protein-protein interactions: Lessons from p53/MDM2. J K Murray, S H Gellman, Biopolymers. 88Murray, J.K.; Gellman, S.H. Targeting protein-protein interactions: Lessons from p53/MDM2. Biopolymers 2007, 88, 657-686.

Targeting oncogenic protein-protein interactions by diversity oriented synthesis and combinatorial chemistry approaches. A G Tzakos, D Fokas, C Johannes, V Moussis, E Hatzimichael, E Briasoulis, Molecules. 16Tzakos, A.G.; Fokas, D.; Johannes, C.; Moussis, V.; Hatzimichael, E.; Briasoulis, E. Targeting oncogenic protein-protein interactions by diversity oriented synthesis and combinatorial chemistry approaches. Molecules 2011, 16, 4408-4427.

Structure of the human Mdmx protein bound to the p53 tumor suppressor transactivation domain. G M Popowicz, A Czarna, T A Holak, Cell Cycle. 7Popowicz, G.M.; Czarna, A.; Holak, T.A. Structure of the human Mdmx protein bound to the p53 tumor suppressor transactivation domain. Cell Cycle 2008, 7, 2441-2443.

PUMA couples the nuclear and cytoplasmic proapoptotic function of p53. J E Chipuk, L Bouchier-Hayes, T Kuwana, D D Newmeyer, D R Green, Science. 309Chipuk, J.E.; Bouchier-Hayes, L.; Kuwana, T.; Newmeyer, D.D.; Green, D.R. PUMA couples the nuclear and cytoplasmic proapoptotic function of p53. Science 2005, 309, 1732-1735.

The Bad guy cooperates with good cop p53: Bad is transcriptionally up-regulated by p53 and forms a Bad/p53 complex at the mitochondria to induce apoptosis. P Jiang, W Du, K Heese, M Wu, Mol. Cell. Biol. 26Jiang, P.; Du, W.; Heese, K.; Wu, M. The Bad guy cooperates with good cop p53: Bad is transcriptionally up-regulated by p53 and forms a Bad/p53 complex at the mitochondria to induce apoptosis. Mol. Cell. Biol. 2006, 26, 9071-9082.

Structure of Bcl-xL-Bak peptide complex: Recognition between regulators of apoptosis. M Sattler, H Liang, D Nettesheim, R P Meadows, J E Harlan, M Eberstadt, H S Yoon, S B Shuker, B S Chang, A J Minn, C B Thompson, S W Fesik, Science. 275Sattler, M.; Liang, H.; Nettesheim, D.; Meadows, R.P.; Harlan, J.E.; Eberstadt, M.; Yoon, H.S.; Shuker, S.B.; Chang, B.S.; Minn, A.J.; Thompson, C.B.; Fesik, S.W. Structure of Bcl-xL-Bak peptide complex: Recognition between regulators of apoptosis. Science 1997, 275, 983-986.

Bacterial redox protein Azurin, tumor suppressor protein p53, and regression of cancer. T Yamada, M Goto, V Punj, O Zaborina, M L Chen, K Kimbara, D Majumdar, E Cunningham, T K D Gupta, A M Chakrabarty, Proc. Natl. Acad. Sci. Natl. Acad. SciUSA99Yamada, T.; Goto, M.; Punj, V.; Zaborina, O.; Chen, M.L.; Kimbara, K.; Majumdar, D.; Cunningham, E.; Gupta, T.K.D.; Chakrabarty, A.M. Bacterial redox protein Azurin, tumor suppressor protein p53, and regression of cancer. Proc. Natl. Acad. Sci. USA 2002, 99, 14098-14103.

Bacterial cupredoxin azurin and its interactions with the tumor suppressor protein p53. V Punj, T K D Gupta, A M Chakrabarty, Biochem. Biophys. Res. Commun. 312Punj, V.; Gupta, T.K.D.; Chakrabarty, A.M. Bacterial cupredoxin azurin and its interactions with the tumor suppressor protein p53. Biochem. Biophys. Res. Commun. 2003, 312, 109-114.

Apoptosis or growth arrest: Modulation of tumor suppressor p53's specificity by bacterial redox protein azurin. T Yamada, Y Hiraoka, M Ikehata, K Kimbara, B S Avner, T K D Gupta, A M Chakrabarty, Proc. Natl. Acad. Sci. 101Yamada, T.; Hiraoka, Y.; Ikehata, M.; Kimbara, K.; Avner, B.S.; Gupta, T.K.D.; Chakrabarty, A.M. Apoptosis or growth arrest: Modulation of tumor suppressor p53's specificity by bacterial redox protein azurin. Proc. Natl. Acad. Sci. USA 2004, 101, 4770-4775.

Single-stranded DNA mimicry in the p53 transactivation domain interaction with replication protein A. E Bochkareva, L Kaustov, A Ayed, G S Yi, Y Lu, A Pineda-Lucena, J C Liao, A L Okorokov, J Milner, C H Arrowsmith, A Bochkarev, Proc. Natl. Acad. Sci. Natl. Acad. SciUSA102Bochkareva, E.; Kaustov, L.; Ayed, A.; Yi, G.S.; Lu, Y.; Pineda-Lucena, A.; Liao, J.C.; Okorokov, A.L.; Milner, J.; Arrowsmith, C.H.; Bochkarev, A. Single-stranded DNA mimicry in the p53 transactivation domain interaction with replication protein A. Proc. Natl. Acad. Sci. USA 2005, 102, 15412-15417.

Structure of the Tfb1/p53 complex: Insights into the interaction between the p62/Tfb1 subunit of TFIIH and the activation domain of p53. P Di Lello, L M Jenkins, T N Jones, B D Nguyen, T Hara, H Yamaguchi, J D Dikeakos, E Appella, P Legault, J G Omichinski, Mol. Cell. 22Di Lello, P.; Jenkins, L.M.; Jones, T.N.; Nguyen, B.D.; Hara, T.; Yamaguchi, H.; Dikeakos, J.D.; Appella, E.; Legault, P.; Omichinski, J.G. Structure of the Tfb1/p53 complex: Insights into the interaction between the p62/Tfb1 subunit of TFIIH and the activation domain of p53. Mol. Cell 2006, 22, 731-740.

Crystal structure of SV40 large T-antigen bound to p53: Interplay between a viral oncoprotein and a cellular tumor suppressor. W Lilyestrom, M G Klein, R Zhang, A Joachimiak, X S Chen, Genes Dev. 20Lilyestrom, W.; Klein, M.G.; Zhang, R.; Joachimiak, A.; Chen, X.S. Crystal structure of SV40 large T-antigen bound to p53: Interplay between a viral oncoprotein and a cellular tumor suppressor. Genes Dev. 2006, 20, 2373-2382.

The C-terminus of p53: The more you learn the less you know. J Ahn, C Prives, Nat. Struct. Biol. 8Ahn, J.; Prives, C. The C-terminus of p53: The more you learn the less you know. Nat. Struct. Biol. 2001, 8, 730-732.

Structure of a Sir2 enzyme bound to an acetylated p53 peptide. J L Avalos, I Celic, S Muhammad, M S Cosgrove, J D Boeke, C Wolberger, Mol. Cell. 10Avalos, J.L.; Celic, I.; Muhammad, S.; Cosgrove, M.S.; Boeke, J.D.; Wolberger, C. Structure of a Sir2 enzyme bound to an acetylated p53 peptide. Mol. Cell 2002, 10, 523-535.

. S Mujtaba, Y He, L Zeng, S Yan, O Plotnikova, Sachchidanand, Mujtaba, S.; He, Y.; Zeng, L.; Yan, S.; Plotnikova, O.; Sachchidanand;

Structural mechanism of the bromodomain of the coactivator CBP in p53 transcriptional activation. R Sanchez, N J Zeleznik-Le, Z Ronai, M M Zhou, Mol. Cell. 13Sanchez, R.; Zeleznik-Le, N.J.; Ronai, Z.; Zhou, M.M. Structural mechanism of the bromodomain of the coactivator CBP in p53 transcriptional activation. Mol. Cell 2004, 13, 251-263.

Regulation of p53 activity through lysine methylation. S Chuikov, J K Kurash, J R Wilson, B Xiao, N Justin, G S Ivanov, K Mckinney, P Tempst, C Prives, S J Gamblin, N A Barlev, D Reinberg, Nature. 432Chuikov, S.; Kurash, J.K.; Wilson, J.R.; Xiao, B.; Justin, N.; Ivanov, G.S.; McKinney, K.; Tempst, P.; Prives, C.; Gamblin, S.J.; Barlev, N.A.; Reinberg, D. Regulation of p53 activity through lysine methylation. Nature 2004, 432, 353-560.

Specificity determinants of recruitment peptides bound to phospho-CDK2/cyclin A. E D Lowe, I Tews, K Y Cheng, N R Brown, S Gul, M E Noble, S J Gamblin, L N Johnson, Biochemistry. 41Lowe, E.D.; Tews, I.; Cheng, K.Y.; Brown, N.R.; Gul, S.; Noble, M.E.; Gamblin, S.J.; Johnson, L.N. Specificity determinants of recruitment peptides bound to phospho-CDK2/cyclin A. Biochemistry 2002, 41, 15625-15634.

. Int. J. Mol. Sci. 201213Int. J. Mol. Sci. 2012, 13

Discovery and cocrystal structure of benzodiazepinedione HDM2 antagonists that activate p53 in cells. B L Grasberger, T Lu, C Schubert, D J Parks, T E Carver, H K Koblish, M D Cummings, L V Lafrance, K L Milkiewicz, R R Calvo, J. Med. Chem. 48Grasberger, B.L.; Lu, T.; Schubert, C.; Parks, D.J.; Carver, T.E.; Koblish, H.K.; Cummings, M.D.; LaFrance, L.V.; Milkiewicz, K.L.; Calvo, R.R.; et al. Discovery and cocrystal structure of benzodiazepinedione HDM2 antagonists that activate p53 in cells. J. Med. Chem. 2005, 48, 909-912.

Structure of the stapled p53 peptide bound to Mdm2. S Baek, P S Kutchukian, G L Verdine, R Huber, T A Holak, K W Lee, G M Popowicz, J. Am. Chem. Soc. 134Baek, S.; Kutchukian, P.S.; Verdine, G.L.; Huber, R.; Holak, T.A.; Lee, K.W.; Popowicz, G.M. Structure of the stapled p53 peptide bound to Mdm2. J. Am. Chem. Soc. 2012, 134, 103-106.

The Ca( 2+ )-dependent interaction of S100B(betabeta) with a peptide derived from p53. R R Rustandi, A C Drohat, D M Baldisseri, P T Wilder, D J Weber, Biochemistry. 37Rustandi, R.R.; Drohat, A.C.; Baldisseri, D.M.; Wilder, P.T.; Weber, D.J. The Ca( 2+ )-dependent interaction of S100B(betabeta) with a peptide derived from p53. Biochemistry 1998, 37, 1951-1960.

Calcium-dependent interaction of S100B with the C-terminal domain of the tumor suppressor p53. C Delphin, M Ronjat, J C Deloulme, G Garin, L Debussche, Y Higashimoto, K Sakaguchi, J Baudier, J. Biol. Chem. 274Delphin, C.; Ronjat, M.; Deloulme, J.C.; Garin, G.; Debussche, L.; Higashimoto, Y.; Sakaguchi, K.; Baudier, J. Calcium-dependent interaction of S100B with the C-terminal domain of the tumor suppressor p53. J. Biol. Chem. 1999, 274, 10539-10544.

A search for inhibitors of S100B, a member of the S100 family of calcium-binding proteins. J Markowitz, A D Mackerell, Jr, D J Weber, Mini Rev. Med. Chem. 7Markowitz, J.; MacKerell, A.D., Jr.; Weber, D.J. A search for inhibitors of S100B, a member of the S100 family of calcium-binding proteins. Mini Rev. Med. Chem. 2007, 7, 609-616.

Insights into drug resistance of mutations D30N and I50V to HIV-1 protease inhibitor TMC-114: Free energy calculation and molecular dynamic simulation. J Chen, S Zhang, X Liu, Q Zhang, Chen, J.; Zhang, S.; Liu, X.; Zhang, Q. Insights into drug resistance of mutations D30N and I50V to HIV-1 protease inhibitor TMC-114: Free energy calculation and molecular dynamic simulation.

. J. Mol. Model. 16J. Mol. Model. 2010, 16, 459-468.

Mdm2 promotes the rapid degradation of p53. Y Haupt, R Maya, A Kazaz, M Oren, Nature. 387Haupt, Y.; Maya, R.; Kazaz, A.; Oren, M. Mdm2 promotes the rapid degradation of p53. Nature 1997, 387, 296-299.

MDM2-arbiter of p53's destruction. D P Lane, P Hall, Trends Biochem. Sci. 22Lane, D.P.; Hall, P.A. MDM2-arbiter of p53's destruction. Trends Biochem. Sci. 1997, 22, 372-374.

The MDM2-p53 interaction. U M Moll, O Petrenko, Mol. Cancer Res. 1Moll, U.M.; Petrenko, O. The MDM2-p53 interaction. Mol. Cancer Res. 2003, 1, 1001-1008.

p53 ubiquitination: Mdm2 and beyond. C L Brooks, W Gu, Mol. Cell. 21Brooks, C.L.; Gu, W. p53 ubiquitination: Mdm2 and beyond. Mol. Cell 2006, 21, 307-315.

MDMX: A novel p53-binding protein with some functional properties of MDM2. A Shvarts, W T Steegenga, N Riteco, T Van Laar, P Dekker, M Bazuine, R C Van Ham, W Van Der Houven Van Oordt, G Hateboer, A J Van Der Eb, A G Jochemsen, EMBO J. 15Shvarts, A.; Steegenga, W.T.; Riteco, N.; van Laar, T.; Dekker, P.; Bazuine, M.; van Ham, R.C.; van der Houven van Oordt, W.; Hateboer, G.; van der Eb, A.J.; Jochemsen, A.G. MDMX: A novel p53-binding protein with some functional properties of MDM2. EMBO J. 1996, 15, 5349-5357.

Mitochondrial MDM4 (MDMX): An unpredicted role in the p53-mediated intrinsic apoptotic pathway. F Mancini, F Moretti, Cell Cycle. 8Mancini, F.; Moretti, F. Mitochondrial MDM4 (MDMX): An unpredicted role in the p53-mediated intrinsic apoptotic pathway. Cell Cycle 2009, 8, 3854-3859.

Regulating the p53 pathway: In vitro hypotheses, in vivo veritas. F Toledo, G M Wahl, Nat. Rev. Cancer. 6Toledo, F.; Wahl, G.M. Regulating the p53 pathway: In vitro hypotheses, in vivo veritas. Nat. Rev. Cancer 2006, 6, 909-923.

Reactivation of p53: From peptides to small molecules. C J Brown, C F Cheok, C S Verma, D P Lane, Trends Pharmacol. Sci. 32Brown, C.J.; Cheok, C.F.; Verma, C.S.; Lane, D.P. Reactivation of p53: From peptides to small molecules. Trends Pharmacol. Sci. 2011, 32, 53-62.

p53 activation by blocking Snail: A novel pharmacological strategy for cancer. S H Lee, B J Park, Curr. Pharm. Des. 17Lee, S.H.; Park, B.J. p53 activation by blocking Snail: A novel pharmacological strategy for cancer. Curr. Pharm. Des. 2011, 17, 610-617.

Study of MDM2 binding to p53-analogues: Affinity, helicity, and applicability to drug design. O Kalid, N Ben-Tal, J. Chem. Inf. Model. 49Kalid, O.; Ben-Tal, N. Study of MDM2 binding to p53-analogues: Affinity, helicity, and applicability to drug design. J. Chem. Inf. Model. 2009, 49, 865-876.

Leucine-rich hydrophobic clusters promote folding of the N-terminus of the intrinsically disordered transactivation domain of p53. L M Espinoza-Fonseca, FEBS Lett. 583Espinoza-Fonseca, L.M. Leucine-rich hydrophobic clusters promote folding of the N-terminus of the intrinsically disordered transactivation domain of p53. FEBS Lett. 2009, 583, 556-560.

Multiple peptide conformations give rise to similar binding affinities: Molecular simulations of p53-MDM2. S G Dastidar, D P Lane, C S Verma, J. Am. Chem. Soc. 130Dastidar, S.G.; Lane, D.P.; Verma, C.S. Multiple peptide conformations give rise to similar binding affinities: Molecular simulations of p53-MDM2. J. Am. Chem. Soc. 2008, 130, 13514-13515.

Targeting the conformational transitions of MDM2 and MDMX: Insights into dissimilarities and similarities of p53 recognition. A Macchiarulo, N Giacchè, A Carotti, M Baroni, G Cruciani, R Pellicciari, J. Chem. Inf. Model. 48Macchiarulo, A.; Giacchè, N.; Carotti, A.; Baroni, M.; Cruciani, G.; Pellicciari, R. Targeting the conformational transitions of MDM2 and MDMX: Insights into dissimilarities and similarities of p53 recognition. J. Chem. Inf. Model. 2008, 48, 1999-2009.

Molecular modeling approaches in the discovery of new drugs for anti-cancer therapy: The investigation of p53-MDM2 interaction and its inhibition by small molecules. A Lauria, M Tutone, M Ippolito, L Pantano, A M Almerico, Curr. Med. Chem. 17Lauria, A.; Tutone, M.; Ippolito, M.; Pantano, L.; Almerico, A.M. Molecular modeling approaches in the discovery of new drugs for anti-cancer therapy: The investigation of p53-MDM2 interaction and its inhibition by small molecules. Curr. Med. Chem. 2010, 17, 3142-3154.

Conformational changes of the p53-binding cleft of MDM2 revealed by molecular dynamics simulations. L M Espinoza-Fonseca, J G Trujillo-Ferrara, Biopolymers. 83Espinoza-Fonseca, L.M.; Trujillo-Ferrara, J.G. Conformational changes of the p53-binding cleft of MDM2 revealed by molecular dynamics simulations. Biopolymers 2006, 83, 365-373.

Multiscale ensemble modeling of intrinsically disordered proteins: P53 N-terminal domain. T Terakawa, S Takada, Biophys. J. 101Terakawa, T.; Takada, S. Multiscale ensemble modeling of intrinsically disordered proteins: P53 N-terminal domain. Biophys. J. 2011, 101, 1450-1458.

Direct observations of conformational distributions of intrinsically disordered p53 peptides using UV Raman and explicit solvent simulations. K Xiong, M C Zwier, N S Myshakina, V M Burger, S A Asher, L T Chong, J. Phys. Chem. A. 115Xiong, K.; Zwier, M.C.; Myshakina, N.S.; Burger, V.M.; Asher, S.A.; Chong, L.T. Direct observations of conformational distributions of intrinsically disordered p53 peptides using UV Raman and explicit solvent simulations. J. Phys. Chem. A 2011, 115, 9520-9527.

Transient stability of the helical pattern of region F19-L22 of the N-terminal domain of p53: A molecular dynamics simulation study. L M Espinoza-Fonseca, . J Trujillo-Ferrara, Biochem. Biophys. Res. Commun. 343Espinoza-Fonseca, L.M.; Trujillo-Ferrara. J.G. Transient stability of the helical pattern of region F19-L22 of the N-terminal domain of p53: A molecular dynamics simulation study. Biochem. Biophys. Res. Commun. 2006, 343, 110-116.

Differences in the transactivation domains of p53 family members: A computational study. J N Mavinahalli, A Madhumalar, R W Beuerman, D P Lane, C Verma, 10.1186/1471-2164-11-S1-S5BMC Genomics. 11Mavinahalli, J.N.; Madhumalar, A.; Beuerman, R.W.; Lane, D.P.; Verma, C. Differences in the transactivation domains of p53 family members: A computational study. BMC Genomics 2010, 11, doi:10.1186/1471-2164-11-S1-S5.

Thermodynamics of p53 binding to hdm2(1-126): Effects of phosphorylation and p53 peptide length. Z Lai, K R Auger, C M Manubay, R A Copeland, Arch. Biochem. Biophys. 381Lai, Z.; Auger, K.R.; Manubay, C.M.; Copeland, R.A. Thermodynamics of p53 binding to hdm2(1-126): Effects of phosphorylation and p53 peptide length. Arch. Biochem. Biophys. 2000, 381, 278-284.

Modulation of the p53-MDM2 interaction by phosphorylation of Thr18: A computational study. H J Lee, D Srinivasan, D Coomber, D P Lane, C S Verma, Cell Cycle. 6Lee, H.J.; Srinivasan, D.; Coomber, D.; Lane, D.P.; Verma, C.S. Modulation of the p53-MDM2 interaction by phosphorylation of Thr18: A computational study. Cell Cycle 2007, 6, 2604-2611.

The electrostatic surface of MDM2 modulates the specificity of its interaction with phosphorylated and unphosphorylated p53 peptides. C J Brown, D Srinivasan, L H Jun, D Coomber, C S Verma, D P Lane, Cell Cycle. 7Brown, C.J.; Srinivasan, D.; Jun, L.H.; Coomber, D.; Verma, C.S.; Lane, D.P. The electrostatic surface of MDM2 modulates the specificity of its interaction with phosphorylated and unphosphorylated p53 peptides. Cell Cycle 2008, 7, 608-610.

Modulation of p53 binding to MDM2: Computational studies reveal important roles of Tyr100. S G Dastidar, D P Lane, C S Verma, 10.1186/1471-2105-10-S15-S6BMC Bioinforma. 10Dastidar, S.G.; Lane, D.P.; Verma, C.S. Modulation of p53 binding to MDM2: Computational studies reveal important roles of Tyr100. BMC Bioinforma. 2009, 10, doi:10.1186/1471-2105- 10-S15-S6.

Forces mediating protein-protein interactions: A computational study of p53 "approaching" MDM2. S G Dastidar, A Madhumalar, G Fuentes, D P Lane, C S Verma, Theor. Chem. Acc. 125Dastidar, S.G.; Madhumalar, A.; Fuentes, G.; Lane, D.P.; Verma, C.S. Forces mediating protein-protein interactions: A computational study of p53 "approaching" MDM2. Theor. Chem. Acc. 2010, 125, 621-635.

Targeting the conformational transitions of MDM2 and MDMX: Insights into key residues affecting p53 recognition. A Carotti, A Macchiarulo, N Giacchè, R Pellicciari, Proteins. 77Carotti, A.; Macchiarulo, A.; Giacchè, N.; Pellicciari, R. Targeting the conformational transitions of MDM2 and MDMX: Insights into key residues affecting p53 recognition. Proteins 2009, 77, 524-535.

C-terminal substitution of MDM2 interacting peptides modulates binding affinity by distinctive mechanisms. C J Brown, S G Dastidar, S T Quah, A Lim, B Chia, C S Verma, 10.1371/journal.pone.0024122PLoS One. 6Brown, C.J.; Dastidar, S.G.; Quah, S.T.; Lim, A.; Chia, B.; Verma, C.S. C-terminal substitution of MDM2 interacting peptides modulates binding affinity by distinctive mechanisms. PLoS One 2011, 6, doi:10.1371/journal.pone.0024122.

Chemical states of the N-terminal "lid" of MDM2 regulate p53 binding: Simulations reveal complexities of modulation. S G Dastidar, D Raghunathan, J Nicholson, T R Hupp, D P Lane, C S Verma, Cell Cycle. 10Dastidar, S.G.; Raghunathan, D.; Nicholson, J.; Hupp, T.R.; Lane, D.P.; Verma, C.S. Chemical states of the N-terminal "lid" of MDM2 regulate p53 binding: Simulations reveal complexities of modulation. Cell Cycle 2011, 10, 82-89.

Molecular Basis of Bcl-XL-p53 Interaction: Insights from Molecular Dynamics Simulations. N Bharatham, S-W Chi, H S Yoon, 10.1371/journal.pone.0026014PLoS One. 6Bharatham, N.; Chi, S-W.; Yoon, H.S. Molecular Basis of Bcl-XL-p53 Interaction: Insights from Molecular Dynamics Simulations. PLoS One 2011, 6, doi:10.1371/journal.pone.0026014.

Modeling the interaction between the N-terminal domain of the tumor suppressor p53 and azurin. M Taranta, A R Bizzarri, S Cannistraro, J. Mol. Recognit. 22Taranta, M.; Bizzarri, A.R.; Cannistraro, S. Modeling the interaction between the N-terminal domain of the tumor suppressor p53 and azurin. J. Mol. Recognit. 2009, 22, 215-222.

Effect of Zn 2+ on DNA recognition and stability of the p53 DNA-binding domain. J Duan, L Nilsson, Biochemistry. 45Duan, J.; Nilsson, L. Effect of Zn 2+ on DNA recognition and stability of the p53 DNA-binding domain. Biochemistry 2006, 45, 7483-7492.

Predicting the coordination geometry for Mg 2+ in the p53 DNA-binding domain: Insights from computational studies. T Wang, X Shao, W Cai, Y Xue, S Wang, X Feng, Phys. Chem. Chem. Phys. 13Wang, T.; Shao, X.; Cai, W.; Xue, Y.; Wang, S.; Feng, X. Predicting the coordination geometry for Mg 2+ in the p53 DNA-binding domain: Insights from computational studies. Phys. Chem. Chem. Phys. 2011, 13, 1140-1151.

Docking study and free energy simulation of the complex between p53 DNA-binding domain and azurin. V De Grandis, A R Bizzarri, S Cannistraro, J. Mol. Recognit. 20De Grandis, V.; Bizzarri, A.R.; Cannistraro, S. Docking study and free energy simulation of the complex between p53 DNA-binding domain and azurin. J. Mol. Recognit. 2007, 20, 215-226.

Modelling the interaction between the p53 DNA-binding domain and the p28 peptide fragment of Azurin. S Santini, A R Bizzarri, S Cannistraro, J. Mol. Recognit. 24Santini, S.; Bizzarri, A.R.; Cannistraro, S. Modelling the interaction between the p53 DNA-binding domain and the p28 peptide fragment of Azurin. J. Mol. Recognit. 2011, 24, 1043-1055.

The role of tetramerization in p53 function. P Chène, Oncogene. 20Chène, P. The role of tetramerization in p53 function. Oncogene 2001, 20, 2611-2617.

Crystal structure of the tetramerization domain of the p53 tumor suppressor at 1.7 angstroms. P D Jeffrey, S Gorina, N P Pavletich, Science. 267Jeffrey, P.D.; Gorina, S.; Pavletich, N.P. Crystal structure of the tetramerization domain of the p53 tumor suppressor at 1.7 angstroms. Science 1995, 267, 1498-1502.

Crystallization and structure solution of p53 (residues 326-356) by molecular replacement using an NMR model as template. P R Mittl, P Chène, M G Grütter, Acta Crystallogr. D. 54Mittl, P.R.; Chène, P.; Grütter, M.G. Crystallization and structure solution of p53 (residues 326-356) by molecular replacement using an NMR model as template. Acta Crystallogr. D 1998, 54, 86-89.

Dimerization of the p53 oligomerization domain: Identification of a folding nucleus by molecular dynamics simulations. L T Chong, C D Snow, Y M Rhee, V S Pande, J. Mol. Biol. 345Chong, L.T.; Snow, C.D.; Rhee, Y.M.; Pande, V.S. Dimerization of the p53 oligomerization domain: Identification of a folding nucleus by molecular dynamics simulations. J. Mol. Biol. 2005, 345, 869-878.

Thermal unfolding simulations of a multimeric protein-transition state and unfolding pathways. J Duan, L Nilsson, Proteins. 59Duan, J.; Nilsson, L. Thermal unfolding simulations of a multimeric protein-transition state and unfolding pathways. Proteins 2005, 59, 170-182.

Disruption of an intermonomer salt bridge in the p53 tetramerization domain results in an increased propensity to form amyloid fibrils. C Galea, P Bowman, R W Kriwacki, Protein Sci. 14Galea, C.; Bowman, P.; Kriwacki, R.W. Disruption of an intermonomer salt bridge in the p53 tetramerization domain results in an increased propensity to form amyloid fibrils. Protein Sci. 2005, 14, 2993-3003.

Folding of tetrameric p53: Oligomerization and tumorigenic mutations induce misfolding and loss of function. D J Lubin, J S Butler, S N Loh, J. Mol. Biol. 395Lubin, D.J.; Butler, J.S.; Loh, S.N. Folding of tetrameric p53: Oligomerization and tumorigenic mutations induce misfolding and loss of function. J. Mol. Biol. 2010, 395, 705-716.

A fluid salt-bridging cluster and the stabilization of p53. T Z Lwin, J J Durant, D Bashford, J. Mol. Biol. 373Lwin, T.Z.; Durant, J.J.; Bashford, D. A fluid salt-bridging cluster and the stabilization of p53. J. Mol. Biol. 2007, 373, 1334-1347.

Characterization of the tumor suppressor protein p53 as a protein kinase C substrate and a S100b-binding protein. J Baudier, C Delphin, D Grunwald, S Khochbin, J J Lawrence, Baudier, J.; Delphin, C.; Grunwald, D.; Khochbin, S.; Lawrence, J.J. Characterization of the tumor suppressor protein p53 as a protein kinase C substrate and a S100b-binding protein.

Proc. Natl. Acad. Sci. Natl. Acad. SciUSA89Proc. Natl. Acad. Sci. USA 1992, 89, 11627-11631.

Proteins of the S100 family regulate the oligomerization of p53 tumor suppressor. M R Fernandez-Fernandez, D B Veprintsev, A R Fersht, Proc. Natl. Acad. Sci. Natl. Acad. SciUSA102Fernandez-Fernandez, M.R.; Veprintsev, D.B.; Fersht, A.R. Proteins of the S100 family regulate the oligomerization of p53 tumor suppressor. Proc. Natl. Acad. Sci. USA 2005, 102, 4735-4740.

Members of the S100 family bind p53 in two distinct ways. M R Fernandez-Fernandez, T J Rutherford, A R Fersht, Protein Sci. 17Fernandez-Fernandez, M.R.; Rutherford, T.J.; Fersht, A.R. Members of the S100 family bind p53 in two distinct ways. Protein Sci. 2008, 17, 1663-1670.

Modulation of the oligomerization state of p53 by differential binding of proteins of the S100 family to p53 monomers and tetramers. J Van Dieck, M R Fernandez-Fernandez, D B Veprintsev, A R Fersht, J. Biol. Chem. 284Van Dieck, J.; Fernandez-Fernandez, M.R.; Veprintsev, D.B.; Fersht, A.R. Modulation of the oligomerization state of p53 by differential binding of proteins of the S100 family to p53 monomers and tetramers. J. Biol. Chem. 2009, 284, 13804-13811.

Structural changes in the C-terminus of Ca 2+ -bound rat S100B (beta beta) upon binding to a peptide derived from the C-terminal regulatory domain of p53. R R Rustandi, D M Baldisseri, A C Drohat, D J Weber, Protein Sci. 8Rustandi, R.R.; Baldisseri, D.M.; Drohat, A.C.; Weber, D.J. Structural changes in the C-terminus of Ca 2+ -bound rat S100B (beta beta) upon binding to a peptide derived from the C-terminal regulatory domain of p53. Protein Sci. 1999, 8, 1743-1751.

Inhibition of p53 transcriptional activity by the S100B calcium-binding protein. J Lin, M Blake, C Tang, D Zimmer, R R Rustandi, D J Weber, F Carrier, J. Biol. Chem. 276Lin, J.; Blake, M.; Tang, C.; Zimmer, D.; Rustandi, R.R.; Weber, D.J.; Carrier, F. Inhibition of p53 transcriptional activity by the S100B calcium-binding protein. J. Biol. Chem. 2001, 276, 35037-35041.

Intrinsically disordered p53 extreme C-terminus binds to S100B(betabeta) through "fly-casting. J Chen, J. Am. Chem. Soc. 131Chen, J. Intrinsically disordered p53 extreme C-terminus binds to S100B(betabeta) through "fly-casting". J. Am. Chem. Soc. 2009, 131, 2088-2089.

Recognition of the tumor suppressor protein p53 and other protein targets by the calcium-binding protein S100B. P T Wilder, J Lin, C L Bair, T H Charpentier, D Yang, M Liriano, K M Varney, A Lee, A B Oppenheim, S Adhya, F Carrier, D J Weber, Biochim. Biophys. Acta. 1763Wilder, P.T.; Lin, J.; Bair, C.L.; Charpentier, T.H.; Yang, D.; Liriano, M.; Varney, K.M.; Lee, A.; Oppenheim, A.B.; Adhya, S.; Carrier, F.; Weber, D.J. Recognition of the tumor suppressor protein p53 and other protein targets by the calcium-binding protein S100B. Biochim. Biophys. Acta 2006, 1763, 1284-1297.

Solution structure of calcium-bound rat S100B(betabeta) as determined by nuclear magnetic resonance spectroscopy. A C Drohat, D M Baldisseri, R R Rustandi, D J Weber, Biochemistry. 37Drohat, A.C.; Baldisseri, D.M.; Rustandi, R.R.; Weber, D.J. Solution structure of calcium-bound rat S100B(betabeta) as determined by nuclear magnetic resonance spectroscopy. Biochemistry 1998, 37, 2729-2740.

S100-Beta, melanoma-inhibiting activity, and lactate dehydrogenase discriminate progressive from nonprogressive American Joint Committee on Cancer stage IV melanoma. M Deichmann, A Benner, M Bock, A Jäckel, K Uhl, V Waldmann, H Näher, J. Clin. Oncol. 17Deichmann, M.; Benner, A.; Bock, M.; Jäckel, A.; Uhl, K.; Waldmann, V.; Näher, H. S100-Beta, melanoma-inhibiting activity, and lactate dehydrogenase discriminate progressive from nonprogressive American Joint Committee on Cancer stage IV melanoma. J. Clin. Oncol. 1999, 17, 1891-1896.

S-100 beta protein in serum, a tumor marker in malignant melanoma-current state of knowledge and clinical experience. A Jäckel, M Deichmann, V Waldmann, M Bock, H Näher, Hautarzt. 50Jäckel, A.; Deichmann, M.; Waldmann, V.; Bock, M.; Näher, H. S-100 beta protein in serum, a tumor marker in malignant melanoma-current state of knowledge and clinical experience. Hautarzt 1999, 50, 250-256.

S100 proteins as cancer biomarkers with focus on S100B in malignant melanoma. R Harpio, R Einarsson, Clin. Biochem. 37Harpio, R.; Einarsson, R. S100 proteins as cancer biomarkers with focus on S100B in malignant melanoma. Clin. Biochem. 2004, 37, 512-518.

Theoretical study on binding of S100B protein. A Gieldon, M Mori, R Del Conte, J. Mol. Model. 13Gieldon, A.; Mori, M.; Del Conte, R. Theoretical study on binding of S100B protein. J. Mol. Model. 2007, 13, 1123-1131.

Modeling the relationship between the p53 C-terminal domain and its binding partners using molecular dynamics. W J Allen, D G Capelluto, C V Finkielstein, D R Bevan, J. Phys. Chem. B. 114Allen, W.J.; Capelluto, D.G.; Finkielstein, C.V.; Bevan, D.R. Modeling the relationship between the p53 C-terminal domain and its binding partners using molecular dynamics. J. Phys. Chem. B 2010, 114, 13201-13213.

Insight into a novel p53 single point mutation (G389E) by molecular dynamics simulations. D Pirolli, C C Alinovi, E Capoluongo, M A Satta, P Concolino, B Giardina, M C De Rosa, Int. J. Mol. Sci. 12Pirolli, D.; Alinovi, C.C.; Capoluongo, E.; Satta, M.A.; Concolino, P.; Giardina, B.; de Rosa, M.C. Insight into a novel p53 single point mutation (G389E) by molecular dynamics simulations. Int. J. Mol. Sci. 2010, 12, 128-140.

Acetylation is indispensable for p53 activation. Y Tang, W Zhao, Y Chen, Y Zhao, W Gu, Cell. 133Tang, Y.; Zhao, W.; Chen, Y.; Zhao, Y.; Gu, W. Acetylation is indispensable for p53 activation. Cell 2008, 133, 612-626.

CREB-binding protein and p300/CBP-associated factor are transcriptional coactivators of the p53 tumor suppressor protein. D M Scolnick, N H Chehab, E S Stavridi, M C Lien, L Caruso, E Moran, S L Berger, T D Halazonetis, Cancer Res. 57Scolnick, D.M.; Chehab, N.H.; Stavridi, E.S.; Lien, M.C.; Caruso, L.; Moran, E.; Berger, S.L.; Halazonetis, T.D. CREB-binding protein and p300/CBP-associated factor are transcriptional coactivators of the p53 tumor suppressor protein. Cancer Res. 1997, 57, 3693-3696.

p53 Transcriptional activity is mediated through the SRC1-interacting domain of CBP/p300. J A Livengood, K E Scoggin, K Van Orden, S J Mcbryant, R S Edayathumangalam, P J Laybourn, J K Nyborg, J. Biol. Chem. 277Livengood, J.A.; Scoggin, K.E.; Van Orden, K.; McBryant, S.J.; Edayathumangalam, R.S.; Laybourn, P.J.; Nyborg, J.K. p53 Transcriptional activity is mediated through the SRC1-interacting domain of CBP/p300. J. Biol. Chem. 2002, 277, 9054-9061.

Transcriptional regulation of tumor suppressor p53 by cAMP-responsive element-binding protein/AMP-activated protein kinase complex in response to glucose deprivation. R Okoshi, K Ando, Y Suenaga, M Sang, N Kubo, H Kizaki, A Nakagawara, T Ozaki, Genes Cells. 14Okoshi, R.; Ando, K.; Suenaga, Y.; Sang, M.; Kubo, N.; Kizaki, H.; Nakagawara, A.; Ozaki, T. Transcriptional regulation of tumor suppressor p53 by cAMP-responsive element-binding protein/AMP-activated protein kinase complex in response to glucose deprivation. Genes Cells 2009, 14, 1429-1440.

Target structure-based discovery of small molecules that block human p53 and CREB binding protein association. ; Sachchidanand, L Resnick-Silverman, S Yan, S Mutjaba, W J Liu, L Zeng, J J Manfredi, M M Zhou, Chem. Biol. 13Sachchidanand; Resnick-Silverman, L.; Yan, S.; Mutjaba, S.; Liu, W.J.; Zeng, L.; Manfredi, J.J.; Zhou, M.M. Target structure-based discovery of small molecules that block human p53 and CREB binding protein association. Chem. Biol. 2006, 13, 81-90.

A small molecule binding to the coactivator CREB-binding protein blocks apoptosis in cardiomyocytes. J C Borah, S Mujtaba, I Karakikes, L Zeng, M Muller, J Patel, N Moshkina, K Morohashi, W Zhang, G Gerona-Navarro, R J Hajjar, M M Zhou, Chem. Biol. 18Borah, J.C.; Mujtaba, S.; Karakikes, I.; Zeng, L.; Muller, M.; Patel, J.; Moshkina, N.; Morohashi, K.; Zhang, W.; Gerona-Navarro, G.; Hajjar, R.J.; Zhou, M.M. A small molecule binding to the coactivator CREB-binding protein blocks apoptosis in cardiomyocytes. Chem. Biol. 2011, 18, 531-541.

Yoel-Rodríguez. G Gerona-Navarro, Gerona-Navarro, G.; Yoel-Rodríguez;

Rational design of cyclic peptide modulators of the transcriptional coactivator CBP: A new class of p53 inhibitors. S Mujtaba, A Frasca, J Patel, L Zeng, A N Plotnikov, R Osman, M M Zhou, J. Am. Chem. Soc. 133Mujtaba, S.; Frasca, A.; Patel, J.; Zeng, L.; Plotnikov, A.N.; Osman, R.; Zhou, M.M. Rational design of cyclic peptide modulators of the transcriptional coactivator CBP: A new class of p53 inhibitors. J. Am. Chem. Soc. 2011, 133, 2040-2043.

The energetics of the acetylation switch in p53-mediated transcriptional activation. K D Eichenbaum, Y Rodríguez, M Mezei, R Osman, Proteins. 78Eichenbaum, K.D.; Rodríguez, Y.; Mezei, M.; Osman, R. The energetics of the acetylation switch in p53-mediated transcriptional activation. Proteins 2010, 78, 447-456.

Polybromo-1-bromodomains bind histone H3 at specific acetyl-lysine positions. R Chandrasekaran, M Thompson, Biochem. Biophys. Res. Commun. 355Chandrasekaran, R.; Thompson, M. Polybromo-1-bromodomains bind histone H3 at specific acetyl-lysine positions. Biochem. Biophys. Res. Commun. 2007, 355, 661-666.

Sirtuins-novel therapeutic targets to treat age-associated diseases. S Lavu, O Boss, P J Elliott, P D Lambert, Nat. Rev. Drug Discov. 7Lavu, S.; Boss, O.; Elliott, P.J.; Lambert, P.D. Sirtuins-novel therapeutic targets to treat age-associated diseases. Nat. Rev. Drug Discov. 2008, 7, 841-853.

The Sirtuin family: Therapeutic targets to treat diseases of aging. J C Milne, J M Denu, Curr. Opin. Chem. Biol. 12Milne, J.C.; Denu, J.M. The Sirtuin family: Therapeutic targets to treat diseases of aging. Curr. Opin. Chem. Biol. 2008, 12, 11-17.

A. hSIR2(SIRT1) functions as an NAD-dependent p53 deacetylase. H Vaziri, S K Dessain, E N Eaton, S I Imai, R A Frye, T K Pandita, L Guarente, R Weinberg, Cell. 107Vaziri, H.; Dessain, S.K.; Eaton, E.N.; Imai, S.I.; Frye, R.A.; Pandita, T.K.; Guarente, L.; Weinberg, R.A. hSIR2(SIRT1) functions as an NAD-dependent p53 deacetylase. Cell 2001, 107, 149-159.

Negative control of p53 by Sir2alpha promotes cell survival under stress. J Luo, A Y Nikolaev, S Imai, D Chen, F Su, A Shiloh, L Guarente, W Gu, Cell. 107Luo, J.; Nikolaev, A.Y.; Imai, S.; Chen, D.; Su, F.; Shiloh, A.; Guarente, L.; Gu, W. Negative control of p53 by Sir2alpha promotes cell survival under stress. Cell 2001, 107, 137-148.

Human SIR2 deacetylates p53 and antagonizes PML/p53-induced cellular senescence. E Langley, M Pearson, M Faretta, U M Bauer, R A Frye, S Minucci, P G Pelicci, T Kouzarides, EMBO J. 21Langley, E.; Pearson, M.; Faretta, M.; Bauer, U.M.; Frye, R.A.; Minucci, S.; Pelicci, P.G.; Kouzarides, T. Human SIR2 deacetylates p53 and antagonizes PML/p53-induced cellular senescence. EMBO J. 2002, 21, 2383-2396.

Human Sir2 and the 'silencing' of p53 activity. J Smith, Trends Cell Biol. 12Smith, J. Human Sir2 and the 'silencing' of p53 activity. Trends Cell Biol. 2002, 12, 404-406.

Methylation-acetylation interplay activates p53 in response to DNA damage. G S Ivanov, T Ivanova, J Kurash, A Ivanov, S Chuikov, F Gizatullin, E M Herrera-Medina, F Rauscher, D Reinberg, N A Barlev, Mol. Cell. Biol. 27Ivanov, G.S.; Ivanova, T.; Kurash, J.; Ivanov, A.; Chuikov, S.; Gizatullin, F.; Herrera-Medina, E.M.; Rauscher, F.; Reinberg, D.; Barlev, N.A. Methylation-acetylation interplay activates p53 in response to DNA damage. Mol. Cell. Biol. 2007, 27, 6756-6769.

Mechanism of product specificity of AdoMet methylation catalyzed by lysine methyltransferases: Transcriptional factor p53 methylation by histone lysine methyltransferase SET7/9. X Zhang, T C Bruice, Biochemistry. 47Zhang, X.; Bruice, T.C. Mechanism of product specificity of AdoMet methylation catalyzed by lysine methyltransferases: Transcriptional factor p53 methylation by histone lysine methyltransferase SET7/9. Biochemistry 2008, 47, 2743-2748.

Modeling a new water channel that allows SET9 to dimethylate p53. Q Bai, Y Shen, X Yao, F Wang, Y Du, Q Wang, N Jin, J Hai, T Hu, J Yang, 10.1371/journal.pone.0019856PLoS One. 6Bai, Q.; Shen, Y.; Yao, X.; Wang, F.; Du, Y.; Wang, Q.; Jin, N.; Hai, J.; Hu, T.; Yang, J. Modeling a new water channel that allows SET9 to dimethylate p53. PLoS One 2011, 6, doi:10.1371/journal.pone.0019856.

Quantum mechanical studies of residue-specific hydrophobic interactions in p53-MDM2 binding. Y Ding, Y Mei, J Z Zhang, J. Phys. Chem. B. 112Ding, Y.; Mei, Y.; Zhang, J.Z. Quantum mechanical studies of residue-specific hydrophobic interactions in p53-MDM2 binding. J. Phys. Chem. B 2008, 112, 11396-11401.

A computational analysis of the binding model of MDM2 with inhibitors. G Hu, D Wang, X Liu, Q Zhang, J. Comput. Aided Mol. Des. 24Hu, G.; Wang, D.; Liu, X.; Zhang, Q. A computational analysis of the binding model of MDM2 with inhibitors. J. Comput. Aided Mol. Des. 2010, 24, 687-697.

Structure-based design of spiro-oxindoles as potent, specific small-molecule inhibitors of the MDM2-p53 interaction. K Ding, Y Lu, Z Nikolovska-Coleska, G Wang, S Qiu, S Shangary, W Gao, D Qin, J Stuckey, K Krajewski, P P Roller, S Wang, J. Med. Chem. 49Ding, K.; Lu, Y.; Nikolovska-Coleska, Z.; Wang, G.; Qiu, S.; Shangary, S.; Gao, W.; Qin, D.; Stuckey, J.; Krajewski, K.; Roller, P.P.; Wang, S. Structure-based design of spiro-oxindoles as potent, specific small-molecule inhibitors of the MDM2-p53 interaction. J. Med. Chem. 2006, 49, 3432-3435.

Insight into mechanism of small molecule inhibitors of the MDM2-p53 interaction: Molecular dynamics simulation and free energy analysis. J Chen, J Wang, B Xu, W Zhu, G Li, J. Mol. Graph. Model. 30Chen, J.; Wang, J.; Xu, B.; Zhu, W.; Li, G. Insight into mechanism of small molecule inhibitors of the MDM2-p53 interaction: Molecular dynamics simulation and free energy analysis. J. Mol. Graph. Model. 2011, 30, 46-53.

Structural determinants of benzodiazepinedione/ peptide-based p53-HDM2 inhibitors using 3D-QSAR, docking and molecular dynamics. F Wang, Y Li, Z Ma, X Wang, Y Wang, J. Mol. Model. 18Wang, F.; Li, Y.; Ma, Z.; Wang, X.; Wang, Y. Structural determinants of benzodiazepinedione/ peptide-based p53-HDM2 inhibitors using 3D-QSAR, docking and molecular dynamics. J. Mol. Model. 2012, 18, 295-306.

Systematic mutational analysis of an ubiquitin ligase (MDM2)-binding peptide: Computational studies. Y Liu, D P Lane, C S Verma, Theor. Chem. Acc. 130Liu, Y.; Lane, D.P.; Verma, C.S. Systematic mutational analysis of an ubiquitin ligase (MDM2)-binding peptide: Computational studies. Theor. Chem. Acc. 2011, 130, 1145-1154.

A computational analysis of interaction mechanisms of peptide and non-peptide inhibitors with MDMX based on molecular dynamics simulation. W Cheng, J Chen, Z Liang, G Li, C Yi, W Wang, K Wang, Comp. Theor. Chem. 984Cheng, W.; Chen, J.; Liang, Z.; Li, G.; Yi, C.; Wang, W.; Wang, K. A computational analysis of interaction mechanisms of peptide and non-peptide inhibitors with MDMX based on molecular dynamics simulation. Comp. Theor. Chem. 2012, 984, 43-50.

Ensemble-based virtual screening reveals dual-inhibitors for the p53-MDM2/MDMX interactions. K Barakat, J Mane, D Friesen, J Tuszynski, J. Mol. Graph. Model. 28Barakat, K.; Mane, J.; Friesen, D.; Tuszynski, J. Ensemble-based virtual screening reveals dual-inhibitors for the p53-MDM2/MDMX interactions. J. Mol. Graph. Model. 2010, 28, 555-568.

Differential binding of p53 and nutlin to MDM2 and MDMX: Computational studies. T L Joseph, A Madhumalar, C J Brown, D P Lane, C S Verma, Cell Cycle. 9Joseph, T.L.; Madhumalar, A.; Brown, C.J.; Lane, D.P.; Verma, C.S. Differential binding of p53 and nutlin to MDM2 and MDMX: Computational studies. Cell Cycle 2010, 9, 1167-1181.

Molecular modeling and molecular dynamics simulation studies on pyrrolopyrimidine-based α-helix mimetic as dual inhibitors of MDM2 and MDMX. S Y Lu, Y J Jiang, J W Zou, T X Wu, J. Mol. Graph. Model. 30Lu, S.Y.; Jiang, Y.J.; Zou, J.W.; Wu, T.X. Molecular modeling and molecular dynamics simulation studies on pyrrolopyrimidine-based α-helix mimetic as dual inhibitors of MDM2 and MDMX. J. Mol. Graph. Model. 2011, 30, 167-178.

Computational studies of difference in binding modes of peptide and non-peptide inhibitors to MDM2/MDMX based on molecular dynamics simulations. J Chen, D Zhang, Y Zhang, G Li, Int. J. Mol. Sci. 13Chen, J.; Zhang, D.; Zhang, Y.; Li, G. Computational studies of difference in binding modes of peptide and non-peptide inhibitors to MDM2/MDMX based on molecular dynamics simulations. Int. J. Mol. Sci. 2012, 13, 2176-2195.

Reactivation of the p53 tumor suppressor pathway by a stapled p53 peptide. F Bernal, A F Tyler, S J Korsmeyer, L D Walensky, G L Verdine, J. Am. Chem. Soc. 129Bernal, F.; Tyler, A.F.; Korsmeyer, S.J.; Walensky, L.D.; Verdine, G.L. Reactivation of the p53 tumor suppressor pathway by a stapled p53 peptide. J. Am. Chem. Soc. 2007, 129, 2456-2457.

Synthesis of cell-permeable stapled peptide dual inhibitors of the p53-Mdm2/Mdmx interactions via photoinduced cycloaddition. M M Madden, A Muppidi, Z Li, X Li, J Chen, Q Lin, Bioorg. Med. Chem. Lett. 21Madden, M.M.; Muppidi, A.; Li, Z.; Li, X.; Chen, J.; Lin, Q. Synthesis of cell-permeable stapled peptide dual inhibitors of the p53-Mdm2/Mdmx interactions via photoinduced cycloaddition. Bioorg. Med. Chem. Lett. 2011, 21, 1472-1475.

Conjugation of spermine enhances cellular uptake of the stapled peptide-based inhibitors of p53-Mdm2 interaction. A Muppidi, X Li, J Chen, Q Lin, Bioorg. Med. Chem. Lett. 21Muppidi, A.; Li, X.; Chen, J.; Lin, Q. Conjugation of spermine enhances cellular uptake of the stapled peptide-based inhibitors of p53-Mdm2 interaction. Bioorg. Med. Chem. Lett. 2011, 21, 7412-7415.

Anticancer drugs: Stapled peptide rescues p53. S Crunkhorn, 10.1038/nrd3342Nat. Rev. Drug Discov. 10Crunkhorn, S. Anticancer drugs: Stapled peptide rescues p53. Nat. Rev. Drug Discov. 2011, 10, doi:10.1038/nrd3342.

Probing the alpha-helical structural stability of stapled p53 peptides: Molecular dynamics simulations and analysis. Z Guo, U Mohanty, J Noehre, T K Sawyer, W Sherman, G Krilov, Chem. Biol. Drug Des. 75Guo, Z.; Mohanty, U.; Noehre, J.; Sawyer, T.K.; Sherman, W.; Krilov, G. Probing the alpha-helical structural stability of stapled p53 peptides: Molecular dynamics simulations and analysis. Chem. Biol. Drug Des. 2010, 75, 348-359.

Stapled peptides in the p53 pathway: Computer simulations reveal novel interactions of the staples with the target protein. T L Joseph, D Lane, C S Verma, Cell Cycle. 9Joseph, T.L.; Lane, D.; Verma, C.S. Stapled peptides in the p53 pathway: Computer simulations reveal novel interactions of the staples with the target protein. Cell Cycle 2010, 9, 4560-4568.

Inhibiting S100B restores p53 levels in primary malignant melanoma cancer cells. J Lin, Q Yang, Z Yan, J Markowitz, P T Wilder, F Carrier, D J Weber, J. Biol. Chem. 279Lin, J.; Yang, Q.; Yan, Z.; Markowitz, J.; Wilder, P.T.; Carrier, F.; Weber, D.J. Inhibiting S100B restores p53 levels in primary malignant melanoma cancer cells. J. Biol. Chem. 2004, 279, 34071-34077.

The calcium-binding protein S100B down-regulates p53 and apoptosis in malignant melanoma. J Lin, Q Yang, P T Wilder, F Carrier, D J Weber, J. Biol. Chem. 285Lin, J.; Yang, Q.; Wilder, P.T.; Carrier, F.; Weber, D.J. The calcium-binding protein S100B down-regulates p53 and apoptosis in malignant melanoma. J. Biol. Chem. 2010, 285, 27487-27498.

Identification and characterization of small molecule inhibitors of the calcium-dependent S100B-p53 tumor suppressor interaction. J Markowitz, I Chen, R Gitti, D M Baldisseri, Y Pan, R Udan, F Carrier, A D Mackerell, Jr, D J Weber, J. Med. Chem. 47Markowitz, J.; Chen, I.; Gitti, R.; Baldisseri, D.M.; Pan, Y.; Udan, R.; Carrier, F.; MacKerell, A.D., Jr.; Weber, D.J. Identification and characterization of small molecule inhibitors of the calcium-dependent S100B-p53 tumor suppressor interaction. J. Med. Chem. 2004, 47, 5085-5093.

Design of inhibitors for S100B. J Markowitz, A D Mackerell, Jr, F Carrier, T H Charpentier, D J Weber, Curr. Top. Med. Chem. 5Markowitz, J.; Mackerell, A.D., Jr.; Carrier, F.; Charpentier, T.H.; Weber, D.J. Design of inhibitors for S100B. Curr. Top. Med. Chem. 2005, 5, 1093-1108.

In vitro screening and structural characterization of inhibitors of the S100B-p53 interaction. P T Wilder, T H Charpentier, M A Liriano, K Gianni, K M Varney, E Pozharski, A Coop, E A Toth, A D Mackerell, D J Weber, Int. J. High Throughput Screen. Wilder, P.T.; Charpentier, T.H.; Liriano, M.A.; Gianni, K.; Varney, K.M.; Pozharski, E.; Coop, A.; Toth, E.A.; Mackerell, A.D.; Weber, D.J. In vitro screening and structural characterization of inhibitors of the S100B-p53 interaction. Int. J. High Throughput Screen. 2010, 2010, 109-126.

Small molecules bound to unique sites in the target protein binding cleft of calcium-bound S100B as characterized by nuclear magnetic resonance and X-ray crystallography. T H Charpentier, P T Wilder, M A Liriano, K M Varney, S Zhong, A Coop, E Pozharski, A D Mackerell, Jr, E A Toth, D J Weber, Biochemistry. 48Charpentier, T.H.; Wilder, P.T.; Liriano, M.A.; Varney, K.M.; Zhong, S.; Coop, A.; Pozharski, E.; MacKerell, A.D., Jr.; Toth, E.A.; Weber, D.J. Small molecules bound to unique sites in the target protein binding cleft of calcium-bound S100B as characterized by nuclear magnetic resonance and X-ray crystallography. Biochemistry 2009, 48, 6202-6212.

Computational screening and design of S100B ligand to block S100B-p53 interaction. J L Whitlow, J F Varughese, Z Zhou, L J Bartolotti, Y Li, J. Mol. Graph. Model. 27Whitlow, J.L.; Varughese, J.F.; Zhou, Z.; Bartolotti, L.J.; Li, Y. Computational screening and design of S100B ligand to block S100B-p53 interaction. J. Mol. Graph. Model. 2009, 27, 969-977.

Molecular dynamics simulation of S100B protein to explore ligand blockage of the interaction with p53 protein. Z Zhou, Y Li, J. Comput. Aided Mol. Des. 23Zhou, Z.; Li, Y. Molecular dynamics simulation of S100B protein to explore ligand blockage of the interaction with p53 protein. J. Comput. Aided Mol. Des. 2009, 23, 705-714.

MDM2 interacts with MDMX through their RING finger domains. S Tanimura, S Ohtsuka, K Mitsui, K Shirouzu, A Yoshimura, M Ohtsubo, FEBS Lett. 447Tanimura, S.; Ohtsuka, S.; Mitsui, K.; Shirouzu, K.; Yoshimura, A.; Ohtsubo, M. MDM2 interacts with MDMX through their RING finger domains. FEBS Lett. 1999, 447, 5-9.

Interplay between MDM2, MDMX, Pirh2 and COP1: The negative regulators of p53. L Wang, G He, P Zhang, X Wang, M Jiang, L Yu, Mol. Biol. Rep. 38Wang, L.; He, G.; Zhang, P.; Wang, X.; Jiang, M.; Yu, L. Interplay between MDM2, MDMX, Pirh2 and COP1: The negative regulators of p53. Mol. Biol. Rep. 2011, 38, 229-236.

S100 proteins interact with the N-terminal domain of MDM2. J Van Dieck, J K Lum, D P Teufel, A R Fersht, FEBS Lett. 584Van Dieck, J.; Lum, J.K.; Teufel, D.P.; Fersht, A.R. S100 proteins interact with the N-terminal domain of MDM2. FEBS Lett. 2010, 584, 3269-3274.

Delaying vascular aging with Chinese medicine: Implications from an overview of the p53 and miR-34s family. M Wang, Y Lei, Chin. J. Integr. Med. 17Wang, M.; Lei, Y. Delaying vascular aging with Chinese medicine: Implications from an overview of the p53 and miR-34s family. Chin. J. Integr. Med. 2011, 17, 635-639.

A triterpenoid from Thalictrum fortunei induces apoptosis in BEL-7402 cells through the P53-induced apoptosis pathway. X Zhang, M Zhao, L Chen, H Jiao, H Liu, L Wang, S Ma, Molecules. 16Zhang, X.; Zhao, M.; Chen, L.; Jiao, H.; Liu, H.; Wang, L.; Ma, S. A triterpenoid from Thalictrum fortunei induces apoptosis in BEL-7402 cells through the P53-induced apoptosis pathway. Molecules 2011, 16, 9505-9519.

Synergistic antitumor effect of β-elemene and etoposide is mediated via induction of cell apoptosis and cell cycle arrest in non-small cell lung carcinoma cells. F Zhang, L Xu, X Qu, M Zhao, B Jin, J Kang, Y Liu, X Hu, Mol. Med. Rep. 4Zhang, F.; Xu, L.; Qu, X.; Zhao, M.; Jin, B.; Kang, J.; Liu, Y.; Hu, X. Synergistic antitumor effect of β-elemene and etoposide is mediated via induction of cell apoptosis and cell cycle arrest in non-small cell lung carcinoma cells. Mol. Med. Rep. 2011, 4, 1189-1193.

Ginsenoside Rb1, Rg1 and three extracts of traditional Chinese medicine attenuate ultraviolet B-induced G1 growth arrest in HaCaT cells and dermal fibroblasts involve down-regulating the expression of p16, p21 and p53. X Y Wang, Y G Wang, Y F Wang, Photodermatol. Photoimmunol. Photomed. 27Wang, X.Y.; Wang, Y.G.; Wang, Y.F. Ginsenoside Rb1, Rg1 and three extracts of traditional Chinese medicine attenuate ultraviolet B-induced G1 growth arrest in HaCaT cells and dermal fibroblasts involve down-regulating the expression of p16, p21 and p53. Photodermatol. Photoimmunol. Photomed. 2011, 27, 203-212.

Subamolide A, a component isolated from Cinnamomum subavenium, induces apoptosis mediated by mitochondria-dependent, p53 and ERK1/2 pathways in human urothelial carcinoma cell line NTUB1. C H Liu, C Y Chen, A M Huang, J H Li, J. Ethnopharmacol. 137Liu, C.H.; Chen, C.Y.; Huang, A.M.; Li, J.H. Subamolide A, a component isolated from Cinnamomum subavenium, induces apoptosis mediated by mitochondria-dependent, p53 and ERK1/2 pathways in human urothelial carcinoma cell line NTUB1. J. Ethnopharmacol. 2011, 137, 503-511.

Ginsenosides as Anticancer Agents: In vitro and in vivo Activities, Structure-Activity Relationships, and Molecular Mechanisms of Action. S A Nag, J J Qin, W Wang, M H Wang, H Wang, R Zhang, 10.3389/fphar.2012.00025Front. Pharmacol. 3Nag, S.A.; Qin, J.J.; Wang, W.; Wang, M.H.; Wang, H.; Zhang, R. Ginsenosides as Anticancer Agents: In vitro and in vivo Activities, Structure-Activity Relationships, and Molecular Mechanisms of Action. Front. Pharmacol. 2012, 3, doi:10.3389/fphar.2012.00025.

This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution license. Basel, Switzerland© 2012 by the authors; licensee MDPI© 2012 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution license (http://