Crystal structures of a complexed and peptide-free membrane protein-binding domain: molecular basis of peptide recognition by PDZ

Unconventional PDZ Recognition Revealed in α7 nAChR-PICK1 Complexes

PDZ domains are modular domains that conventionally bind to C terminal or internal motifs of target proteins to control cellular functions through the regulation of protein complex assemblies. Almost all reported structures of PDZ-target protein complexes rely on fragments or peptides as target proteins. No intact target protein complexed with PDZ was structurally characterized. In this study, we used NMR spectroscopy and other biochemistry and biophysics tools to uncover insights into structural coupling between the PDZ domain of protein interacting with C-kinase 1 (PICK1) and α7 nicotinic acetylcholine receptors (α7 nAChR). Notably, the intracellular domains of both α7 nAChR and PICK1 PDZ exhibit a high degree of plasticity in their coupling. Specifically, the MA helix of α7 nAChR interacts with residues lining the canonical binding site of the PICK1 PDZ, while flexible loops also engage in protein-protein interactions. Both hydrophobic and electrostatic interactions mediate the coupling. Overall, the resulting structure of the α7 nAChR-PICK1 complex reveals an unconventional PDZ binding mode, significantly expanding the repertoire of functionally important PDZ interactions.

Voltage-induced calcium release in Caenorhabditis elegans body muscles

Type 1 voltage-activated calcium channels (CaV1) in the plasma membrane trigger calcium release from the sarcoplasmic reticulum (SR) by two mechanisms. In voltage-induced calcium release (VICR), CaV1 voltage sensing domains are directly coupled to ryanodine receptors (RYRs), an SR calcium channel. In calcium-induced calcium release (CICR), calcium ions flowing through activated CaV1 channels bind and activate RYR channels. VICR is thought to occur exclusively in vertebrate skeletal muscle while CICR occurs in all other muscles (including all invertebrate muscles). Here, we use calcium-activated SLO-2 potassium channels to analyze CaV1-SR coupling in Caenorhabditis elegans body muscles. SLO-2 channels were activated by both VICR and external calcium. VICR-mediated SLO-2 activation requires two SR calcium channels (RYRs and IP3 Receptors), JPH-1/Junctophilin, a PDZ (PSD95, Dlg1, ZO-1 domain) binding domain (PBD) at EGL-19/CaV1's carboxy-terminus, and SHN-1/Shank (a scaffolding protein that binds EGL-19's PBD). Thus, VICR occurs in invertebrate muscles.

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.

Functional protein dynamics in a crystal

Proteins are molecular machines and to understand how they work, we need to understand how they move. New pump-probe time-resolved X-ray diffraction methods open up ways to initiate and observe protein motions with atomistic detail in crystals on biologically relevant timescales. However, practical limitations of these experiments demands parallel development of effective molecular dynamics approaches to accelerate progress and extract meaning. Here, we establish robust and accurate methods for simulating dynamics in protein crystals, a nontrivial process requiring careful attention to equilibration, environmental composition, and choice of force fields. With more than seven milliseconds of sampling of a single chain, we identify critical factors controlling agreement between simulation and experiments and show that simulated motions recapitulate ligand-induced conformational changes. This work enables a virtuous cycle between simulation and experiments for visualizing and understanding the basic functional motions of proteins.

Functional Protein Dynamics in a Crystal

Proteins are molecular machines and to understand how they work, we need to understand how they move. New pump-probe time-resolved X-ray diffraction methods open up ways to initiate and observe protein motions with atomistic detail in crystals on biologically relevant timescales. However, practical limitations of these experiments demands parallel development of effective molecular dynamics approaches to accelerate progress and extract meaning. Here, we establish robust and accurate methods for simulating dynamics in protein crystals, a nontrivial process requiring careful attention to equilibration, environmental composition, and choice of force fields. With more than seven milliseconds of sampling of a single chain, we identify critical factors controlling agreement between simulation and experiments and show that simulated motions recapitulate ligand-induced conformational changes. This work enables a virtuous cycle between simulation and experiments for visualizing and understanding the basic functional motions of proteins.

The microvillar protocadherin CDHR5 associates with EBP50 to promote brush border assembly

Transporting epithelial cells of the gut and kidney interact with their luminal environment through a densely packed collection of apical microvilli known as a brush border (BB). Proper brush border assembly depends on the intermicrovillar adhesion complex (IMAC), a protocadherin-based adhesion complex found at the distal tips of microvilli that mediates adhesion between neighboring protrusions to promote their organized packing. Loss of the IMAC adhesion molecule Cadherin-related family member 5 (CDHR5) results in significant brush border defects, though the functional properties of this protocadherin have not been thoroughly explored. Here, we show that the cytoplasmic tail of CDHR5 contributes to its correct apical targeting and functional properties in an isoform-specific manner. Library screening identified the Ezrin-associated scaffolds EBP50 and E3KARP as cytoplasmic binding partners for CDHR5. Consistent with this, loss of EBP50 disrupted proper brush border assembly with cells exhibiting markedly reduced apical IMAC levels. Together, our results shed light on the apical targeting determinants of CDHR5 and further define the interactome of the IMAC involved in brush border assembly.

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.

Structure-based design of peptidomimetic inhibitors of PSD-95 with improved affinity for the PDZ3 domain

Aberrant brain-derived neurotrophic factor (BDNF) signaling has been proposed to contribute to the pathophysiology of depression and other neurological disorders such as Angelman syndrome. We have previously shown that targeting the tropomyosin receptor kinase B/postsynaptic density protein-95 (PSD-95) nexus in the BDNF signaling pathway by peptidomimetic inhibitors is a promising approach for therapeutic intervention. Here, we used structure-based knowledge to develop a new Syn3 peptidomimetic compound series that fuses peptides derived from the PSD-95-binding protein SynGAP to our prototype compound CN2097. The new compounds target the PSD-95 PDZ3 domain and adjoining αC helix to achieve bivalent binding that results in up to 7-fold stronger affinity compared to CN2097. These compounds were designed to improve CN2097 specificity for the PSD-95 PDZ3 domain, and structure-activity relationship studies were performed to improve their resistance to proteolysis.

