Structural divergence is more extensive than sequence divergence for a family of intrinsically disordered proteins

The p53 transactivation domain (p53TAD) is an intrinsically disordered protein (IDP) domain that undergoes coupled folding and binding when interacting with partner proteins like the E3 ligase, MDM2, and the 70 kDa subunit of replication protein A, RPA70. The secondary structure and dynamics of six closely related mammalian homologues of p53TAD were investigated using nuclear magnetic resonance (NMR) spectroscopy. Differences in both transient secondary structure and backbone dynamics were observed for the homologues. Many of these differences were localized to the binding sites for MDM2 and RPA70. The amount of transient helical secondary structure observed for the MDM2 binding site was lower for the dog and mouse homologues, compared with human, and the amount of transient helical secondary structure observed for the RPA70 binding site was higher for guinea pig and rabbit, compared with human. Differences in the amount of transient helical secondary structure observed for the MDM2 binding site were directly related to amino acid substitutions occurring on the solvent exposed side of the amphipathic helix that forms during the p53TAD/MDM2 interaction. Differences in the amount of transient helical secondary structure were not as easily explained for the RPA70 binding site because of its extensive sequence divergence. Clustering analysis shows that the divergence in the transient secondary structure of the p53TAD homologues exceeds the amino acid sequence divergence. In contrast, strong correlations were observed between the backbone dynamics of the homologues and the sequence identity matrix, suggesting that the dynamic behavior of IDPs is a conserved evolutionary feature. Proteins 2013; 81:1686–1698. © 2013 Wiley Periodicals, Inc.     

Elucidation of Ligand-Dependent Modulation of Disorder-Order Transitions in the Oncoprotein MDM2

Numerous biomolecular interactions involve unstructured protein regions, but how to exploit such interactions to enhance the affinity of a lead molecule in the context of rational drug design remains uncertain. Here clarification was sought for cases where interactions of different ligands with the same disordered protein region yield qualitatively different results. Specifically, conformational ensembles for the disordered lid region of the N-terminal domain of the oncoprotein MDM2 in the presence of different ligands were computed by means of a novel combination of accelerated molecular dynamics, umbrella sampling, and variational free energy profile methodologies. The resulting conformational ensembles for MDM2, free and bound to p53 TAD (17-29) peptide identify lid states compatible with previous NMR measurements. Remarkably, the MDM2 lid region is shown to adopt distinct conformational states in the presence of different small-molecule ligands. Detailed analyses of small-molecule bound ensembles reveal that the ca. 25-fold affinity improvement of the piperidinone family of inhibitors for MDM2 constructs that include the full lid correlates with interactions between ligand hydrophobic groups and the C-terminal lid region that is already partially ordered in apo MDM2. By contrast, Nutlin or benzodiazepinedione inhibitors, that bind with similar affinity to full lid and lid-truncated MDM2 constructs, interact additionally through their solubilizing groups with N-terminal lid residues that are more disordered in apo MDM2.     

Posttranslational Modifications Affect the Interaction of S100 Proteins with Tumor Suppressor p53

Proteins of the S100 family bind to the intrinsically disordered transactivation domain (TAD; residues 1-57) and C-terminus (residues 293-393) of the tumor suppressor p53. Both regions provide sites that are subject to posttranslational modifications, such as phosphorylation and acetylation, that can alter the affinity for interacting proteins such as p300 and MDM2. Here, we found that S100A1, S100A2, S100A4, S100A6, and S100B bound to two subdomains of the TAD (TAD1 and TAD2). Both subdomains were mandatory for high-affinity binding to S100 proteins. Phosphorylation of Ser and Thr residues increased the affinity for the p53 TAD. Conversely, acetylation and phosphorylation of the C-terminus of p53 decreased the affinity for S100A2 and S100B. In contrast, we found that nitrosylation of S100B caused a minor increase in binding to the p53 C-terminus, whereas binding to the TAD remained unaffected. As activation of p53 is usually accompanied by phosphorylation and acetylation at several sites, our results suggest that a shift in binding from the C-terminus in favor of the N-terminus occurs upon the modification of p53. We propose that binding to the p53 TAD might be involved in the stimulation of p53 activity by S100 proteins.     

Molecular mechanism of the interaction between MDM2 and p53

