Disruptor: Computational identification of oncogenic mutants disrupting protein-protein and protein-DNA interactions

We report an Osprey-based computational protocol to prospectively identify oncogenic mutations that act via disruption of molecular interactions. It is applicable to analyse both protein-protein and protein-DNA interfaces and it is validated on a dataset of clinically relevant mutations. In addition, it is used to predict previously uncharacterised patient mutations in CDK6 and p16 genes, which are experimentally confirmed to impair complex formation.

Targeting p53 pathways: mechanisms, structures, and advances in therapy

The TP53 tumor suppressor is the most frequently altered gene in human cancers, and has been a major focus of oncology research. The p53 protein is a transcription factor that can activate the expression of multiple target genes and plays critical roles in regulating cell cycle, apoptosis, and genomic stability, and is widely regarded as the "guardian of the genome". Accumulating evidence has shown that p53 also regulates cell metabolism, ferroptosis, tumor microenvironment, autophagy and so on, all of which contribute to tumor suppression. Mutations in TP53 not only impair its tumor suppressor function, but also confer oncogenic properties to p53 mutants. Since p53 is mutated and inactivated in most malignant tumors, it has been a very attractive target for developing new anti-cancer drugs. However, until recently, p53 was considered an "undruggable" target and little progress has been made with p53-targeted therapies. Here, we provide a systematic review of the diverse molecular mechanisms of the p53 signaling pathway and how TP53 mutations impact tumor progression. We also discuss key structural features of the p53 protein and its inactivation by oncogenic mutations. In addition, we review the efforts that have been made in p53-targeted therapies, and discuss the challenges that have been encountered in clinical development.

Irinotecan and vandetanib create synergies for treatment of pancreatic cancer patients with concomitant TP53 and KRAS mutations

BACKGROUND

The most frequently mutated gene pairs in pancreatic adenocarcinoma (PAAD) are KRAS and TP53, and our goal is to illustrate the multiomics and molecular dynamics landscapes of KRAS/TP53 mutation and also to obtain prospective novel drugs for KRAS- and TP53-mutated PAAD patients. Moreover, we also made an attempt to discover the probable link amid KRAS and TP53 on the basis of the abovementioned multiomics data.

METHOD

We utilized TCGA & Cancer Cell Line Encyclopedia data for the analysis of KRAS/TP53 mutation in a multiomics manner. In addition to that, we performed molecular dynamics analysis of KRAS and TP53 to produce mechanistic descriptions of particular mutations and carcinogenesis.

RESULT

We discover that there is a significant difference in the genomics, transcriptomics, methylomics, and molecular dynamics pattern of KRAS and TP53 mutation from the matching wild type in PAAD, and the prognosis of pancreatic cancer is directly linked with a particular mutation of KRAS and protein stability. Screened drugs are potentially effective in PAAD patients.

CONCLUSIONS

KRAS and TP53 prognosis of PAAD is directly associated with a specific mutation of KRAS. Irinotecan and vandetanib are prospective drugs for PAAD patients with KRASG12Dmutation and TP53 mutation.

Roles of computational modelling in understanding p53 structure, biology, and its therapeutic targeting

The transcription factor p53 plays pivotal roles in numerous biological processes, including the suppression of tumours. The rich availability of biophysical data aimed at understanding its structure–function relationships since the 1990s has enabled the application of a variety of computational modelling techniques towards the establishment of mechanistic models. Together they have provided deep insights into the structure, mechanics, energetics, and dynamics of p53. In parallel, the observation that mutations in p53 or changes in its associated pathways characterize several human cancers has resulted in a race to develop therapeutic modulators of p53, some of which have entered clinical trials. This review describes how computational modelling has played key roles in understanding structural-dynamic aspects of p53, formulating hypotheses about domains that are beyond current experimental investigations, and the development of therapeutic molecules that target the p53 pathway.

Variable Mutations at the p53-R273 Oncogenic Hotspot Position Leads to Altered Properties

