Single-molecule analysis of interaction between p53TAD and MDM2 using aerolysin nanopores

Protein–protein interactions (PPIs) are regarded as important, but undruggable targets. Intrinsically disordered p53 transactivation domain (p53TAD) mediates PPI with mouse double minute 2 (MDM2), which is an attractive anticancer target for therapeutic intervention. Here, using aerolysin nanopores, we probed the p53TAD peptide/MDM2 interaction and its modulation by small-molecule PPI inhibitors or p53TAD phosphorylation. Although the p53TAD peptide showed short-lived (<100 ms) translocation, the protein complex induced the characteristic extraordinarily long-lived (0.1 s ∼ tens of min) current blockage, indicating that the MDM2 recruitment by p53TAD peptide almost fully occludes the pore. Simultaneously, the protein complex formation substantially reduced the event frequency of short-lived peptide translocation. Notably, the addition of small-molecule PPI inhibitors, Nutlin-3 and AMG232, or Thr18 phosphorylation of p53TAD peptide, were able to diminish the extraordinarily long-lived events and restore the short-lived translocation of the peptide rescued from the complex. Taken together, our results elucidate a novel mechanism of single-molecule sensing for analyzing PPIs and their inhibitors using aerolysin nanopores. This novel methodology may contribute to remarkable improvements in drug discovery targeted against undruggable PPIs.

Using aerolysin nanopores, we probed protein–protein interaction (PPI) between p53TAD and MDM2 and its modulation by small-molecule PPI inhibitors and p53TAD phosphorylation.     

Interplay between intrinsically disordered proteins inside membraneless protein liquid droplets

Intrinsically disordered proteins (IDPs) in cells phase separate to form diverse membraneless organelles, which have condensed liquid droplet-like properties and often contain multiple IDPs. However, how potential interactions between different IDPs affect the dynamic behavior of these protein droplets is largely unknown. Here, we develop a rapid IDP clustering system to generate protein droplets with varied residue compositions and examine diverse interacting IDPs inside droplets. Three different IDP droplets actively recruited other diverse IDPs inside droplets with extremely varied enrichment (inside/outside) degrees (over 100-fold variation) under highly crowded conditions. The recruited IDPs were mostly mobile even inside highly immobile droplets. Among the five tested IDPs, the disordered region of Ddx4 helicase with its unique multiple charged residue blocks was noticeably influenced by droplet mobility. We also discovered that droplets of different IDPs could rapidly fuse to each other. Interestingly, some droplets were heterogeneously fused with segregated subcompartments, and this segregation was enhanced by droplet maturation and was more apparent for specific IDP pairs, in which the polar and charged residue compositions are highly different. The present study not only reports multiple peculiar behaviors of interacting IDP pairs inside droplets but also provides valuable information on generating membraneless organelle models with controllable droplet properties.

Membraneless droplets of intrinsically disordered proteins (IDPs) with varied residue compositions uniquely interact with each other as droplets and clients.     

Phase separation of p53 precedes aggregation and is affected by oncogenic mutations and ligands

Mutant p53 tends to form aggregates with amyloid properties, especially amyloid oligomers inside the nucleus, which are believed to cause oncogenic gain-of-function (GoF). The mechanism of the formation of the aggregates in the nucleus remains uncertain. The present study demonstrated that the DNA-binding domain of p53 (p53C) underwent phase separation (PS) on the pathway to aggregation under various conditions. p53C phase separated in the presence of the crowding agent polyethylene glycol (PEG). Similarly, mutant p53C (M237I and R249S) underwent PS; however, the process evolved to a solid-like phase transition faster than that in the case of wild-type p53C. The data obtained by microscopy of live cells indicated that transfection of mutant full-length p53 into the cells tended to result in PS and phase transition (PT) in the nuclear compartments, which are likely the cause of the GoF effects. Fluorescence recovery after photobleaching (FRAP) experiments revealed liquid characteristics of the condensates in the nucleus. Mutant p53 tended to undergo gel- and solid-like phase transitions in the nucleus and in nuclear bodies demonstrated by slow and incomplete recovery of fluorescence after photobleaching. Polyanions, such as heparin and RNA, were able to modulate PS and PT in vitro. Heparin apparently stabilized the condensates in a gel-like state, and RNA apparently induced a solid-like state of the protein even in the absence of PEG. Conditions that destabilize p53C into a molten globule conformation also produced liquid droplets in the absence of crowding. The disordered transactivation domain (TAD) modulated both phase separation and amyloid aggregation. In summary, our data provide mechanistic insight into the formation of p53 condensates and conditions that may result in the formation of aggregated structures, such as mutant amyloid oligomers, in cancer. The pathway of mutant p53 from liquid droplets to gel-like and solid-like (amyloid) species may be a suitable target for anticancer therapy.