We have investigated the kinetic and thermodynamic basis of the p53-MDM2 interaction using a set of peptides based on residues 15-29 of p53. Wild-type p53 peptide bound MDM2 with a dissociation constant of 580nM. Phosphorylation of S15 and S20 did not affect binding, but T18 phosphorylation weakened binding tenfold, indicating that phosphorylation of only T18 is responsible for abrogating p53-MDM2 binding. Truncation to residues 17-26 increased affinity 13-fold, but further truncation to 19-26 abolished binding. NMR studies of the binding of the p53-derived peptides revealed global conformational changes of the overall structure of MDM2, stretching far beyond the binding cleft, indicating significant changes in the domain dynamics of MDM2 upon ligand binding.     

The conformationally flexible S9-S10 linker region in the core domain of p53 contains a novel MDM2 binding site whose mutation increases ubiquitination of p53 in vivo

Although the N-terminal BOX-I domain of the tumor suppressor protein p53 contains the primary docking site for MDM2, previous studies demonstrated that RNA stabilizes the MDM2.p53 complex using a p53 mutant lacking the BOX-I motif. In vitro assays measuring the specific activity of MDM2 in the ligand-free and RNA-bound state identified a novel MDM2 interaction site in the core domain of p53. As defined using phage-peptide display, the RNA.MDM2 isoform exhibited a notable switch in peptide binding specificity, with enhanced affinity for novel peptide sequences in either p53 or small nuclear ribonucleoprotein-U (snRNP-U) and substantially reduced affinity for the primary p53 binding site in the BOX-I domain. The consensus binding site for the RNA.MDM2 complex within p53 is SGXLLGESXF, which links the S9-S10 beta-sheets flanking the BOX-IV and BOX-V motifs in the core domain and which is a site of reversible conformational flexibility in p53. Mutation of conserved amino acids in the linker at Ser(261) and Leu(264), which bridges the S9-S10 beta-sheets, stimulated p53 activity from reporter templates and increased MDM2-dependent ubiquitination of p53. Furthermore, mutation of the conserved Phe(270) within the S10 beta-sheet resulted in a mutant p53, which binds more stably to RNA.MDM2 complexes in vitro and which is strikingly hyper-ubiquitinated in vivo. Introducing an Ala(19) mutation into the p53(F270A) protein abolished both RNA.MDM2 complex binding and hyper-ubiquitination in vivo, thus indicating that p53(F270A) protein hyper-ubiquitination depends upon MDM2 binding to its primary site in the BOX-I domain. Together, these data identify a novel MDM2 binding interface within the S9-S10 beta-sheet region of p53 that plays a regulatory role in modulating the rate of MDM2-dependent ubiquitination of p53 in cells.     

Modeling the accessible conformations of the intrinsically unstructured transactivation domain of p53

Internuclear distances derived from paramagnetic relaxation enhancement (PRE) data were used to restrain molecular dynamics simulations of the intrinsically unstructured transactivation domain of the tumor suppressor protein, p53. About 1000 structures were simulated using ensemble averaging of replicate molecules to compensate for the inherent bias in the PRE-derived distances. Gyration radii measurements on these structures show that the p53 transactivation domain (p53TAD) is statistically predominantly in a partially collapsed state that is unlike the open structure that is found for p53TAD bound to either the E3 ubiquitin ligase, MDM2, or the 70 kDa subunit of replication protein A, RPA70. Contact regions that potentially mediate the collapse were identified and found to consist of mostly hydrophobic residues. The identified contact regions preferentially place the MDM2 and RPA70 binding regions in close proximity. We show that our simulations thoroughly sample the available range of conformations and that a fraction of the molecules are in an open state that would be competent for binding either MDM2 or RPA70. We also show that the Stokes radius estimated from the average gyration radius of the ensemble is in good agreement with the value determined using size exclusion chromatography. Finally, the presence of a persistent loop localized to a PXP motif was identified. Serine residues flanking the PXP motif become phosphorylated in response to DNA damage, and we postulate that this will perturb the equilibrium population to more open conformations.     

On the origin of the stereoselective affinity of Nutlin-3 geometrical isomers for the MDM2 protein

