The structural flexibility of the shank1 PDZ domain is important for its binding to different ligands

Biophysical and structural analyses of the interaction between the SHANK1 PDZ domain and an internal SLiM

The PDZ (Postsynaptic density protein-95[PSD-95]/Discs-large) domain, prevalent as a recognition module, has attracted significant attention given its ability to specifically recognize ligands with consensus motifs (also termed PDZ binding motifs [PBMs]). PBMs typically bear a C-terminal carboxylate as a recognition handle and have been extensively characterized, whilst internal ligands are less well known. Here we characterize a short linear motif (SLiM) - EESTSFQGP - as an internal PBM based on its strong binding affinity towards the SHANK1 PDZ domain (SHANK1656-762 hereafter referred to as SHANK1). Using the acetylated analogue Ac-EESTSFQGP-CONH2 as a competitor for the interaction of SHANK1 with FAM-Ahx-EESTSFQGP-CONH2 or a typical fluorophore-labelled C-terminal PBM - GKAP - FITC-Ahx-EAQTRL-COOH - the internal SLiM was demonstrated to show comparable low-micromolar IC50 by competition fluorescent anisotropy. To gain further insight into the internal ligand interaction at the molecular level, we obtained the X-ray co-crystal structure of the Ac-EESTSFQGP-CONH2/SHANK1 complex and compared this to the Ac-EAQTRL-COOH/SHANK1 complex. The crystallographic studies reveal that the SHANK1 backbones for the two interactions overlap significantly. The main structural differences were shown to result from the flexible loops which reorganize to accommodate the two PBMs with distinct lengths and terminal groups. In addition, the two C-terminal residues Gly and Pro in Ac-EESTSFQGP-CONH2 were shown not to participate in interaction with the target protein, implying further truncation and structural modification using peptidomimetic approaches on this sequence may be feasible. Taken together, the SLiM Ac-EESTSFQGP-CONH2 holds potential as an internal ligand for targeting SHANK1.

Simulated complexes formed from a set of postsynaptic proteins suggest a localised effect of a hypomorphic Shank mutation

BACKGROUND

The postsynaptic density is an elaborate protein network beneath the postsynaptic membrane involved in the molecular processes underlying learning and memory. The postsynaptic density is built up from the same major proteins but its exact composition and organization differs between synapses. Mutations perturbing protein: protein interactions generally occurring in this network might lead to effects specific for cell types or processes, the understanding of which can be especially challenging.

RESULTS

In this work we use systems biology-based modeling of protein complex distributions in a simplified set of major postsynaptic proteins to investigate the effect of a hypomorphic Shank mutation perturbing a single well-defined interaction. We use data sets with widely variable abundances of the constituent proteins. Our results suggest that the effect of the mutation is heavily dependent on the overall availability of all the protein components of the whole network and no trivial correspondence between the expression level of the directly affected proteins and overall complex distribution can be observed.

CONCLUSIONS

Our results stress the importance of context-dependent interpretation of mutations. Even the weakening of a generally occurring protein: protein interaction might have well-defined effects, and these can not easily be predicted based only on the abundance of the proteins directly affected. Our results provide insight on how cell-specific effects can be exerted by a mutation perturbing a generally occurring interaction even when the wider interaction network is largely similar.

Shank promotes action potential repolarization by recruiting BK channels to calcium microdomains

Mutations altering the scaffolding protein Shank are linked to several psychiatric disorders, and to synaptic and behavioral defects in mice. Among its many binding partners, Shank directly binds CaV1 voltage activated calcium channels. Here we show that the C. elegans SHN-1/Shank promotes CaV1 coupling to calcium activated potassium channels. Mutations inactivating SHN-1, and those preventing SHN-1 binding to EGL-19/CaV1 all increase action potential durations in body muscles. Action potential repolarization is mediated by two classes of potassium channels: SHK-1/KCNA and SLO-1 and SLO-2 BK channels. BK channels are calcium-dependent, and their activation requires tight coupling to EGL-19/CaV1 channels. SHN-1's effects on AP duration are mediated by changes in BK channels. In shn-1 mutants, SLO-2 currents and channel clustering are significantly decreased in both body muscles and neurons. Finally, increased and decreased shn-1 gene copy number produce similar changes in AP width and SLO-2 current. Collectively, these results suggest that an important function of Shank is to promote microdomain coupling of BK with CaV1.

Query-guided protein-protein interaction inhibitor discovery