Mutations in p53 protein, especially in the DNA-binding domain, is one of the major hallmarks of cancer. The R273 position is a DNA-contact position and has several oncogenic variants. Surprisingly, cancer patients carrying different mutant variants of R273 in p53 have different survival rates, indicating that the DNA-contact inhibition may not be the sole reason for reduced survival with R273 variants. Here, we probed the properties of three major oncogenic variants of the wild-type (WT) p53: [R273H]p53, [R273C]p53, and [R273L]p53. Using a series of biophysical, biochemical, and theoretical simulation studies, we observe that these oncogenic variants of the p53 not only suffer a loss in DNA binding, but they also show distinct structural stability, aggregation, and toxicity profiles. The WTp53 and the [R273H]p53 show the least destabilization and aggregation propensity. [R273C]p53 aggregation is disulfide mediated, leading to cross-β, thioflavin-T-positive aggregates, whereas hydrophobic interactions dominate self-assembly in [R273L]p53, leading to a mixture of amyloid and amorphous aggregates. Molecular dynamics simulations indicate different contact maps and secondary structures for the different variants along the course of the simulations. Our study indicates that each of the R273 variants has its own distinct property of stability and self-assembly, the molecular basis of which may lead to different types of cancer pathogenesis in vivo. These studies will aid the design of therapeutic strategies for cancer using residue-specific or process-specific protein aggregation as a target.

The cancer-associated, gain-of-function TP53 variant P152Lp53 activates multiple signaling pathways implicated in tumorigenesis

TP53 is the most frequently mutated tumor suppressor gene in many cancers, yet biochemical characterization of several of its reported mutations with probable biological significance have not been accomplished enough. Specifically, missense mutations in TP53 can contribute to tumorigenesis through gain-of-function of biochemical and biological properties that stimulate tumor growth. Here, we identified a relatively rare mutation leading to a proline to leucine substitution (P152L) in TP53 at the very end of its DNA-binding domain (DBD) in a sample from an Indian oral cancer patient. Although the P152Lp53 DBD alone bound to DNA, the full-length protein completely lacked binding ability at its cognate DNA motifs. Interestingly, P152Lp53 could efficiently tetramerize, and the mutation had only a limited impact on the structure and stability of full-length p53. Significantly, when we expressed this variant in a TP53-null cell line, it induced cell motility, proliferation, and invasion compared with a vector-only control. Also, enhanced tumorigenic potential was observed when P152Lp53-expressing cells were xenografted into nude mice. Investigating the effects of P152Lp53 expression on cellular pathways, we found that it is associated with up-regulation of several pathways, including cell-cell and cell-extracellular matrix signaling, epidermal growth factor receptor signaling, and Rho-GTPase signaling, commonly active in tumorigenesis and metastasis. Taken together, our findings provide a detailed account of the biochemical and cellular alterations associated with the cancer-associated P152Lp53 variant and establish it as a gain-of-function TP53 variant.

Simulations of mutant p53 DNA binding domains reveal a novel druggable pocket

The DNA binding domain (DBD) of the tumor suppressor p53 is the site of several oncogenic mutations. A subset of these mutations lowers the unfolding temperature of the DBD. Unfolding leads to the exposure of a hydrophobic β-strand and nucleates aggregation which results in pathologies through loss of function and dominant negative/gain of function effects. Inspired by the hypothesis that structural changes that are associated with events initiating unfolding in DBD are likely to present opportunities for inhibition, we investigate the dynamics of the wild type (WT) and some aggregating mutants through extensive all atom explicit solvent MD simulations. Simulations reveal differential conformational sampling between the WT and the mutants of a turn region (S6-S7) that is contiguous to a known aggregation-prone region (APR). The conformational properties of the S6-S7 turn appear to be modulated by a network of interacting residues. We speculate that changes that take place in this network as a result of the mutational stress result in the events that destabilize the DBD and initiate unfolding. These perturbations also result in the emergence of a novel pocket that appears to have druggable characteristics. FDA approved drugs are computationally screened against this pocket.

In silico identification of rescue sites by double force scanning

Motivation