Mutant p53 tends to form aggregates with amyloid properties, especially amyloid oligomers inside the nucleus, which are believed to cause oncogenic gain-of-function (GoF).     

Determinants for intrinsically disordered protein recruitment into phase-separated protein condensates

Multivalent interactions between amino acid residues of intrinsically disordered proteins (IDPs) drive phase separation of these proteins into liquid condensates, forming various membrane-less organelles in cells. These interactions between often biased residues of IDPs are also likely involved in selective recruitment of many other IDPs into condensates. However, determining factors for this IDP recruitment into protein condensates are not understood yet. Here, we quantitatively examined recruitment tendencies of various IDPs with different sequence compositions into IDP-clustered condensates both in vitro as well as in cells. Condensate-forming IDP scaffolds, recruited IDP clients, and phase separation conditions were carefully varied to find key factors for selective IDP partitioning in protein condensates. Regardless of scaffold sequences, charged residues in client IDPs assured potent IDP recruitment, likely via strong electrostatic interactions, where positive residues could further enhance recruitment, possibly with cation–pi interactions. Notably, poly-ethylene glycol, a widely used crowding reagent for in vitro phase separation, abnormally increased IDP recruitment, indicating the need for careful use of crowding conditions. Tyrosines of IDP clients also strongly participated in recruitment both in vitro and in cells. Lastly, we measured recruitment degrees by more conventional interactions between folded proteins instead of disordered proteins. Surprisingly, recruitment forces by an even moderate protein interaction (K_d ∼ 5 μM) were substantially stronger than those by natural IDP–IDP interactions. The present data offer valuable information on how cells might organize protein partitioning on various protein condensates.

Diverse interactions between folded and disordered proteins collectively dictate selective protein recruitment into bimolecular condensates.     

How well does molecular simulation reproduce environment-specific conformations of the intrinsically disordered peptides PLP,TP2 and ONEG?

Understanding the conformational ensembles of intrinsically disordered proteins and peptides (IDPs) in their various biological environments is essential for understanding their mechanisms and functional roles in the proteome, leading to a greater knowledge of, and potential treatments for, a broad range of diseases. To determine whether molecular simulation is able to generate accurate conformational ensembles of IDPs, we explore the structural landscape of the PLP peptide (an intrinsically disordered region of the proteolipid membrane protein) in aqueous and membrane-mimicking solvents, using replica exchange with solute scaling (REST2), and examine the ability of four force fields (ff14SB, ff14IDPSFF, CHARMM36 and CHARMM36m) to reproduce literature circular dichroism (CD) data. Results from variable temperature (VT) ¹H and Rotating frame Overhauser Effect SpectroscopY (ROESY) nuclear magnetic resonance (NMR) experiments are also presented and are consistent with the structural observations obtained from the simulations and CD. We also apply the optimum simulation protocol to TP2 and ONEG (a cell-penetrating peptide (CPP) and a negative control peptide, respectively) to gain insight into the structural differences that may account for the observed difference in their membrane-penetrating abilities. Of the tested force fields, we find that CHARMM36 and CHARMM36m are best suited to the study of IDPs, and accurately predict a disordered to helical conformational transition of the PLP peptide accompanying the change from aqueous to membrane-mimicking solvents. We also identify an α-helical structure of TP2 in the membrane-mimicking solvents and provide a discussion of the mechanistic implications of this observation with reference to the previous literature on the peptide. From these results, we recommend the use of CHARMM36m with the REST2 protocol for the study of environment-specific IDP conformations. We believe that the simulation protocol will allow the study of a broad range of IDPs that undergo conformational transitions in different biological environments.

