Modulation of p53 N-terminal transactivation domain 2 conformation ensemble and kinetics by phosphorylation

Impact of Nonnative Interactions on the Binding Kinetics of Intrinsically Disordered p53 with MDM2: Insights from All-Atom Simulation and Markov State Model Analysis

Intrinsically disordered proteins (IDPs) lack a well-defined tertiary structure but are essential players in various biological processes. Their ability to undergo a disorder-to-order transition upon binding to their partners, known as the folding-upon-binding process, is crucial for their function. One classical example is the intrinsically disordered transactivation domain (TAD) of the tumor suppressor protein p53, which quickly forms a structured α-helix after binding to its partner MDM2, with clinical significance for cancer treatment. However, the contribution of nonnative interactions between the IDP and its partner to the rapid binding kinetics, as well as their interplay with native interactions, is not well understood at the atomic level. Here, we used molecular dynamics simulation and Markov state model (MSM) analysis to study the folding-upon-binding mechanism between p53-TAD and MDM2. Our results suggest that the system progresses from the nascent encounter complex to the well-structured encounter complex and finally reaches the native complex, following an induced-fit mechanism. We found that nonnative hydrophobic and hydrogen bond interactions, combined with native interactions, effectively stabilize the nascent and well-structured encounter complexes. Among the nonnative interactions, Leu25p53-Leu54MDM2 and Leu25p53-Phe55MDM2 are particularly noteworthy, as their interaction strength is close to the optimum. Evidently, strengthening or weakening these interactions could both adversely affect the binding kinetics. Overall, our findings suggest that nonnative interactions are evolutionarily optimized to accelerate the binding kinetics of IDPs in conjunction with native interactions.

The Role of Hydrogen Sulfide in the Localization and Expression of p53 and Cell Death in the Nervous Tissue in Traumatic Brain Injury and Axotomy

Traumatic brain injury (TBI) is one of the leading causes of disability and death worldwide. It is characterized by various molecular-cellular events, with the main ones being apoptosis and damage to axons. To date, there are no clinically effective neuroprotective drugs. In this study, we examined the role of hydrogen sulfide (H2S) in the localization and expression of the key pro-apoptotic protein p53, as well as cell death in the nervous tissue in TBI and axotomy. We used a fast donor (sodium sulphide, Na2S) H2S and a classic inhibitor (aminooxyacetic acid, AOAA) of cystathionine β-synthase (CBS), which is a key enzyme in H2S synthesis. These studies were carried out on three models of neurotrauma in vertebrates and invertebrates. As a result, it was found that Na2S exhibits a pronounced neuroprotective effect that reduces the number of TUNEL-positive neurons and glial cells in TBI and apoptotic glia in axotomy. This effect could be realized through the Na2S-dependent decrease in the level of p53 in the cells of the nervous tissue of vertebrates and invertebrates, which we observed in our study. We also observed the opposite effect when using AOAA, which indicates the important role of CBS in the regulation of p53 expression and death of neurons and glial cells in TBI and axotomy.

Phosphorylation and specific DNA improved the incorporation ability of p53 into functional condensates

The transcription factor p53 acted as a critical tumor suppressor by activating the expression of various target genes to regulate diverse cellular responses. The phosphorylation of p53 influenced the binding of p53 to promotor-specific DNA and the choice of cell fate. In this study, we found that full-length wild-type p53 and pol II CTD could form heterotypic phase separation condensates in vitro. The heterotypic condensates of p53 and pol II CTD were mediated by electrostatic and hydrophobic interactions between pol II CTD and multiple domains of p53. The mobility of heterotypic p53 and pol II CTD droplets was significantly higher than that of p53 droplet. The phosphorylation promoted p53 to be recruited into pol II CTD droplets and transcription condensates. The specific DNA could further enhance the incorporation ability of p53 into functional condensates. Therefore, we proposed that the p53 droplet might be in a mediate state, the mutations resulting in p53 mutants with gain-of-function impelled the aggregate of p53, while the phosphorylation promoted p53 to be recruited into functional condensates as a client molecule to exert its function. This study might provide insights into the regulation mechanism that the phosphorylation and nuclei acid affected the phase behavior of p53.

Dynamic structures of intrinsically disordered proteins related to the general transcription factor TFIIH, nucleosomes, and histone chaperones

Advances in structural analysis by cryogenic electron microscopy (cryo-EM) and X-ray crystallography have revealed the tertiary structures of various chromatin-related proteins, including transcription factors, RNA polymerases, nucleosomes, and histone chaperones; however, the dynamic structures of intrinsically disordered regions (IDRs) in these proteins remain elusive. Recent studies using nuclear magnetic resonance (NMR), together with molecular dynamics (MD) simulations, are beginning to reveal dynamic structures of the general transcription factor TFIIH complexed with target proteins including the general transcription factor TFIIE, the tumor suppressor p53, the cell cycle protein DP1, the DNA repair factors XPC and UVSSA, and three RNA polymerases, in addition to the dynamics of histone tails in nucleosomes and histone chaperones. In complexes of TFIIH, the PH domain of the p62 subunit binds to an acidic string formed by the IDR in TFIIE, p53, XPC, UVSSA, DP1, and the RPB6 subunit of three RNA polymerases by a common interaction mode, namely extended string-like binding of the IDR on the positively charged surface of the PH domain. In the nucleosome, the dynamic conformations of the N-tails of histones H2A and H2B are correlated, while the dynamic conformations of the N-tails of H3 and H4 form a histone tail network dependent on their modifications and linker DNA. The acidic IDRs of the histone chaperones of FACT and NAP1 play important roles in regulating the accessibility to histone proteins in the nucleosome.