Autoregulation of the LIM kinases by their PDZ domain

LIM domain kinases (LIMK) are important regulators of actin cytoskeletal remodeling. These protein kinases phosphorylate the actin depolymerizing factor cofilin to suppress filament severing, and are key nodes between Rho GTPase cascades and actin. The two mammalian LIMKs, LIMK1 and LIMK2, contain consecutive LIM domains and a PDZ domain upstream of the C-terminal kinase domain. The roles of the N-terminal regions are not fully understood, and the function of the PDZ domain remains elusive. Here, we determine the 2.0 Å crystal structure of the PDZ domain of LIMK2 and reveal features not previously observed in PDZ domains including a core-facing arginine residue located at the second position of the 'x-Φ-G-Φ' motif, and that the expected peptide binding cleft is shallow and poorly conserved. We find a distal extended surface to be highly conserved, and when LIMK1 was ectopically expressed in yeast we find targeted mutagenesis of this surface decreases growth, implying increased LIMK activity. PDZ domain LIMK1 mutants expressed in yeast are hyperphosphorylated and show elevated activity in vitro. This surface in both LIMK1 and LIMK2 is critical for autoregulation independent of activation loop phosphorylation. Overall, our study demonstrates the functional importance of the PDZ domain to autoregulation of LIMKs.

Sec8 specifically interacts with the PDZ2 domain of synapse associated protein 102 (SAP102)

The exocyst is an evolutionarily conserved protein complex tethering secretory vesicles before their docking and fusion with the plasma membrane. The complex also plays important roles in cell migration, synaptogenesis, and neurite outgrowth. One of its subunits, Sec8, was reported to interact with two major synaptic scaffolding proteins SAP102 and PSD-95 that share high sequence homology and contain three PDZ domains at their N-terminal region. The interaction is via the binding of the C-terminal ITTV motif in Sec8 to the PDZ domains of the two synaptic proteins. However, it remains elusive to which PDZ domain(s) Sec8 binds and how their interaction occurs. Here we reported a 2.5 Å resolution crystal structure of the C-terminal half of rat Sec8 containing the ITTV motif. The structure shows that Sec8 contains an enormously long helix at its C-terminus, which bears a unique long "spacer" of 14 residues to bridge the ITTV motif to the compact core of Sec8. We found that Sec8 preferentially binds PDZ2 over PDZ1 and PDZ3 of SAP102. Deletion of the spacer completely abolished the binding of Sec8 to SAP102. Overall, our structural studies, biochemical data and modeling analyses altogether provide an explanation for how Sec8 interacts with SAP102.

nNOS and Neurological, Neuropsychiatric Disorders: A 20-Year Story

In the central nervous system, nitric oxide (NO), a free gas with multitudinous bioactivities, is mainly produced from the oxidation of L-arginine by neuronal nitric oxide synthase (nNOS). In the past 20 years, the studies in our group and other laboratories have suggested a significant involvement of nNOS in a variety of neurological and neuropsychiatric disorders. In particular, the interactions between the PDZ domain of nNOS and its adaptor proteins, including post-synaptic density 95, the carboxy-terminal PDZ ligand of nNOS, and the serotonin transporter, significantly influence the subcellular localization and functions of nNOS in the brain. The nNOS-mediated protein-protein interactions provide new attractive targets and guide the discovery of therapeutic drugs for neurological and neuropsychiatric disorders. Here, we summarize the work on the roles of nNOS and its association with multiple adaptor proteins on neurological and neuropsychiatric disorders.

The Human PDZome 2.0: Characterization of a New Resource to Test for PDZ Interactions by Yeast Two-Hybrid

PSD95-disc large-zonula occludens (PDZ) domains are globular modules of 80-90 amino acids that co-evolved with multicellularity. They commonly bind to carboxy-terminal sequences of a plethora of membrane-associated proteins and influence their trafficking and signaling. We previously built a PDZ resource (PDZome) allowing us to unveil human PDZ interactions by Yeast two-hybrid. Yet, this resource is incomplete according to the current knowledge on the human PDZ proteome. Here we built the PDZome 2.0 library for Yeast two-hybrid, based on a PDZ library manually curated from online resources. The PDZome2.0 contains 305 individual clones (266 PDZ domains in isolation and 39 tandems), for which all boundaries were designed based on available PDZ structures. Using as bait the E6 oncoprotein from HPV16, a known promiscuous PDZ interactor, we show that PDZome 2.0 outperforms the previous resource.

Structure-based development of new cyclic compounds targeting PSD-95 PDZ3 domain

Aberrant BDNF signaling has been proposed to contribute to the pathophysiology of depression and other neurological disorders such as Angelman syndrome. We have previously shown that targeting the TrkB / PSD-95 nexus by peptidomimetic inhibitors is a promising approach for therapeutic intervention. Here we used structure-based knowledge to develop a new peptidomimetic compound series that fuses SynGAP-derived peptides to our prototype compound CN2097. These compounds target the PSD-95 PDZ3 domain and adjoining αC helix to achieve bivalent binding that results in up to 7-fold stronger affinity compared to CN2097. These compounds were designed to improve CN2097 specificity for the PDZ3 domain and limited SAR studies have been performed to improve their resistance to proteolysis.

Membrane recruitment of the polarity protein Scribble by the cell adhesion receptor TMIGD1