A protocol for simulating intrinsically disordered peptides in aqueous and hydrophobic solvents is proposed. Results from four force fields are compared with experiment. CHARMM36m performs the best for the simulated IDPs in all environments.     

Markov state models and NMR uncover an overlooked allosteric loop in p53

The tumor suppressor p53 is the most frequently mutated gene in human cancer, and thus reactivation of mutated p53 is a promising avenue for cancer therapy. Analysis of wildtype p53 and the Y220C cancer mutant long-timescale molecular dynamics simulations with Markov state models and validation by NMR relaxation studies has uncovered the involvement of loop L6 in the slowest motions of the protein. Due to its distant location from the DNA-binding surface, the conformational dynamics of this loop has so far remained largely unexplored. We observe mutation-induced stabilization of alternate L6 conformations, distinct from all experimentally-determined structures, in which the loop is both extended and located further away from the DNA-interacting surface. Additionally, the effect of the L6-adjacent Y220C mutation on the conformational landscape of the functionally-important loop L1 suggests an allosteric role to this dynamic loop and the inactivation mechanism of the mutation. Finally, the simulations reveal a novel Y220C cryptic pocket that can be targeted for p53 rescue efforts. Our approach exemplifies the power of the MSM methodology for uncovering intrinsic dynamic and kinetic differences among distinct protein ensembles, such as for the investigation of mutation effects on protein function.

Wildtype and Y220C L1 and L6 loops conformational landscape, with MSM-identified L6 states highlighted on the right.     

Molecular biosensing system based on intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) that undergo structural transition upon binding their target molecules are becoming increasingly known. IDPs, because of their binding specificity and induced folding properties, can serve as biological recognition elements for sensing applications. In this paper, BRCA1, an IDP, was utilized as the biological recognition element to detect tumor suppressor protein p53 through the BRCA1/p53 binding interaction to serve as a proof-of-concept for the use of IDPs as recognition elements. The binding resulted in a disordered-to-ordered BRCA1 conformational change, as seen in our circular dichroism (CD) measurements. This conformational change in BRCA1 (residues 219-498) was utilized in the detection of p53 (residues 311-393) via both intrinsic and extrinsic fluorescent probes. Intrinsic tryptophan residues within the BRCA1 sequence detected p53 (311-393) with a detection limit of 0.559 nM (0.112 pmol). Two environmentally sensitive fluorophores, tetramethylrhodamine-5-maleimide (TMR) and 6-((5-dimethylaminonaphthalene-1-sulfonyl)amino)hexanoic acid, succinimidyl ester (dansyl-X, SE) were conjugated to BRCA1 (219-498). Dansyl-X, SE-conjugated BRCA1 (219-498) detected p53 (311-393) with a detection limit of 1.50 nM (0.300 pmol). The sensitivities for TMR and dansyl-X, SE-conjugated BRCA1 for the detection of p53 were nearly threefold and twofold higher, respectively, than the sensitivity reported using intrinsic BRCA1 tryptophan fluorescence. CD measurements did not reveal a disruption of p53/dye-conjugated BRCA1 binding, thus validating the applicability of environmentally sensitive fluorophores as transduction moieties to detect molecules which bind to IDPs and induce a structural change.     

An Integrated Mass Spectrometry Based Approach to Probe the Structure of the Full‐Length Wild‐Type Tetrameric p53 Tumor Suppressor

