Molecular dynamics study on the inhibition mechanisms of ReACp53 peptide for p53-R175H mutant aggregation

Multiple-omics analysis of aggrephagy-related cellular patterns and development of an aggrephagy-related signature for hepatocellular carcinoma

BACKGROUND

Protein aggrephagy, a selected autophagy process response for degrading protein aggregates, plays a critical role in various cancers. However, its regulatory mechanisms and clinical implications in hepatocellular carcinoma (HCC) remain largely unexplored.

METHODS

We integrated bulk RNA-seq data from TCGA and single-cell RNA sequencing (scRNA-seq) data from GEO databases to systematically analyze aggrephagy-related genes (AGGRGs) in HCC. Prognostic aggrephagy-related genes (AGGRGs) were identified through univariate Cox and LASSO regression analyses, followed by the construction of a risk prediction model. Patients were stratified into high- and low-risk groups based on the median risk score. Comparative analyses were performed to assess clinical outcomes, pathway enrichment, and drug sensitivity. Independent risk factors were incorporated a nomogram using univariate and multivariate Cox regression. At the single-cell level, the AGG scores were calculated using AUCell algorithm, and cell interactions and pseudotime trajectory analyses were conducted. Finally, protein levels of key AGGRG was assessed via tissue microarray.

RESULTS

Eight AGGRGs (PFKP, TPX2, UBE2S, GOT2, ST6GALNAC4, ADAM15, G6PD, and KPNA2) were identified as prognostic markers for HCC. The high-risk group exhibited significantly worse survival outcomes, heightened drug resistance, and enrichment in cell cycle, mTORC1 signaling, and reactive oxygen species pathways. Single-cell transcriptomic analysis revealed 11 distinct cell types within the HCC tumor microenvironment (TME), including hepatocytes, T cells, NK cells, macrophages, monocytes, dendritic cells, plasma B cells, mature B cells, mast cells, endothelial cells, and fibroblasts. Hepatocytes exhibited the highest AGGRG scores and were associated with metabolic reprograming, proliferation, and immune evasion. Further subclustering of malignant hepatocytes using inferCNV revealed eight functionally heterogeneous subpopulations with extensive intercellular crosstalk. Trajectory analysis showed G6PD- and CCNB1-expressing subpopulations in early-to-intermediate differentiation states, whereas C3 and ARGs marked terminal differentiation. Notably, G6PD was predominantly expressed in early and mid-stages, while KPNA2, PFKP, and TPX2 were upregulated in advanced tumor states. Immunohistochemical (IHC) validation confirmed significant overexpression of G6PD in HCC tissues compared to adjacent normal tissues.

CONCLUSION

These findings provide a molecular framework for targeting aggrephagy pathways in HCC treatment strategies.

Experimental kinetic mechanism of P53 condensation-amyloid aggregation

The tumor suppressor p53 modulates the transcription of a variety of genes, constituting a protective barrier against anomalous cellular proliferation. High-frequency "hotspot" mutations result in loss of function by the formation of amyloid-like aggregates that correlate with cancerous progression. We show that full-length p53 undergoes spontaneous homotypic condensation at submicromolar concentrations and in the absence of crowders to yield dynamic coacervates that are stoichiometrically dissolved by DNA. These coacervates fuse and evolve into hydrogel-like clusters with strong thioflavin T binding capacity, which further evolve into fibrillar species with a clearcut branching growth pattern. The amyloid-like coacervates can be rescued by the human papillomavirus master regulator E2 protein to yield large regular droplets. Furthermore, we kinetically dissected an overall condensation mechanism, which consists of a nucleation-growth process by the sequential addition of p53 tetramers, leading to discretely sized and monodisperse early condensates followed by coalescence into bead-like coacervates that slowly evolve to the fibrillar species. Our results suggest strong similarities to condensation-to-amyloid transitions observed in neurological aggregopathies. Mechanistic insights uncover novel key early and intermediate stages of condensation that can be targeted for p53 rescuing drug discovery.

