Origins of PDZ domain ligand specificity. Structure determination and mutagenesis of the Erbin PDZ domain

Dishevelled2 activates WGEF via its interaction with a unique internal peptide motif of the GEF

The Wnt-planar cell polarity (Wnt-PCP) pathway is crucial in establishing cell polarity during development and tissue homoeostasis. This pathway is found to be dysregulated in many pathological conditions, including cancer and autoimmune disorders. The central event in Wnt-PCP pathway is the activation of Weak-similarity guanine nucleotide exchange factor (WGEF) by the adapter protein Dishevelled (Dvl). The PDZ domain of Dishevelled2 (Dvl2PDZ) binds and activates WGEF by releasing it from its autoinhibitory state. However, the actual Dvl2PDZ binding site of WGEF and the consequent activation mechanism of the GEF have remained elusive. Using biochemical and molecular dynamics studies, we show that a unique "internal-PDZ binding motif" (IPM) of WGEF mediates the WGEF-Dvl2PDZ interaction to activate the GEF. The residues at P2, P0, P-2 and P-3 positions of IPM play an important role in stabilizing the WGEFpep-Dvl2PDZ interaction. Furthermore, MD simulations of modelled Dvl2PDZ-WGEFIPM peptide complexes suggest that WGEF-Dvl2PDZ interaction may differ from the reported Dvl2PDZ-IPM interactions. Additionally, the apo structure of human Dvl2PDZ shows conformational dynamics different from its IPM peptide bound state, suggesting an induced fit mechanism for the Dvl2PDZ-peptide interaction. The current study provides a model for Dvl2 induced activation of WGEF.

A complete map of specificity encoding for a partially fuzzy protein interaction

Thousands of human proteins function by binding short linear motifs embedded in intrinsically disordered regions. How affinity and specificity are encoded in these binding domains and the motifs themselves is not well understood. The evolvability of binding specificity - how rapidly and extensively it can change upon mutation - is also largely unexplored, as is the contribution of 'fuzzy' dynamic residues to affinity and specificity in protein-protein interactions. Here we report the first complete map of specificity encoding for a globular protein domain. Quantifying >200,000 energetic interactions between a PDZ domain and its ligand identifies 20 major energetically coupled pairs of sites that control specificity. These are organized into six modules, with most mutations in each module reprogramming specificity for a single position in the ligand. Nine of the major energetic couplings controlling specificity are between structural contacts and 11 have an allosteric mechanism of action. The dynamic tail of the ligand is more robust to mutation than the structured residues but contributes additively to binding affinity and communicates with structured residues to enable changes in specificity. Our results quantify the binding specificities of >1,800 globular proteins to reveal how specificity is encoded and provide a direct comparison of the encoding of affinity and specificity in structured and dynamic molecular recognition.

Disulfide-constrained peptide scaffolds enable a robust peptide-therapeutic discovery platform

Peptides present an alternative modality to immunoglobulin domains or small molecules for developing therapeutics to either agonize or antagonize cellular pathways associated with diseases. However, peptides often suffer from poor chemical and physical stability, limiting their therapeutic potential. Disulfide-constrained peptides (DCP) are naturally occurring and possess numerous desirable properties, such as high stability, that qualify them as drug-like scaffolds for peptide therapeutics. DCPs contain loop regions protruding from the core of the molecule that are amenable to peptide engineering via direct evolution by use of phage display technology. In this study, we have established a robust platform for the discovery of peptide therapeutics using various DCPs as scaffolds. We created diverse libraries comprising seven different DCP scaffolds, resulting in an overall diversity of 2 x 1011. The effectiveness of this platform for functional hit discovery has been extensively evaluated, demonstrating a hit rate comparable to that of synthetic antibody libraries. By utilizing chemically synthesized and in vitro folded peptides derived from selections of phage displayed DCP libraries, we have successfully generated functional inhibitors targeting the HtrA1 protease. Through affinity maturation strategies, we have transformed initially weak binders against Notch2 with micromolar Kd values to high-affinity ligands in the nanomolar range. This process highlights a viable hit-to-lead progression. Overall, our platform holds significant potential to greatly enhance the discovery of peptide therapeutics.

The binding of Mint/X11 PDZ domains to Ca V 2 calcium channels predates bilaterian animals

PDZ domain mediated interactions with voltage-gated calcium (Ca V ) channel C-termini play important roles in localizing and compartmentalizing membrane Ca 2+ signaling. The first such interaction discovered was between the neuronal multi-domain protein Mint-1, and the presynaptc calcium channel Ca V 2.2 in mammals. Although the physiological significance of this interaction is unclear, its occurrence in vertebrates and bilaterian invertebrates suggests important and conserved functions. In this study, we explore the evolutionary origins of Mint and its interaction with Ca V 2 channels. Phylogenetic and structural in silico analyses revealed that Mint is an animal-specific gene, like Ca V 2 channels, which bears a highly divergent N-terminus but strongly conserved C-terminus comprised of a phosphotyrosine binding domain, two tandem PDZ domains (PDZ-1 and PDZ-2), and a C-terminal auto-inhibitory element that binds and inhibits PDZ-1. Also deeply conserved are other Mint interacting proteins, namely amyloid precursor and related proteins, presenilins, neurexin, as well as CASK and Veli which form a tripartite complex with Mint in bilaterians. Through yeast 2-hybrid and bacterial 2-hybrid experiments, we show that Mint and Ca V 2 channels from cnidarians and placozoans interact in vitro , and in situ hybridization revealed co-expression of corresponding transcripts in dissociated neurons from the cnidarian Nematostella vectensis . Unexpectedly, the Mint orthologue from the ctenophore Hormiphora californiensis was able to strongly bind the divergent C-terminal ligands of cnidarian and placozoan Ca V 2 channels, despite neither the ctenophore Mint, nor the placozoan and cnidarian orthologues, binding the ctenophore Ca V 2 channel C-terminus. Altogether, our analyses provide a model for the emergence of this interaction in early animals first via adoption of a PDZ ligand by Ca V 2 channels, followed by sequence changes in the ligand that caused a modality switch for binding to Mint.