Protein-protein interactions (PPIs) are central to biological mechanisms, and can serve as compelling targets for drug discovery. Yet, the discovery of small molecule inhibitors of PPIs remains challenging given the large and typically shallow topography of the interacting protein surfaces. Here, we describe a general approach to the discovery of orthosteric PPI inhibitors that mimic specific secondary protein structures. Initially, hot residues at protein-protein interfaces are identified in silico or from experimental data, and incorporated into secondary structure-based queries. Virtual libraries of small molecules are then shape-matched against the queries, and promising ligands docked to target proteins. The approach is exemplified experimentally using two unrelated PPIs that are mediated by an α-helix (p53/hDM2) and a β-strand (GKAP/SHANK1-PDZ). In each case, selective PPI inhibitors are discovered with low μM activity as determined by a combination of fluorescence anisotropy and 1H-15N HSQC experiments. In addition, hit expansion yields a series of PPI inhibitors with defined structure-activity relationships. It is envisaged that the generality of the approach will enable discovery of inhibitors of a wide range of unrelated secondary structure-mediated PPIs.

Identification of β-strand mediated protein-protein interaction inhibitors using ligand-directed fragment ligation

β-Strand mediated protein-protein interactions (PPIs) represent underexploited targets for chemical probe development despite representing a significant proportion of known and therapeutically relevant PPI targets. β-Strand mimicry is challenging given that both amino acid side-chains and backbone hydrogen-bonds are typically required for molecular recognition, yet these are oriented along perpendicular vectors. This paper describes an alternative approach, using GKAP/SHANK1 PDZ as a model and dynamic ligation screening to identify small-molecule replacements for tranches of peptide sequence. A peptide truncation of GKAP functionalized at the N- and C-termini with acylhydrazone groups was used as an anchor. Reversible acylhydrazone bond exchange with a library of aldehyde fragments in the presence of the protein as template and in situ screening using a fluorescence anisotropy (FA) assay identified peptide hybrid hits with comparable affinity to the GKAP peptide binding sequence. Identified hits were validated using FA, ITC, NMR and X-ray crystallography to confirm selective inhibition of the target PDZ-mediated PPI and mode of binding. These analyses together with molecular dynamics simulations demonstrated the ligands make transient interactions with an unoccupied basic patch through electrostatic interactions, establishing proof-of-concept that this unbiased approach to ligand discovery represents a powerful addition to the armory of tools that can be used to identify PPI modulators.

Integrated analysis of Shank1 PDZ interactions with C-terminal and internal binding motifs

PDZ domains constitute a large family of modular domains that are well-known for binding C-terminal motifs of target proteins. Some of them also bind to internal PDZ binding motifs (PDZbms), but this aspect of the PDZ interactome is poorly studied. Here we explored internal PDZbm-mediated interactions using the PDZ domain of Shank1 as a model. We identified a series of human Shank1 ligands with C-terminal or internal PDZbms using proteomic peptide-phage display, and established that while the consensus sequence of C-terminal ligands is x-T-x-(L/F)-COOH, the consensus of internal PDZbm is exclusively x-T-x-F-x, where x is any amino acid. We found that the affinities of PDZbm interactions are in the low micromolar range. The crystal structure of the complex between Shank1 PDZ and an internal PDZbm revealed that the binding mode of internal PDZbms was similar to that of C-terminal ligands. Pull-down experiments confirmed that both C-terminal and internal PDZbm interactions can occur in the context of full-length proteins. Our study expands the interactome of Shank1 and hints at a largely unexplored interaction space of PDZ domains.

Structural characterization and computational analysis of PDZ domains in Monosiga brevicollis

Identification of the molecular networks that facilitated the evolution of multicellular animals from their unicellular ancestors is a fundamental problem in evolutionary cellular biology. Choanoflagellates are recognized as the closest extant nonmetazoan ancestors to animals. These unicellular eukaryotes can adopt a multicellular-like "rosette" state. Therefore, they are compelling models for the study of early multicellularity. Comparative studies revealed that a number of putative human orthologs are present in choanoflagellate genomes, suggesting that a subset of these genes were necessary for the emergence of multicellularity. However, previous work is largely based on sequence alignments alone, which does not confirm structural nor functional similarity. Here, we focus on the PDZ domain, a peptide-binding domain which plays critical roles in myriad cellular signaling networks and which underwent a gene family expansion in metazoan lineages. Using a customized sequence similarity search algorithm, we identified 178 PDZ domains in the Monosiga brevicollis proteome. This includes 11 previously unidentified sequences, which we analyzed using Rosetta and homology modeling. To assess conservation of protein structure, we solved high-resolution crystal structures of representative M. brevicollis PDZ domains that are homologous to human Dlg1 PDZ2, Dlg1 PDZ3, GIPC, and SHANK1 PDZ domains. To assess functional conservation, we calculated binding affinities for mbGIPC, mbSHANK1, mbSNX27, and mbDLG-3 PDZ domains from M. brevicollis. Overall, we find that peptide selectivity is generally conserved between these two disparate organisms, with one possible exception, mbDLG-3. Overall, our results provide novel insight into signaling pathways in a choanoflagellate model of primitive multicellularity.

