A Novel Determinant of PSMD9 PDZ Binding Guides the Evolution of the First Generation of Super Binding Peptides

The interaction network of the proteasome assembly chaperone PSMD9 regulates proteostasis

Functional networks in cells are created by physical, genetic, and regulatory interactions. Mapping them and annotating their functions by available methods remains a challenge. We use affinity purification mass spectrometry (AP-MS) coupled with SLiMFinder to discern such a network involving 26S proteasome non-ATPase regulatory subunit 9 (PSMD9), a chaperone of proteasome assembly. Approximately 20% of proteins within the PSMD9 interactome carry a short linear motif (SLiM) of the type 'EXKK'. The binding of purified PSMD9 with the peptide sequence ERKK, proteins heterogeneous nuclear ribonucleoproteins A2/B1 (hnRNPA2B1; containing ERKK), and peroxiredoxin-6 (PRDX6; containing EAKK) provided proof of principle for this motif-driven network. The EXKK motif in the peptide primarily interacts with the coiled-coil N domain of PSMD9, a unique interaction not reported for any coiled-coil domain. PSMD9 knockout (KO) HEK293 cells experience endoplasmic reticulum (ER) stress and respond by increasing the unfolded protein response (UPR) and reducing the formation of aggresomes and lipid droplets. Trans-expression of PSMD9 in the KO cells rescues lipid droplet formation. Overexpression of PSMD9 in HEK293 cells results in reduced UPR, and increased lipid droplet and aggresome formation. The outcome argues for the prominent role of PSMD9 in maintaining proteostasis. Probable mechanisms involve the binding of PSMD9 to binding immunoglobulin protein (BIP/GRP78; containing EDKK), an endoplasmic reticulum chaperone and key regulator of the UPR, and fatty acid synthase (FASN; containing ELKK), involved in fatty acid synthesis/lipid biogenesis. We propose that PSMD9 acts as a buffer in the cellular milieu by moderating the UPR and enhancing aggresome formation to reduce stress-induced proteotoxicity. Akin to waves created in ponds that perpetuate to a distance, perturbing the levels of PSMD9 would cause ripples down the networks, affecting final reactions in the pathway, one of which is altered proteostasis.

How Peptides Bind to PSD-95/Discs-Large/ZO-1 Domains

PSD-95/discs-large/ZO-1 (PDZ) domains form a large family of adaptor proteins that bind to the C-terminal tails of their binding partner proteins. Via extensive molecular dynamics simulations and alchemical free energy calculations, we characterized the binding modi of phosphorylated and unphosphorylated EQVSAV peptides and of a EQVEAV phosphate mimic to the hPTP1E PDZ2 and MAGI1 PDZ1 domains. The simulations reproduced the well-known binding characteristics such as tight coordination of the peptidic carboxyl tail and pronounced hydrogen bonding between the peptide backbone and the backbone atoms of a β-sheet in PDZ. Overall, coordination by hPTP1E PDZ2 appeared tighter than by MAGI1 PDZ1. Simulations of wild-type PDZ and arginine mutants suggest that contacts with Arg79/85 in hPTP1E/MAGI1 are more important for the EQVEAV peptide than for EQVSAV. Alchemical free energy calculations and PaCS-MD simulations could well reproduce the difference in binding free energy between unphosphorylated EQVSAV and EQVEAV peptides and the absolute binding free energy of EQVSAV. However, likely due to small force field inaccuracies, the simulations erroneously favored binding of the phosphorylated peptide instead of its unphosphorylated counterpart, which is in contrast to the experiment.

A druggable pocket on PSMD10Gankyrin that can accommodate an interface peptide and doxorubicin

BACKGROUND

PSMD10Gankyrin, a proteasomal chaperone is also an oncoprotein. Overexpression of PSMD10Gankyrin is associated with poor prognosis and survival in many cancers. Therefore, PSMD10Gankyrin is a sought-after drug target in many hard-to-treat cancers. However, its surface appears flat and undruggable. Here, we build on our earlier discovery of a common hot spot region that defined the interface of multiple interacting partners of PSMD10Gankyrin to expose vulnerable spots for a peptide and a small molecule inhibitor.