Scribble (Scrib) is a multidomain polarity protein and member of the leucine-rich repeat and PDZ domain (LAP) protein family. A loss of Scrib expression is associated with disturbed apical-basal polarity and tumor formation. The tumor-suppressive activity of Scrib correlates with its membrane localization. Despite the identification of numerous Scrib-interacting proteins, the mechanisms regulating its membrane recruitment are not fully understood. Here, we identify the cell adhesion receptor TMIGD1 as a membrane anchor of Scrib. TMIGD1 directly interacts with Scrib through a PDZ domain-mediated interaction and recruits Scrib to the lateral membrane domain in epithelial cells. We characterize the association of TMIGD1 with each Scrib PDZ domain and describe the crystal structure of the TMIGD1 C-terminal peptide complexed with PDZ domain 1 of Scrib. Our findings describe a mechanism of Scrib membrane localization and contribute to the understanding of the tumor-suppressive activity of Scrib.

Additive energetic contributions of multiple peptide positions determine the relative promiscuity of viral and human sequences for PDZ domain targets

Protein-protein interactions that involve recognition of short peptides are critical in cellular processes. Protein-peptide interaction surface areas are relatively small and shallow, and there are often overlapping specificities in families of peptide-binding domains. Therefore, dissecting selectivity determinants can be challenging. PDZ domains are a family of peptide-binding domains located in several intracellular signaling and trafficking pathways. These domains are also directly targeted by pathogens, and a hallmark of many oncogenic viral proteins is a PDZ-binding motif. However, amidst sequences that target PDZ domains, there is a wide spectrum in relative promiscuity. For example, the viral HPV16 E6 oncoprotein recognizes over double the number of PDZ domain-containing proteins as the cystic fibrosis transmembrane conductance regulator (CFTR) in the cell, despite similar PDZ targeting-sequences and identical motif residues. Here, we determine binding affinities for PDZ domains known to bind either HPV16 E6 alone or both CFTR and HPV16 E6, using peptides matching WT and hybrid sequences. We also use energy minimization to model PDZ-peptide complexes and use sequence analyses to investigate this difference. We find that while the majority of single mutations had marginal effects on overall affinity, the additive effect on the free energy of binding accurately describes the selectivity observed. Taken together, our results describe how complex and differing PDZ interactomes can be programmed in the cell.

Biological Role and Aberrant Overexpression of Syntenin-1 in Cancer: Potential Role as a Biomarker and Therapeutic Target

Syntenin-1 is a 298 amino acid protein codified by the melanoma differentiation-associated gene-9 (MDA-9). Structurally, it is composed of four domains: N-terminal, PDZ1, PDZ2, and C-terminal. The PDZ domains of syntenin-1 are involved in the stability and interaction with other molecules such as proteins, glycoproteins, and lipids. Domains are also associated with several biological functions such as the activation of signaling pathways related to cell-to-cell adhesion, signaling translation, and the traffic of intracellular lipids, among others. The overexpression of syntenin-1 has been reported in glioblastoma, colorectal, melanoma, lung, prostate, and breast cancer, which promotes tumorigenesis by regulating cell migration, invasion, proliferation, angiogenesis, apoptosis, and immune response evasion, and metastasis. The overexpression of syntenin-1 in samples has been associated with worst prognostic and recurrence, whereas the use of inhibitors such as shRNA, siRNA, and PDZli showed a diminution of the tumor size and reduction in metastasis and invasion. Syntenin-1 has been suggested as a potential biomarker and therapeutic target in cancer for developing more effective diagnostic/prognostic tests or passive/active immunotherapies.

A Significant Difference in Core PDZ Interactivity of SARS-CoV, SARS-CoV2 and MERS-CoV Protein E Peptide PDZ Motifs In Vitro

Small structural E protein of coronaviruses uses its C-terminal PDZ motif to compromise the cellular PDZ interactome. In this work we compared core PDZ interactivity of small (seven amino acids) peptide PDZ motifs, originating from the envelope proteins of recently transmitted coronaviruses SARS-CoV, SARS-CoV2, and MERS-CoV. As the interaction targets we used 23 domains of the largest PDZ proteins MUPP1/MPDZ and PATJ/INAD. Results revealed exceptional affinity and interaction promiscuity of MERS-CoV PDZ motif in vitro, suggesting an increased probability of potential PDZ targets in vivo. We hypothesize that together with its known ability to enter the cells from both apical and basolateral sites, this might further contribute to its elevated disruption of cellular PDZ pathways and higher virulence.

Allostery Frustrates the Experimentalist

Proteins interact with other proteins, with nucleic acids, lipids, carbohydrates and various small molecules in the living cell. These interactions have been quantified and structurally characterized in numerous studies such that we today have a comprehensive picture of protein structure and function. However, proteins are dynamic and even folded proteins are likely more heterogeneous than they appear in most descriptions. One property of proteins that relies on dynamics and heterogeneity is allostery, the ability of a protein to change structure and function upon ligand binding to an allosteric site. Over the last decades the concept of allostery was broadened to embrace all types of long-range interactions across a protein including purely entropic changes without a conformational change in single protein domains. But with this re-definition came a problem: How do we measure allostery? In this opinion, we discuss some caveats arising from the quantitative description of single-domain allostery from an experimental perspective and how the limitations cannot be separated from the definition of allostery per se. Furthermore, we attempt to tie together allostery with the concept of frustration in an effort to investigate the links between these two complex, and yet general, properties of proteins. We arrive at the conclusion that the sensitivity to perturbation of allosteric networks in single protein domains is too large for the networks to be of significant biological relevance.

Coevolution-based prediction of key allosteric residues for protein function regulation

Allostery is fundamental to many biological processes. Due to the distant regulation nature, how allosteric mutations, modifications, and effector binding impact protein function is difficult to forecast. In protein engineering, remote mutations cannot be rationally designed without large-scale experimental screening. Allosteric drugs have raised much attention due to their high specificity and possibility of overcoming existing drug-resistant mutations. However, optimization of allosteric compounds remains challenging. Here, we developed a novel computational method KeyAlloSite to predict allosteric site and to identify key allosteric residues (allo-residues) based on the evolutionary coupling model. We found that protein allosteric sites are strongly coupled to orthosteric site compared to non-functional sites. We further inferred key allo-residues by pairwise comparing the difference of evolutionary coupling scores of each residue in the allosteric pocket with the functional site. Our predicted key allo-residues are in accordance with previous experimental studies for typical allosteric proteins like BCR-ABL1, Tar, and PDZ3, as well as key cancer mutations. We also showed that KeyAlloSite can be used to predict key allosteric residues distant from the catalytic site that are important for enzyme catalysis. Our study demonstrates that weak coevolutionary couplings contain important information of protein allosteric regulation function. KeyAlloSite can be applied in studying the evolution of protein allosteric regulation, designing and optimizing allosteric drugs, and performing functional protein design and enzyme engineering.