An integrated view of p53 dynamics, function, and reactivation

The tumor suppressor p53 plays a vital role in responding to cell stressors such as DNA damage, hypoxia, and tumor formation by inducing cell-cycle arrest, senescence, or apoptosis. Expression level alterations and mutational frequency implicates p53 in most human cancers. In this review, we show how both computational and experimental methods have been used to provide an integrated view of p53 dynamics, function, and reactivation potential. We argue that p53 serves as an exceptional case study for developing methods in modeling intrinsically disordered proteins. We describe how these methods can be leveraged to improve p53 reactivation molecule design and other novel therapeutic modalities, such as PROteolysis TARgeting Chimeras (PROTACs).

Computational Investigation on the p53-MDM2 Interaction Using the Potential of Mean Force Study

Murine double minute 2 (MDM2) proteins are found to be overproduced by many human tumors in order to inhibit the functioning of p53 molecules, a tumor suppressor protein. Thus, reactivating p53 functioning in cancer cells by disrupting p53-MDM2 interactions may offer a significant approach in cancer treatment. However, the structural characterization of the p53-MDM2 complex at the atomistic level and the mechanism of binding/unbinding of the p53-MDM2 complex still remain unclear. Therefore, we demonstrate here the probable binding (unbinding) pathway of transactivation domain 1 of p53 during the formation (dissociation) of the p53-MDM2 complex in terms of free energy as a function of reaction coordinate from the potential of mean force (PMF) study using two different force fields: ff99SB and ff99SB-ILDN. From the PMF plot, we noticed the PMF to have a minimum value at a p53-MDM2 separation of 12 Å, with a dissociation energy of 30 kcal mol^-1. We also analyzed the conformational dynamics and stability of p53 as a function of its distance of separation from MDM2. The secondary structure content (helix and turns) in p53 was found to vary with its distance of separation from MDM2. The p53-MDM2 complex structure with lowest potential energy was isolated from the ensemble at the reaction coordinate corresponding to the minimum PMF value and subjected to molecular dynamics simulation to identify the interface surface area, interacting residues at the interface, and the stability of the complex. The simulation results highlight the importance of hydrogen bonds and the salt bridge between Lys94 of MDM2 and Glu17 of p53 in the stability of the p53-MDM2 complex. We also carried out the binding free energy calculations and the per residue energy decomposition analyses of the interface residues of the p53-MDM2 complex. We found that the binding affinity between MDM2 and p53 is indeed high [ΔG _bind = -7.29 kcal mol^-1 from molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) and ΔG _bind = -53.29 kcal mol^-1 from molecular mechanics/generalized borne surface area]. The total binding energy obtained using the MM/PBSA method was noticed to be closer to the experimental values (-6.4 to -9.0 kcal mol^-1). The p53-MDM2 complex binding profile was observed to follow the same trend even in the duplicate simulation run and also in the simulation carried out with different force fields. We found that Lys51, Leu54, Tyr100, and Tyr104 from MDM2 and the residues Phe19, Trp23, and Leu26 from p53 provide the highest energy contributions for the p53-MDM2 interaction. Our findings highlight the prominent structural and binding characteristics of the p53-MDM2 complex that may be useful in designing potential inhibitors to disrupt the p53-MDM2 interactions.

Features of molecular recognition of intrinsically disordered proteins via coupled folding and binding

The sequence-structure-function paradigm of proteins has been revolutionized by the discovery of intrinsically disordered proteins (IDPs) or intrinsically disordered regions (IDRs). In contrast to traditional ordered proteins, IDPs/IDRs are unstructured under physiological conditions. The absence of well-defined three-dimensional structures in the free state of IDPs/IDRs is fundamental to their function. Folding upon binding is an important mode of molecular recognition for IDPs/IDRs. While great efforts have been devoted to investigating the complex structures and binding kinetics and affinities, our knowledge on the binding mechanisms of IDPs/IDRs remains very limited. Here, we review recent advances on the binding mechanisms of IDPs/IDRs. The structures and kinetic parameters of IDPs/IDRs can vary greatly, and the binding mechanisms can be highly dependent on the structural properties of IDPs/IDRs. IDPs/IDRs can employ various combinations of conformational selection and induced fit in a binding process, which can be templated by the target and/or encoded by the IDP/IDR. Further studies should provide deeper insights into the molecular recognition of IDPs/IDRs and enable the rational design of IDP/IDR binding mechanisms in the future.