We present an integrated approach for investigating the topology of proteins through native mass spectrometry (MS) and cross‐linking/MS, which we applied to the full‐length wild‐type p53 tetramer. For the first time, the two techniques were combined in one workflow to obtain not only structural insight in the p53 tetramer, but also information on the cross‐linking efficiency and the impact of cross‐linker modification on the conformation of an intrinsically disordered protein (IDP). P53 cross‐linking was monitored by native MS and as such, our strategy serves as a quality control for different cross‐linking reagents. Our approach can be applied to the structural investigation of various protein systems, including IDPs and large protein assemblies, which are challenging to study by the conventional methods used for protein structure characterization.     

The importance of the compact disordered state in the fuzzy interactions between intrinsically disordered proteins

The intrinsically disordered C-terminal domain (CTD) of protein 4.1G is able to specifically bind a 26-residue intrinsically disordered region of NuMA, forming a dynamic fuzzy complex. As one of a few cases of extremely fuzzy interactions between two intrinsically disordered proteins/regions (IDPs/IDRs) without induced folding, the principle of the binding is unknown. Here, we combined experimental and computational methods to explore the detailed mechanism of the interaction between 4.1G-CTD and NuMA. MD simulations suggest that the kinetic hub states in the structure ensemble of 4.1G-CTD are favorable in the fuzzy complex. The feature of these hub states is that the binding ‘hot spot’ motifs βA and βB exhibit β strand propensities and are well packed to each other. The binding between 4.1G-CTD and NuMA is disrupted at low pH, which changes the intramolecular packing of 4.1G-CTD and weakens the packing between βA and βB motifs. Low pH conditions also lead to increased hydrodynamic radius and acceleration of backbone dynamics of 4.1G-CTD. All these results underscore the importance of tertiary structural arrangements and overall compactness of 4.1G-CTD in its binding to NuMA, i.e. the compact disordered state of 4.1G-CTD is crucial for binding. Different from the short linear motifs (SLiMs) that are often found to mediate IDP interactions, 4.1G-CTD functions as an intrinsically disordered domain (IDD), which is a functional and structural unit similar to conventional protein domains. This work sheds light on the molecular recognition mechanism of IDPs/IDRs and expands the conventional structure-function paradigm in protein biochemistry.

Tertiary structural arrangements and overall compactness is important for interactions between IDPs.     

Semisynthetic ‘designer’ p53 sheds light on a phosphorylation–acetylation relay

The tumor suppressor protein p53 is a master regulator of cell fate. The activity of p53 is controlled by a plethora of posttranslational modifications (PTMs). However, despite extensive research, the mechanisms of this regulation are still poorly understood due to a paucity of biochemical studies with p53 carrying defined PTMs. Here, we report a protein semi-synthesis approach to access site-specifically modified p53. We synthesized a set of chemically homogeneous full-length p53 carrying one (Ser20ph and Ser15ph) or two (Ser15,20ph) naturally occurring, damage-associated phosphoryl marks. Refolding and biochemical characterization of semisynthetic p53 variants confirmed their structural and functional integrity. Furthermore, we show that phosphorylation within the N-terminal domain directly enhances p300-dependent acetylation approximately twofold, consistent with the role of these marks in p53 activation. Given that the p53 N-terminus is a hotspot for PTMs, we believe that our approach will contribute greatly to a mechanistic understanding of how p53 is controlled by PTMs.

‘Designer’ p53: semi-synthesis of the tumor suppressor protein p53 via native chemical ligation enables in vitro structure–activity studies to reveal how this master regulator of cell fate is itself regulated by phosphorylation.     

Anticancer gold(iii)-bisphosphine complex alters the mitochondrial electron transport chain to induce in vivo tumor inhibition