Allosteric Signaling in PDZ Energetic Networks: Embedding Error Analysis

Allosteric signaling in proteins has been known for some half a century, yet how the signal traverses the protein remains an active area of research. Recently, the importance of electrostatics to achieve long-range signaling has become increasingly appreciated. Our laboratory has been working on developing network approaches to capture such interactions. In this study, we turn our attention to the well-studied allosteric model protein, PDZ. We study the allosteric dynamics on a per-residue basis in key constructs involving the PDZ domain, its allosteric effector, and its peptide ligand. We utilize molecular dynamics trajectories to create the networks for the constructs to explore the allosteric effect by plotting the heat kernel results onto axes defined by principal components. We introduce a new metric to quantitate the volume sampled by a residue in the latent space. We relate our findings to PDZ and the greater field of allostery.

Additive energetic contributions of multiple peptide positions determine the relative promiscuity of viral and human sequences for PDZ domain targets

Protein-protein interactions that include recognition of short sequences of amino acids, or peptides, are critical in cellular processes. Protein-peptide interaction surface areas are relatively small and shallow, and there are often overlapping specificities in families of peptide-binding domains. Therefore, dissecting selectivity determinants can be challenging. PDZ domains are an example of a peptide-binding domain located in several intracellular signaling and trafficking pathways, which form interactions critical for the regulation of receptor endocytic trafficking, tight junction formation, organization of supramolecular complexes in neurons, and other biological systems. These domains are also directly targeted by pathogens, and a hallmark of many oncogenic viral proteins is a PDZ-binding motif. However, amidst sequences that target PDZ domains, there is a wide spectrum in relative promiscuity. For example, the viral HPV16 E6 oncoprotein recognizes over double the number of PDZ domain-containing proteins as the cystic fibrosis transmembrane conductance regulator (CFTR) in the cell, despite similar PDZ targeting-sequences and identical motif residues. Here, we determine binding affinities for PDZ domains known to bind either HPV16 E6 alone or both CFTR and HPV16 E6, using peptides matching WT and hybrid sequences. We also use energy minimization to model PDZ-peptide complexes and use sequence analyses to investigate this difference. We find that while the majority of single mutations had a marginal effect on overall affinity, the additive effect on the free energy of binding accurately describes the selectivity observed. Taken together, our results describe how complex and differing PDZ interactomes can be programmed in the cell.

Investigating the allosteric response of the PICK1 PDZ domain to different ligands with all-atom simulations

The PDZ family is comprised of small modular domains that play critical roles in the allosteric modulation of many cellular signaling processes by binding to the C-terminal tail of different proteins. As dominant modular proteins that interact with a diverse set of peptides, it is of particular interest to explore how different binding partners induce different allosteric effects on the same PDZ domain. Because the PICK1 PDZ domain can bind different types of ligands, it is an ideal test case to answer this question and explore the network of interactions that give rise to dynamic allostery. Here, we use all-atom molecular dynamics simulations to explore dynamic allostery in the PICK1 PDZ domain by modeling two PICK1 PDZ systems: PICK1 PDZ-DAT and PICK1 PDZ-GluR2. Our results suggest that ligand binding to the PICK1 PDZ domain induces dynamic allostery at the αA helix that is similar to what has been observed in other PDZ domains. We found that the PICK1 PDZ-ligand distance is directly correlated with both dynamic changes of the αA helix and the distance between the αA helix and βB strand. Furthermore, our work identifies a hydrophobic core between DAT/GluR2 and I35 as a key interaction in inducing such dynamic allostery. Finally, the unique interaction patterns between different binding partners and the PICK1 PDZ domain can induce unique dynamic changes to the PICK1 PDZ domain. We suspect that unique allosteric coupling patterns with different ligands may play a critical role in how PICK1 performs its biological functions in various signaling networks.

Structural deficits in key domains of Shank2 lead to alterations in postsynaptic nanoclusters and to a neurodevelopmental disorder in humans

Postsynaptic scaffold proteins such as Shank, PSD-95, Homer and SAPAP/GKAP family members establish the postsynaptic density of glutamatergic synapses through a dense network of molecular interactions. Mutations in SHANK genes are associated with neurodevelopmental disorders including autism and intellectual disability. However, no SHANK missense mutations have been described which interfere with the key functions of Shank proteins believed to be central for synapse formation, such as GKAP binding via the PDZ domain, or Zn2+-dependent multimerization of the SAM domain. We identify two individuals with a neurodevelopmental disorder carrying de novo missense mutations in SHANK2. The p.G643R variant distorts the binding pocket for GKAP in the Shank2 PDZ domain and prevents interaction with Thr(-2) in the canonical PDZ ligand motif of GKAP. The p.L1800W variant severely delays the kinetics of Zn2+-dependent polymerization of the Shank2-SAM domain. Structural analysis shows that Trp1800 dislodges one histidine crucial for Zn2+ binding. The resulting conformational changes block the stacking of helical polymers of SAM domains into sheets through side-by-side contacts, which is a hallmark of Shank proteins, thereby disrupting the highly cooperative assembly process induced by Zn2+. Both variants reduce the postsynaptic targeting of Shank2 in primary cultured neurons and alter glutamatergic synaptic transmission. Super-resolution microscopy shows that both mutants interfere with the formation of postsynaptic nanoclusters. Our data indicate that both the PDZ- and the SAM-mediated interactions of Shank2 contribute to the compaction of postsynaptic protein complexes into nanoclusters, and that deficiencies in this process interfere with normal brain development in humans.