A deleterious amino acid change in a protein can be compensated by a second-site rescue mutation. These compensatory mechanisms can be mimicked by drugs. In particular, the location of rescue mutations can be used to identify protein regions that can be targeted by small molecules to reactivate a damaged mutant.

Results

We present the first general computational method to detect rescue sites. By mimicking the effect of mutations through the application of forces, the double force scanning (DFS) method identifies the second-site residues that make the protein structure most resilient to the effect of pathogenic mutations. We tested DFS predictions against two datasets containing experimentally validated and putative evolutionary-related rescue sites. A remarkably good agreement was found between predictions and experimental data. Indeed, almost half of the rescue sites in p53 was correctly predicted by DFS, with 65% of remaining sites in contact with DFS predictions. Similar results were found for other proteins in the evolutionary dataset.

Availability and implementation

The DFS code is available under GPL at https://fornililab.github.io/dfs/.

Supplementary information

Supplementary data are available at Bioinformatics online.

R248Q mutation--Beyond p53-DNA binding

R248 in the DNA binding domain (DBD) of p53 interacts directly with the minor groove of DNA. Earlier nuclear magnetic resonance (NMR) studies indicated that the R248Q mutation resulted in conformation changes in parts of DBD far from the mutation site. However, how information propagates from the mutation site to the rest of the DBD is still not well understood. We performed a series of all-atom molecular dynamics (MD) simulations to dissect sterics and charge effects of R248 on p53-DBD conformation: (i) wild-type p53 DBD; (ii) p53 DBD with an electrically neutral arginine side-chain; (iii) p53 DBD with R248A; (iv) p53 DBD with R248W; and (v) p53 DBD with R248Q. Our results agree well with experimental observations of global conformational changes induced by the R248Q mutation. Our simulations suggest that both charge- and sterics are important in the dynamics of the loop (L3) where the mutation resides. We show that helix 2 (H2) dynamics is altered as a result of a change in the hydrogen bonding partner of D281. In turn, neighboring L1 dynamics is altered: in mutants, L1 predominantly adopts the recessed conformation and is unable to interact with the major groove of DNA. We focused our attention the R248Q mutant that is commonly found in a wide range of cancer and observed changes at the zinc-binding pocket that might account for the dominant negative effects of R248Q. Furthermore, in our simulations, the S6/S7 turn was more frequently solvent exposed in R248Q, suggesting that there is a greater tendency of R248Q to partially unfold and possibly lead to an increased aggregation propensity. Finally, based on the observations made in our simulations, we propose strategies for the rescue of R248Q mutants.

Effects of p21-activated kinase 1 inhibition on 11q13-amplified ovarian cancer cells

p21-activated kinases (PAKs) are Cdc42/Rac–activated serine-threonine protein kinases that regulate of several key cancer-relevant signaling pathways, such as the Mek/Erk, PI3K/Akt, and Wnt/b-catenin signaling pathways. Pak1 is frequently overexpressed and/or hyperactivated in different human cancers, including human breast, ovary, prostate, and brain cancer, due to amplification of the PAK1 gene in an 11q13 amplicon. Genetic or pharmacological inactivation of Pak1 has been shown to reduce proliferation of different cancer cells in vitro and reduce tumor progression in vivo.

In this work, we examined the roles of Pak1 in cellular and animal models of PAK1-amplified ovarian cancer. We found that inhibition of Pak1 leads to decreased proliferation and migration in PAK1 amplified/overexpressed ovarian cancer cells, and has no effect in cell that lack such amplification/overexpression. Further, we observed that loss of Pak1 function causes 11q13 amplified ovarian cancer cells to arrest in the G2/M phase of the cell cycle. This arrest correlates with activation of p53 and p21^Cip and decreased expression of cyclin B1. These findings suggest that small molecule inhibitors of Pak1 may play a therapeutic role in the ~25% of ovarian cancers characterized by PAK1 gene amplification.

Structural and sequential context of p53: A review of experimental and theoretical evidence