Expanding the chemical diversity of metal complexes provides a robust platform to generate functional bioactive reagents. To access an excellent repository of metal-based compounds for probe/drug discovery, we capitalized on the rich chemistry of gold to create organometallic gold(iii) compounds by ligand tuning. We obtained novel organogold(iii) compounds bearing a 1,2-bis(diphenylphosphino)benzene ligand, providing structural diversity with optimal physiological stability. Biological evaluation of the lead compound AuPhos-89 demonstrates mitochondrial complex I-mediated alteration of the mitochondrial electron transport chain (ETC) to drive respiration and diminish cellular energy in the form of adenosine triphosphate (ATP). Mechanism-of-action efforts, RNA-Seq, quantitative proteomics, and NCI-60 screening reveal a highly potent anticancer agent that modulates mitochondrial ETC. AuPhos-89 inhibits the tumor growth of metastatic triple negative breast cancer and represents a new strategy to study the modulation of mitochondrial respiration for the treatment of aggressive cancer and other disease states where mitochondria play a pivotal role in the pathobiology.

Expanding the chemical diversity of metal complexes provides a robust platform to generate functional bioactive reagents.     

Chemical control of peptide material phase transitions

Progressive solute-rich polymer phase transitions provide pathways for achieving ordered supramolecular assemblies. Intrinsically disordered protein domains specifically regulate information in biological networks via conformational ordering. Here we consider a molecular tagging strategy to control ordering transitions in polymeric materials and provide a proof-of-principle minimal peptide phase network captured with a dynamic chemical network.

Substrate initiated assembly of a dynamic chemical network.     

Probing conformational hotspots for the recognition and intervention of protein complexes by lysine reactivity profiling

Probing the conformational and functional hotspot sites within aqueous native protein complexes is still a challenging task. Herein, a mass spectrometry (MS)-based two-step isotope labeling-lysine reactivity profiling (TILLRP) strategy is developed to quantify the reactivities of lysine residues and probe the molecular details of protein–protein interactions as well as evaluate the conformational interventions by small-molecule active compounds. The hotspot lysine sites that are crucial to the SARS-CoV-2 S1–ACE2 combination could be successfully probed, such as S1 Lys⁴¹⁷ and Lys⁴⁴⁴. Significant alteration of the reactivities of lysine residues at the interaction interface of S1-RBD Lys³⁸⁶–Lys⁴⁶² was observed during the formation of complexes, which might be utilized as indicators for investigating the S1-ACE2 dynamic recognition and intervention at the molecular level in high throughput.

A mass spectrometry-based two-step isotope labeling-lysine reactivity profiling strategy is developed to probe the molecular details of protein–protein interactions and evaluate the conformational interventions by small-molecule active compounds.     

Accelerating high-throughput virtual screening through molecular pool-based active learning

Structure-based virtual screening is an important tool in early stage drug discovery that scores the interactions between a target protein and candidate ligands. As virtual libraries continue to grow (in excess of 10⁸ molecules), so too do the resources necessary to conduct exhaustive virtual screening campaigns on these libraries. However, Bayesian optimization techniques, previously employed in other scientific discovery problems, can aid in their exploration: a surrogate structure–property relationship model trained on the predicted affinities of a subset of the library can be applied to the remaining library members, allowing the least promising compounds to be excluded from evaluation. In this study, we explore the application of these techniques to computational docking datasets and assess the impact of surrogate model architecture, acquisition function, and acquisition batch size on optimization performance. We observe significant reductions in computational costs; for example, using a directed-message passing neural network we can identify 94.8% or 89.3% of the top-50 000 ligands in a 100M member library after testing only 2.4% of candidate ligands using an upper confidence bound or greedy acquisition strategy, respectively. Such model-guided searches mitigate the increasing computational costs of screening increasingly large virtual libraries and can accelerate high-throughput virtual screening campaigns with applications beyond docking.

Bayesian optimization can accelerate structure-based virtual screening campaigns by minimizing the total number of simulations performed while still identifying the vast majority of computational hits.     

An integrated in silico screening strategy for identifying promising disruptors of p53-MDM2 interaction