A short structural extension dictates the early stages of folding of a PDZ domain

PDZ domains are highly abundant protein-protein interaction modules in human. One of the most extensively characterized PDZ domain, the third PDZ domain from PSD-95 (PDZ3), contains an α-helical C-terminal extension that has a key role in the function of the domain. Here we compared the folding of PDZ3 with a truncated variant (PDZ3Δα3), lacking the additional helix, by means of the so-called Φ-value analysis, an experimental technique that allows inferring the structure of folding transition states. Experiments reveal subtle but detectable differences in the folding of PDZ3Δα3 versus PDZ3, as probed by structural characterization of the folding transition states. These differences appear more remarkable in the early stages of folding, with a detectable shift of the folding nucleus. The presented results allow demonstrating that the native state exerts a weak bias at the early stages of folding, which appear to be characterized by alternative pathways.

ProxyBind: a Compendium of Binding Sites for Proximity-Induced Pharmacology

Proximity-induced pharmacology (ProxPharm) is a novel paradigm in drug discovery where a small molecule brings two proteins in close proximity to elicit a signal, generally from one protein onto another. The potential of ProxPharm compounds as a new therapeutic modality is firmly established by proteolysis targeting chimeras (PROTACs) that bring an E3 ubiquitin ligase in proximity to a target protein to induce ubiquitination and subsequent degradation of the target. The concept can be expanded to induce other post-translational modifications via the recruitment of different types of protein-modifying enzymes. To survey the human proteome for opportunities in proximity pharmacology, we systematically mapped non-catalytic drug binding pockets on the structure of protein-modifying enzymes available from the Protein Databank. In addition to binding sites exploited by previously reported ProxPharm compounds, we identified putative ligandable non-catalytic pockets in 236 kinases, 45 phosphatases, 37 deubiquitinases, 14 methyltransferases, 11 acetyltransferases, 13 glycosyltransferases, 4 deacetylases, 7 demethylases and 2 glycosidases, including cavities occupied by chemical matter that may serve as starting points for future ProxPharm compounds. This systematic survey confirms that proximity pharmacology is a versatile modality with largely unexplored and promising potential and reveals novel opportunities to pharmacologically rewire molecular circuitries. All data is available from the ProxyBind database at https://polymorph.sgc.utoronto.ca/proxybind/index.php.

GCAT: A network model of mutational influences between amino acid positions in PSD95pdz3

Proteins exist for more than 3 billion years: proof of a sustainable design. They have mechanisms coping with internal perturbations (e.g., amino acid mutations), which tie genetic backgrounds to diseases or drug therapy failure. One difficulty to grasp these mechanisms is the asymmetry of amino acid mutational impact: a mutation at position i in the sequence, which impact a position j does not imply that the mutation at position j impacts the position i. Thus, to distinguish the influence of the mutation of i on j from the influence of the mutation of j on i, position mutational influences must be represented with directions. Using the X ray structure of the third PDZ domain of PDS-95 (Protein Data Bank 1BE9) and in silico mutations, we build a directed network called GCAT that models position mutational influences. In the GCAT, a position is a node with edges that leave the node (out-edges) for the influences of the mutation of the position on other positions and edges that enter the position (in-edges) for the influences of the mutation of other positions on the position. 1BE9 positions split into four influence categories called G, C, A and T going from positions influencing on average less other positions and influenced on average by less other positions (category C) to positions influencing on average more others positions and influenced on average by more other positions (category T). The four categories depict position neighborhoods in the protein structure with different tolerance to mutations.

Nonequilibrium Modeling of the Elementary Step in PDZ3 Allosteric Communication

While allostery is of paramount importance for protein signaling and regulation, the underlying dynamical process of allosteric communication is not well understood. The PDZ3 domain represents a prime example of an allosteric single-domain protein, as it features a well-established long-range coupling between the C-terminal α₃-helix and ligand binding. In an intriguing experiment, Hamm and co-workers employed photoswitching of the α₃-helix to initiate a conformational change of PDZ3 that propagates from the C-terminus to the bound ligand within 200 ns. Performing extensive nonequilibrium molecular dynamics simulations, the modeling of the experiment reproduces the measured time scales and reveals a detailed picture of the allosteric communication in PDZ3. In particular, a correlation analysis identifies a network of contacts connecting the α₃-helix and the core of the protein, which move in a concerted manner. Representing a one-step process and involving direct α₃-ligand contacts, this cooperative transition is considered as the elementary step in the propagation of conformational change.

Atomic resolution protein allostery from the multi-state structure of a PDZ domain

Recent methodological advances in solution NMR allow the determination of multi-state protein structures and provide insights into structurally and dynamically correlated protein sites at atomic resolution. This is demonstrated in the present work for the well-studied PDZ2 domain of protein human tyrosine phosphatase 1E for which protein allostery had been predicted. Two-state protein structures were calculated for both the free form and in complex with the RA-GEF2 peptide using the exact nuclear Overhauser effect (eNOE) method. In the apo protein, an allosteric conformational selection step comprising almost 60% of the domain was detected with an "open" ligand welcoming state and a "closed" state that obstructs the binding site by changing the distance between the β-sheet 2, α-helix 2, and sidechains of residues Lys38 and Lys72. The observed induced fit-type apo-holo structural rearrangements are in line with the previously published evolution-based analysis covering ~25% of the domain with only a partial overlap with the protein allostery of the open form. These presented structural studies highlight the presence of a dedicated highly optimized and complex dynamic interplay of the PDZ2 domain owed by the structure-dynamics landscape.

Intestinal brush border formation requires a TMIGD1-based intermicrovillar adhesion complex