METHODS

High throughput virtual screening was used to screen compounds against PSMD10Gankyrin. Interaction of PSMD10Gankyrin with the drug or protein (CLIC1) or peptide was studied using any one or more of these techniques; Microscale Thermophoresis, limited trypsinolysis, SPR and ITC. Cytotoxic effect of doxorubicin was evaluated using MTT assay.

RESULTS

We identified doxorubicin as the first-generation small molecule inhibitor of PSMD10Gankyrin. K116 and to a lesser extent R41 on PSMD10Gankyrin contribute to the bulk of binding energy for the peptide EEVD, CLIC1 and doxorubicin. We further demonstrate that PSMD10Gankyrin is an intended target for doxorubicin in cells.

GENERAL SIGNIFICANCE

Drug design against protein interactions in general and PSMD10Gankyrin in particular, remains a challenge. We provide consolidated biophysical evidence for the use of a shared interface motif EEVD as a possible inhibitor of interaction network in cancers driven by PSMD10Gankyrin. We identify a chemical scaffold for designing novel inhibitors of PSMD10Gankyrin. These findings will impact the field of protein interactions in the context of disease biology/drug discovery.

Recent Developments in Free Energy Calculations for Drug Discovery

The grand challenge in structure-based drug design is achieving accurate prediction of binding free energies. Molecular dynamics (MD) simulations enable modeling of conformational changes critical to the binding process, leading to calculation of thermodynamic quantities involved in estimation of binding affinities. With recent advancements in computing capability and predictive accuracy, MD based virtual screening has progressed from the domain of theoretical attempts to real application in drug development. Approaches including the Molecular Mechanics Poisson Boltzmann Surface Area (MM-PBSA), Linear Interaction Energy (LIE), and alchemical methods have been broadly applied to model molecular recognition for drug discovery and lead optimization. Here we review the varied methodology of these approaches, developments enhancing simulation efficiency and reliability, remaining challenges hindering predictive performance, and applications to problems in the fields of medicine and biochemistry.

PSMD9 ribosomal protein network maintains nucleolar architecture and WT p53 levels

Capitalizing on an unexpected observation that multiple free ribosomal proteins co-purify/pull-down with PSMD9, we report here for the first time that PSMD9 is necessary to maintain the morphology and integrity of the nucleolus. As seen by NPM1 immunofluorescence and electron microscopy, the nucleolar structure is clearly disrupted in PSMD9 null MCF7 breast cancer cells. The resultant stress is pronounced leading to the accumulation of WT p53 and slow growth. A dual insult with Actinomycin D exasperates the nucleolar stress in these cells which fail to recover in stipulated time. This double insult in the WT cells enhances the interaction of PSMD9 with ribosomal subunits. Our data also reveals that in PSMD9 null cells, ribosomal proteins RPS25 and RPL15 fail to localise in the nucleolus. We speculate that the interaction of PSMD9 with multiple free ribosome subunits has at least two important implications: a) PSMD9 plays a role in trafficking of ribosomal proteins into the nucleolus, therefore contributing to the maintenance of structural and morphological organization of the membrane-less nucleolar compartment; b) under conditions that induce nucleolar stress, PSMD9-Ribosomal Protein interaction protects WT MCF7 breast cancer cells from slow growth and eventual death. This possibility renders the domains of PSMD9 to be attractive drug targets in the context of cancer and other multiple ribosome-associated disorders.

Evolution of biophysical tools for quantitative protein interactions and drug discovery