Computational Analysis of Energy Landscapes Reveals Dynamic Features That Contribute to Binding of Inhibitors to CFTR-Associated Ligand

The CFTR-associated ligand PDZ domain (CALP) binds to the cystic fibrosis transmembrane conductance regulator (CFTR) and mediates lysosomal degradation of mature CFTR. Inhibition of this interaction has been explored as a therapeutic avenue for cystic fibrosis. Previously, we reported the ensemble-based computational design of a novel peptide inhibitor of CALP, which resulted in the most binding-efficient inhibitor to date. This inhibitor, kCAL01, was designed using osprey and evinced significant biological activity in in vitro cell-based assays. Here, we report a crystal structure of kCAL01 bound to CALP and compare structural features against iCAL36, a previously developed inhibitor of CALP. We compute side-chain energy landscapes for each structure to not only enable approximation of binding thermodynamics but also reveal ensemble features that contribute to the comparatively efficient binding of kCAL01. Finally, we compare the previously reported design ensemble for kCAL01 vs the new crystal structure and show that, despite small differences between the design model and crystal structure, significant biophysical features that enhance inhibitor binding are captured in the design ensemble. This suggests not only that ensemble-based design captured thermodynamically significant features observed in vitro, but also that a design eschewing ensembles would miss the kCAL01 sequence entirely.

Molecular Dynamics Simulations of Wild Type and Mutants of SAPAP in Complexed with Shank3

Specific interactions between scaffold protein SH3 and multiple ankyrin repeat domains protein 3 (Shank3) and synapse-associated protein 90/postsynaptic density-95–associated protein (SAPAP) are essential for excitatory synapse development and plasticity. In a bunch of human neurological diseases, mutations on Shank3 or SAPAP are detected. To investigate the dynamical and thermodynamic properties of the specific binding between the N-terminal extended PDZ (Post-synaptic density-95/Discs large/Zonaoccludens-1) domain (N-PDZ) of Shank3 and the extended PDZ binding motif (E-PBM) of SAPAP, molecular dynamics simulation approaches were used to study the complex of N-PDZ with wild type and mutated E-PBM peptides. To compare with the experimental data, ₉₇₄QTRL₉₇₇ and ₉₆₆IEIYI₉₇₀ of E-PBM peptide were mutated to prolines to obtain the M4P and M5P system, respectively. Conformational analysis shows that the canonical PDZ domain is stable while the βN extension presents high flexibility in all systems, especially for M5P. The high flexibility of βN extension seems to set up a barrier for the non-specific binding in this area and provide the basis for specific molecular recognition between Shank3 and SAPAP. The wild type E-PBM tightly binds to N-PDZ during the simulation while loss of binding is observed in different segments of the mutated E-PBM peptides. Energy decomposition and hydrogen bonds analysis show that M4P mutations only disrupt the interactions with canonical PDZ domain, but the interactions with βN1′ remain. In M5P system, although the interactions with βN1′ are abolished, the binding between peptide and the canonical PDZ domain is not affected. The results indicate that the interactions in the two-binding site, the canonical PDZ domain and the βN1′ extension, contribute to the binding between E-PBM and N-PDZ independently. The binding free energies calculated by MM/GBSA (Molecular Mechanics/Generalized Born Surface Area) are in agreement with the experimental binding affinities. Most of the residues on E-PBM contribute considerably favorable energies to the binding except A963 and D964 in the N-terminal. The study provides information to understand the molecular basis of specific binding between Shank3 and SAPAP, as well as clues for design of peptide inhibitors.

Shank is a dose-dependent regulator of Cav1 calcium current and CREB target expression

Shank is a post-synaptic scaffolding protein that has many binding partners. Shank mutations and copy number variations (CNVs) are linked to several psychiatric disorders, and to synaptic and behavioral defects in mice. It is not known which Shank binding partners are responsible for these defects. Here we show that the C. elegans SHN-1/Shank binds L-type calcium channels and that increased and decreased shn-1 gene dosage alter L-channel current and activity-induced expression of a CRH-1/CREB transcriptional target (gem-4 Copine), which parallels the effects of human Shank copy number variations (CNVs) on Autism spectrum disorders and schizophrenia. These results suggest that an important function of Shank proteins is to regulate L-channel current and activity induced gene expression.