DexDesign: A new OSPREY-based algorithm for designing de novo D-peptide inhibitors

With over 270 unique occurrences in the human genome, peptide-recognizing PDZ domains play a central role in modulating polarization, signaling, and trafficking pathways. Mutations in PDZ domains lead to diseases such as cancer and cystic fibrosis, making PDZ domains attractive targets for therapeutic intervention. D-peptide inhibitors offer unique advantages as therapeutics, including increased metabolic stability and low immunogenicity. Here, we introduce DexDesign, a novel OSPREY-based algorithm for computationally designing de novo D-peptide inhibitors. DexDesign leverages three novel techniques that are broadly applicable to computational protein design: the Minimum Flexible Set, K*-based Mutational Scan, and Inverse Alanine Scan, which enable exponential reductions in the size of the peptide sequence search space. We apply these techniques and DexDesign to generate novel D-peptide inhibitors of two biomedically important PDZ domain targets: CAL and MAST2. We introduce a new framework for analyzing de novo peptides-evaluation along a replication/restitution axis-and apply it to the DexDesign-generated D-peptides. Notably, the peptides we generated are predicted to bind their targets tighter than their targets' endogenous ligands, validating the peptides' potential as lead therapeutic candidates. We provide an implementation of DexDesign in the free and open source computational protein design software OSPREY.

DexDesign: an OSPREY-based algorithm for designing de novo D-peptide inhibitors

With over 270 unique occurrences in the human genome, peptide-recognizing PDZ domains play a central role in modulating polarization, signaling, and trafficking pathways. Mutations in PDZ domains lead to diseases such as cancer and cystic fibrosis, making PDZ domains attractive targets for therapeutic intervention. D-peptide inhibitors offer unique advantages as therapeutics, including increased metabolic stability and low immunogenicity. Here, we introduce DexDesign, a novel OSPREY-based algorithm for computationally designing de novo D-peptide inhibitors. DexDesign leverages three novel techniques that are broadly applicable to computational protein design: the Minimum Flexible Set, K*-based Mutational Scan, and Inverse Alanine Scan. We apply these techniques and DexDesign to generate novel D-peptide inhibitors of two biomedically important PDZ domain targets: CAL and MAST2. We introduce a framework for analyzing de novo peptides-evaluation along a replication/restitution axis-and apply it to the DexDesign-generated D-peptides. Notably, the peptides we generated are predicted to bind their targets tighter than their targets' endogenous ligands, validating the peptides' potential as lead inhibitors. We also provide an implementation of DexDesign in the free and open source computational protein design software OSPREY.

Development of a Potent Cyclic Peptide Inhibitor of the nNOS/PSD-95 Interaction

The complex between the N-methyl-d-aspartate receptor (NMDAR), neuronal nitric oxide synthase (nNOS), and the postsynaptic density protein-95 (PSD-95) is an attractive therapeutic target for the treatment of acute ischemic stroke. The complex is formed via the PDZ protein domains of PSD-95, and efforts to disrupt the complex have generally been based on C-terminal peptides derived from the NMDAR. However, nNOS binds PSD-95 through a β-hairpin motif, providing an alternative starting point for developing PSD-95 inhibitors. Here, we designed a cyclic nNOS β-hairpin mimetic peptide and generated cyclic nNOS β-hairpin peptide arrays with natural and unnatural amino acids (AAs), which provided molecular insights into this interaction. We then optimized cyclic peptides and identified a potent inhibitor of the nNOS/PSD-95 interaction, with the highest affinity reported thus far for a peptide macrocycle inhibitor of PDZ domains, which serves as a template for the development of treatment for acute ischemic stroke.

The Repeating, Modular Architecture of the HtrA Proteases

A conserved, 26-residue sequence [AA(X₂)[A/G][G/L](X₂)GDV[I/L](X₂)[V/L]NGE(X₁)V(X₆)] and corresponding structure repeating module were identified within the HtrA protease family using a non-redundant set (N = 20) of publicly available structures. While the repeats themselves were far from sequence perfect, they had notable conservation to a statistically significant level. Three or more repetitions were identified within each protein despite being statistically expected to randomly occur only once per 1031 residues. This sequence repeat was associated with a six stranded antiparallel β-barrel module, two of which are present in the core of the structures of the PA clan of serine proteases, while a modified version of this module could be identified in the PDZ-like domains. Automated structural alignment methods had difficulties in superimposing these β-barrels, but the use of a target human HtrA2 structure showed that these modules had an average RMSD across the set of structures of less than 2 Å (mean and median). Our findings support Dayhoff’s hypothesis that complex proteins arose through duplication of simpler peptide motifs and domains.

Domain Analysis and Motif Matcher (DAMM): A Program to Predict Selectivity Determinants in Monosiga brevicollis PDZ Domains Using Human PDZ Data

Choanoflagellates are single-celled eukaryotes with complex signaling pathways. They are considered the closest non-metazoan ancestors to mammals and other metazoans and form multicellular-like states called rosettes. The choanoflagellate Monosiga brevicollis contains over 150 PDZ domains, an important peptide-binding domain in all three domains of life (Archaea, Bacteria, and Eukarya). Therefore, an understanding of PDZ domain signaling pathways in choanoflagellates may provide insight into the origins of multicellularity. PDZ domains recognize the C-terminus of target proteins and regulate signaling and trafficking pathways, as well as cellular adhesion. Here, we developed a computational software suite, Domain Analysis and Motif Matcher (DAMM), that analyzes peptide-binding cleft sequence identity as compared with human PDZ domains and that can be used in combination with literature searches of known human PDZ-interacting sequences to predict target specificity in choanoflagellate PDZ domains. We used this program, protein biochemistry, fluorescence polarization, and structural analyses to characterize the specificity of A9UPE9_MONBE, a M. brevicollis PDZ domain-containing protein with no homology to any metazoan protein, finding that its PDZ domain is most similar to those of the DLG family. We then identified two endogenous sequences that bind A9UPE9 PDZ with <100 μM affinity, a value commonly considered the threshold for cellular PDZ-peptide interactions. Taken together, this approach can be used to predict cellular targets of previously uncharacterized PDZ domains in choanoflagellates and other organisms. Our data contribute to investigations into choanoflagellate signaling and how it informs metazoan evolution.