Approximately 27 million people are suffering from cancer that contains either an inactivating missense mutation of TP53 gene or partially abrogated p53 signaling pathway. Concerted action of folded and intrinsically disordered domains accounts for multi-faceted role of p53. The intricacy of dynamic p53 structure is believed to shed light on its cellular activity for developing new cancer therapies. In this review, insights into structural details of p53, diverse single point mutations affecting its core domain, thermodynamic understanding and therapeutic strategies for pharmacological rescue of p53 function has been illustrated. An effort has been made here to bridge the structural and sequential evidence of p53 from experimental to computational studies. First, we focused on the individual domains and the crucial protein-protein or DNA-protein contacts that determine conformation and dynamic behavior of p53. Next, the oncogenic mutations associated with cancer and its contribution to thermodynamic fluctuation has been discussed. Thus the emerging anti-cancer strategies include targeting of destabilized cancer mutants with selective inhibition of its negative regulators. Recent advances in development of small molecule inhibitors and peptides exploiting p53-MDM2 interaction has been included. In a nutshell, this review attempts to describe structural biology of p53 which provide new openings for structure-guided rescue.

Mapping the structural and dynamical features of multiple p53 DNA binding domains: insights into loop 1 intrinsic dynamics

The transcription factor p53 regulates cellular integrity in response to stress. p53 is mutated in more than half of cancerous cells, with a majority of the mutations localized to the DNA binding domain (DBD). In order to map the structural and dynamical features of the DBD, we carried out multiple copy molecular dynamics simulations (totaling 0.8 μs). Simulations show the loop 1 to be the most dynamic element among the DNA-contacting loops (loops 1-3). Loop 1 occupies two major conformational states: extended and recessed; the former but not the latter displays correlations in atomic fluctuations with those of loop 2 (~24 Å apart). Since loop 1 binds to the major groove whereas loop 2 binds to the minor groove of DNA, our results begin to provide some insight into the possible mechanism underpinning the cooperative nature of DBD binding to DNA. We propose (1) a novel mechanism underlying the dynamics of loop 1 and the possible tread-milling of p53 on DNA and (2) possible mutations on loop 1 residues to restore the transcriptional activity of an oncogenic mutation at a distant site.

Anti-apoptotic roles for the mutant p53R248Q through suppression of p53-regulated apoptosis-inducing protein 1 in the RA-derived fibroblast-like synoviocyte cell line MH7A

We previously reported that somatic mutations in the p53 gene accumulated at a higher frequency in AID(activation induced cytidine deaminase)(+) RA-FLS, which may result in the malfunction of p53, causing the tumor-like properties of RA-FLS. Among the p53 mutations identified from 3 sources of AID(+) RA-FLS, we focused on the p53R248Q mutation because it was reported to enhance the invasiveness of lung cancer cells and to have dominant-negative activity for pro-apoptotic molecules. We obtained cDNA encoding the p53R248Q mutant and introduced it into the MH7A RA-FLS cell line. P53R248Q dramatically suppressed the expression of the pro-apoptotic molecule p53AIP1 even under oxidative stress, which normally upregulates p53AIP1, leading to apoptosis. Moreover, overexpression of p53AIP1 increased apoptosis, whereas p53AIP1 knockdown rescued the cells from apoptosis. Together, these studies indicate the critical role of p53AIP1 under DNA damaging stresses for cell fate determination in RA-FLS containing the p53R248Q mutation.

Structure and stability insights into tumour suppressor p53 evolutionary related proteins

The p53 family of genes and their protein products, namely, p53, p63 and p73, have over one billion years of evolutionary history. Advances in computational biology and genomics are enabling studies of the complexities of the molecular evolution of p53 protein family to decipher the underpinnings of key biological conditions spanning from cancer through to various metabolic and developmental disorders and facilitate the design of personalised medicines. However, a complete understanding of the inherent nature of the thermodynamic and structural stability of the p53 protein family is still lacking. This is due, to a degree, to the lack of comprehensive structural information for a large number of homologous proteins and to an incomplete knowledge of the intrinsic factors responsible for their stability and how these might influence function. Here we investigate the thermal stability, secondary structure and folding properties of the DNA-binding domains (DBDs) of a range of proteins from the p53 family using biophysical methods. While the N- and the C-terminal domains of the p53 family show sequence diversity and are normally targets for post-translational modifications and alternative splicing, the central DBD is highly conserved. Together with data obtained from Molecular Dynamics simulations in solution and with structure based homology modelling, our results provide further insights into the molecular properties of evolutionary related p53 proteins. We identify some marked structural differences within the p53 family, which could account for the divergence in biological functions as well as the subtleties manifested in the oligomerization properties of this family.