Mechanistic insights into the allosteric inactivation mechanism of ZAP-70 induced by the hot spot W165C mutation

Zeta chain-associated protein kinase 70 (ZAP-70) is a non-receptor tyrosine kinase that interacts with the activated T-cell receptor to transduce downstream signals, and thus plays an important role in the adaptive immune system. The biphosphorylated immunotyrosine-based activation motifs (ITAM-Y2P) binds to the N-SH2 and C-SH2 domains of ZAP-70 to promote the activation of ZAP-70. The present study explores molecular mechanisms of allosteric inactivation of ZAP-70 induced by the hot spot W165C mutation through atomically detailed molecular dynamics simulation approaches. We report microsecond-length simulations of two states of the tandem SH2 domains of ZAP-70 in complex with the ITAM-Y2P motif, including the wild-type and W165C mutant. Extensive analysis of local flexibility and dynamical correlated motions show that W165C mutation changes coupled motions of protein domains and community networks. The binding affinities of the ITAM-Y2P motif to the wild-type and W165C mutant of ZAP-70 are predicted using binding free energy calculations. The results suggest that the driving force to decrease the binding affinity in the W165C mutant derives from the difference in the protein-protein electrostatic interactions. Moreover, the per-residue free energy decomposition unravels that the contributions from residues in the phosphorylated Tyr315 (pY315) binding site, in particular pY315 of ITAM-Y2P, and Arg43, Tyr240 of ZAP-70, are the key determinants for the loss of binding affinity. This study may insights into our understanding of the pathological mechanism of ZAP-70.Communicated by Ramaswamy H. Sarma.

Oncogenic R248W mutation induced conformational perturbation of the p53 core domain and the structural protection by proteomimetic amyloid inhibitor ADH-6

The involvement of p53 aggregation in cancer pathogenesis emphasizes the importance of unraveling the mechanisms underlying mutation-induced p53 destabilization. And understanding how small molecule inhibitors prevent the conversion of p53 into aggregation-primed conformations is pivotal for the development of therapeutics targeting p53-aggregation-associated cancers. A recent experimental study highlights the efficacy of the proteomimetic amyloid inhibitor ADH-6 in stabilizing R248W p53 and inhibiting its aggregation in cancer cells by interacting with the p53 core domain (p53C). However, it remains mostly unclear how R248W mutation induces destabilization of p53C and how ADH-6 stabilizes this p53C mutant and inhibits its aggregation. Herein, we conducted all-atom molecular dynamics simulations of R248W p53C in the absence and presence of ADH-6, as well as that of wild-type (WT) p53C. Our simulations reveal that the R248W mutation results in a shift of helix H2 and β-hairpin S2-S2' towards the mutation site, leading to the destruction of their neighboring β-sheet structure. This further facilitates the formation of a cavity in the hydrophobic core, and reduces the stability of the β-sandwich. Importantly, two crucial aggregation-prone regions (APRs) S9 and S10 are disturbed and more exposed to solvent in R248W p53C, which is conducive to p53C aggregation. Intriguingly, ADH-6 dynamically binds to the mutation site and multiple destabilized regions in R248W p53C, partially inhibiting the shift of helix H2 and β-hairpin S2-S2', thus preventing the disruption of the β-sheets and the formation of the cavity. ADH-6 also reduces the solvent exposure of APRs S9 and S10, which disfavors the aggregation of R248W p53C. Moreover, ADH-6 can preserve the WT-like dynamical network of R248W p53C. Our study elucidates the mechanisms underlying the oncogenic R248W mutation induced p53C destabilization and the structural protection of p53C by ADH-6.

Advanced computational approaches to understand protein aggregation