The stereoselective affinity of small-molecule binding to proteins is typically broadly explained in terms of the thermodynamics of the final bound complex. Using Brownian dynamics simulations, we show that the preferential binding of the MDM2 protein to the geometrical isomers of Nutlin-3, an effective anticancer lead that works by inhibiting the interaction between the proteins p53 and MDM2, can be explained by kinetic arguments related to the formation of the MDM2:Nutlin-3 encounter complex. This is a diffusively bound state that forms prior to the final bound complex. We find that the MDM2 protein stereoselectivity for the Nutlin-3a enantiomer stems largely from the destabilization of the encounter complex of its mirror image enantiomer Nutlin-3b, by the K70 residue that is located away from the binding site. On the other hand, the trans-Nutlin-3a diastereoisomer exhibits a shorter residence time in the vicinity of MDM2 compared with Nutlin-3a due to destabilization of its encounter complex by the collective interaction of pairs of charged residues on either side of the binding site: Glu25 and Lys51 on one side, and Lys94 and Arg97 on the other side. This destabilization is largely due to the electrostatic potential of the trans-Nutlin-3a isomer being largely positive over extended continuous regions around its structure, which are otherwise well-identified into positive and negative regions in the case of the Nutlin-3a isomer. Such rich insight into the binding processes underlying biological selectivity complements the static view derived from the traditional thermodynamic analysis of the final bound complex. This approach, based on an explicit consideration of the dynamics of molecular association, suggests new avenues for kinetics-based anticancer drug development and discovery.     

Computational studies and peptidomimetic design for the human p53-MDM2 complex

The interaction between human p53 and MDM2 is a key event in controlling cell growth. Many studies have suggested that a p53 mimic would be sufficient to inhibit MDM2 to reduce cell growth in cancerous tissue. In order to design a potent p53 mimic, molecular dynamics (MD) simulations were used to examine the binding interface and the effect of mutating key residues in the human p53-MDM2 complex. The Generalized Born surface area (GBSA) method was used to estimate free energies of binding, and a computational alanine-scanning approach was used to calculate the relative effects in the free energy of binding for key mutations. Our calculations determine the free energy of binding for a model p53-MDM2 complex to be -7.4 kcal/mol, which is in very good agreement with the experimentally determined values (-6.6--8.8 kcal/mol). The alanine-scanning results are in good agreement with experimental data and calculations by other groups. We have used the information from our studies of human p53-MDM2 to design a beta-peptide mimic of p53. MD simulations of the mimic bound to MDM2 estimate a free energy of binding of -8.8 kcal/mol. We have also applied alanine scanning to the mimic-MDM2 complex and reveal which mutations are most likely to alter the binding affinity, possibly giving rise to escape mutants. The mimic was compared to nutlins, a new class of inhibitors that block the formation of the p53-MDM2 complex. There are interesting similarities between the nutlins and our mimic, and the differences point to ways that both inhibitors may be improved. Finally, an additional hydrophobic pocket is noted in the interior of MDM2. It may be possible to design new inhibitors to take advantage of that pocket.     

A spatiotemporal characterization of the effect of p53 phosphorylation on its interaction with MDM2

The interaction of p53 and MDM2 is modulated by the phosphorylation of p53. This mechanism is key to activating p53, yet its molecular determinants are not fully understood. To study the spatiotemporal characteristics of this molecular process we carried out Brownian dynamics simulations of the interactions of the MDM2 protein with a p53 peptide in its wild type state and when phosphorylated at Thr18 (pThr18) and Ser20 (pSer20). We found that p53 phosphorylation results in concerted changes in the topology of the interaction landscape in the diffusively bound encounter complex domain. These changes hinder phosphorylated p53 peptides from binding to MDM2 well before reaching the binding site. The underlying mechanism appears to involve shift of the peptide away from the vicinity of the MDM2 protein, peptide reorientation, and reduction in peptide residence time relative to wild-type p53 peptide. pThr18 and pSr20 p53 peptides experience reduction in residence times by factors of 13.6 and 37.5 respectively relative to the wild-type p53 peptide, indicating a greater role for Ser20 phosphorylation in abrogating p53 MDM2 interactions. These detailed insights into the effect of phosphorylation on molecular interactions are not available from conventional experimental and theoretical approaches and open up new avenues that incorporate molecular interaction dynamics, for stabilizing p53 against MDM2, which is a major focus of anticancer drug lead development.     

A spatiotemporal characterization of the effect of p53 phosphorylation on its interaction with MDM2