Rescue of deleterious mutations by the compensatory Y30F mutation in ketosteroid isomerase

Proteins have evolved to compensate for detrimental mutations. However, compensatory mechanisms for protein defects are not well understood. Using ketosteroid isomerase (KSI), we investigated how second-site mutations could recover defective mutant function and stability. Previous results revealed that the Y30F mutation rescued the Y14F, Y55F and Y14F/Y55F mutants by increasing the catalytic activity by 23-, 3- and 1.3-fold, respectively, and the Y55F mutant by increasing the stability by 3.3 kcal/mol. To better understand these observations, we systematically investigated detailed structural and thermodynamic effects of the Y30F mutation on these mutants. Crystal structures of the Y14F/Y30F and Y14F/Y55F mutants were solved at 2.0 and 1.8 previoulsy solved structures of wild-type and other mutant KSIs. Structural analyses revealed that the Y30F mutation partially restored the active-site cleft of these mutant KSIs. The Y30F mutation also increased Y14F and Y14F/Y55F mutant stability by 3.2 and 4.3 kcal/mol, respectively, and the melting temperatures of the Y14F, Y55F and Y14F/Y55F mutants by 6.4°C, 5.1°C and 10.0°C, respectively. Compensatory effects of the Y30F mutation on stability might be due to improved hydrophobic interactions because removal of a hydroxyl group from Tyr30 induced local compaction by neighboring residue movement and enhanced interactions with surrounding hydrophobic residues in the active site. Taken together, our results suggest that perturbed active-site geometry recovery and favorable hydrophobic interactions mediate the role of Y30F as a secondsite suppressor.

A novel p53 mutant found in iatrogenic urothelial cancers is dysfunctional and can be rescued by a second-site global suppressor mutation

Exposure to herbal remedies containing the carcinogen aristolochic acid (AA) has been widespread in some regions of the world. Rare A→T TP53 mutations were recently discovered in AA-associated urothelial cancers. The near absence of these mutations among all other sequenced human tumors suggests that they could be biologically silent. There are no cell banks with established lines derived from human tumors with which to explore the influence of the novel mutants on p53 function and cellular behavior. To investigate their impact, we generated isogenic mutant clones by integrase-mediated cassette exchange at the p53 locus of platform (null) murine embryonic fibroblasts and kidney epithelial cells. Common tumor mutants (R248W, R273C) were compared with the AA-associated mutants N131Y, R249W, and Q104L. Assays of cell proliferation, migration, growth in soft agar, apoptosis, senescence, and gene expression revealed contrasting outcomes on cellular behavior following introduction of N131Y or Q104L. The N131Y mutant demonstrated a phenotype akin to common tumor mutants, whereas Q104L clone behavior resembled that of cells with wild-type p53. Wild-type p53 responses were restored in double-mutant cells harboring N131Y and N239Y, a second-site rescue mutation, suggesting that pharmaceutical reactivation of p53 function in tumors expressing N131Y could have therapeutic benefit. N131Y is likely to contribute directly to tumor phenotype and is a promising candidate biomarker of AA exposure and disease. Rare mutations thus do not necessarily point to sites where amino acid exchanges are phenotypically neutral. Encounter with mutagenic insults targeting cryptic sites can reveal specific signature hotspots.

The substitutions G245C and G245D in the Zn(2+)-binding pocket of the p53 protein result in differences of conformational flexibility of the DNA-binding domain