A Thermodynamic Analysis of the Binding Specificity between Four Human PDZ Domains and Eight Host, Viral and Designed Ligands

PDZ domains are binding modules mostly involved in cell signaling and cell–cell junctions. These domains are able to recognize a wide variety of natural targets and, among the PDZ partners, viruses have been discovered to interact with their host via a PDZ domain. With such an array of relevant and diverse interactions, PDZ binding specificity has been thoroughly studied and a traditional classification has grouped PDZ domains in three major specificity classes. In this work, we have selected four human PDZ domains covering the three canonical specificity-class binding mode and a set of their corresponding binders, including host/natural, viral and designed PDZ motifs. Through calorimetric techniques, we have covered the entire cross interactions between the selected PDZ domains and partners. The results indicate a rather basic specificity in each PDZ domain, with two of the domains that bind their cognate and some non-cognate ligands and the two other domains that basically bind their cognate partners. On the other hand, the host partners mostly bind their corresponding PDZ domain and, interestingly, the viral ligands are able to bind most of the studied PDZ domains, even those not previously described. Some viruses may have evolved to use of the ability of the PDZ fold to bind multiple targets, with resulting affinities for the virus–host interactions that are, in some cases, higher than for host–host interactions.

Comprehensive Assessment of the Relationship Between Site-2 Specificity and Helix α2 in the Erbin PDZ Domain

PDZ domains are key players in signalling pathways. These modular domains generally recognize short linear C-terminal stretches of sequences in proteins that organize the formation of complex multi-component assemblies. The development of new methodologies for the characterization of the molecular principles governing these interactions is critical to fully understand the functional diversity of the family and to elucidate biological functions for family members. Here, we applied an in vitro evolution strategy to explore comprehensively the capacity of PDZ domains for specific recognition of different amino acids at a key position in C-terminal peptide ligands. We constructed a phage-displayed library of the Erbin PDZ domain by randomizing the binding site^-2 and adjacent residues, which are all contained in helix α2, and we selected for variants binding to a panel of peptides representing all possible position^-2 residues. This approach generated insights into the basis for the common natural class I and II specificities, demonstrated an alternative basis for a rare natural class III specificity for Asp^-2, and revealed a novel specificity for Arg^-2 that has not been reported in natural PDZ domains. A structure of a PDZ-peptide complex explained the minimum requirement for switching specificity from class I ligands containing Thr/Ser^-2 to class II ligands containing hydrophobic residues at position^-2. A second structure explained the molecular basis for the specificity for ligands containing Arg^-2. Overall, the evolved PDZ variants greatly expand our understanding of site^-2 specificities and the variants themselves may prove useful as building blocks for synthetic biology.

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.

The leucine-rich repeat signaling scaffolds Shoc2 and Erbin: cellular mechanism and role in disease

Leucine-rich repeat-containing proteins (LRR proteins) are involved in supporting a large number of cellular functions. In this review, we summarize recent advancements in understanding functions of the LRR proteins as signaling scaffolds. In particular, we explore what we have learned about the mechanisms of action of the LRR scaffolds Shoc2 and Erbin and their roles in normal development and disease. We discuss Shoc2 and Erbin in the context of their multiple known interacting partners in various cellular processes and summarize often unexpected functions of these proteins through analysis of their roles in human pathologies. We also review these LRR scaffold proteins as promising therapeutic targets and biomarkers with potential application across various pathologies.

Specificity in PDZ-peptide interaction networks: Computational analysis and review

Globular PDZ domains typically serve as protein-protein interaction modules that regulate a wide variety of cellular functions via recognition of short linear motifs (SLiMs). Often, PDZ mediated-interactions are essential components of macromolecular complexes, and disruption affects the entire scaffold. Due to their roles as linchpins in trafficking and signaling pathways, PDZ domains are attractive targets: both for controlling viral pathogens, which bind PDZ domains and hijack cellular machinery, as well as for developing therapies to combat human disease. However, successful therapeutic interventions that avoid off-target effects are a challenge, because each PDZ domain interacts with a number of cellular targets, and specific binding preferences can be difficult to decipher. Over twenty-five years of research has produced a wealth of data on the stereochemical preferences of individual PDZ proteins and their binding partners. Currently the field lacks a central repository for this information. Here, we provide this important resource and provide a manually curated, comprehensive list of the 271 human PDZ domains. We use individual domain, as well as recent genomic and proteomic, data in order to gain a holistic view of PDZ domains and interaction networks, arguing this knowledge is critical to optimize targeting selectivity and to benefit human health.

Comprehensive analysis of all evolutionary paths between two divergent PDZ domain specificities

To understand the molecular evolution of functional diversity in protein families, we comprehensively investigated the consequences of all possible mutation combinations separating two peptide-binding domains with highly divergent specificities. We analyzed the Erbin PDZ domain (Erbin-PDZ), which exhibits canonical type I specificity, and a synthetic Erbin-PDZ variant (E-14) that differs at six positions and exhibits an atypical specificity that closely resembles that of the natural Pdlim4 PDZ domain (Pdlim4-PDZ). We constructed a panel of 64 PDZ domains covering all possible transitions between Erbin-PDZ and E-14 (i.e., the panel contained variants with all possible combinations of either the Erbin-PDZ or E-14 sequence at the six differing positions). We assessed the specificity profiles of the 64 PDZ domains using a C-terminal phage-displayed peptide library containing all possible genetically encoded heptapeptides. The specificity profiles clustered into six distinct groups, showing that intermediate domains can be nodes for the evolution of divergent functions. Remarkably, three substitutions were sufficient to convert the specificity of Erbin-PDZ to that of Pdlim4-PDZ, whereas Pdlim4-PDZ contains 71 differences relative to Erbin-PDZ. X-ray crystallography revealed the structural basis for specificity transition: a single substitution in the center of the binding site, supported by contributions from auxiliary substitutions, altered the main chain conformation of the peptide ligand to resemble that of ligands bound to Pdlim4-PDZ. Our results show that a very small set of mutations can dramatically alter protein specificity, and these findings support the hypothesis whereby complex protein functions evolve by gene duplication followed by cumulative mutations.