Intestinal epithelial cells absorb nutrients through the brush border, composed of dense arrays of highly ordered microvilli at their apical membranes. A protocadherin-based intermicrovillar adhesion complex localized at microvilli tips mediates microvilli packing and organization. Here, we identified a second adhesion complex localized at the proximal base region of microvilli. This complex contained the immunoglobulin superfamily member TMIGD1, which directly interacted with the microvillar scaffolding proteins EBP50 and E3KARP. Complex formation with EBP50 required the activation of EBP50 by the actin-binding protein ezrin and was enhanced by the dephosphorylation of Ser¹⁶² in the PDZ2 domain of EBP50 by the phosphatase PP1α. Binding of the EBP50-ezrin complex to TMIGD1 enhanced the dynamic turnover of EBP50 at microvilli. Enterocyte-specific inactivation of Tmigd1 in mice resulted in microvillar blebbing, loss of intermicrovillar adhesion, and perturbed brush border formation. Thus, we identified a second adhesion complex in microvilli and propose a mechanism that promotes microvillar formation and dynamics.

Three Binding Conformations of BIO124 in the Pocket of the PICK1 PDZ Domain

The PDZ family has drawn attention as possible drug targets because of the domains’ wide ranges of function and highly conserved binding pockets. The PICK1 PDZ domain has been proposed as a possible drug target because the interactions between the PICK1 PDZ domain and the GluA2 subunit of the AMPA receptor have been shown to progress neurodegenerative diseases. BIO124 has been identified as a sub µM inhibitor of the PICK1–GluA2 interaction. Here, we use all-atom molecular dynamics simulations to reveal the atomic-level interaction pattern between the PICK1 PDZ domain and BIO124. Our simulations reveal three unique binding conformations of BIO124 in the PICK1 PDZ binding pocket, referred to here as state 0, state 1, and state 2. Each conformation is defined by a unique hydrogen bonding network and a unique pattern of hydrophobic interactions between BIO124 and the PICK1 PDZ domain. Interestingly, each conformation of BIO124 results in different dynamic changes to the PICK1 PDZ domain. Unlike states 1 and 2, state 0 induces dynamic coupling between BIO124 and the αA helix. Notably, this dynamic coupling with the αA helix is similar to what has been observed in other PDZ–ligand complexes. Our analysis indicates that the interactions formed between BIO124 and I35 may be the key to inducing dynamic coupling with the αA helix. Lastly, we suspect that the conformational shifts observed in our simulations may affect the stability and thus the overall effectiveness of BIO124. We propose that a physically larger inhibitor may be necessary to ensure sufficient interactions that permit stable binding between a drug and the PICK1 PDZ domain.

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.

Versatile Vibrational Energy Sensors for Proteins

Vibrational energy transfer (VET) is emerging as key mechanism for protein functions, possibly playing an important role for energy dissipation, allosteric regulation, and enzyme catalysis. A deep understanding of VET is required to elucidate its role in such processes. Ultrafast VIS-pump/IR-probe spectroscopy can detect pathways of VET in proteins. However, the requirement of having a VET donor and a VET sensor installed simultaneously limits the possible target proteins and sites; to increase their number we compare six IR labels regarding their utility as VET sensors. We compare these labels in terms of their FTIR, and VET signature in VET donor-sensor dipeptides in different solvents. Furthermore, we incorporated four of these labels in PDZ3 to assess their capabilities in more complex systems. Our results show that different IR labels can be used interchangeably, allowing for free choice of the right label depending on the system under investigation and the methods available.

Human Serine Racemase Weakly Binds the Third PDZ Domain of PSD-95

Human serine racemase (hSR) is a pyridoxal-5'-phosphate (PLP)-dependent dimer that catalyzes the formation of D-serine from L-serine, as well as the dehydration of both L- and D-serine to pyruvate and ammonia. As D-serine is a co-agonist of N-methyl-D-aspartate receptors (NMDARs), hSR is a key enzyme in glutamatergic neurotransmission. hSR activity is finely regulated by Mg²⁺, ATP, post-translational modifications, and the interaction with protein partners. In particular, the C-terminus of murine SR binds the third PDZ domain (PDZ3) of postsynaptic density protein 95 (PSD-95), a member of the membrane-associated guanylate kinase (MAGUK) family involved in the trafficking and localization of glutamate receptors. The structural details of the interaction and the stability of the complex have not been elucidated yet. We evaluated the binding of recombinant human PSD-95 PDZ3 to hSR by glutaraldehyde cross-linking, pull-down assays, isothermal titration calorimetry, nuclear magnetic resonance, and enzymatic assays. Overall, a weak interaction was observed, confirming the binding for the human orthologs but supporting the hypothesis that a third protein partner (i.e., stargazin) is required for the regulation of hSR activity by PSD-95 and to stabilize their interaction.

The novel human HtrA2 ortholog in zebrafish: New molecular insight and challenges into the imbalance of homeostasis

High temperature requirement A2 (HtrA2) contributes to regulating mitochondrial quality control and maintaining the balance between the death and survival of cells and living organisms. However, the molecular mechanism of HtrA2 in physiological and pathophysiological processes remains unclear. HtrA2 exhibits multifaceted characteristics according to the expression levels and acts opposite functions depending on its subcellular localization. Thus, innovative technologies and systems that can be freely manipulated at the quantitative, biochemical, molecular and cellular levels are needed to address not only the challenges faced by HtrA2 research but also the general obstacles to protein research. Here, we are the first to identify zebrafish HtrA2 (zHtrA2) as the true ortholog of human HtrA2 (hHtrA2), by in silico sequence analysis of genomic DNA and molecular biological techniques, which is highly conserved structurally and functionally as a serine protease and cell death regulator. The zHtrA2 protein is primarily localized in the mitochondria, where alanine-exposed mature zHtrA2 ((A)-zHtrA2) is generated by removing 111 residues at the N-terminus of pro-zHtrA2. The (A)-zHtrA2 released from the mitochondria into the cytosol induces the caspase cascade by binding to and inhibiting hXIAP, a cognate partner of hHtrA2. Notably, zHtrA2 has well conserved properties of serine protease that specifically cleaves hParkin, a cognate substrate of hHtrA2. Interestingly, cytosolic (M)-zHtrA2, which does not bind hXIAP, induces atypical cell death in a serine protease-dependent manner, as occurs in hHtrA2. Thus, the zebrafish-zHtrA2 system can be used to clarify the crucial role of HtrA2 in maintaining the survival of living organisms and provide an opportunity to develop novel therapeutics for HtrA2-associated diseases, such as neurodegenerative diseases and cancer, which are caused by dysregulation of HtrA2.