Transcription activation of the proapoptotic target genes is a means by which the p53 protein implements its function of tumor suppression. Zn(2+) is a known regulator of p53 binding to the target genes. We have previously obtained an evidence that amino acid substitutions in the p53 Zn(2+)-binding pocket can presumably exert an influence on Zn(2+) position in the Zn(2+)-p53 complex and thereby affect p53 binding to DNA. With these background considerations, our aim was to estimate the effect of the putative changes in the Zn(2+) position in its binding pocket due to the G245C and G245D substitutions on the conformation of the p53 DNA-binding motif. Statistical analysis of the molecular dynamics (MD) trajectories of the mutant p53-Zn(2+) complexes was used to detect significant deviations in conformation of the mutant p53 forms. MD simulations demonstrated that (1) the two substitutions in the Zn(2+)-binding pocket caused changes in the conformation of the p53 DNA-binding motif, as compared with the wild-type (WT) p53; (2) binding of Zn(2+) to the p53 mutant forms reduced the effect of the substitutions on conformational change; and (3) Zn(2+) binding in the normal position compensated the effect of the mutations on the conformation in comparison to the altered Zn(2+) position.

Structural effects of the L145Q, V157F, and R282W cancer-associated mutations in the p53 DNA-binding core domain

The p53 tumor suppressor is a transcription factor involved in many important signaling pathways, such as apoptosis and cell-cycle arrest. In over half of human cancers, p53 function is compromised by a mutation in its gene. Mutations in the p53 DNA-binding core domain destabilize the structure and reduce DNA-binding activity. We performed molecular dynamics simulations at physiological temperature to study the structural and dynamic effects of the L145Q, V157F, and R282W cancer-associated mutations in comparison to the wild-type protein. While there were common regions of destabilization in the mutant simulations, structural changes particular to individual mutations were also observed. Significant backbone deviations of the H2 helix and S7-S8 loop were observed in all mutant simulations; the H2 helix binds to DNA. In addition, the L145Q and V157F mutations, which are located in the β-sandwich core of the domain, disrupted the β-sheet structure and the loop-sheet-helix motif. The R282W mutation caused distortion of the loop-sheet-helix motif, but otherwise this mutant was similar to the wild-type structure. The introduction of these mutations caused rearrangement of the DNA-binding surface, consistent with their reduced DNA-binding activity. The simulations reveal detailed effects of the mutations on the stability and dynamics of p53 that may provide insight for therapeutic approaches.

Mutant p53 R248Q but not R248W enhances in vitro invasiveness of human lung cancer NCI-H1299 cells

More than half of all human cancers are associated with mutations of the TP53 gene. In regard to the functional interaction with the remaining wild-type (WT) p53 allele, p53 mutations are classified into two types, recessive and dominant-negative (DN) mutations. The latter mutant protein has a DN activity over the remaining WT allele. We previously showed that the DN p53 mutant was useful as a predictor of poor outcome or a risk factor for metastatic recurrence in patients with some types of cancers, regardless of the presence or absence of loss of heterozygosity (LOH) of WT p53, suggesting that the DN p53 had 'gain-of-function (GOF)' activity besides the transdominance function. In this study, we investigated GOF activity of two DN p53 mutants which had a point mutation at codon 248 (R248Q and R248W), one of the hot spots, by transfecting them respectively into H1299 cells which originally expressed no p53 protein. Growth activity of the transfectants with the two mutants was not different from that of parent or Mock transfectants. Meanwhile, in vitro invasions of Matrigel and type I collagen gel by R248Q-transfectants were significantly higher than those by R248W-transfectants or the control cells. However, there were no differences in cell motile activities, expressions of extracellular matrix-degradative enzymes such as matrix metalloproteinases, urokinase-type plasminogen activator and heparanase, and their inhibitors, between R248Q- and R248W-transfectants. These findings indicate that the p53 mutants have a different quality in GOF activities even if the mutations occurred at the same codon. And detailed information of the status of p53, including transdominancy and GOF activity, is expected to be useful for diagnosis and therapeutic strategy fitting the individual patients.