Scribble, Erbin, and Lano redundantly regulate epithelial polarity and apical adhesion complex

The basolateral protein Scribble (Scrib), a member of the LAP protein family, is essential for epithelial apicobasal polarity (ABP) in Drosophila However, a conserved function for this protein in mammals is unclear. Here we show that the crucial role for Scrib in ABP has remained obscure due to the compensatory function of two other LAP proteins, Erbin and Lano. A combined Scrib/Erbin/Lano knockout disorganizes the cell-cell junctions and the cytoskeleton. It also results in mislocalization of several apical (Par6, aPKC, and Pals1) and basolateral (Llgl1 and Llgl2) identity proteins. These defects can be rescued by the conserved "LU" region of these LAP proteins. Structure-function analysis of this region determined that the so-called LAPSDb domain is essential for basolateral targeting of these proteins, while the LAPSDa domain is essential for supporting the membrane basolateral identity and binding to Llgl. In contrast to the key role in Drosophila, mislocalization of Llgl proteins does not appear to be critical in the scrib ABP phenotype.

A Multireporter Bacterial 2-Hybrid Assay for the High-Throughput and Dynamic Assay of PDZ Domain-Peptide Interactions

The accurate determination of protein-protein interactions has been an important focus of molecular biology toward which much progress has been made due to the continuous development of existing and new technologies. However, current methods can have limitations, including scale and restriction to high affinity interactions, limiting our understanding of a large subset of these interactions. Here, we describe a modified bacterial-hybrid assay that employs combined selectable and scalable reporters that enable the sensitive screening of large peptide libraries followed by the sorting of positive interactions by the level of reporter output. We have applied this tool to characterize a set of human and E. coli PDZ domains. Our results are consistent with prior characterization of these proteins, and the improved sensitivity increases our ability to predict known and novel in vivo binding partners. This approach allows for the recovery of a wide range of affinities with a high throughput method that does not sacrifice the scale of the screen.

Epitope Mapping Using Yeast Display and Next Generation Sequencing

Monoclonal antibodies are the largest class of therapeutic proteins due in part to their ability to bind an antigen with a high degree of affinity and specificity. A precise determination of their epitope is important for gaining insights into their therapeutic mechanism of action and to help differentiate antibodies that bind the same antigen. Here, we describe a method to precisely and efficiently map the epitopes of multiple antibodies in parallel over the course of just several weeks. This approach is based on a combination of rational library design, yeast surface display, and next generation DNA sequencing and provides quantitative insights into the epitope residues most critical for the antibody-antigen interaction. As an example, we will use this method to map the epitopes of several antibodies that neutralize alpha toxin from Staphylococcus aureus.

Drosophila melanogaster Guk-holder interacts with the Scribbled PDZ1 domain and regulates epithelial development with Scribbled and Discs Large

Epithelial cell polarity is controlled by components of the Scribble polarity module, and its regulation is critical for tissue architecture and cell proliferation and migration. In Drosophila melanogaster, the adaptor protein Guk-holder (Gukh) binds to the Scribbled (Scrib) and Discs Large (Dlg) components of the Scribble polarity module and plays an important role in the formation of neuromuscular junctions. However, Gukh's role in epithelial tissue formation and the molecular basis for the Scrib-Gukh interaction remain to be defined. We now show using isothermal titration calorimetry that the Scrib PDZ1 domain is the major site for an interaction with Gukh. Furthermore, we defined the structural basis of this interaction by determining the crystal structure of the Scrib PDZ1-Gukh complex. The C-terminal PDZ-binding motif of Gukh is located in the canonical ligand-binding groove of Scrib PDZ1 and utilizes an unusually extensive network of hydrogen bonds and ionic interactions to enable binding to PDZ1 with high affinity. We next examined the role of Gukh along with those of Scrib and Dlg in Drosophila epithelial tissues and found that Gukh is expressed in larval-wing and eye-epithelial tissues and co-localizes with Scrib and Dlg at the apical cell cortex. Importantly, we show that Gukh functions with Scrib and Dlg in the development of Drosophila epithelial tissues, with depletion of Gukh enhancing the eye- and wing-tissue defects caused by Scrib or Dlg depletion. Overall, our findings reveal that Scrib's PDZ1 domain functions in the interaction with Gukh and that the Scrib-Gukh interaction has a key role in epithelial tissue development in Drosophila.

Structural basis for the differential interaction of Scribble PDZ domains with the guanine nucleotide exchange factor β-PIX

Scribble is a highly conserved protein regulator of cell polarity that has been demonstrated to function as a tumor suppressor or, conversely, as an oncogene in a context-dependent manner, and it also controls many physiological processes ranging from immunity to memory. Scribble consists of a leucine-rich repeat domain and four PDZ domains, with the latter being responsible for most of Scribble's complex formation with other proteins. Given the similarities of the Scribble PDZ domain sequences in their binding grooves, it is common for these domains to show overlapping preferences for the same ligand. Yet, Scribble PDZ domains can still exhibit unique binding profiles toward other ligands. This raises the fundamental question as to how these PDZ domains discriminate ligands and exert specificities in Scribble complex formation. To better understand how Scribble PDZ domains direct cell polarity signaling, we investigated here their interactions with the well-characterized Scribble binding partner β-PIX, a guanine nucleotide exchange factor. We report the interaction profiles of all isolated Scribble PDZ domains with a β-PIX peptide. We show that Scribble PDZ1 and PDZ3 are the major interactors with β-PIX and reveal a distinct binding hierarchy in the interactions between the individual Scribble PDZ domains and β-PIX. Furthermore, using crystal structures of PDZ1 and PDZ3 bound to β-PIX, we define the structural basis for Scribble's ability to specifically engage β-PIX via its PDZ domains and provide a mechanistic platform for understanding Scribble-β-PIX-coordinated cellular functions such as directional cell migration.