DOI:http://dx.doi.org/10.7554/eLife.18931.001

Discovery of short linear motif-mediated interactions through phage display of intrinsically disordered regions of the human proteome

The intrinsically disordered regions of eukaryotic proteomes are enriched in short linear motifs (SLiMs), which are of crucial relevance for cellular signaling and protein regulation; many mediate interactions by providing binding sites for peptide-binding domains. The vast majority of SLiMs remain to be discovered highlighting the need for experimental methods for their large-scale identification. We present a novel proteomic peptide phage display (ProP-PD) library that displays peptides representing the disordered regions of the human proteome, allowing direct large-scale interrogation of most potential binding SLiMs in the proteome. The performance of the ProP-PD library was validated through selections against SLiM-binding bait domains with distinct folds and binding preferences. The vast majority of identified binding peptides contained sequences that matched the known SLiM-binding specificities of the bait proteins. For SHANK1 PDZ, we establish a novel consensus TxF motif for its non-C-terminal ligands. The binding peptides mostly represented novel target proteins, however, several previously validated protein-protein interactions (PPIs) were also discovered. We determined the affinities between the VHS domain of GGA1 and three identified ligands to 40-130 μm through isothermal titration calorimetry, and confirmed interactions through coimmunoprecipitation using full-length proteins. Taken together, we outline a general pipeline for the design and construction of ProP-PD libraries and the analysis of ProP-PD-derived, SLiM-based PPIs. We demonstrated the methods potential to identify low affinity motif-mediated interactions for modular domains with distinct binding preferences. The approach is a highly useful complement to the current toolbox of methods for PPI discovery.

Crystal structure of the PDZ domain of mouse Dishevelled 1 and its interaction with CXXC5

Dishevelled (Dvl) plays a crucial role in Wnt signaling by interacting with membrane-bound receptors and downstream molecules through its PDZ domain. CXXC5 is one of the key molecules that interacts with Dvl and negatively regulates the Wnt/β-catenin pathway in osteoblast differentiation. Recently, the Dvl-CXXC5 interaction has been identified as an excellent target for osteoporosis treatment. Therefore, it is desirable to have detailed structural information for the Dvl-CXXC5 interaction. Although solution structures of the Dvl1 PDZ domain have been reported, a high-resolution crystal structure would provide detailed sidechain information that is essential for drug development. Here, we determined the first crystal structure of the Dvl-1 PDZ domain at a resolution of 1.76 Å, and compared it with its previously reported solution structure. The Dvl1 PDZ domain crystal belonged to the space group H32 with unit-cell parameters a = b = 72.837, c = 120.616, α = β = 90.00, γ = 120.00. The crystal structure of Dvl1 PDZ shared its topology with the previously reported structure determined by nuclear magnetic resonance (NMR); however, the crystal structure was quite different from the solution structure in both the secondary structural region and the ligand-binding pocket. Molecular modeling based on NMR and X-ray crystallographic data yielded detailed information about the Dvl1/CXXC5 interaction, which will be useful for designing inhibitors.

Structure and function of POTRA domains of Omp85/TPS superfamily

The Omp85/TPS (outer-membrane protein of 85 kDa/two-partner secretion) superfamily is a ubiquitous and major class of β-barrel proteins. This superfamily is restricted to the outer membranes of gram-negative bacteria, mitochondria, and chloroplasts. The common architecture, with an N-terminus consisting of repeats of soluble polypeptide-transport-associated (POTRA) domains and a C-terminal β-barrel pore is highly conserved. The structures of multiple POTRA domains and one full-length TPS protein have been solved, yet discovering roles of individual POTRA domains has been difficult. This review focuses on similarities and differences between POTRA structures, emphasizing POTRA domains in autotrophic organisms including plants and cyanobacteria. Unique roles, specific for certain POTRA domains, are examined in the context of POTRA location with respect to their attachment to the β-barrel pore, and their degree of biological dispensability. Finally, because many POTRA domains may have the ability to interact with thousands of partner proteins, possible modes of these interactions are also explored.

Plasticity of PDZ domains in ligand recognition and signaling

The PDZ domain is a protein-protein interacting module that plays an important role in the organization of signaling complexes. The recognition of short intrinsically disordered C-terminal peptide motifs is the archetypical PDZ function, but the functional repertoire of this versatile module also includes recognition of internal peptide sequences, dimerization and phospholipid binding. The PDZ function can be tuned by various means such as allosteric effects, changes of physiological buffer conditions and phosphorylation of PDZ domains and/or ligands, which poses PDZ domains as dynamic regulators of cell signaling. This review is focused on the plasticity of the PDZ interactions.