The interaction of p53 and MDM2 is modulated by the phosphorylation of p53. This mechanism is key to activating p53, yet its molecular determinants are not fully understood. To study the spatiotemporal characteristics of this molecular process we carried out Brownian dynamics simulations of the interactions of the MDM2 protein with a p53 peptide in its wild type state and when phosphorylated at Thr18 (pThr18) and Ser20 (pSer20). We found that p53 phosphorylation results in concerted changes in the topology of the interaction landscape in the diffusively bound encounter complex domain. These changes hinder phosphorylated p53 peptides from binding to MDM2 well before reaching the binding site. The underlying mechanism appears to involve shift of the peptide away from the vicinity of the MDM2 protein, peptide reorientation, and reduction in peptide residence time relative to wild-type p53 peptide. pThr18 and pSr20 p53 peptides experience reduction in residence times by factors of 13.6 and 37.5 respectively relative to the wild-type p53 peptide, indicating a greater role for Ser20 phosphorylation in abrogating p53 MDM2 interactions. These detailed insights into the effect of phosphorylation on molecular interactions are not available from conventional experimental and theoretical approaches and open up new avenues that incorporate molecular interaction dynamics, for stabilizing p53 against MDM2, which is a major focus of anticancer drug lead development.     

MDM2 Inhibits Axin-Induced p53 Activation Independently of its E3 Ligase Activity

MDM2 plays a crucial role in negatively regulating the functions of tumor suppressor p53. Here we show that MDM2 can inhibit Axin-stimulated p53-dependent apoptosis by suppressing p53 phosphorylation at Ser 46 and apoptosis-related p53 transactivational activity. Interestingly, the ubiquitin E3 ligase activity of MDM2 is not required for this inhibitory effect. Mechanically, either wildtype MDM2 or its E3-dead mutant, disrupts the Axin-based HIPK2/p53 complex formation by blocking the binding of p53 and HIPK2 to Axin. MDM2Δp53, a deletion mutant that lacks p53 binding domain fails to exert the inhibitory effect, demonstrating that the interaction of MDM2 and p53, but not its E3 ligase activity toward p53 plays key role in suppressing Axin-stimulated p53 activation. Our results thus have revealed a novel aspect of the mechanism by which MDM2 regulates p53 activities.     

Why is F19Ap53 unable to bind MDM2? Simulations suggest crack propagation modulates binding

Why doesn’t the F19A mutant of p53 bind to MDM2? Binding thermodynamics have suggested that the loss of packing interactions upon mutating Phe into Ala sidechain results in destabilizing the binding free energy between p53 and MDM2. Does this mutation also modulate the initial recognition between p53 and MDM2? We look at atomistic computer simulations of the process of the initial encounter between wild type p53 peptide and its F19A mutant with the N-terminal domain of MDM2. These simulations show that binding is characterized by a complex multistep process. It starts with the capture of F19 of wild type p53 by certain residues in the MDM2 binding pocket. This initial step anchors the peptide onto the surface of MDM2, and with the consequent reduction in the search space of the peptide, the peptide docks into the partially occluded surface of MDM2. This is similar to a crack forming in an otherwise occluded hydrophobic cavity in MDM2, and the peptide, docked through F19, modulates the propagation of this crack, which subsequently results in the stepwise docking of the rest of the peptide through insertions of W23 and L26. The lack of the bulky sidechain of F in the F19A mutant results in the absence of the initial “grasp” complex, and hence the mutant peptide diffuses randomly on the surface of MDM2 without binding. This is the first such demonstration of the possibility that a “kinetic” effect may partly underlie the destabilized thermodynamics of binding of F19A and is a feature that appears to be conserved in evolution. The observations by Wallace et al. (Mol Cell 2006; 23:251–63) that despite the inability of F19A to bind at the N-terminal domain of MDM2, it gets ubiquitinated, can now be partly understood based on a mechanism whereby the occupation of the binding pocket by ligands/peptides induces, via crack propagation and the dynamics of gatekeeper Y100, the ubiquitination signal for interactions between the acidic domain of MDM2 and the DNA binding domain of p53.     

NMR chemical shift and relaxation measurements provide evidence for the coupled folding and binding of the p53 transactivation domain