MFPred: Rapid and accurate prediction of protein-peptide recognition multispecificity using self-consistent mean field theory

Multispecificity-the ability of a single receptor protein molecule to interact with multiple substrates-is a hallmark of molecular recognition at protein-protein and protein-peptide interfaces, including enzyme-substrate complexes. The ability to perform structure-based prediction of multispecificity would aid in the identification of novel enzyme substrates, protein interaction partners, and enable design of novel enzymes targeted towards alternative substrates. The relatively slow speed of current biophysical, structure-based methods limits their use for prediction and, especially, design of multispecificity. Here, we develop a rapid, flexible-backbone self-consistent mean field theory-based technique, MFPred, for multispecificity modeling at protein-peptide interfaces. We benchmark our method by predicting experimentally determined peptide specificity profiles for a range of receptors: protease and kinase enzymes, and protein recognition modules including SH2, SH3, MHC Class I and PDZ domains. We observe robust recapitulation of known specificities for all receptor-peptide complexes, and comparison with other methods shows that MFPred results in equivalent or better prediction accuracy with a ~10-1000-fold decrease in computational expense. We find that modeling bound peptide backbone flexibility is key to the observed accuracy of the method. We used MFPred for predicting with high accuracy the impact of receptor-side mutations on experimentally determined multispecificity of a protease enzyme. Our approach should enable the design of a wide range of altered receptor proteins with programmed multispecificities.

How to Train a Cell-Cutting-Edge Molecular Tools

In biological systems, the formation of molecular complexes is the currency for all cellular processes. Traditionally, functional experimentation was targeted to single molecular players in order to understand its effects in a cell or animal phenotype. In the last few years, we have been experiencing rapid progress in the development of ground-breaking molecular biology tools that affect the metabolic, structural, morphological, and (epi)genetic instructions of cells by chemical, optical (optogenetic) and mechanical inputs. Such precise dissection of cellular processes is not only essential for a better understanding of biological systems, but will also allow us to better diagnose and fix common dysfunctions. Here, we present several of these emerging and innovative techniques by providing the reader with elegant examples on how these tools have been implemented in cells, and, in some cases, organisms, to unravel molecular processes in minute detail. We also discuss their advantages and disadvantages with particular focus on their translation to multicellular organisms for in vivo spatiotemporal regulation. We envision that further developments of these tools will not only help solve the processes of life, but will give rise to novel clinical and industrial applications.

Investigating Protein-Peptide Interactions Using the Schrödinger Computational Suite

The Schrödinger software suite contains a broad array of computational chemistry and molecular modeling tools that can be used to study the interaction of peptides with proteins. These include molecular docking using Glide and Piper, relative binding free energy predictions with FEP+, conformational searches using MacroModel and Desmond, and structural refinement using Prime and PrimeX. In this review we provide a comprehensive overview of these tools and describe their potential application in the identification and optimization of peptide ligands for proteins.

Fluorescence-based ATG8 sensors monitor localization and function of LC3/GABARAP proteins

Autophagy is a cellular surveillance pathway that balances metabolic and energy resources and transports specific cargos, including damaged mitochondria, other broken organelles, or pathogens for degradation to the lysosome. Central components of autophagosomal biogenesis are six members of the LC3 and GABARAP family of ubiquitin-like proteins (mATG8s). We used phage display to isolate peptides that possess bona fide LIR (LC3-interacting region) properties and are selective for individual mATG8 isoforms. Sensitivity of the developed sensors was optimized by multiplication, charge distribution, and fusion with a membrane recruitment (FYVE) or an oligomerization (PB1) domain. We demonstrate the use of the engineered peptides as intracellular sensors that recognize specifically GABARAP, GABL1, GABL2, and LC3C, as well as a bispecific sensor for LC3A and LC3B. By using an LC3C-specific sensor, we were able to monitor recruitment of endogenous LC3C to Salmonella during xenophagy, as well as to mitochondria during mitophagy. The sensors are general tools to monitor the fate of mATG8s and will be valuable in decoding the biological functions of the individual LC3/GABARAPs.

Exploring Multiple Binding Modes Using Confined Replica Exchange Molecular Dynamics

Molecular docking is extensively applied to determine the position of a ligand on its receptor despite the rather poor correspondence between docking scores and experimental binding affinities found in several studies, especially for systems structurally unrelated with those used in the scoring functions' training sets. Here, we present a method for the prediction of binding modes and binding free energies, which uses replica exchange molecular dynamics in combination with a receptor-shaped piecewise potential, confining the ligand in the proximity of the receptor surface and limiting the accessible conformational space of interest. We assess our methodology with a set of protein receptor-ligand test cases. In every case studied, the method is able to locate the ligand on the experimentally known receptor binding site, and it gives as output the binding free energy. The added value of our approach with respect to other available methods is that it quickly performs a conformational space search, providing a set of bound (or unbound) configurations, which can be used to determine phenomenological structural and energetic properties of an experimental binding state as a result of contributions provided by diversified multiple binding poses.

Homogeneous and Nonradioactive High-Throughput Screening Platform for the Characterization of Kinase Inhibitors in Cell Lysates