Structure-mechanics statistical learning uncovers mechanical relay in proteins

A protein's adaptive response to its substrates is one of the key questions driving molecular physics and physical chemistry. This work employs the recently developed structure-mechanics statistical learning method to establish a mechanical perspective. Specifically, by mapping all-atom molecular dynamics simulations onto the spring parameters of a backbone-side-chain elastic network model, the chemical moiety specific force constants (or mechanical rigidity) are used to assemble the rigidity graph, which is the matrix of inter-residue coupling strength. Using the S1A protease and the PDZ3 signaling domain as examples, chains of spatially contiguous residues are found to exhibit prominent changes in their mechanical rigidity upon substrate binding or dissociation. Such a mechanical-relay picture thus provides a mechanistic underpinning for conformational changes, long-range communication, and inter-domain allostery in both proteins, where the responsive mechanical hotspots are mostly residues having important biological functions or significant mutation sensitivity.

Protein residues exhibit specific routes of mechanical relay as the adaptive responses to substrate binding or dissociation. On such physically contiguous connections, residues experience prominent changes in their coupling strengths.

Efficient Unnatural Protein Production by Pyrrolysyl-tRNA Synthetase With Genetically Fused Solubility Tags

Introducing non-canonical amino acids (ncAAs) by engineered orthogonal pairs of aminoacyl-tRNA synthetases and tRNAs has proven to be a highly useful tool for the expansion of the genetic code. Pyrrolysyl-tRNA synthetase (PylRS) from methanogenic archaeal and bacterial species is particularly attractive due to its natural orthogonal reactivity in bacterial and eukaryotic cells. However, the scope of such a reprogrammed translation is often limited, due to low yields of chemically modified target protein. This can be the result of substrate specificity engineering, which decreases the aminoacyl-tRNA synthetase stability and reduces the abundance of active enzyme. We show that the solubility and folding of these engineered enzymes can become a bottleneck for the production of ncAA-containing proteins in vivo. Solubility tags derived from various species provide a strategy to remedy this issue. We find the N-terminal fusion of the small metal binding protein from Nitrosomonas europaea to the PylRS sequence to improve enzyme solubility and to boost orthogonal translation efficiency. Our strategy enhances the production of site-specifically labelled proteins with a variety of engineered PylRS variants by 200–540%, and further allows triple labeling. Even the wild-type enzyme gains up to 245% efficiency for established ncAA substrates.

A scalable strategy to solve structures of PDZ domains and their complexes

The human PDZome represents one of the largest globular domain families in the human proteome, with 266 instances. These globular domains typically interact with C-terminal peptide motifs found in thousands of human proteins. Despite previous efforts, not all PDZ domains have experimentally solved structures and most of their complexes remain to be solved. Here, a simple and cost-effective strategy is proposed for the crystallization of PDZ domains and their complexes. A human annexin A2 fusion tag was used as a crystallization chaperone and the structures of nine PDZ domains were solved, including five domains that had not yet been solved. Finally, these novel experimental structures were compared with AlphaFold predictions and it is speculated how predictions and experimental methods could cooperate in order to investigate the structural landscapes of entire domain families and interactomes.

Allosterism in the PDZ Family

Dynamic allosterism allows the propagation of signal throughout a protein. The PDZ (PSD-95/Dlg1/ZO-1) family has been named as a classic example of dynamic allostery in small modular domains. While the PDZ family consists of more than 200 domains, previous efforts have primarily focused on a few well-studied PDZ domains, including PTP-BL PDZ2, PSD-95 PDZ3, and Par6 PDZ. Taken together, experimental and computational studies have identified regions of these domains that are dynamically coupled to ligand binding. These regions include the αA helix, the αB lower-loop, and the αC helix. In this review, we summarize the specific residues on the αA helix, the αB lower-loop, and the αC helix of PTP-BL PDZ2, PSD-95 PDZ3, and Par6 PDZ that have been identified as participants in dynamic allostery by either experimental or computational approaches. This review can serve as an index for researchers to look back on the previously identified allostery in the PDZ family. Interestingly, our summary of previous work reveals clear consistencies between the domains. While the PDZ family has a low sequence identity, we show that some of the most consistently identified allosteric residues within PTP-BL PDZ2 and PSD-95 PDZ3 domains are evolutionarily conserved. These residues include A46/A347, V61/V362, and L66/L367 on PTP-BL PDZ2 and PSD-95 PDZ3, respectively. Finally, we expose a need for future work to explore dynamic allostery within (1) PDZ domains with multiple binding partners and (2) multidomain constructs containing a PDZ domain.

PSNtools for standalone and web-based structure network analyses of conformational ensembles

Structure graphs, in which interacting amino acids/nucleotides correspond to linked nodes, represent cutting-edge tools to investigate macromolecular function.

The graph-based approach defined as Protein Structure Network (PSN) was initially implemented in the Wordom software and subsequently in the webPSN server. PSNs are computed either on a molecular dynamics (MD) trajectory (PSN-MD) or on a single structure. In the latter case, information on atomic fluctuations is inferred from the Elastic Network Model-Normal Mode Analysis (ENM-NMA) (PSN-ENM). While Wordom performs both PSN-ENM and PSN-MD analyses but without output post-processing, the webPSN server performs only single-structure PSN-EMN but assisting the user in input setup and output analysis.