The interaction between the acidic transactivation domain of the human tumor suppressor protein p53 (p53TAD) and the 70 kDa subunit of human replication protein A (hRPA70) was investigated using heteronuclear magnetic resonance spectroscopy. A ¹H–¹⁵N heteronuclear single quantum coherence (HSQC) titration experiment was performed on a ¹⁵N-labeled fragment of hRPA70, containing the N-terminal 168 residues (hRPA70_1–168) and p53TAD. HRPA70_1–168 residues important for binding were identified and found to be localized to a prominent basic cleft. This binding site overlapped with a previously identified single-stranded DNA-binding site, suggesting that a competitive binding mechanism may regulate the formation of p53TAD–hRPA70 complex. The amide ¹H and ¹⁵N chemical shifts of an uniformly ¹⁵N-labeled sample of p53TAD were also monitored before and after the addition of unlabeled hRPA70_1–168. In the presence of unlabeled hRPA70_1–168, resonance lineshapes increased and corresponding intensity reductions were observed for specific p53TAD residues. The largest intensity reductions were observed for p53TAD residues 42–56. Minimal binding was observed between p53TAD and a mutant form of hRPA70_1–168, where the basic cleft residue R41 was changed to a glutamic acid (R41E), demonstrating that ionic interactions play an important role in specifying the binding interface. The region of p53TAD most affected by binding hRPA70_1–168 was found to have some residual alpha helical and beta strand structure; however, this structure was not stabilized by binding hRPA70_1–168. ¹⁵N relaxation experiments were performed to monitor changes in backbone dynamics of p53TAD when bound to hRPA70_1–168. Large changes in both the transverse (R₂) and rotating frame (R_1ρ) relaxation rates were observed for a subset of the p53TAD residues that had ¹H–¹⁵N HSQC resonance intensity reductions during the complex formation. The folding of p53TAD upon complex formation is suggested by the pattern of changes observed for both R₂ and R_1ρ. A model that couples the formation of a weak encounter complex between p53TAD and hRPA70_1–168 to the folding of p53TAD is discussed in the context of a functional role for the p53–hRPA70 complex in DNA repair.     

Structure of free MDM2 N-terminal domain reveals conformational adjustments that accompany p53-binding

Critical to the inhibitory action of the oncogene product, MDM2, on the tumour suppressor, p53, is association of the N-terminal domain of MDM2 (MDM2N) with the transactivation domain of p53. The structure of MDM2N was previously solved with a p53-derived peptide, or small-molecule ligands, occupying its binding cleft, but no structure of the non-liganded MDM2N (i.e. the apo-form) has been reported. Here, we describe the solution structure and dynamics of apo-MDM2N and thus reveal the nature of the conformational changes in MDM2N that accompany binding of p53. The new structure suggests that p53 effects displacement of an N-terminal segment of apo-MDM2N that occludes access to the shallow end of the p53-binding cleft. MDM2N must also undergo an expansion upon binding, achieved through a rearrangement of its two pseudosymetrically related sub-domains resulting in outward displacements of the secondary structural elements that comprise the walls and floor of the p53-binding cleft. MDM2N becomes more rigid and stable upon binding p53. Conformational plasticity of the binding cleft of apo-MDM2N could allow the parent protein to bind specifically to several different partners, although, to date, all the known liganded structures of MDM2N are highly similar to one another. The results indicate that the more open conformation of the binding cleft of MDM2N observed in structures of complexes with small molecules and peptides is a more suitable one for ligand discovery and optimisation.     

Structure of the p53 transactivation domain in complex with the nuclear receptor coactivator binding domain of CREB binding protein

The activity and stability of the tumor suppressor p53 are regulated by interactions with key cellular proteins such as MDM2 and CBP/p300. The transactivation domain (TAD) of p53 contains two subdomains (AD1 and AD2) and interacts directly with the N-terminal domain of MDM2 and with several domains of CBP/p300. Here we report the NMR structure of the full-length p53 TAD in complex with the nuclear coactivator binding domain (NCBD) of CBP. Both the p53 TAD and NCBD are intrinsically disordered and fold synergistically upon binding, as evidenced by the observed increase in helicity and increased level of dispersion of the amide proton resonances. The p53 TAD folds to form a pair of helices (denoted Pα1 and Pα2), which extend from Phe19 to Leu25 and from Pro47 to Trp53, respectively. In the complex, the NCBD forms a bundle of three helices (Cα1, residues 2066-2075; Cα2, residues 2081-2092; and Cα3, residues 2095-2105) with a hydrophobic groove into which p53 helices Pα1 and Pα2 dock. The polypeptide chain between the p53 helices remains flexible and makes no detectable intermolecular contacts with the NCBD. Complex formation is driven largely by hydrophobic contacts that form a stable intermolecular hydrophobic core. A salt bridge between D49 of p53 and R2105 of NCBD may contribute to the binding specificity. The structure provides the first insights into simultaneous binding of the AD1 and AD2 motifs to a target protein.