Protein kinases are directly implicated in many human diseases; therefore, kinase inhibitors show great promises as new therapeutic drugs. In an effort to facilitate the screening and the characterization of kinase inhibitors, a novel application of the AlphaScreen technology was developed to monitor JNK activity from (1) purified kinase preparations and (2) endogenous kinase from cell lysates preactivated with different cytokines. The authors confirmed that both adenosine triphosphate (ATP) competitive as well as peptide-based JNK inhibitors were able to block the activity of both recombinant and HepG2 endogenous JNK activity. Using the same luminescence technique adapted for binding studies, the authors characterized peptide inhibitor mechanisms by measuring the binding affinity of the inhibitors for JNK. Because of the versatility of the technology, this cell-based JNK kinase assay could be adapted to other kinases and would represent a powerful tool to evaluate endogenous kinase activity and test a large number of potential inhibitors in a more physiologically relevant environment.

Local RhoA activation induces cytokinetic furrows independent of spindle position and cell cycle stage

The GTPase RhoA promotes contractile ring assembly and furrow ingression during cytokinesis. Although many factors that regulate RhoA during cytokinesis have been characterized, the spatiotemporal regulatory logic remains undefined. We have developed an optogenetic probe to gain tight spatial and temporal control of RhoA activity in mammalian cells and demonstrate that cytokinetic furrowing is primarily regulated at the level of RhoA activation. Light-mediated recruitment of a RhoGEF domain to the plasma membrane leads to rapid induction of RhoA activity, leading to assembly of cytokinetic furrows that partially ingress. Furthermore, furrow formation in response to RhoA activation is not temporally or spatially restricted. RhoA activation is sufficient to generate furrows at both the cell equator and cell poles, in both metaphase and anaphase. Remarkably, furrow formation can be initiated in rounded interphase cells, but not adherent cells. These results indicate that RhoA activation is sufficient to induce assembly of functional contractile rings and that cell rounding facilitates furrow formation.

Ligand binding to the PDZ domains of postsynaptic density protein 95

Cellular scaffolding and signalling is generally governed by multidomain proteins, where each domain has a particular function. Postsynaptic density protein 95 (PSD-95) is involved in synapse formation and is a typical example of such a multidomain protein. Protein-protein interactions of PSD-95 are well studied and include the following three protein ligands: (i)N-methyl-d-aspartate-type ionotropic glutamate receptor subunit GluN2B, (ii) neuronal nitric oxide synthase and (iii) cysteine-rich protein (CRIPT), all of which bind to one or more of the three PDZ domains in PSD-95. While interactions for individual PDZ domains of PSD-95 have been well studied, less is known about the influence of neighbouring domains on the function of the respective individual domain. We therefore performed a systematic study on the ligand-binding kinetics of PSD-95 using constructs of different size for PSD-95 and its ligands. Regarding the canonical peptide-binding pocket and relatively short peptides (up to 15-mer), the PDZ domains in PSD-95 by and large work as individual binding modules. However, in agreement with previous studies, residues outside of the canonical binding pocket modulate the affinity of the ligands. In particular, the dissociation of the 101 amino acid CRIPT from PSD-95 is slowed down at least 10-fold for full-length PSD-95 when compared with the individual PDZ3 domain.

Adaptable hydrogel networks with reversible linkages for tissue engineering

Adaptable hydrogels have recently emerged as a promising platform for three-dimensional (3D) cell encapsulation and culture. In conventional, covalently crosslinked hydrogels, degradation is typically required to allow complex cellular functions to occur, leading to bulk material degradation. In contrast, adaptable hydrogels are formed by reversible crosslinks. Through breaking and re-formation of the reversible linkages, adaptable hydrogels can be locally modified to permit complex cellular functions while maintaining their long-term integrity. In addition, these adaptable materials can have biomimetic viscoelastic properties that make them well suited for several biotechnology and medical applications. In this review, an overview of adaptable-hydrogel design considerations and linkage selections is presented, with a focus on various cell-compatible crosslinking mechanisms that can be exploited to form adaptable hydrogels for tissue engineering.

Individual protomers of a G protein-coupled receptor dimer integrate distinct functional modules

Recent advances in proteomic technology reveal G-protein-coupled receptors (GPCRs) are organized as large, macromolecular protein complexes in cell membranes, adding a new layer of intricacy to GPCR signaling. We previously reported the α_1D-adrenergic receptor (ADRA1D)—a key regulator of cardiovascular, urinary and CNS function—binds the syntrophin family of PDZ domain proteins (SNTA, SNTB1, and SNTB2) through a C-terminal PDZ ligand interaction, ensuring receptor plasma membrane localization and G-protein coupling. To assess the uniqueness of this novel GPCR complex, 23 human GPCRs containing Type I PDZ ligands were subjected to TAP/MS proteomic analysis. Syntrophins did not interact with any other GPCRs. Unexpectedly, a second PDZ domain protein, scribble (SCRIB), was detected in ADRA1D complexes. Biochemical, proteomic, and dynamic mass redistribution analyses indicate syntrophins and SCRIB compete for the PDZ ligand, simultaneously exist within an ADRA1D multimer, and impart divergent pharmacological properties to the complex. Our results reveal an unprecedented modular dimeric architecture for the ADRA1D in the cell membrane, providing unexpected opportunities for fine-tuning receptor function through novel protein interactions in vivo, and for intervening in signal transduction with small molecules that can stabilize or disrupt unique GPCR:PDZ protein interfaces.

Multiple paramagnetic effects through a tagged reporter protein

Paramagnetic effects provide unique information about the structure and dynamics of biomolecules. We developed a method in which the lanthanoid tag is not directly attached to the protein of interest, but instead to a "reporter" protein, which binds and then transmits paramagnetic information to the target. The designed method allows access to a large number of paramagnetic restraints and residual dipolar couplings produced from independent molecular alignments in high-molecular-weight proteins with unknown 3D structure.

Predicting affinity- and specificity-enhancing mutations at protein-protein interfaces