Here we release for the first time the standalone software PSNtools, which allows calculation and post-processing of PSN analyses carried out either on single structures or on conformational ensembles. Relevant unique and novel features of PSNtools are either comparisons of two networks or computations of consensus networks on sets of homologous/analogous macromolecular structures or conformational ensembles. Network comparisons and consensus serve to infer differences in functionally different states of the same system or network-based signatures in groups of bio-macromolecules sharing either the same functionality or the same fold.

In addition to the new software, here we release also an updated version of the webPSN server, which allows performing an interactive graphical analysis of PSN-MD, following the upload of the PSNtools output.

PSNtools, the auxiliary binary version of Wordom software, and the WebPSN server are freely available at http://webpsn.hpc.unimo.it/wpsn3.php.

Using azobenzene photocontrol to set proteins in motion

Controlling the activity of proteins with azobenzene photoswitches is a potent tool for manipulating their biological function. With the help of light, it is possible to change binding affinities, control allostery or manipulate complex biological processes, for example. Additionally, owing to their intrinsically fast photoisomerization, azobenzene photoswitches can serve as triggers that initiate out-of-equilibrium processes. Such switching of the activity initiates a cascade of conformational events that can be accessed with time-resolved methods. In this Review, we show how the potency of azobenzene photoswitching can be combined with transient spectroscopic techniques to disclose the order of events and experimentally observe biomolecular interactions in real time. This strategy will further our understanding of how a protein can accommodate, adapt and readjust its structure to answer an incoming signal, revealing more of the dynamical character of proteins.

Site-specific phosphorylation of PSD-95 dynamically regulates the postsynaptic density as observed by phase separation

Postsynaptic density protein 95 is a key scaffolding protein in the postsynaptic density of excitatory glutamatergic neurons, organizing signaling complexes primarily via its three PSD-95/Discs-large/Zona occludens domains. PSD-95 is regulated by phosphorylation, but technical challenges have limited studies of the molecular details. Here, we genetically introduced site-specific phosphorylations in single, tandem, and full-length PSD-95 and generated a total of 11 phosphorylated protein variants. We examined how these phosphorylations affected binding to known interaction partners and the impact on phase separation of PSD-95 complexes and identified two new phosphorylation sites with opposing effects. Phosphorylation of Ser78 inhibited phase separation with the glutamate receptor subunit GluN2B and the auxiliary protein stargazin, whereas phosphorylation of Ser116 induced phase separation with stargazin only. Thus, by genetically introducing phosphoserine site-specifically and exploring the impact on phase separation, we have provided new insights into the regulation of PSD-95 by phosphorylation and the dynamics of the PSD.

Mechanical couplings of protein backbone and side chains exhibit scale-free network properties and specific hotspots for function

•

Statistical learning from protein dynamics unravels rigidities in interaction network.

•

Backbone and side-chain mechanical couplings exhibit scale-free network properties.

•

Graphical depiction of network rigidities captures sequence co-evolution patterns.

•

Functional sites at secondary structure peripheries are mechanical hotspots.

•

Our rigidity scores are compelling metrics for residue biological significance.

A backbone-side-chain elastic network model (bsENM) is devised in this contribution to decipher the network of molecular interactions during protein dynamics. The chemical details in 5 μs all-atom molecular dynamics (MD) simulation are mapped onto the bsENM spring constants by self-consistent iterations. The elastic parameters obtained by this structure-mechanics statistical learning are then used to construct inter-residue rigidity graphs for the chemical components in protein amino acids. A key discovery is that the mechanical coupling strengths of both backbone and side chains exhibit heavy-tailed distributions and scale-free network properties. In both rat trypsin and PDZ3 proteins, the statistically prominent modes of rigidity graphs uncover the sequence-specific coupling patterns and mechanical hotspots. Based on the contributions to graphical modes, our residue rigidity scores in backbone and side chains are found to be very useful metrics for the biological significance. Most functional sites have high residue rigidity scores in side chains while the biologically important glycines are generally next to mechanical hotspots. Furthermore, prominent modes in the rigidity graphs involving side chains oftentimes coincide with the co-evolution patterns due to evolutionary restraints. The bsENM specifically devised to resolve the protein chemical character thus provides useful means for extracting functional information from all-atom MD.

The biology of papillomavirus PDZ associations: what do they offer papillomaviruses?

The high-risk α-type papillomaviruses have a C-terminal PDZ-binding motif (PBM) on one of the two major oncoproteins E6 or E7; the vast majority on E6. The PBM is essential for the high-risk HPV life cycle, for episomal maintenance of the virus genome, and for maintaining the mitotic stability of the infected cell. The question is why only these viruses have PBMs - are there specific constraints imposed by the mucosal epithelium in which these viruses replicate? However the low-risk α-HPVs, such as HPV-6 and HPV-11 replicate extremely efficiently without a PBM, while viruses of the alpha8 group, such as HPV-40, replicate well with a very primitive PBM. So what does PDZ-binding capacity contribute to the fitness of the virus?

Targeting PDZ domains as potential treatment for viral infections, neurodegeneration and cancer

The interaction between proteins is a fundamental event for cellular life that is generally mediated by specialized protein domains or modules. PDZ domains are the largest class of protein-protein interaction modules, involved in several cellular pathways such as signal transduction, cell-cell junctions, cell polarity and adhesion, and protein trafficking. Because of that, dysregulation of PDZ domain function often causes the onset of pathologies, thus making this family of domains an interesting pharmaceutical target. In this review article we provide an overview of the structural and functional features of PDZ domains and their involvement in the cellular and molecular pathways at the basis of different human pathologies. We also discuss some of the strategies that have been developed with the final goal to hijack or inhibit the interaction of PDZ domains with their ligands. Because of the generally low binding selectivity of PDZ domain and the scarce efficiency of small molecules in inhibiting PDZ binding, this task resulted particularly difficult to pursue and still demands increasing experimental efforts in order to become completely feasible and successful in vivo.

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.