The p53 protein, also called guardian of the genome, plays a critical role in the cell cycle regulation and apoptosis. This protein is frequently inactivated in several types of human cancer by abnormally high levels of its negative regulator, mouse double minute 2 (MDM2). As a result, restoration of p53 function by inhibiting p53-MDM2 protein–protein interaction has been pursued as a compelling strategy for cancer therapy. To date, a limited number of small-molecules have been reported as effective p53−MDM2 inhibitors. X-ray structures of MDM2 in complex with some ligands are available in Protein Data Bank and herein, these data have been exploited to efficiently identify new p53-MDM2 interaction antagonists through a hierarchical virtual screening strategy. For this purpose, the first step was aimed at compiling a focused library of 686,630 structurally suitable compounds, from PubChem database, similar to two known effective inhibitors, Nutlin-3a and DP222669. These compounds were subjected to the subsequent structure-based approaches (quantum polarized ligand docking and molecular dynamics simulation) to select potential compounds with highest binding affinity for MDM2 protein. Additionally, ligand binding energy, ADMET properties and PAINS analysis were also considered as filtering criteria for selecting the most promising drug-like molecules. On the basis of these analyses, three top-ranked hit molecules, CID_118439641, CID_60452010 and CID_3106907, were found to have acceptable pharmacokinetics properties along with superior in silico inhibitory ability towards the p53-MDM2 interaction compared to known inhibitors. Molecular docking and molecular dynamics results well confirmed the interactions of the final selected compounds with critical residues within p53 binding site on the MDM2 hydrophobic clefts with satisfactory thermodynamics stability. Consequently, the new final scaffolds identified by the presented computational approach could offer a set of guidelines for designing promising anti-cancer agents targeting p53-MDM2 interaction.     

Molecular dynamics simulation of S100B protein to explore ligand blockage of the interaction with p53 protein

As a tumor suppressor, p53 plays an important role in cancer suppression. The biological function of p53 as a tumor suppressor is disabled when it binds to S100B. Developing the ligands to block the S100B-p53 interaction has been proposed as one of the most important approaches to the development of anti-cancer agents. We screened a small compound library against the binding interface of S100B and p53 to identify potential compounds to interfere with the interaction. The ligand-binding effect on the S100B-p53 interaction was explored by molecular dynamics at the atomic level. The results show that the ligand bound between S100B and p53 propels the two proteins apart by about 2 Å compared to the unligated S100B-p53 complex. The binding affinity of S100B and p53 decreases by ~8.5–14.6 kcal/mol after a ligand binds to the interface from the original unligated state of the S100B-p53 complex. Ligand-binding interferes with the interaction of S100B and p53. Such interference could impact the association of S100B and p53, which would free more p53 protein from the pairing with S100B and restore the biological function of p53 as a tumor suppressor. The analysis of the binding mode and ligand structural features would facilitate our effort to identify and design ligands to block S100B-p53 interaction effectively. The results from the work suggest that developing ligands targeting the interface of S100B and p53 could be a promising approach to recover the normal function of p53 as a tumor suppressor.     

Isothermal kinase-triggered supramolecular assemblies as drug sensitizers

Protein kinases, the main regulators of a vast map of cellular processes, are the most attractive targets in drug discovery. Despite a few successful examples of protein kinase inhibitors, the drug discovery strategy of downregulating protein kinase activity has been quite limited and often fails even in animal models. Here, we utilize protein kinase A (PKA) activity to design PKA-triggered supramolecular assemblies with anticancer activities. Grafting a suitable peptide to PNIPAM raises the critical temperature of the LCST polymer above body temperature. Interestingly, the corresponding phosphorylated polymer has a critical temperature below body temperature, making this peptide-appended PNIPAM a suitable polymer for the PKA-triggered supramolecular assembly process. PKA-triggered assembly occurs selectively in PKA-upregulated MCF-7 cells, which disturbs the cytoskeleton and sensitizes cancer cells against doxorubicin. The chemosensitization is also observed in vivo to identify effective tumor inhibitors with satisfactory biocompatibility. Overall, this phosphorylation-induced (in principle, PKA-catalyzed) supramolecular assembly opens up a promising chemotherapy strategy for combating kinase-upregulated cancer.

A nonapeptide grafted LCST polymer undergoes enzymatic phosphorylation to assemble, which selectively disrupts PKA overexpressing cancer cells via kinetics targeting.