Manipulations of PPIs (protein-protein interactions) are important for many biological applications such as synthetic biology and drug design. Combinatorial methods have been traditionally used for such manipulations, failing, however, to explain the effects achieved. We developed a computational method for prediction of changes in free energy of binding due to mutation that bring about deeper understanding of the molecular forces underlying binding interactions. Our method could be used for computational scanning of binding interfaces and subsequent analysis of the interfacial sequence optimality. The computational method was validated in two biological systems. Computational saturated mutagenesis of a high-affinity complex between an enzyme AChE (acetylcholinesterase) and a snake toxin Fas (fasciculin) revealed the optimal nature of this interface with only a few predicted affinity-enhancing mutations. Binding measurements confirmed high optimality of this interface and identified a few mutations that could further improve interaction fitness. Computational interface scanning of a medium-affinity complex between TIMP-2 (tissue inhibitor of metalloproteinases-2) and MMP (matrix metalloproteinase) 14 revealed a non-optimal nature of the binding interface with multiple mutations predicted to stabilize the complex. Experimental results corroborated our computational predictions, identifying a large number of mutations that improve the binding affinity for this interaction and some mutations that enhance binding specificity. Overall, our computational protocol greatly facilitates the discovery of affinity- and specificity-enhancing mutations and thus could be applied for design of potent and highly specific inhibitors of any PPI.

Mapping of the binding landscape for a picomolar protein-protein complex through computation and experiment

Our understanding of protein evolution would greatly benefit from mapping of binding landscapes, i.e., changes in protein-protein binding affinity due to all single mutations. However, experimental generation of such landscapes is a tedious task due to a large number of possible mutations. Here, we use a simple computational protocol to map the binding landscape for two homologous high-affinity complexes, involving a snake toxin fasciculin and acetylcholinesterase from two different species. To verify our computational predictions, we experimentally measure binding between 25 Fas mutants and the 2 enzymes. Both computational and experimental results demonstrate that the Fas sequence is close to the optimum when interacting with its targets, yet a few mutations could further improve Kd, kon, and koff. Our computational predictions agree well with experimental results and generate distributions similar to those observed in other high-affinity PPIs, demonstrating the potential of simple computational protocols in capturing realistic binding landscapes.

C-terminal COOH of integrin β1 is necessary for β1 association with the kindlin-2 adapter protein

Protein-protein interactions are driving forces in cellular processes. As a prime example, transmembrane integrins link extracellular matrix and intracellular proteins, resulting in bidirectional signaling that regulates cell migration, proliferation, differentiation, and survival. Here we provide the first evidence that interaction between the integrin β1 cytoplasmic tail and kindlin-2, a member of a family of adapters implicated in human disease pathogenesis, is mainly governed by the β1 C-terminal carboxylate moiety and is required for laterality organ development in zebrafish. Affinity measurements indicate that this unusual protein-protein interaction mode is coordinated by a putative carboxylate-binding motif in the kindlin-2 FERM subdomain F3. Contrary to the C terminus of proteins that engage PDZ domains, the C-terminal three residues of β1, per se, do not contribute to kindlin-2 binding or to laterality organ development. Thus, by employing zebrafish as an in situ physiological tool to correlate protein structure and function, we have discovered an unexpected association chemistry between an integrin and a key adapter involved in integrin signaling.

Stereochemical preferences modulate affinity and selectivity among five PDZ domains that bind CFTR: comparative structural and sequence analyses

PDZ domain interactions are involved in signaling and trafficking pathways that coordinate crucial cellular processes. Alignment-based PDZ binding motifs identify the few most favorable residues at certain positions along the peptide backbone. However, sequences that bind the CAL (CFTR-associated ligand) PDZ domain reveal only a degenerate motif that overpredicts the true number of high-affinity interactors. Here, we combine extended peptide-array motif analysis with biochemical techniques to show that non-motif "modulator" residues influence CAL binding. The crystallographic structures of 13 CAL:peptide complexes reveal defined, but accommodating stereochemical environments at non-motif positions, which are reflected in modulator preferences uncovered by multisequence substitutional arrays. These preferences facilitate the identification of high-affinity CAL binding sequences and differentially affect CAL and NHERF PDZ binding. As a result, they also help determine the specificity of a PDZ domain network that regulates the trafficking of CFTR at the apical membrane.

Inhibitory mechanism of an allosteric antibody targeting the glucagon receptor

Elevated glucagon levels and increased hepatic glucagon receptor (GCGR) signaling contribute to hyperglycemia in type 2 diabetes. We have identified a monoclonal antibody that inhibits GCGR, a class B G-protein coupled receptor (GPCR), through a unique allosteric mechanism. Receptor inhibition is mediated by the binding of this antibody to two distinct sites that lie outside of the glucagon binding cleft. One site consists of a patch of residues that are surface-exposed on the face of the extracellular domain (ECD) opposite the ligand-binding cleft, whereas the second binding site consists of residues in the αA helix of the ECD. A docking model suggests that the antibody does not occlude the ligand-binding cleft. We solved the crystal structure of GCGR ECD containing a naturally occurring G40S mutation and found a shift in the register of the αA helix that prevents antibody binding. We also found that alterations in the αA helix impact the normal function of GCGR. We present a model for the allosteric inhibition of GCGR by a monoclonal antibody that may form the basis for the development of allosteric modulators for the treatment of diabetes and other class B GPCR-related diseases.

Structure-based design of PDZ ligands as inhibitors of 5-HT(2A) receptor/PSD-95 PDZ1 domain interaction possessing anti-hyperalgesic activity

