- Open Access
Shaping the regulation of the p53 mRNA tumour suppressor: the co-evolution of genetic signatures
BMC Cancer volume 19, Article number: 915 (2019)
Structured RNA regulatory motifs exist from the prebiotic stages of the RNA world to the more complex eukaryotic systems. In cases where a functional RNA structure is within the coding sequence a selective pressure drives a parallel co-evolution of the RNA structure and the encoded peptide domain. The p53-MDM2 axis, describing the interactions between the p53 tumor suppressor and the MDM2 E3 ubiquitin ligase, serves as particularly useful model revealing how secondary RNA structures have co-evolved along with corresponding interacting protein motifs, thus having an impact on protein – RNA and protein – protein interactions; and how such structures developed signal-dependent regulation in mammalian systems. The p53(BOX-I) RNA sequence binds the C-terminus of MDM2 and controls p53 synthesis while the encoded peptide domain binds MDM2 and controls p53 degradation. The BOX-I peptide domain is also located within p53 transcription activation domain. The folding of the p53 mRNA structure has evolved from temperature-regulated in pre-vertebrates to an ATM kinase signal-dependent pathway in mammalian cells. The protein – protein interaction evolved in vertebrates and became regulated by the same signaling pathway. At the same time the protein - RNA and protein - protein interactions evolved, the p53 trans-activation domain progressed to become integrated into a range of cellular pathways. We discuss how a single synonymous mutation in the BOX-1, the p53(L22 L), observed in a chronic lymphocyte leukaemia patient, prevents the activation of p53 following DNA damage. The concepts analysed and discussed in this review may serve as a conceptual mechanistic paradigm of the co-evolution and function of molecules having roles in cellular regulation, or the aetiology of genetic diseases and how synonymous mutations can affect the encoded protein.
The p53 tumor suppressor and its main regulator, the MDM2 E3 ubiquitin ligase, constitute a fine model to understand molecular co-evolution, conservation and adaptability, as well as the molecular basis of cancer or other genetic diseases [1,2,3,4,5]. It was recently shown that the BOX-I motif of the transactivation domain of p53 has co-evolved with its regulator MDM2, both at the RNA and the protein levels, leading to the evolution of an intimate regulatory script adopting distinct roles in various stress-induced signalling pathways . Under normal conditions, MDM2 promotes the degradation of p53 via protein-protein interactions while following genotoxic stress, the stress sensor ATM kinase is activated by double-stranded DNA breaks and phosphorylates MDM2(S395) inducing a conformational change which dramatically increases the affinity of MDM2 for the p53 mRNA [6,7,8,9,10,11,12,13,14]. The stress-induced MDM2-p53, protein-RNA interaction leads to the stabilisation of p53 via a mechanism whereby MDM2 becomes a positive regulator of p53 [10, 15, 16]. The p53-MDM2 axis contributes a few very important implications and may serve as a paradigm, both mechanistically and conceptually, to understand mechanisms of cellular signalling, the role of intrinsically disordered domains, the role of molecular signatures and interacting motifs as well as their co-evolution, deriving from selective pressure.
This review further highlights novel discoveries on functional interactions between molecular partners, both at the protein - protein and the protein - RNA interaction levels and how p53 evolved from an ancient p53/p63/p73 protein having roles in development, to become a tumor suppressor with numerous interacting partners and functions [17,18,19,20]. Findings from in vitro studies on co-evolutionary structural modifications on the interacting motifs and the stereochemically flanking domains on p53 and MDM2 regulating the expression and stabilisation both at the RNA and protein levels from pre-vertebrates, are presented and discussed. These results are set into context with previous evidences supporting a model whereby RNA structures interacting with peptidic motifs may have co-evolved from early prebiotic environments of the RNA world hypothesis to adopt an intimate biochemical relationship with various molecular and cellular functions. The concepts discussed here thus give insights on the nature of the guiding force of the evolution and on a strategy to identify molecular profiling signatures within key players regulating the cellular processes or the development of genetic diseases.
Transition from an RNA world
Life-forms require at least three biopolymers (DNA, RNA and proteins) that mediate the biochemical processes of DNA replication, transcription and RNA translation. This well-orchestrated complex machinery strongly implies its evolution from a more simple system . Recently, a chemical reaction network building up 9 of the 11 intermediates of the biological Krebs (or tricarboxylic acid) cycle, was observed. Such chemical reactions may represent prebiotic precursors to core metabolic pathways . The RNA World hypothesis, describing an intermediate stage of life  is a proposed model of ancient biochemistry where structured RNA acquires catalytic properties [24, 25]. Close geological settings and environmental conditions undergoing specific changes (e.g. in the pH) and mixtures of simple chemical compounds could form the required precursors for the prebiotic RNA synthesis. Indeed, a mixture of hydrogen cyanide (HCN) and hydrogen sulphide (H2S) activated by ultraviolet light was shown to effectively form the required precursors of nucleotides, amino acids and lipids [26, 27]. The first polynucleotides are suggested to be small oligomers formed randomly or by non-enzymatic template-copying, via such conditions that promote a feedback between molecular activity and fitness, whereby certain sequences gained a competitive advantage. Chemical properties such as the charge and the hydrophobicity [21, 28] and an early achieved biopolymeric chirality are suggested to have adequately promoted the selection of certain RNA sequences from a vast heterogenous pool of chemical precursors which may catalyse the formation of amide bonds [29, 30] and enforce an enantiomeric selection of peptidic products .
As such, early metabolic processes which evolved with a selective preference for certain amino acids, associated with cognate adapter RNAs (tRNA) and their pairing to a linear genetic molecule (mRNA) leading to ordered peptide synthesis (ribosome), is possible . A selection of higher activity sequences can be enforced within lipid vesicle compartments, resembling to encapsulated cell-like systems forming nucleocapsids and protein assemblies that acquire virus-like genome packaging and protection [33, 34]. Such assemblies are currently developed for biomedical applications generating artificial protocells and translational regulatory systems . The transition from an RNA world to self-replicating systems employing DNA, could explain how DNA genomes have evolved via molecular interactions leading to increased biological complexity [36, 37] and to the formation of an early RNA enzyme (ribozyme) capable of copying RNA molecules [23, 38, 39] or even assembled DNA genomes . Additionally, it was recently shown that 3D structural motifs promoting viroid RNA folding, constitute a critical constraint in the RNA–RNA and RNA–protein interactions regulating the viroid genome evolution. Accordingly, only mutations which did not disrupt the structure and function were retained in the population . Such examples illustrate how specific secondary RNA structures and motifs may resemble to nucleation points acting as molecular scaffolds on which interfaces co-evolve along with cognate sites on partnering molecules, to expand the functional cues of their interplay.
In conclusion, the transition form a system as described in the RNA world hypothesis, consisting of small peptides and small RNAs encapsulated in enclosed settings, could be mediated via a co-evolutionary simultaneous lengthening of precursor molecules which selectively reinforced the formation and survival of structured co-evolving motifs .
Shaping the evolutionary path: the interplay between RNA and protein
It has become evident that secondary RNA structures and specific motifs play a crucial role in the interaction with proteins  and the underlying mechanisms of protein–RNA interaction have started to unravel [44, 45]. Mutations taking place on structurally sensitive sites on the mRNA that have a disruptive effect on the secondary structure cause aberrant gene expression and their adaptation depends on the expression levels . A recent study, has presented a comprehensive resource database for synonymous mutations from human cancers, emphasising how such single point mutations may have an impact on the expression as well as on the mRNA secondary structure, splicing, RNA stability, RNA folding, and co-translational protein folding . It becomes clear that synonymous mutations affecting certain mRNA secondary structures, may have a strong impact on vital cellular processes. Indeed, the structure of an RNA may favor the accumulation of genetic variation in proteins  or regulate the efficiency of translation via several mechanisms. Examples include the 5′ UTR–mediated initiation and stabilization , the group of riboswitches  and the double-stranded RNA-activated protein kinase . RNA structural features also modulate the dynamics of protein folding during protein synthesis . Specific RNA motifs have been proposed to drive the high mutation rate and the genetic viral variation leading to the formation of quasispecies  in the human immunodeficiency virus 1 (HIV-1) [52, 54] as a result of a reciprocal interference between selection at the RNA and protein levels. A fine example of RNA adaptation in a co-evolutionary setting is how the HIV-1 uses the APOBEC3-Vif interaction (host) to modulate its own mutation rate in harsh or variable environments . In eukaryotes, the biological specificity of RNA binding proteins (RBPs) for RNA is affected by the RNA structure and its concentration [56, 57]. Indeed, post-transcriptional gene regulation and RNA processing mediated by proteins adapted to splice introns (RNA splicing) results from their co-evolution with partner RNAs . RBPs play critical roles in post-transcriptional regulation of gene expression in all domains of life , having roles in development, cell cycle and signalling mechanisms.
Proteins with RNA-related functions constitute the 7.5% of the human proteome (a set of 1542 proteins) [60, 61] and the occurrence of intrinsically disordered regions (IDRs) within RBPs appears to be conserved and expanded from yeast to humans, often in the form of repeats that co-evolved with the increasing complexity of eukaryotic transcriptomes . IDRs are abundant in RBPs and approximately 20% of mammalian RBPs contain disordered regions which correspond to over 80% of their complete structures/sequences . IDRs are subjected to strong sequence constraints , implying crucial functional roles in the intermolecular and intramolecular interactions. A fine example is the RBP-dependent regulation of the p53 tumor suppressor, whereby RBPs, such as RBM38, interact with the 5′ or 3′ UTR of p53 mRNA thus playing pivotal roles in both maintaining genomic integrity and tumor suppression . The p53 itself comprises IDRs dramatically increasing the potential of the p53 interactome  and an extensive description of p53 mRNA interacting factors in context of various cellular conditions and stress responses, was presented recently . A computational model describing how proteins associate with extended RNA elements may aid the discrimination between regulated and non-regulated RNAs and reveal potential regulatory motifs, by giving insights on the consequences of mutational events towards the binding activity . In another computational approach, about 80 cellular processes that can be regulated at the post-transcriptional level were identified by analysis of protein–mRNA interaction networks from more than 800 human RBPs; and mechanistic regulatory hypotheses were proposed . These studies show that there is growing interest in identifying cellular roles of RNA - protein interfaces in the context of human diseases and it would be of great interest and usefulness to include the well-studied p53 mRNA - MDM2 axis.
Other factors influencing the RNA-protein interplay thus shaping their co-evolutionary path by balancing the intrinsic conflict between stability and aggregation, include the chaperone Hsp90 . Hsp90 modulates trade-offs among protein stability and increased aggregation and hydrophobicity, acting both at the polypeptidic level and the RNA level. Acting as a chaperone, it reduces the aggregation of intermediate folding states of proteins while regulating the translation speed, thus leading to more stable and versatile proteins. In the eukaryotic cell, Hsp90 promotes evolutionary changes in otherwise entrenched developmental processes, by functioning as a buffer allowing the accumulation of silent sequence variants under neutral conditions, which become selected upon conditions while Hsp90 activity is compromised. In this way, certain variants gain the prospect to get enriched by selection and even rapidly become independent of the Hsp90 mutation .
Additionally, it has been shown that post-translational modifications (PTMs) are known to regulate RNA-binding proteins  as well as proteins that are regulated via complex mechanisms employing RNA. It has indeed become evident that the interaction of secondary RNA structures with respective RBDs in proteins is favoured in sequences including IDRs and in substrates where modifications take place. These factors, along with the activity of chaperones and of post-transcriptional and post-translational mechanisms, regulate the morphological co-evolution of such RNA-protein partners. A particularly interesting and useful example is the interaction of the p53 tumor suppressor with MDM2, which is employed in this review as a model to describe the co-evolution of protein-mRNA structures with functional outcomes; and it is discussed in detail in the following sections.
Τhe evolution of the p53 family from an ancestral molecule
In mammals, the p53 family consists of the p53, p63, and p73 proteins. P63 and p73 have important roles in development [72, 73]. While the p63 protein has a role in skin and epithelial development, p73 has an important role in neuronal development and differentiation . Activated by cellular stresses such as genotoxic stress [8, 74] and endoplasmic reticulum stress , p53 has evolved to play a key role in cellular homeostasis as a tumor suppressor and as a signal response factor preventing oncogenesis (Fig. 1, a).
The proteins of the p53 family have evolved from whole genome duplications of a single ancestral p53/p63/p73 protein in the vertebrates while invertebrate p53 superfamily members appear to have a p63-like domain structure [77, 78]. Evolutionary studies have identified various homologues of the ancestral protein in the placozoan Trichoplax adhaerens , in Nematoda roundworms (Caenorhabditis) , in the pre-vertebrate Ciona intestinalis , in the invertebrate Mytilus trossulus , in the early vertebrate cartilaginous fish, as well as in Cnidaria like the starlet sea anemone (Nematostella vectensis), and in flies (Drosophila) . The phylogenetic relationships of metazoan genes containing p53/p63/p73 transactivation domains (TADs) and corresponding p53/p63/p73BD on MDM2, in context of their co-evolution, have been thoroughly investigated and structural models comparing the FxxxWxxL motifs (box-I) were presented . A detailed phylogenetic analysis of the evolution of the p53 family with a particular focus on the TAD and its co-evolution with the binding domain on MDM2 showed that these sequences are significantly conserved from early metazoan time . Such predictive and experimental studies focusing on the interfaces of structured RNA and protein binding epitopes, are highly encouraged. The function of this ancestral gene in the early metazoan sea anemone was suggested to be related with the protection of the germ-line gametes from DNA damage while the first functional change coincided with the development of stem cells and progenitor cells, designed to regenerate somatic tissues over the life time of adult organisms. As with germ cells, this has allowed the possibility of unlimited cell growth and the development of cancer . However, the function of the ancestral p53/p63/p73 and how it later evolved into three functionally diverse proteins that share a certain homology, has only recently started to unravel. Phylogenetic analyses have shown that the p53/p63/p73 family has been duplicated multiple times during the evolution, from the invertebrates . In the vertebrate lineage leading to fishes, reptiles and mammals, duplications gave rise to three distinct genes (p53, p63 and p73) that are retained in the majority of species . Such p53 gene duplications are suggested to contribute to an enhanced induction of apoptosis via a hyperactive TP53 signaling pathway activated in response to DNA damage and they have been linked with the evolution of large body sizes and with the resolution of Peto′s paradox that correlates the body size with the risk of cancer [2, 82]. A phylogenetic analysis of the evolution of the ancestral p53/p63/p73 protein from the pre-vertebrate Ciona intestinalis and the p53, p63 and p73 proteins from human, shows that the full-length sequence of the ancestral protein is more closely related to human p53 and p73, as compared to p63 (Fig. 1, b I). Similarly, the aligned sequence of the N-terminal within the ancestral transactivation domain (TAD) (aa 1–33) shows a higher similarity with the corresponding sequences on p53 (aa 1–36) and p73 (aa 1–32) as compared to p63 (aa 1–30) (Fig. 1, b II). During evolution, a distinguished accelerated expansion of various gene families coding for transcription factors has been observed and it is often predated by modular domain rearrangements, forming new sub-families in terms of protein–protein interactions. This separation allows for radical shifts in the functional spectrum of duplicated transcription factors . Such domain rearrangements can take place on modular regions involved in molecular binding and the modular organization of binding sites is shown to confer robustness and functional flexibility, facilitating the evolution of protein interactions . Gene evolution at the domain level has been an interesting point of study . An example of modular evolution is the BOX-I motif of p53, analysed here.
A key evolutionary signature: the BOX-I motif of the transactivation domain of p53
All members of the mammalian p53 family are transcription factors containing an N-terminal transactivation domain (TAD), a DNA-binding domain and a C-terminal oligomerization domain. Even though they exhibit a high overall sequence and structural similarity, their N-terminal TADs are not well conserved [4, 86, 87]. The TAD sequence of the mammalian p53 is well-conserved and forms an amphipathic a-helix domain (aa. 19–26, BOX-I) [7, 88, 89] which exhibits steric complementarity with the N-terminal hydrophobic cleft (aa 26–108) of the E3 ubiquitin ligase MDM2 oncoprotein. The TAD domain of p53 confers the p53 transcriptional transactivation in mammals and is thus under selective pressure in the evolution. The p53 transactivation activity is mediated via two adjacent but functionally specialized domains (TAD1 and TAD2), transactivating different target genes and effector pathways [90, 91]. Yet, TAD1 plays a predominant role over TAD2, and is required for DNA damage-induced G1 arrest and apoptosis, but not for RAS-induced senescence in fibroblasts, explaining the evolutionary pressure to split the transactivation function into two domains with different panels of co-factors and modifying complexes that confer a more robust and dynamic (context-dependent) transactivation activity . The BOX-I motif of the TAD1 of p53 (aa. 19–26) is the most conserved region of p53 and was shown to play a dual role on its interaction and regulation of p53 by MDM2; one at the RNA level and one at the protein level [14, 93]. These interactions evolved from the pre-vertebrate ancestral p53(ΒΟΧ-Ι) as a result of co-evolution with the MDM2 protein . In other studies structural models comparing FxxxWxxL peptidic motifs (BOX-I) have also shown that the p53(BOX-I) - MDM2 interplay dates back to early metazoan time . However, to our knowledge, there are no additional studies focusing on the molecular evolution of the p53(box-I) mRNA structure and on the interactions with MDM2 homologues from various species. Such experimental and predictive studies with an emphasis on the functional roles of such molecular interfaces on co-evolved molecular partners, are highly encouraged and anticipated.
The co-evolution of p53 and MDM2 proteins
The TAD domain of p53 comprises the co-evolved epitopes mediating the p53 - MDM2 interaction. P53 is tightly regulated by MDM2. Under normal conditions, MDM2 negatively regulates p53 by promoting p53 ubiquitination and degradation via the 26S proteasomal pathway and by blocking its TA activity. However, during genotoxic stress, MDM2 switches to become a positive regulator of p53. The underlying mechanism also involves MDMX (HDMX), a close homolog to MDM2 which has no E3 ubiquitin ligase activity [93, 94]. Phosphorylation at serine MDM2(S395) and MDMX(S403) by the ATM kinase induces conformational changes both on MDM2 and MDMX [94, 95], allowing their RING domains to bind the p53 mRNA sequence that encodes the BOX-I motif , resulting to an increase in p53’s rate of synthesis and to the suppression of the MDM2’s E3 ubiquitin ligase activity [10, 89, 93, 96]. Hence, two MDM2-interacting motifs with opposing functions towards p53 expression have evolved from the same genomic sequence of p53: one at the mRNA level and one peptidic [4, 15] and these peptide- and RNA- motifs interacting with MDM2 are encoded by the same conserved BOX-I sequence. In an attempt to describe how this molecular co-evolution took place, it has been shown that the evolution of the p53 peptide- and RNA- interactions with MDM2 have been influenced and selected by different cues from pre-vertebrates to vertebrates . Short engineered peptides of the C. intestinalis or mammalian BOX-I sequences bind similarly well to both either C. intestinalis or mammalian MDM2 proteins. However, when the full length p53 proteins were tested, Ci-p53wt showed an interaction with either Hu-MDM2 and Ci-MDM2 protein, but the Hu-p53wt showed poor affinity for Ci-MDM2 in ELISA experiments . These results strongly indicated that the allosteric interference imposed by the C-terminal flanking region of the BOX-I domain prevents the p53-MDM2 protein–protein interaction. Indeed, it has been experimentally shown in the pre-vertebrate C. intestinalis that the region encompassing the residues Q41 to F56, prevents the interaction and deletion of these flanking residues (Q41 to F56) not only restored the interaction between C. intestinalis p53 protein and MDM2 but it also induced the binding to the DO-I mAB which binds the BOX-I domain . This was confirmed in Proximity Ligation Assay (PLA) experiments . The PLA microscopy is a state-of-the-art microscopy-based technique that provides an in situ semi-quantitative estimation of endogenous molecular interactions/associations that occur in a low frequency and/or mediated by molecules of a low concentration; and are thus considered to be of a low abundance, in comparison to numerous other interactions taking place in the cell [97,98,99]. It offers unique information on the sub-cellular localisation of molecular interactions between proteins or between RNA and protein. These findings strongly indicated that p53 has evolved from an ancestral p53/p63/p73 gene to interact with MDM2 by elimination of the encoding flanking region of BOX-I which paved the way for MDM2 to take a negative regulatory role on p53. These findings highlighted that computational studies and molecular evolutionary models based on sequence alignment, should also take into account 3D structural evidences or models. In line with this notion, an evolutionary study on a molecular level of the p53/p63/p73 TAD domain with regard to protein disorder and regulatory properties, showed similarities in the phosphorylation pattern of vertebrate p53 and mollusk and annelid p53/p63/p73, implying that functional properties of regulation via phosphorylation were already present in p53/p63/p73 from deuterostomes (e.g. Chordata) and protostomes (e.g. Mollusca and Arthropoda) . Additionally, a phylogenetic analysis on the p53/p63/p73 and MDM proteins from phyla that retain the interaction domains TAD and p53/p63/p73 BD (binding domain), based on both vertebrate and invertebrate species showed that the signaling pathway of the TAD and p53/p63/p73BD has co-evolved or disappeared in distinct lineages . In line with this, evolutionary studies have suggested that the MDM2 – p53 interaction is present in early metazoans and was later lost in some species, but it is not clear whether this event occurred in the pre-vertebrates or earlier in the evolution [1, 5].
In conclusion, alternations in the intrinsically disordered p53 N-terminal sequence of pre-vertebrates shaped the conformation and the presentation of the conserved BOX-I peptide motif (TAD) to interact with MDM2 in vertebrates, allowing MDM2 to become a negative regulator of p53.
The co-evolution of the p53 mRNA - MDM2 protein interaction
In mammals the p53 mRNA - MDM2 protein interaction is facilitated via MDMX [12, 100] which is found in mammals but not in the pre-vertebrate C.intestinalis. Following phosphorylation by the ATM kinase during the DNA damage response, MDMX binds the nascent p53 mRNA forming an RNA platform/structure on which MDM2 can bind and stimulate p53 synthesis. This mRNA structure is found on the 5′ of the p53 mRNA coding sequence (within the first 120 nt) and consists of three stem loops [12, 101]. In the pre-vertebrates the catalytic role of MDMX for the p53 mRNA - MDM2 interaction is influenced by temperature. Temperature was shown to govern the folding of the ancestral Ci-p53 mRNA and its binding to Ci-MDM2, by inducing a stem-loop structure within the box-I RNA motif. This structure interestingly resembles to the mammalian p53 homolog that has a high affinity for MDMX and MDM2, following DNA damage [4, 12, 15], (Fig. 2, a), indicating a putative role of the ancestral MDM2 as a positive regulator of p53. The RNA structures were solved experimentally by the in vivo DMS Footprint Assay and the in vitro DMS modification assay  that use dimethyl sulphate (DMS)-modified RNA and sequencing to reveal the RNA structure. The box-I-coding Ci-p53 mRNA sequence adopts a temperature-dependent structure that governs the interaction with Ci-MDM2. This interaction was experimentally confirmed by RNA co-immunoprecipitation coupled with qPCR and ELISA and it was visualized on fixed embryos using a modified version of the PLA, the RNA PLA [4, 10, 102]. Similar findings of temperature-dependent structural adaptations of transcripts were identified in Protozoa, Nematoda, Cnidaria, and Tunicata  while temperature was shown to influence mRNA structures and translation in various organisms [104,105,106], including bacteria , yeast ; corals , and plants . It has been particularly interesting to track in the evolution how the temperature-dependent structured pre-vertebrate p53(box-I) mRNA becomes folded in an MDMX–dependent fashion in vertebrates and mammals while at the same time the encoded BOX-I motif and its flanking region within the encoded TAD, evolved to facilitate the protein - protein interaction  (Fig. 2, b).
In conclusion, the steric complementarity of the MDM2 - p53 mRNA interaction requires a specific p53(box-I) mRNA structure, which has evolved within the p53 mRNA sequence. This temperature-induced RNA structure in the pre-vertebrates evolved to become MDM2 and MDMX-dependent, indicating a putative role of the ancestral MDM2 as a positive regulator of p53.
The regulation of p53 by MDM2 at the protein and RNA levels
The p53 tumor suppressor and the p53 regulatory axis is a fine model constituting a complex ATM kinase - dependent regulation system [10, 14, 111, 112] that involves: (1) the employment of post-transcriptional regulation by RBPs ; (2) the employment of IDRs conferring multi-functionality ; (3) the employment of post-translational modifications following different stress responses and (4) the formation of a multifactorial mechanism, employing various additional proteins, as it has evolved from a more simple system [4, 14]. Recent findings showing the functional consequences of a single synonymous cancer mutation on the p53 mRNA (L22 L) that abrogates the interaction of the p53 mRNA with the MDM2 E3 ubiquitin ligase protein , highlight the notion that certain structured RNA motifs constitute a signature of regulation and urge the attention for more detailed studies on similar molecular partners and networks. One of the key points of this review is to analyse and discuss the role of the p53 mRNA in the regulation of p53 and to give insights on the notion that certain structured motifs, such as the box-I motif of p53, may be encoded by sequences which are prone to adaptive mutation (sequence-dependent ‘hotspot’ ) thus constituting co-evolved structures allowing molecular interactions between RNA and protein.
The conservation of the mammalian MDM2 - p53(box-I), protein-protein and protein-RNA interactions has become evident and their functional consequences have started to unravel. Under normal conditions MDM2 acts as a negative regulator of p53. MDM2 binds the p53(BOX-I) protein and targets it for degradation via the 26S proteasomal pathway. However, following DNA damage, the ATM kinase is autophosphorylated and phosphorylates p53(S15) preventing the p53-MDM2 protein-protein interaction. Additionally, ATM phosphorylates MDM2(S395) which in turn develops an affinity for the p53 mRNA. This event also leads to the prevention of the binding of MDM2 on the p53 protein. These phosphorylation events lead to the stabilisation of p53 [3, 114,115,116,117]. Thus, there are at least two phosphorylation events leading to an MDM2-dependent regulation of p53: (a) phosphorylation of p53(S15) abrogating the binding of MDM2 and preventing the ubiquitination and the proteasomal targeting of p53; and (b) phosphorylation of MDM2(S395) promoting the interaction of MDM2 with the p53 mRNA thus leading to enhanced p53 synthesis. A recent study experimentally showed how these events are orchestrated via a mechanism whereby phosphorylated MDM2(S395) acts as a p53 mRNA-dependent carrier that mediates the formation of a complex at the cytoplasm. This complex, consists of the MDM2, the p53 mRNA, the ribosomal proteins RPL5 and RPL11 and the ATM kinase and was shown to promote the synthesis and activation of p53 protein while preventing its degradation . To analyse these interactions, the p53(L22 L) synonymous mutation was used. L22 L at the codon in position 22, CUA to CUG, is located in the apical loop of the hairpin U180-A218, and was initially observed in a chronic lymphocyte leukaemia patient [101, 118]. The p53(L22 L) synonymous mutation is located at the TAD domain of p53 which constitutes the MDM2 binding site and it prevents the interaction of the p53 mRNA with MDM2 by disrupting the p53 hairpin-MDM2 interaction . It was also shown to impair p53 and Δ40p53 synthesis [118, 119] and to lead to a poor stabilization of the encoded p53 protein following genotoxic stress [14, 15]. P53(L22 L) showed poor affinity for MDM2 and this resulted to a poor MDM2-dependent positive regulation effect towards p53, following genotoxic stress. Additionally, the MDM2(C305F) and MDM2(C308Y) mutants that prevent the interactions with the ribosomal proteins RPL5 and RPL11 respectively, failed to stimulate the over-expression of p53 following DNA damage. These mutants remained localised in the nucleolus and did not reach the cytoplasm to form the ribosome including the RPL5, RPL11 and the p53 mRNA (as shown by immunofluorescence). The work of Karakostis et al,  presents direct evidences showing that the interactions of MDM2 with (a) the p53 mRNA and with (b) the ribosomal proteins RPL5 and RPL11 are required for the translocation of the complex to the cytoplasm, supporting the concept that MDM2 acts as a carrier of RP proteins and of the p53 mRNA at the cytoplasm which promotes the formation of the ribosome that translates the p53 mRNA (Fig. 3).
Consequently and in parallel with the release of MDM2 from the translating p53 mRNA, the under-synthesis nascent p53 peptide is phosphorylated on S15 by the ATM kinase which has hitch-hiked at the p53-translating ribosome by MDM2 and it is thus localised in close proximity with the nascent p53 peptide. This phosphorylation prevents the interaction of the newly synthesized p53 peptide with free (not binding the p53 mRNA) MDM2, leading to its stabilization. Evidently, as p53 synthesis continues, MDM2 will be displaced from the translating p53 mRNA. In human, the binding platform consists of three stem loops within the first 120 nt of the encoding sequence of the p53 mRNA. Due stereochemical constrains, ATM may bind the nascent p53 protein and phosphorylate S15 only at the point when about 45 aa have been synthesised. This is qualitatively consistent with MDM2 becoming displaced from the p53 mRNA. In parallel, MDM2 interacts with RPL5 and RPL11 as well with other ribosomal proteins, such as the RPL37, RPS15 and RPS20 independently of p53, as it has been shown in H1299 cells (p53 null) [120, 121]. Conclusively, MDM2 remains bound on the translating ribosome until p53 is synthesised and phosphorylated by ATM at S15. This was further confirmed by PLA and PLEA on purified polysomes . The Proximity Ligation ELISA (PLEA) is a novel assay that combines the PLA with the ELISA in order to obtain quantitative affinity values of the interactions under study [14, 122]. PLEA offers the possibility to study the interaction of three molecules, by using three primary antibodies; the capture antibody used to capture the complexed set of proteins on a well plate and a set of two additional antibodies is used for the amplification of the PLA signal . As a result, this model describes how ATM targets the nascent p53 peptide and phosphorylates the S15, thus serving as a rapid mechanism whereby an activated p53 pool is synthesised in response to DNA damage . This comes in line with the notion that post-translational modifications on p53, including phosphorylation, acetylation, ubiquitination, sumoylation, neddylation and methylation may not only promote p53 stabilization but also confer p53 specificity to certain target genes , as well as regulate the specificity of post-translational modifications on p53, in response to other stresses, such as oxidative stress  or ER stress  and such studies are highly anticipated.
As a general conclusion from the described model, mutations on the mRNA can signal post-translational modifications of the encoded protein and may involve a fine-tuning activity by which the stabilization and the activity of the encoded protein towards certain stresses (in this case DNA damage), is decided and regulated at the cognate RNA level. This novel notion of functional aspects of silent mutations has only recently started to unravel . The description of the profound effects of synonymous mutations in the coding region of p53 regulating the interaction between the p53 mRNA and the MDM2 protein, opens up new processes by which p53 translation is controlled and emphasizes the role of the RNA itself in the regulation of p53 . Additional roles of p53 mRNA as a target for signalling pathways are listed in a recent review  and another example of RNA/protein coevolution comes for the ribosomal protein L22 and neighbouring 23S RNA. This work, based on a theoretic computational method, provided with direct evidence of RNA/protein coevolution driven by biophysical constraints on their interacting RNA and protein chains [128, 129]. It was shown that there is a co-evolutionary pattern (an increased co-evolution likelihood) on residue triplets, and not on duplets, for which mutations in an RNA nucleotide results in a change in residue distribution for proximal amino acids .
Functional consequences of co-evolved molecules
The p53-MDM2 co-evolutionary pattern described here, gives insights on the notion that certain genetically conserved signatures such as the p53(box-I) motif and the corresponding structural domains regulating its presentation to the interacting partner, such as the flanking sequence of the p53(TAD), evolve allosterically and intramolecularly. Extreme interchanging environments limiting the functional barriers while increasing the selective pressure, may favor strict cellular regulation and dedicated molecular networking. As described here, the temperature-regulated ancestral p53/63/73(TAD) mRNA structure, managing its interaction with the ancestral MDM2/X in the pre-vertebrate Ciona intestinalis, has co-evolved with the MDM family, to become regulated by MDMX in post-vertebrate species. This is an example of co-evolution between interacting and interdependent partners which is favored over the molecular evolution influenced by environmental inducers which may be interchangeable and constantly variable, increasing the risk of lethal phenotypes. The information gained by exploring the sequences of such molecular partners and how they co-evolved to become interdependent at a molecular (structural), cellular (expression, regulation) and functional level, may help in identifying the roles of specific genetic variations and could be considered as a major component of the driving force of the evolution.
In order to gain insights into the functional consequences of transcription factors and targets that modulate the expression of networks to lead the phenotypic diversification among species as well as being implicated in diseases, future studies should focus on mutations occurring both in exons and introns and aim to characterize such variations such as single nucleotide polymorphisms (SNPs) . The transcription factor p53 plays a key role in cancer suppression by preventing the proliferation of cells carrying potentially cancer-prone mutations. This function results from the co-evolution of an ancestral p53/p63/p73 molecule with MDM2 which shaped p53 superfamily to become integrated into a range of cellular pathways and molecular interactions, acquiring hundreds of distinct functions in metabolism, genome integrity, ageing, immune cell response and stress responses; it thus is not surprising that cancer cells require its inactivation. One important point highlighted here is that the co-evolution of p53 with its regulator, MDM2, lies on the evolutionary pressure on the BOX-I and the TA domain, both at the RNA and the protein level (Fig. 2, b).
An MDM2-dependent mechanism leading to the stabilisation of p53 was presented in a model by Karakostis et al,  describing downstream effects of the p53 mRNA - MDM2 interaction. The experimental evidences were derived by co-immunoprecipitation, sandwich ELISA and PLA experiments and were further confirmed by PLEA. PLEA confirmed that MDM2 and ATM are co-localised on purified ribosomes, actively translating p53 and that the nascent p53 peptide is phosphorylated on Ser15 epitope by ATM . This is in line with previously described ribonucleoprotein complexes composed of 5S RNA, RPL5 protein, MDM2 proteins, p53 protein  and with studies showing the roles of RPL11 and other ribosomal proteins on the regulation of MDM2 [131,132,133,134]. MDM2 is also implicated in sensing dysfunctional ribosomal biogenesis leading to the activation of p53 [135,136,137,138] and it will be of high interest to test the concepts of this model to study the mechanism whereby the ribosomal precursor complex is formed and how p53 is activated following ribosomal stress. This model explains previous puzzling observations using the artificial silent mutant p53TriM mRNA, which carries synonymous mutations in codons 17, 18 and 19 and has an increased affinity for MDM2 under normal conditions. p53TriM exhibits a higher rate of MDM2-dependent translation but also a higher rate of MDM2-dependent p53 degradation [14, 15]. Hence, according to the described mechanism, an increase in the MDM2 - p53 mRNA interaction facilitates the consequent increase in MDM2-dependent p53 synthesis. This is in line with p53TriM mRNA exhibiting an affinity for MDM2 which is not regulated by the DNA damage response and the ATM phosphorylations, explaining why the nascent p53 in the absence of ATM is efficiently degraded by MDM2 . The observed increased degradation can also be explained by the fact that once the phosphorylated MDM2(S395) is released from the ribosome and from the p53 mRNA which is being translated, it efficiently ubiquitinates the p53 protein that is not phosphorylated by ATM at S15 .
Overall, molecular cell biology and evolutionary theories are interlinked with the understanding of mechanisms involved in the development of genetic variations with functional outcomes. Such variations are promoted during the evolution as a result of adaptation, either directly towards environmental factors; or via interactions between co-evolved interacting molecular partners, as described for p53 and MDM2; or even driven by co-evolutionary interactions between species, as described in viruses and bacteria [139,140,141]. Molecular, biochemical and immunochemical tools for identifying and characterising the interaction of co-evolved variants, are currently being developed, improved and expanded. Techniques such as the DMS footprint for studying RNA structures when combined with the PLA and the PLEA, for studying dynamic molecular interactions of low abundance in situ at the sub-cellular level, open new horizons in studying in vitro and in cellulo molecular interactions and establish the possibility to functionally characterise the effect of genetic variations (mutations or polymorphisms) on the molecular interactions; and the co-evolution of the molecular partners. Such functional genetic variations represent key targets for understanding and treating genetic diseases such as cancer  as well as for understanding how synonymous mutations, which have largely been ignored, urge for a closer examination . Targeting these genetic variations may be extremely useful in sequencing-based approaches employed for genetic profiling in clinical diagnoses for personalized medicine. Considering variants of co-evolved partners, which may be located both on protein-coding regions (exons) of the genome as well as on intronic and non-coding regions, may be of crucial importance in developing efficient profiling molecular networks. Additionally, it is important for the employed bioinformatic analyses to accurately predict and interpret the effects of such variants involved in gene-expression or RNA splicing . Such studies focusing on synonymous variants and on the structural variations both on the mRNA and the encoded protein motifs, induced either by mutations or via evolutionary cues of co-evolved functional remnants from ancient interactions, may help in answering the question of how cancer shapes evolution and how evolution shapes cancer, and may be extremely useful in developing effective diagnoses and therapies by identifying key genetic signatures involved in genetic diseases (Fig. 4). The p53-MDM2 axis constitutes an ideal model with significant impact on cancer therapeutics  that may facilitate the understanding of various systems and interaction networks employing protein-mRNA interactions and refine predictive studies aiming to identify biomarkers of genetic diseases for translation research [47, 68, 145].
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
- Ci; C. intestinalis :
Enzyme-linked immunosorbent assay
Human immunodeficiency virus 1
Intrinsically disordered regions
Proximity ligation assay
Proximity ligation ELISA
Serine aa 395
Serine aa 15
Single nucleotide polymorphism
Lane DP, Cheok CF, Brown C, Madhumalar A, Ghadessy FJ, Verma C. Mdm2 and p53 are highly conserved from placozoans to man. Cell Cycle. 2010;9(3):540–7.
Abegglen LM, Caulin AF, Chan A, Lee K, Robinson R, Campbell MS, Kiso WK, Schmitt DL, Waddell PJ, Bhaskara S, et al. Potential mechanisms for Cancer resistance in elephants and comparative cellular response to DNA damage in humans. Jama. 2015;314(17):1850–60.
Coffill CR, Lee AP, Siau JW, Chee SM, Joseph TL, Tan YS, Madhumalar A, Tay BH, Brenner S, Verma CS, et al. The p53-Mdm2 interaction and the E3 ligase activity of Mdm2/Mdm4 are conserved from lampreys to humans. Genes Dev. 2016;30(3):281–92.
Karakostis K, Ponnuswamy A, Fusee LT, Bailly X, Laguerre L, Worall E, Vojtesek B, Nylander K, Fahraeus R. p53 mRNA and p53 protein structures have evolved independently to interact with MDM2. Mol Biol Evol. 2016;33(5):1280–92.
Aberg E, Saccoccia F, Grabherr M, Ore WYJ, Jemth P, Hultqvist G. Evolution of the p53-MDM2 pathway. BMC Evol Biol. 2017;17(1):177.
Canman CE, Lim DS, Cimprich KA, Taya Y, Tamai K, Sakaguchi K, Appella E, Kastan MB, Siliciano JD. Activation of the ATM kinase by ionizing radiation and phosphorylation of p53. Science. 1998;281(5383):1677–9.
Wallace M, Worrall E, Pettersson S, Hupp TR, Ball KL. Dual-site regulation of MDM2 E3-ubiquitin ligase activity. Mol Cell. 2006;23(2):251–63.
Maclaine NJ, Hupp TR. The regulation of p53 by phosphorylation: a model for how distinct signals integrate into the p53 pathway. Aging. 2009;1(5):490–502.
MacLaine NJ, Hupp TR. How phosphorylation controls p53. Cell Cycle. 2011;10(6):916–21.
Gajjar M, Candeias MM, Malbert-Colas L, Mazars A, Fujita J, Olivares-Illana V, Fahraeus R. The p53 mRNA-Mdm2 interaction controls Mdm2 nuclear trafficking and is required for p53 activation following DNA damage. Cancer Cell. 2012;21(1):25–35.
Ponnuswamy A, Hupp T, Fahraeus R. Concepts in MDM2 signaling: allosteric regulation and feedback loops. Genes Cancer. 2012;3(3–4):291–7.
Malbert-Colas L, Ponnuswamy A, Olivares-Illana V, Tournillon AS, Naski N, Fahraeus R. HDMX folds the nascent p53 mRNA following activation by the ATM kinase. Mol Cell. 2014;54(3):500–11.
Aymard F, Legube G. A TAD closer to ATM. Mol Cell Oncol. 2016;3(3):e1134411.
Karakostis K, Vadivel Gnanasundram S, Lopez I, Thermou A, Wang L, Nylander K, Olivares-Illana V, Fahraeus R. A single synonymous mutation determines the phosphorylation and stability of the nascent protein. J Mol Cell Biol. 2019;1;11(3):187-199.
Candeias MM, Malbert-Colas L, Powell DJ, Daskalogianni C, Maslon MM, Naski N, Bourougaa K, Calvo F, Fahraeus R. P53 mRNA controls p53 activity by managing Mdm2 functions. Nat Cell Biol. 2008;10(9):1098–105.
Medina-Medina I, Garcia-Beltran P, de la Mora-de la Mora I, Oria-Hernandez J, Millot G, Fahraeus R, Reyes-Vivas H, Sampedro JG, Olivares-Illana V. Allosteric interactions by p53 mRNA governs HDM2 E3 ubiquitin ligase specificity under different conditions. Mol Cell Biol. 2016;36(16):2195-205.
Dotsch V, Bernassola F, Coutandin D, Candi E, Melino G. p63 and p73, the ancestors of p53. Cold Spring Harb Perspect Biol. 2010;2(9):a004887.
Cheok CF, Verma CS, Baselga J, Lane DP. Translating p53 into the clinic. Nat Rev Clin Oncol. 2011;8(1):25–37.
Golomb L, Volarevic S, Oren M. p53 and ribosome biogenesis stress: the essentials. FEBS Lett. 2014;588(16):2571–9.
Kruiswijk F, Labuschagne CF, Vousden KH. p53 in survival, death and metabolic health: a lifeguard with a licence to kill. Nat Rev Mol Cell Biol. 2015;16(7):393–405.
Pressman A, Blanco C, Chen IA. The RNA world as a model system to study the origin of life. Curr Biol. 2015;25(19):R953–63.
Muchowska KB, Varma SJ, Moran J. Synthesis and breakdown of universal metabolic precursors promoted by iron. Nature. 2019;569(7754):104–7.
Gilbert W. Origin of life: the RNA world. Nature. 1986;319:618.
Woese CR, Dugre DH, Dugre SA, Kondo M, Saxinger WC. On the fundamental nature and evolution of the genetic code. Cold Spring Harb Symp Quant Biol. 1966;31:723–36.
Crick FH. The origin of the genetic code. J Mol Biol. 1968;38(3):367–79.
Powner MW, Gerland B, Sutherland JD. Synthesis of activated pyrimidine ribonucleotides in prebiotically plausible conditions. Nature. 2009;459(7244):239–42.
Patel BH, Percivalle C, Ritson DJ, Duffy CD, Sutherland JD. Common origins of RNA, protein and lipid precursors in a cyanosulfidic protometabolism. Nat Chem. 2015;7(4):301–7.
Blanco C, Bayas M, Yan F, Chen IA. Analysis of evolutionarily independent protein-RNA complexes yields a criterion to evaluate the relevance of prebiotic scenarios. Curr Biol. 2018;28(4):526–37 e525.
Wiegand TW, Janssen RC, Eaton BE. Selection of RNA amide synthases. Chem Biol. 1997;4(9):675–83.
Zhang B, Cech TR. Peptide bond formation by in vitro selected ribozymes. Nature. 1997;390(6655):96–100.
Tamura K. Origin of amino acid homochirality: relationship with the RNA world and origin of tRNA aminoacylation. Bio Systems. 2008;92(1):91–8.
Penny D. Evolutionary biology: relativity for molecular clocks. Nature. 2005;436(7048):183–4.
Szathmary E, Smith JM. The major evolutionary transitions. Nature. 1995;374(6519):227–32.
Butterfield GL, Lajoie MJ, Gustafson HH, Sellers DL, Nattermann U, Ellis D, Bale JB, Ke S, Lenz GH, Yehdego A, et al. Evolution of a designed protein assembly encapsulating its own RNA genome. Nature. 2017;552(7685):415–20.
Saito H, Yamada A, Ohmori R, Kato Y, Yamanaka T, Yoshikawa K, Inoue T: Towards constructing synthetic cells: RNA/RNP evolution and cell-free translational systems in giant liposomes. In: 2007 International Symposium on Micro-NanoMechatronics and Human Science: 11–14 Nov. 2007 2007; 2007: 286–291.
Yarus M, Caporaso JG, Knight R. Origins of the genetic code: the escaped triplet theory. Annu Rev Biochem. 2005;74:179–98.
Cojocaru R, Unrau PJ. Transitioning to DNA genomes in an RNA world. eLife. 2017;6:e32330.
Cech TR. The RNA worlds in context. Cold Spring Harb Perspect Biol. 2012;4(7):a006742.
Sloan KE, Bohnsack MT, Watkins NJ. The 5S RNP couples p53 homeostasis to ribosome biogenesis and nucleolar stress. Cell Rep. 2013;5(1):237–47.
Martin DP, Murrell B, Golden M, Khoosal A, Muhire B. RDP4: Detection and analysis of recombination patterns in virus genomes. Virus evolution. 2015;1(1):vev003.
Wang Y, Zirbel CL, Leontis NB, Ding B. RNA 3-dimensional structural motifs as a critical constraint of viroid RNA evolution. PLoS Pathog. 2018;14(2):e1006801.
van der Gulik PT, Speijer D. How amino acids and peptides shaped the RNA world. Life. 2015;5(1):230–46.
Kruger DM, Neubacher S, Grossmann TN. Protein-RNA interactions: structural characteristics and hotspot amino acids. Rna. 2018;24(11):1457–65.
Dominguez D, Freese P, Alexis MS, Su A, Hochman M, Palden T, Bazile C, Lambert NJ, Van Nostrand EL, Pratt GA, et al. Sequence, structure, and context preferences of human RNA binding proteins. Mol Cell. 2018;70(5):854–67 e859.
Hu W, Qin L, Li M, Pu X, Guo Y. A structural dissection of protein–RNA interactions based on different RNA base areas of interfaces. RSC Adv. 2018;8(19):10582–92.
Gu W, Li M, Xu Y, Wang T, Ko JH, Zhou T. The impact of RNA structure on coding sequence evolution in both bacteria and eukaryotes. BMC Evol Biol. 2014;14:87.
Sharma Y, Miladi M, Dukare S, Boulay K, Caudron-Herger M, Gross M, Backofen R, Diederichs S. A pan-cancer analysis of synonymous mutations. Nat Commun. 2019;10(1):2569.
Sanjuan R, Borderia AV. Interplay between RNA structure and protein evolution in HIV-1. Mol Biol Evol. 2011;28(4):1333–8.
Kudla G, Murray AW, Tollervey D, Plotkin JB. Coding-sequence determinants of gene expression in Escherichia coli. Science. 2009;324(5924):255–8.
Tucker BJ, Breaker RR. Riboswitches as versatile gene control elements. Curr Opin Struct Biol. 2005;15(3):342–8.
Kaempfer R. RNA sensors: novel regulators of gene expression. EMBO Rep. 2003;4(11):1043–7.
Watts JM, Dang KK, Gorelick RJ, Leonard CW, Bess JW Jr, Swanstrom R, Burch CL, Weeks KM. Architecture and secondary structure of an entire HIV-1 RNA genome. Nature. 2009;460(7256):711–6.
Berkhout B. HIV-1 as RNA evolution machine. RNA Biol. 2011;8(2):225–9.
Nagata S, Imai J, Makino G, Tomita M, Kanai A. Evolutionary analysis of HIV-1 pol proteins reveals representative residues for viral subtype differentiation. Front Microbiol. 2017;8:2151.
Munk C, Jensen BE, Zielonka J, Haussinger D, Kamp C. Running loose or getting lost: how HIV-1 counters and capitalizes on APOBEC3-induced mutagenesis through its Vif protein. Viruses. 2012;4(11):3132–61.
Jankowsky E, Harris ME. Specificity and nonspecificity in RNA-protein interactions. Nat Rev Mol Cell Biol. 2015;16(9):533–44.
Yang S, Wang J, Ng RT. Inferring RNA sequence preferences for poorly studied RNA-binding proteins based on co-evolution. BMC bioinformatics. 2018;19(1):96.
Lamech LT, Mallam AL, Lambowitz AM. Evolution of RNA-protein interactions: non-specific binding led to RNA splicing activity of fungal mitochondrial tyrosyl-tRNA synthetases. PLoS Biol. 2014;12(12):e1002028.
Hall KB. RNA-protein interactions. Curr Opin Struct Biol. 2002;12(3):283–8.
Gerstberger S, Hafner M, Tuschl T. A census of human RNA-binding proteins. Nat Rev Genet. 2014;15(12):829–45.
Hentze MW, Castello A, Schwarzl T, Preiss T. A brave new world of RNA-binding proteins. Nat Rev Mol Cell Biol. 2018;19(5):327–41.
Beckmann BM, Horos R, Fischer B, Castello A, Eichelbaum K, Alleaume AM, Schwarzl T, Curk T, Foehr S, Huber W, et al. The RNA-binding proteomes from yeast to man harbour conserved enigmRBPs. Nat Commun. 2015;6:10127.
Coppin L, Leclerc J, Vincent A, Porchet N, Pigny P. Messenger RNA Life-Cycle in Cancer Cells: Emerging Role of Conventional and Non-Conventional RNA-Binding Proteins? Int J Mol Sci. 2018;19(3):650.
Lucchesi C, Zhang J, Chen X. Modulation of the p53 family network by RNA-binding proteins. Transl Cancer Res. 2016;5(6):676–84.
Uversky VN. p53 Proteoforms and Intrinsic Disorder: An Illustration of the Protein Structure-Function Continuum Concept. Int J Mol Sci. 2016;17(11):1874.
Haronikova L, Olivares-Illana V, Wang L, Karakostis K, Chen S, Fahraeus R. The p53 mRNA: an integral part of the cellular stress response. Nucleic Acids Res. 2019;47(7):3257–71.
Zhou Q, Kunder N, De la Paz JA, Lasley AE, Bhat VD, Morcos F, Campbell ZT. Global pairwise RNA interaction landscapes reveal core features of protein recognition. Nat Commun. 2018;9(1):2511.
Zanzoni A, Spinelli L, Ribeiro DM, Tartaglia GG, Brun C. Post-transcriptional regulatory patterns revealed by protein-RNA interactions. Sci Rep. 2019;9(1):4302.
Geller R, Pechmann S, Acevedo A, Andino R, Frydman J. Hsp90 shapes protein and RNA evolution to balance trade-offs between protein stability and aggregation. Nat Commun. 2018;9(1):1781.
Rutherford SL, Lindquist S. Hsp90 as a capacitor for morphological evolution. Nature. 1998;396(6709):336–42.
Lovci MT, Bengtson MH, Massirer KB. Post-translational modifications and RNA-binding proteins. Adv Exp Med Biol. 2016;907:297–317.
Wallingford JB, Seufert DW, Virta VC, Vize PD. p53 activity is essential for normal development in Xenopus. Curr Biol. 1997;7(10):747–57.
Levrero M, De Laurenzi V, Costanzo A, Gong J, Wang JY, Melino G. The p53/p63/p73 family of transcription factors: overlapping and distinct functions. J Cell Sci. 2000;113(Pt 10):1661–70.
Moll UM, Petrenko O. The MDM2-p53 interaction. Mol Cancer Res. 2003;1(14):1001–8.
Lopez I, Tournillon AS, Prado Martins R, Karakostis K, Malbert-Colas L, Nylander K, Fahraeus R. p53-mediated suppression of BiP triggers BIK-induced apoptosis during prolonged endoplasmic reticulum stress. Cell Death Differ. 2017;24(10):1717–29.
Dereeper A, Audic S, Claverie JM, Blanc G. BLAST-EXPLORER helps you building datasets for phylogenetic analysis. BMC Evol Biol. 2010;10:8.
Harms KL, Chen X. The functional domains in p53 family proteins exhibit both common and distinct properties. Cell Death Differ. 2006;13(6):890–7.
Rutkowski R, Hofmann K, Gartner A. Phylogeny and function of the invertebrate p53 superfamily. Cold Spring Harb Perspect Biol. 2010;2(7):a001131.
Belyi VA, Ak P, Markert E, Wang H, Hu W, Puzio-Kuter A, Levine AJ. The origins and evolution of the p53 family of genes. Cold Spring Harb Perspect Biol. 2010;2(6):a001198.
Muttray AF, O'Toole TF, Morrill W, Van Beneden RJ, Baldwin SA. An invertebrate mdm homolog interacts with p53 and is differentially expressed together with p53 and ras in neoplastic Mytilus trossulus haemocytes. Comp Biochem Physiol B, Biochem Mol Biol. 2010;156(4):298–308.
Dos Santos HG, Nunez-Castilla J, Siltberg-Liberles J. Functional diversification after gene duplication: paralog specific regions of structural disorder and phosphorylation in p53, p63, and p73. PLoS One. 2016;11(3):e0151961.
Sulak M, Fong L, Mika K, Chigurupati S, Yon L, Mongan NP, Emes RD, Lynch VJ. Correction: TP53 copy number expansion is associated with the evolution of increased body size and an enhanced DNA damage response in elephants. eLife. 2016;5:e24307.
Schmitz JF, Zimmer F, Bornberg-Bauer E. Mechanisms of transcription factor evolution in Metazoa. Nucleic Acids Res. 2016;44(13):6287–97.
del Sol A, Carbonell P. The modular organization of domain structures: insights into protein-protein binding. PLoS Comput Biol. 2007;3(12):e239.
Wu YC, Rasmussen MD, Kellis M. Evolution at the subgene level: domain rearrangements in the Drosophila phylogeny. Mol Biol Evol. 2012;29(2):689–705.
Dohn M, Zhang S, Chen X. p63alpha and DeltaNp63alpha can induce cell cycle arrest and apoptosis and differentially regulate p53 target genes. Oncogene. 2001;20(25):3193–205.
Melino G, De Laurenzi V, Vousden KH. p73: Friend or foe in tumorigenesis. Nat Rev Cancer. 2002;2(8):605–15.
Kussie PH, Gorina S, Marechal V, Elenbaas B, Moreau J, Levine AJ, Pavletich NP. Structure of the MDM2 oncoprotein bound to the p53 tumor suppressor transactivation domain. Science. 1996;274(5289):948–53.
Honda R, Tanaka H, Yasuda H. Oncoprotein MDM2 is a ubiquitin ligase E3 for tumor suppressor p53. FEBS Lett. 1997;420(1):25–7.
Brady CA, Jiang D, Mello SS, Johnson TM, Jarvis LA, Kozak MM, Kenzelmann Broz D, Basak S, Park EJ, McLaughlin ME, et al. Distinct p53 transcriptional programs dictate acute DNA-damage responses and tumor suppression. Cell. 2011;145(4):571–83.
Sullivan KD, Galbraith MD, Andrysik Z, Espinosa JM. Mechanisms of transcriptional regulation by p53. Cell Death Differ. 2018;25(1):133–43.
Raj N, Attardi LD. The Transactivation Domains of the p53 Protein. Cold Spring Harb Perspect Med. 2017;7(1):a026047.
Medina-Medina I, Garcia-Beltran P, de la Mora-de la Mora I, Oria-Hernandez J, Millot G, Fahraeus R, Reyes-Vivas H, Sampedro JG, Olivares-Illana V. Allosteric interactions by p53 mRNA govern HDM2 E3 ubiquitin ligase specificity under different conditions. Mol Cell Biol. 2016;36(16):2195–205.
Hernandez-Monge J, Rousset-Roman AB, Medina-Medina I, Olivares-Illana V. Dual function of MDM2 and MDMX toward the tumor suppressors p53 and RB. Genes & cancer. 2016;7(9–10):278–87.
Banin S, Moyal L, Shieh S, Taya Y, Anderson CW, Chessa L, Smorodinsky NI, Prives C, Reiss Y, Shiloh Y, et al. Enhanced phosphorylation of p53 by ATM in response to DNA damage. Science. 1998;281(5383):1674–7.
Gannon HS, Woda BA, Jones SN. ATM phosphorylation of Mdm2 Ser394 regulates the amplitude and duration of the DNA damage response in mice. Cancer Cell. 2012;21(5):668–79.
Soderberg O, Gullberg M, Jarvius M, Ridderstrale K, Leuchowius KJ, Jarvius J, Wester K, Hydbring P, Bahram F, Larsson LG, et al. Direct observation of individual endogenous protein complexes in situ by proximity ligation. Nat Methods. 2006;3(12):995–1000.
Gullberg M, Goransson C, Fredriksson S. Duolink-"In-cell Co-IP" for visualization of protein interactions in situ. Nat Meth. 2011;8:982.
Koos B, Andersson L, Clausson CM, Grannas K, Klaesson A, Cane G, Soderberg O. Analysis of protein interactions in situ by proximity ligation assays. Curr Top Microbiol Immunol. 2014;377:111–26.
Tournillon AS, Lopez I, Malbert-Colas L, Findakly S, Naski N, Olivares-Illana V, Karakostis K, Vojtesek B, Nylander K, Fahraeus R. p53 binds the mdmx mRNA and controls its translation. Oncogene. 2017;2;36(5):723-730.
Oscier DG, Gardiner AC, Mould SJ, Glide S, Davis ZA, Ibbotson RE, Corcoran MM, Chapman RM, Thomas PW, Copplestone JA, et al. Multivariate analysis of prognostic factors in CLL: clinical stage, IGVH gene mutational status, and loss or mutation of the p53 gene are independent prognostic factors. Blood. 2002;100(4):1177–84.
Weibrecht I, Lundin E, Kiflemariam S, Mignardi M, Grundberg I, Larsson C, Koos B, Nilsson M, Soderberg O. In situ detection of individual mRNA molecules and protein complexes or post-translational modifications using padlock probes combined with the in situ proximity ligation assay. Nat Protoc. 2013;8(2):355–72.
Marz M, Vanzo N, Stadler PF. Temperature-dependent structural variability of RNAs: spliced leader RNAs and their evolutionary history. J Bioinforma Comput Biol. 2010;8(1):1–17.
Zilberstein D, Shapira M. The role of pH and temperature in the development of Leishmania parasites. Annu Rev Microbiol. 1994;48:449–70.
Close A. Development by degrees. Nurs Stand. 1999;13(41):59.
Bellas J, Beiras R, Vazquez E. A standardisation of Ciona intestinalis (Chordata, Ascidiacea) embryo-larval bioassay for ecotoxicological studies. Water Res. 2003;37(19):4613–22.
Kortmann J, Narberhaus F. Bacterial RNA thermometers: molecular zippers and switches. Nat Rev Microbiol. 2012;10(4):255–65.
Meyer M, Plass M, Perez-Valle J, Eyras E, Vilardell J. Deciphering 3'ss selection in the yeast genome reveals an RNA thermosensor that mediates alternative splicing. Mol Cell. 2011;43(6):1033–9.
Meyer E, Aglyamova GV, Matz MV. Profiling gene expression responses of coral larvae (Acropora millepora) to elevated temperature and settlement inducers using a novel RNA-Seq procedure. Mol Ecol. 2011;20(17):3599–616.
Vandivier L, Li F, Zheng Q, Willmann M, Chen Y, Gregory B. Arabidopsis mRNA secondary structure correlates with protein function and domains. Plant Signal Behav. 2013;8(6):e24301.
DeHart CJ, Chahal JS, Flint SJ, Perlman DH. Extensive post-translational modification of active and inactivated forms of endogenous p53. Mol Cell Proteomics. 2014;13(1):1–17.
Fan QD, Wu G, Liu ZR. Dynamics of posttranslational modifications of p53. Comput Math Methods Med. 2014;2014:245610.
Baugh EH, Ke H, Levine AJ, Bonneau RA, Chan CS. Why are there hotspot mutations in the TP53 gene in human cancers? Cell Death Differ. 2018;25(1):154–60.
Haupt Y, Maya R, Kazaz A, Oren M. Mdm2 promotes the rapid degradation of p53. Nature. 1997;387(6630):296–9.
Kubbutat MH, Jones SN, Vousden KH. Regulation of p53 stability by Mdm2. Nature. 1997;387(6630):299–303.
Maya R, Balass M, Kim ST, Shkedy D, Leal JF, Shifman O, Moas M, Buschmann T, Ronai Z, Shiloh Y, et al. ATM-dependent phosphorylation of Mdm2 on serine 395: role in p53 activation by DNA damage. Genes Dev. 2001;15(9):1067–77.
Naski N, Gajjar M, Bourougaa K, Malbert-Colas L, Fahraeus R, Candeias MM. The p53 mRNA-Mdm2 interaction. Cell Cycle. 2009;8(1):31–4.
Grover R, Sharathchandra A, Ponnuswamy A, Khan D, Das S. Effect of mutations on the p53 IRES RNA structure: implications for de-regulation of the synthesis of p53 isoforms. RNA Biol. 2011;8(1):132–42.
Swiatkowska A, Zydowicz P, Gorska A, Suchacka J, Dutkiewicz M, Ciesiolka J. The role of structural elements of the 5′-terminal region of p53 mRNA in translation under stress conditions assayed by the antisense oligonucleotide approach. PLoS One. 2015;10(10):e0141676.
Marechal V, Elenbaas B, Piette J, Nicolas JC, Levine AJ. The ribosomal L5 protein is associated with mdm-2 and mdm-2-p53 complexes. Mol Cell Biol. 1994;14(11):7414–20.
Daftuar L, Zhu Y, Jacq X, Prives C. Ribosomal proteins RPL37, RPS15 and RPS20 regulate the Mdm2-p53-MdmX network. PLoS One. 2013;8(7):e68667.
Tong QH, Tao T, Xie LQ, Lu HJ. ELISA-PLA: a novel hybrid platform for the rapid, highly sensitive and specific quantification of proteins and post-translational modifications. Biosens Bioelectron. 2016;80:385–91.
Marcel V, Nguyen Van Long F, Diaz JJ. 40 Years of Research Put p53 in Translation. Cancers. 2018;10(5):152.
Guo Z, Kozlov S, Lavin MF, Person MD, Paull TT. ATM activation by oxidative stress. Science. 2010;330(6003):517–21.
Stavridi ES, Halazonetis TD. p53 and stress in the ER. Genes Dev. 2004;18(3):241–4.
Fahraeus R, Marin M, Olivares-Illana V. Whisper mutations: cryptic messages within the genetic code. Oncogene. 2016;35(29):3753–9.
Lane DP, Verma CS. p53: updates on mechanisms, biology and therapy (I). J Mol Cell Biol. 2019;11(3):185–6.
Fares MA, Travers SA. A novel method for detecting intramolecular coevolution: adding a further dimension to selective constraints analyses. Genetics. 2006;173(1):9–23.
Brandman R, Brandman Y, Pande VS. Sequence coevolution between RNA and protein characterized by mutual information between residue triplets. PLoS One. 2012;7(1):e30022.
Jegga AG, Inga A, Menendez D, Aronow BJ, Resnick MA. Functional evolution of the p53 regulatory network through its target response elements. Proc Natl Acad Sci U S A. 2008;105(3):944–9.
Lohrum MA, Ludwig RL, Kubbutat MH, Hanlon M, Vousden KH. Regulation of HDM2 activity by the ribosomal protein L11. Cancer Cell. 2003;3(6):577–87.
Zheng J, Lang Y, Zhang Q, Cui D, Sun H, Jiang L, Chen Z, Zhang R, Gao Y, Tian W, et al. Structure of human MDM2 complexed with RPL11 reveals the molecular basis of p53 activation. Genes Dev. 2015;29(14):1524–34.
Liu Y, Deisenroth C, Zhang Y. RP-MDM2-p53 pathway: linking ribosomal biogenesis and tumor surveillance. Trends in cancer. 2016;2(4):191–204.
Sagar V, Caldarola S, Aria V, Monteleone V, Fuoco C, Gargioli C, Cannata S, Loreni F. PIM1 destabilization activates a p53-dependent response to ribosomal stress in cancer cells. Oncotarget. 2016;7(17):23837–49.
Lindstrom MS, Deisenroth C, Zhang Y. Putting a finger on growth surveillance: insight into MDM2 zinc finger-ribosomal protein interactions. Cell Cycle. 2007;6(4):434–7.
Lindstrom MS, Jin A, Deisenroth C, White Wolf G, Zhang Y. Cancer-associated mutations in the MDM2 zinc finger domain disrupt ribosomal protein interaction and attenuate MDM2-induced p53 degradation. Mol Cell Biol. 2007;27(3):1056–68.
Donati G, Peddigari S, Mercer CA, Thomas G. 5S ribosomal RNA is an essential component of a nascent ribosomal precursor complex that regulates the Hdm2-p53 checkpoint. Cell Rep. 2013;4(1):87–98.
Bursac S, Brdovcak MC, Donati G, Volarevic S. Activation of the tumor suppressor p53 upon impairment of ribosome biogenesis. Biochim Biophys Acta. 2014;1842(6):817–30.
Bottery MJ, Wood AJ, Brockhurst MA. Adaptive modulation of antibiotic resistance through intragenomic coevolution. Nature ecology & evolution. 2017;1(9):1364–9.
Huang W, Traulsen A, Werner B, Hiltunen T, Becks L. Dynamical trade-offs arise from antagonistic coevolution and decrease intraspecific diversity. Nat Commun. 2017;8(1):2059.
Bottery MJ, Wood AJ, Brockhurst MA. Temporal dynamics of bacteria-plasmid coevolution under antibiotic selection. ISME J. 2019;13:559–562.
Casas-Selves M, Degregori J. How cancer shapes evolution, and how evolution shapes cancer. Evolution. 2011;4(4):624–34.
Kamps R, Brandao RD, Bosch BJ, Paulussen AD, Xanthoulea S, Blok MJ, Romano A. Next-Generation Sequencing in Oncology: Genetic Diagnosis, Risk Prediction and Cancer Classification. Int J Mol Sci. 2017;18(2):308.
Hientz K, Mohr A, Bhakta-Guha D, Efferth T. The role of p53 in cancer drug resistance and targeted chemotherapy. Oncotarget. 2017;8(5):8921–46.
Franco M, Jeggari A, Peuget S, Bottger F, Selivanova G, Alexeyenko A. Prediction of response to anti-cancer drugs becomes robust via network integration of molecular data. Sci Rep. 2019;9(1):2379.
Cancerforskningsfonden Norr [AMP 16-823]; Cancerfonden (180296), Vetenskapsradet and MEYS-NPS [I-L01413].
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Karakostis, K., Fåhraeus, R. Shaping the regulation of the p53 mRNA tumour suppressor: the co-evolution of genetic signatures. BMC Cancer 19, 915 (2019). https://doi.org/10.1186/s12885-019-6118-y
- RNA world
- Ciona intestinalis
- Transcription factor
- Intrinsically disordered proteins
- Protein-RNA interactions
- mRNA translation
- Molecular basis of disease