Protein aggregation is a widespread phenomenon implicated in debilitating diseases like Alzheimer's, Parkinson's, and cataracts, presenting complex hurdles for the field of molecular biology. In this review, we explore the evolving realm of computational methods and bioinformatics tools that have revolutionized our comprehension of protein aggregation. Beginning with a discussion of the multifaceted challenges associated with understanding this process and emphasizing the critical need for precise predictive tools, we highlight how computational techniques have become indispensable for understanding protein aggregation. We focus on molecular simulations, notably molecular dynamics (MD) simulations, spanning from atomistic to coarse-grained levels, which have emerged as pivotal tools in unraveling the complex dynamics governing protein aggregation in diseases such as cataracts, Alzheimer's, and Parkinson's. MD simulations provide microscopic insights into protein interactions and the subtleties of aggregation pathways, with advanced techniques like replica exchange molecular dynamics, Metadynamics (MetaD), and umbrella sampling enhancing our understanding by probing intricate energy landscapes and transition states. We delve into specific applications of MD simulations, elucidating the chaperone mechanism underlying cataract formation using Markov state modeling and the intricate pathways and interactions driving the toxic aggregate formation in Alzheimer's and Parkinson's disease. Transitioning we highlight how computational techniques, including bioinformatics, sequence analysis, structural data, machine learning algorithms, and artificial intelligence have become indispensable for predicting protein aggregation propensity and locating aggregation-prone regions within protein sequences. Throughout our exploration, we underscore the symbiotic relationship between computational approaches and empirical data, which has paved the way for potential therapeutic strategies against protein aggregation-related diseases. In conclusion, this review offers a comprehensive overview of advanced computational methodologies and bioinformatics tools that have catalyzed breakthroughs in unraveling the molecular basis of protein aggregation, with significant implications for clinical interventions, standing at the intersection of computational biology and experimental research.

Elucidating the Mechanisms of R248Q Mutation-Enhanced p53 Aggregation and Its Inhibition by Resveratrol

Aggregation of p53 mutants can result in loss-of-function, gain-of-function, and dominant-negative effects that contribute to tumor growth. Revealing the mechanisms underlying mutation-enhanced p53 aggregation and dissecting how small molecule inhibitors prevent the conversion of p53 into aggregation-primed conformations are fundamentally important for the development of novel therapeutics for p53 aggregation-associated cancers. A recent experimental study shows that resveratrol (RSV) has an inhibitory effect on the aggregation of hot-spot R248Q mutant of the p53 core domain (p53C), while pterostilbene (PT) exhibits a relatively poor inhibitory efficacy. However, the conformational properties of the R248Q mutant leading to its enhanced aggregation propensity and the inhibitory mechanism of RSV against p53C aggregation are not well understood. Herein, we performed extensive all-atom molecular dynamics simulations on R248Q p53C in the absence and presence of RSV/PT, as well as wild-type (WT) p53C. Our simulations reveal that loop L3, where the mutation resides, remains compact in WT p53C, while it becomes extended in the R248Q mutant. The extension of loop L3 weakens the interactions between loop L3 and two crucial aggregation-prone regions (APRs) of p53C, leading to impaired interactions within the APRs and their structural destabilization as well as p53C. The destabilized APRs in the R248Q mutant are more exposed than in WT p53C, which is conducive to p53C aggregation. RSV has a higher preference to bind to R248Q p53C than PT. This binding not only stabilizes loop L3 of R248Q mutant to its WT-like conformation, preventing L3-extension-caused APRs' destabilization but also reduces APRs' solvent exposure, thereby inhibiting R248Q p53C aggregation. However, PT exhibits a lower hydrogen-bonding capability and a higher self-association propensity, which would lead to a reduced p53C binding and a weakened inhibitory effect on R248Q mutant aggregation. Our study provides mechanistic insights into the R248Q mutation-enhanced aggregation propensity and RSV's potent inhibition against R248Q p53C aggregation.

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.

Insights into Allosteric Mechanisms of the Lung-Enriched p53 Mutants V157F and R158L