Disrupting the interaction between the PDZ protein PSD-95 and the C-terminal domain of the 5-HT2A serotonin receptor has been shown to reduce hyperalgesia in a rodent model of neuropathic pain. Here, we designed and synthesized PDZ ligands capable of binding to the first PDZ domain (PDZ1) of the PSD-95 protein and evaluated their biological activity in vitro and in vivo. A series of substituted indoles was identified by docking simulations, and six novel analogues were synthesized. Three analogues displayed strong interactions with the first PDZ domain (PDZ1) of PDZ-95 in (1)H-(15)N heteronuclear single-quantum coherence (HSQC) experiments and two of them were able to inhibit the interaction between PSD-95 and the 5-HT2A receptor in vitro. We identified compound 8b as the analogue able to significantly suppress mechanical hyperalgesia in an experimental model of traumatic neuropathic pain in the rat. This effect was suppressed by the coadministration of the 5-HT2A receptor antagonist M100907, consistent with an inhibitory effect upon 5-HT2A receptor/PSD-95 interaction. Finally, we determined an NMR-restraint driven model structure for the PSD95 PDZ1/8b complex, which confirms that indole 8b binds to the putative PDZ-ligand binding site.

Probing the role of backbone hydrogen bonds in protein-peptide interactions by amide-to-ester mutations

One of the most frequent protein-protein interaction modules in mammalian cells is the postsynaptic density 95/discs large/zonula occludens 1 (PDZ) domain, involved in scaffolding and signaling and emerging as an important drug target for several diseases. Like many other protein-protein interactions, those of the PDZ domain family involve formation of intermolecular hydrogen bonds: C-termini or internal linear motifs of proteins bind as β-strands to form an extended antiparallel β-sheet with the PDZ domain. Whereas extensive work has focused on the importance of the amino acid side chains of the protein ligand, the role of the backbone hydrogen bonds in the binding reaction is not known. Using amide-to-ester substitutions to perturb the backbone hydrogen-bonding pattern, we have systematically probed putative backbone hydrogen bonds between four different PDZ domains and peptides corresponding to natural protein ligands. Amide-to-ester mutations of the three C-terminal amides of the peptide ligand severely affected the affinity with the PDZ domain, demonstrating that hydrogen bonds contribute significantly to ligand binding (apparent changes in binding energy, ΔΔG = 1.3 to >3.8 kcal mol(-1)). This decrease in affinity was mainly due to an increase in the dissociation rate constant, but a significant decrease in the association rate constant was found for some amide-to-ester mutations suggesting that native hydrogen bonds have begun to form in the transition state of the binding reaction. This study provides a general framework for studying the role of backbone hydrogen bonds in protein-peptide interactions and for the first time specifically addresses these for PDZ domain-peptide interactions.

Predicting PDZ domain mediated protein interactions from structure

Background

PDZ domains are structural protein domains that recognize simple linear amino acid motifs, often at protein C-termini, and mediate protein-protein interactions (PPIs) in important biological processes, such as ion channel regulation, cell polarity and neural development. PDZ domain-peptide interaction predictors have been developed based on domain and peptide sequence information. Since domain structure is known to influence binding specificity, we hypothesized that structural information could be used to predict new interactions compared to sequence-based predictors.

Results

We developed a novel computational predictor of PDZ domain and C-terminal peptide interactions using a support vector machine trained with PDZ domain structure and peptide sequence information. Performance was estimated using extensive cross validation testing. We used the structure-based predictor to scan the human proteome for ligands of 218 PDZ domains and show that the predictions correspond to known PDZ domain-peptide interactions and PPIs in curated databases. The structure-based predictor is complementary to the sequence-based predictor, finding unique known and novel PPIs, and is less dependent on training–testing domain sequence similarity. We used a functional enrichment analysis of our hits to create a predicted map of PDZ domain biology. This map highlights PDZ domain involvement in diverse biological processes, some only found by the structure-based predictor. Based on this analysis, we predict novel PDZ domain involvement in xenobiotic metabolism and suggest new interactions for other processes including wound healing and Wnt signalling.

Conclusions

We built a structure-based predictor of PDZ domain-peptide interactions, which can be used to scan C-terminal proteomes for PDZ interactions. We also show that the structure-based predictor finds many known PDZ mediated PPIs in human that were not found by our previous sequence-based predictor and is less dependent on training–testing domain sequence similarity. Using both predictors, we defined a functional map of human PDZ domain biology and predict novel PDZ domain function. Users may access our structure-based and previous sequence-based predictors at http://webservice.baderlab.org/domains/POW.

Generating thermal stable variants of protein domains through phage display

Often in protein design research, one desires to generate thermally stable variants of a protein or domain. One route to identifying mutations that yield domains that remain folded and active at a higher temperature is through the use of directed evolution. A library of protein domain variants can be generated by mutagenic PCR, expressed on the surface of bacteriophage M13, and subjected to heat, such that the unfolded forms of the domain, showing reduced or no binding activity, are lost during subsequent affinity selection, whereas variants that still retain binding to their target are selected and enriched with each subsequent round of affinity selection. This approach takes advantage of the fact that bacteriophage M13 particles are heat stable and resistant to many proteases and protein denaturants. We present the application of this general approach to generating thermally stable variants of a eukaryotic peptide-binding domain. The benefits of producing such variants are that they typically express at high levels in Escherichia coli (30-60 mg/L shake flask) and remain soluble in solution at higher concentrations for longer periods of time than the wild-type form of the domain. The process of library generation and screening generally requires about one month of effort, and yields variants with >10 °C increase in thermal stability, as measured in a simple fluorescence-based thermal shift assay. It is anticipated that thermally stable variants will serve as excellent scaffolds for generating affinity reagents to a variety of targets of interest.

Ligand binding by PDZ domains

The postsynaptic density protein-95/disks large/zonula occludens-1 (PDZ) protein domain family is one of the most common protein-protein interaction modules in mammalian cells, with paralogs present in several hundred human proteins. PDZ domains are found in most cell types, but neuronal proteins, for example, are particularly rich in these domains. The general function of PDZ domains is to bring proteins together within the appropriate cellular compartment, thereby facilitating scaffolding, signaling, and trafficking events. The many functions of PDZ domains under normal physiological as well as pathological conditions have been reviewed recently. In this review, we focus on the molecular details of how PDZ domains bind their protein ligands and their potential as drug targets in this context.