With millions of signalling events occurring simultaneously, cells process a continuous flux of information. The genesis, processing, and regulation of information are dictated by a huge network of protein interactions. This is proven by the fact that alterations in the levels of proteins, single amino acid changes, post-translational modifications, protein products arising out of gene fusions alter the interaction landscape leading to diseases such as congenital disorders, deleterious syndromes like cancer, and crippling diseases like the neurodegenerative disorders which are often fatal. Needless to say, there is an immense effort to understand the biophysical basis of such direct interactions between any two proteins, the structure, domains, and sequence motifs involved in tethering them, their spatio-temporal regulation in cells, the structure of the network, and their eventual manipulation for intervention in diseases. In this chapter, we will deliberate on a few techniques that allow us to dissect the thermodynamic and kinetic aspects of protein interaction, how innovation has rendered some of the traditional techniques applicable for rapid analysis of multiple samples using small amounts of material. These advances coupled with automation are catching up with the genome-wide or proteome-wide studies aimed at identifying new therapeutic targets. The chapter will also summarize how some of these techniques are suited either in the standalone mode or in combination with other biophysical techniques for the drug discovery process.

Straightforward Incorporation of Multiple Ligand Types into Molecular Dynamics Simulations for Efficient Binding Site Detection and Characterization

Binding site identification and characterization is an important initial step in structure-based drug design. To account for the effects of protein flexibility and solvation, several cosolvent molecular dynamics (MD) simulation methods that incorporate small organic molecules into the protein's solvent box to probe for binding sites have been developed. However, most of these methods are highly inefficient, as they allow for the use of only one probe type at a time, which means that multiple sets of simulations have to be performed to map different types of binding sites. The high probe concentrations used in some of these methods also necessitate the use of artificial repulsive forces to prevent the probes from aggregating. Here, we present multiple-ligand-mapping MD (mLMMD), a method that incorporates multiple types of probes for simultaneous and efficient mapping of different types of binding sites without the need for introduction of artificial forces that may cause unintended mapping artifacts. We validate the method on a diverse set of 10 proteins and show that the mLMMD probes are able to reliably identify hydrophobic, hydrogen-bonding, charged, and cryptic binding sites in all of the test cases. Our results also highlight the potential utility of mLMMD for virtual screening and rational drug design.

PSMD9 expression correlates with recurrence after radiotherapy in patients with cervical cancer

In the current retrospective cohort study, the expression of the Proteasome 26S non-ATPase Subunit 9 (PSMD9) was investigated in 102 patients with cervical cancer. The rat homologue of PSMD9, Bridge-1, was identified as a binding protein of the transcription factors PDX-1 and E-12 via its PDZ-domain. The aim of the current study was to evaluate the prognostic or predictive value of PSMD9 expression as a biomarker for patients with cervical cancer. Tissue microarrays were constructed from formalin-fixed paraffin-embedded tissue specimens of cervical cancer and peritumoral stroma after hysterectomy and a Bridge-1 antibody was used to perform immunohistochemistry. The immunoreactions were analyzed using an immunoreactive score, which evaluated the number of positive cells as well as their intensity of PSMD9 expression. A misinterpretation of statistically significant results after multiple testing was controlled by the false discovery rate correction using the algorithm of Benjamini and Hochberg. All tumor tissues and almost all peritumoral stroma tissues expressed PSMD9. The PSMD9 expression in tumor tissues was significantly higher compared with the peritumoral stroma. PSMD9 expression correlated significantly with the expression of the proliferation marker MIB-1. Patients with stronger PSMD9 expression tended to exhibit a higher odds ratio for the recurrence of the disease in all patients (n=102) as well as in the subgroup of 47 patients having received a combined chemoradiotherapy following hysterectomy. In the group of 62 patients having that received radiotherapy following hysterectomy, which included the chemoradiotherapy patients, a higher PSMD9 expression significantly increased the odds for a recurrence to 1.983-fold even after FDR correction (P=0.0304). In conclusion, PSMD9 was indicated to be overexpressed in tumor tissues and associated with tumor cell proliferation. Therefore, PSMD9 may be useful as a tumor marker. Furthermore, increased PSMD9 overexpression may be used to predict resistance against radiation.