Lung cancer is a leading fatal malignancy in humans. p53 mutants exhibit not only loss of tumor suppressor capability but also oncogenic gain-of-function, contributing to lung cancer initiation, progression and therapeutic resistance. Research shows that p53 mutants V157F and R158L occur with high frequency in lung squamous cell carcinomas. Revealing their conformational dynamics is critical for developing novel lung therapies. Here, we used all-atom molecular dynamics (MD) simulations to investigate the effect of V157F and R158L substitutions on the structural properties of the p53 core domain (p53C). Compared to wild-type (WT) p53C, both V157F and R158L mutants display slightly lesser β-sheet structure, larger radius of gyration, larger volume and larger exposed surface area, showing aggregation-prone structural characteristics. The aggregation-prone fragments (residues 249–267 and 268–282) of two mutants are more exposed to water solution than that of WT p53C. V157F and R158L mutation sites can affect the conformation switch of loop 1 through long-range associations. Simulations also reveal that the local structure and conformation around the V157F and R158L mutation sites are in a dynamic equilibrium between the misfolded and properly folded conformations. These results provide molecular mechanistic insights into allosteric mechanisms of the lung-enriched p53 mutants.

Regulation of p53 Function by Formation of Non-Nuclear Heterologous Protein Complexes

A transcription factor p53 is activated upon cellular exposure to endogenous and exogenous stresses, triggering either homeostatic correction or cell death. Depending on the stress level, often measurable as DNA damage, the dual outcome is supported by p53 binding to a number of regulatory and metabolic proteins. Apart from the nucleus, p53 localizes to mitochondria, endoplasmic reticulum and cytosol. We consider non-nuclear heterologous protein complexes of p53, their structural determinants, regulatory post-translational modifications and the role in intricate p53 functions. The p53 heterologous complexes regulate the folding, trafficking and/or action of interacting partners in cellular compartments. Some of them mainly sequester p53 (HSP proteins, G6PD, LONP1) or its partners (RRM2B, PRKN) in specific locations. Formation of other complexes (with ATP2A2, ATP5PO, BAX, BCL2L1, CHCHD4, PPIF, POLG, SOD2, SSBP1, TFAM) depends on p53 upregulation according to the stress level. The p53 complexes with SIRT2, MUL1, USP7, TXN, PIN1 and PPIF control regulation of p53 function through post-translational modifications, such as lysine acetylation or ubiquitination, cysteine/cystine redox transformation and peptidyl-prolyl cis-trans isomerization. Redox sensitivity of p53 functions is supported by (i) thioredoxin-dependent reduction of p53 disulfides, (ii) inhibition of the thioredoxin-dependent deoxyribonucleotide synthesis by p53 binding to RRM2B and (iii) changed intracellular distribution of p53 through its oxidation by CHCHD4 in the mitochondrial intermembrane space. Increasing knowledge on the structure, function and (patho)physiological significance of the p53 heterologous complexes will enable a fine tuning of the settings-dependent p53 programs, using small molecule regulators of specific protein–protein interactions of p53.

Combining ReACp53 with Carboplatin to Target High-Grade Serous Ovarian Cancers

Ovarian malignancies are a leading cause of cancer-related death for US women. High-grade serous ovarian carcinomas (HGSOCs), the most common ovarian cancer subtype, are aggressive tumors with poor outcomes. Mutations in TP53 are common in HGSOCs, with a subset resulting in p53 aggregation and misregulation. ReACp53 is a peptide designed to inhibit mutant p53 aggregation and has been shown efficacious in targeting cancer cells in vitro and in vivo. As p53 regulates apoptosis, combining ReACp53 with carboplatin represents a logical therapeutic strategy. The efficacy of this combinatorial approach was tested in eight ovarian cancer cell lines and 10 patient HGSOC samples using an in vitro organoid drug assay, with the SynergyFinder tool utilized for calculating drug interactions. Results demonstrate that the addition of ReACp53 to carboplatin enhanced tumor cell targeting in the majority of samples tested, with synergistic effects measured in 2 samples, additivity measured in 14 samples, and antagonism measured in 1 sample. This combination was found to be synergistic in OVCAR3 ovarian cancer cells in vitro through enhanced apoptosis, and survival of mice bearing OVCAR3 intraperitoneal xenografts was extended when treated with the addition of ReACp53 to carboplatin versus carboplatin alone. Results suggest that carboplatin and ReACp53 may be a potential strategy in targeting a subset of HGSOCs.