Interaction prediction and classification of PDZ domains

E4orf1: The triple agent of adenovirus - Unraveling its roles in oncogenesis, infectious obesity and immune responses in virus replication and vector therapy

Human Adenoviruses (HAdV) are nearly ubiquitous pathogens comprising numerous sub-types that infect various tissues and organs. Among many encoded proteins that facilitate viral replication and subversion of host cellular processes, the viral E4orf1 protein has emerged as an intriguing yet under-investigated player in the complex interplay between the virus and its host. E4orf1 has gained attention as a metabolism activator and oncogenic agent, while recent research is showing that E4orf1 may play a more important role in modulating cellular pathways such as PI3K-Akt-mTOR, Ras, the immune response and further HAdV replication stages than previously anticipated. In this review, we aim to explore the structure, molecular mechanisms, and biological functions of E4orf1, shedding light on its potentially multifaceted roles during HAdV infection, including metabolic diseases and oncogenesis. Furthermore, we discuss the role of functional E4orf1 in biotechnological applications such as Adenovirus (AdV) vaccine vectors and oncolytic AdV. By dissecting the intricate relationships between HAdV types and E4orf1 proteins, this review provides valuable insights into viral pathogenesis and points to promising areas of future research.

PDZ Peptide of the ZO-1 Protein Significantly Increases UTP-Induced MUC8 Anti-Inflammatory Mucin Overproduction in Human Airway Epithelial Cells

Mucus hyperproduction and hypersecretion are observed often in respiratory diseases. MUC8 is a glycoprotein synthesized by epithelial cells and generally expressed in the respiratory track. However, the physiological mechanism by which extracellular nucleotides induce MUC8 gene expression in human airway epithelial cells is unclear. Here, we show that UTP could induce MUC8 gene expression through P2Y2-PLCβ3-Ca2+ activation. Because the full-length cDNA sequence of MUC8 has not been identified, a specific siRNA-MUC8 was designed based on the partial cDNA sequence of MUC8. siRNA-MUC8 significantly increased TNF-α production and decreased IL-1Ra production, suggesting that MUC8 may downregulate UTP/P2Y2-induced airway inflammation. Interestingly, the PDZ peptide of ZO-1 protein strongly abolished UTP-induced TNF-α production and increased IL-1Ra production and MUC8 gene expression. In addition, the PDZ peptide dramatically increased the levels of UTP-induced ZO proteins and TEER (trans-epithelial electrical resistance). These results show that the anti-inflammatory mucin MUC8 may contribute to homeostasis, and the PDZ peptide can be a novel therapeutic candidate for UTP-induced airway inflammation.

Claudin Barriers on the Brink: How Conflicting Tissue and Cellular Priorities Drive IBD Pathogenesis

Inflammatory bowel diseases (IBDs) are characterized by acute or chronic recurring inflammation of the intestinal mucosa, often with increasing severity over time. Life-long morbidities and diminishing quality of life for IBD patients compel a search for a better understanding of the molecular contributors to disease progression. One unifying feature of IBDs is the failure of the gut to form an effective barrier, a core role for intercellular complexes called tight junctions. In this review, the claudin family of tight junction proteins are discussed as they are a fundamental component of intestinal barriers. Importantly, claudin expression and/or protein localization is altered in IBD, leading to the supposition that intestinal barrier dysfunction exacerbates immune hyperactivity and disease. Claudins are a large family of transmembrane structural proteins that constrain the passage of ions, water, or substances between cells. However, growing evidence suggests non-canonical claudin functions during mucosal homeostasis and healing after injury. Therefore, whether claudins participate in adaptive or pathological IBD responses remains an open question. By reviewing current studies, the possibility is assessed that with claudins, a jack-of-all-trades is master of none. Potentially, a robust claudin barrier and wound restitution involve conflicting biophysical phenomena, exposing barrier vulnerabilities and a tissue-wide frailty during healing in IBD.

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.

Interactions of the protein tyrosine phosphatase PTPN3 with viral and cellular partners through its PDZ domain: insights into structural determinants and phosphatase activity

The human protein tyrosine phosphatase non-receptor type 3 (PTPN3) is a phosphatase containing a PDZ (PSD-95/Dlg/ZO-1) domain that has been found to play both tumor-suppressive and tumor-promoting roles in various cancers, despite limited knowledge of its cellular partners and signaling functions. Notably, the high-risk genital human papillomavirus (HPV) types 16 and 18 and the hepatitis B virus (HBV) target the PDZ domain of PTPN3 through PDZ-binding motifs (PBMs) in their E6 and HBc proteins respectively. This study focuses on the interactions between the PTPN3 PDZ domain (PTPN3-PDZ) and PBMs of viral and cellular protein partners. We solved the X-ray structures of complexes between PTPN3-PDZ and PBMs of E6 of HPV18 and the tumor necrosis factor-alpha converting enzyme (TACE). We provide new insights into key structural determinants of PBM recognition by PTPN3 by screening the selectivity of PTPN3-PDZ recognition of PBMs, and by comparing the PDZome binding profiles of PTPN3-recognized PBMs and the interactome of PTPN3-PDZ. The PDZ domain of PTPN3 was known to auto-inhibit the protein's phosphatase activity. We discovered that the linker connecting the PDZ and phosphatase domains is involved in this inhibition, and that the binding of PBMs does not impact this catalytic regulation. Overall, the study sheds light on the interactions and structural determinants of PTPN3 with its cellular and viral partners, as well as on the inhibitory role of its PDZ domain on its phosphatase activity.

The unexpected versatility of ALP/Enigma family proteins

One of the most intriguing features of multicellular animals is their ability to move. On a cellular level, this is accomplished by the rearrangement and reorganization of the cytoskeleton, a dynamic network of filamentous proteins which provides stability and structure in a stationary context, but also facilitates directed movement by contracting. The ALP/Enigma family proteins are a diverse group of docking proteins found in numerous cellular milieus and facilitate these processes among others. In vertebrates, they are characterized by having a PDZ domain in combination with one or three LIM domains. The family is comprised of CLP-36 (PDLIM1), Mystique (PDLIM2), ALP (PDLIM3), RIL (PDLIM4), ENH (PDLIM5), ZASP (PDLIM6), and Enigma (PDLIM7). In this review, we will outline the evolution and function of their protein domains which confers their versatility. Additionally, we highlight their role in different cellular environments, focusing specifically on recent advances in muscle research using Drosophila as a model organism. Finally, we show the relevance of this protein family to human myopathies and the development of muscle-related diseases.

CRISPR/Cas9-based genome-wide screening for deubiquitinase subfamily identifies USP1 regulating MAST1-driven cisplatin-resistance in cancer cells

Background: Cisplatin is one of the frontline anticancer agents. However, development of cisplatin-resistance limits the therapeutic efficacy of cisplatin-based treatment. The expression of microtubule-associated serine/threonine kinase 1 (MAST1) is a primary factor driving cisplatin-resistance in cancers by rewiring the MEK pathway. However, the mechanisms responsible for MAST1 regulation in conferring drug resistance is unknown.

Methods: We implemented a CRISPR/Cas9-based, genome-wide, dual screening system to identify deubiquitinating enzymes (DUBs) that govern cisplatin resistance and regulate MAST1 protein level. We analyzed K48- and K63-linked polyubiquitination of MAST1 protein and mapped the interacting domain between USP1 and MAST1 by immunoprecipitation assay. The deubiquitinating effect of USP1 on MAST1 protein was validated using rescue experiments, in vitro deubiquitination assay, immunoprecipitation assays, and half-life analysis. Furthermore, USP1-knockout A549 lung cancer cells were generated to validate the deubiquitinating activity of USP1 on MAST1 abundance. The USP1-MAST1 correlation was evaluated using bioinformatics tool and in different human clinical tissues. The potential role of USP1 in regulating MAST1-mediated cisplatin resistance was confirmed using a series of in vitro and in vivo experiments. Finally, the clinical relevance of the USP1-MAST1 axis was validated by application of small-molecule inhibitors in a lung cancer xenograft model in NSG mice.

Results: The CRISPR/Cas9-based dual screening system identified USP1 as a novel deubiquitinase that interacts, stabilizes, and extends the half-life of MAST1 by preventing its K48-linked polyubiquitination. The expression analysis across human clinical tissues revealed a positive correlation between USP1 and MAST1. USP1 promotes MAST1-mediated MEK1 activation as an underlying mechanism that contributes to cisplatin-resistance in cancers. Loss of USP1 led to attenuation of MAST1-mediated cisplatin-resistance both in vitro and in vivo. The combined pharmacological inhibition of USP1 and MAST1 using small-molecule inhibitors further abrogated MAST1 level and synergistically enhanced cisplatin efficacy in a mouse xenograft model.

Conclusions: Overall, our study highlights the role of USP1 in the development of cisplatin resistance and uncovers the regulatory mechanism of MAST1-mediated cisplatin resistance in cancers. Co-treatment with USP1 and MAST1 inhibitors abrogated tumor growth and synergistically enhanced cisplatin efficacy, suggesting a novel alternative combinatorial therapeutic strategy that could further improve MAST1-based therapy in patients with cisplatin-resistant tumors.

Sparse group selection and analysis of function-related residue for protein-state recognition

Machine learning methods have helped to advance wide range of scientific and technological field in recent years, including computational chemistry. As the chemical systems could become complex with high dimension, feature selection could be critical but challenging to develop reliable machine learning based prediction models, especially for proteins as bio-macromolecules. In this study, we applied sparse group lasso (SGL) method as a general feature selection method to develop classification model for an allosteric protein in different functional states. This results into a much improved model with comparable accuracy (Acc) and only 28 selected features comparing to 289 selected features from a previous study. The Acc achieves 91.50% with 1936 selected feature, which is far higher than that of baseline methods. In addition, grouping protein amino acids into secondary structures provides additional interpretability of the selected features. The selected features are verified as associated with key allosteric residues through comparison with both experimental and computational works about the model protein, and demonstrate the effectiveness and necessity of applying rigorous feature selection and evaluation methods on complex chemical systems.

Molecular basis of the interaction of the human tyrosine phosphatase PTPN3 with the hepatitis B virus core protein

Interactions between the hepatitis B virus core protein (HBc) and host cell proteins are poorly understood, although they may be essential for the propagation of the virus and its pathogenicity. HBc has a C-terminal PDZ (PSD-95, Dlg1, ZO-1)-binding motif (PBM) that is responsible for interactions with host PDZ domain-containing proteins. In this work, we focused on the human protein tyrosine phosphatase non-receptor type 3 (PTPN3) and its interaction with HBc. We solved the crystal structure of the PDZ domain of PTPN3 in complex with the PBM of HBc, revealing a network of interactions specific to class I PDZ domains despite the presence of a C-terminal cysteine in this atypical PBM. We further showed that PTPN3 binds the HBc protein within capsids or as a homodimer. We demonstrate that overexpression of PTPN3 significantly affects HBV infection in HepG2 NTCP cells. Finally, we performed proteomics studies on both sides by pull-down assays and screening of a human PDZ domain library. We identified a pool of human PBM-containing proteins that might interact with PTPN3 in cells and that could be in competition with the HBc PBM during infection, and we also identified potential cellular partners of HBc through PDZ-PBM interactions. This study opens up many avenues of future investigations into the pathophysiology of HBV.

Advances in aquaporin-2 trafficking mechanisms and their implications for treatment of water balance disorders

In mammals, conservation of body water is critical for survival and is dependent on the kidneys' ability to minimize water loss in the urine during periods of water deprivation. The collecting duct water channel aquaporin-2 (AQP2) plays an essential role in this homeostatic response by facilitating water reabsorption along osmotic gradients. The ability to increase the levels of AQP2 in the apical plasma membrane following an increase in plasma osmolality is a rate-limiting step in water reabsorption, a process that is tightly regulated by the antidiuretic hormone arginine vasopressin (AVP). In this review, the focus is on the role of the carboxyl-terminus of AQP2 as a key regulatory point for AQP2 trafficking. We provide an overview of AQP2 structure, disease-causing mutations in the AQP2 carboxyl-terminus, the role of posttranslational modifications such as phosphorylation and ubiquitylation in the tail domain, and their implications for balanced trafficking of AQP2. Finally, we discuss how various modifications of the AQP2 tail facilitate selective protein-protein interactions that modulate the AQP2 trafficking mechanism.

A Novel Heterozygous Missense Variant (c.667G>T;p.Gly223Cys) in USH1C That Interferes With Cadherin-Related 23 and Harmonin Interaction Causes Autosomal Dominant Nonsyndromic Hearing Loss

Background

Pathogenic variants of USH1C, encoding a PDZ-domain-containing protein called harmonin, have been known to cause autosomal recessive syndromic or nonsyndromic hearing loss (NSHL). We identified a causative gene in a large Korean family with NSHL showing a typical pattern of autosomal dominant (AD) inheritance.

Methods

Exome sequencing was performed for five affected and three unaffected individuals in this family. Following identification of a candidate gene variant, segregation analysis and functional studies, including circular dichroism and biolayer interferometry experiments, were performed.

Results

A novel USH1C heterozygous missense variant (c.667G>T;p.Gly223Cys) was shown to segregate with the NSHL phenotype in this family. This variant affects an amino acid residue located in the highly conserved carboxylate-binding loop of the harmonin PDZ2 domain and is predicted to disturb the interaction with cadherin-related 23 (cdh23). The affinity of the variant PDZ2 domain for a biotinylated synthetic peptide containing the PDZ-binding motif of cdh23 was approximately 16-fold lower than that of the wild-type PDZ2 domain and that this inaccessibility of the binding site was caused by a conformational change in the variant PDZ2 domain.

Conclusions

A heterozygous variant of USH1C that interferes with the interaction between cdh23 and harmonin causes novel AD-NSHL.

Multifaceted Regulation of PTEN Subcellular Distributions and Biological Functions

Phosphatase and tensin homolog deleted on chromosome 10 (PTEN) is a tumor suppressor gene frequently found to be inactivated in over 30% of human cancers. PTEN encodes a 54-kDa lipid phosphatase that serves as a gatekeeper of the phosphoinositide 3-kinase pathway involved in the promotion of multiple pro-tumorigenic phenotypes. Although the PTEN protein plays a pivotal role in carcinogenesis, cumulative evidence has implicated it as a key signaling molecule in several other diseases as well, such as diabetes, Alzheimer's disease, and autism spectrum disorders. This finding suggests that diverse cell types, especially differentiated cells, express PTEN. At the cellular level, PTEN is widely distributed in all subcellular compartments and organelles. Surprisingly, the cytoplasmic compartment, not the plasma membrane, is the predominant subcellular location of PTEN. More recently, the finding of a secreted 'long' isoform of PTEN and the presence of PTEN in the cell nucleus further revealed unexpected biological functions of this multifaceted molecule. At the regulatory level, PTEN activity, stability, and subcellular distribution are modulated by a fascinating array of post-translational modification events, including phosphorylation, ubiquitination, and sumoylation. Dysregulation of these regulatory mechanisms has been observed in various human diseases. In this review, we provide an up-to-date overview of the knowledge gained in the last decade on how different functional domains of PTEN regulate its biological functions, with special emphasis on its subcellular distribution. This review also highlights the findings of published studies that have reported how mutational alterations in specific PTEN domains can lead to pathogenesis in humans.

MDA-9/Syntenin: An emerging global molecular target regulating cancer invasion and metastasis

With few exceptions, metastasis is the terminal stage of cancer with limited therapeutic options. Metastasis consists of numerous phenotypic and genotypic alterations of cells that are directly and indirectly induced by multiple intrinsic (cellular) and extrinsic (micro-environmental) factors. To metastasize, a cancer cell often transitions from an epithelial to mesenchymal morphology (EMT), modifies the extracellular matrix, forms emboli and survives in the circulation, escapes immune surveillance, adheres to sites distant from the initial tumor and finally develops a blood supply (angiogenesis) and colonizes in a secondary niche (a micrometastasis). Scientific advances have greatly enhanced our understanding of the precise molecular and genetic changes, operating independently or collectively, that lead to metastasis. This review focuses on a unique gene, melanoma differentiation associated gene-9 (also known as Syntenin-1; Syndecan Binding Protein (sdcbp); mda-9/syntenin), initially cloned and characterized from metastatic human melanoma and shown to be a pro-metastatic gene. In the last two decades, our comprehension of the diversity of actions of MDA-9/Syntenin on cellular phenotype has emerged. MDA-9/Sytenin plays pivotal regulatory roles in multiple signaling cascades and orchestrates both metastatic and non-metastatic events. Considering the relevance of this gene in controlling cancer invasion and metastasis, approaches have been developed to uniquely and selectively target this gene. We also provide recent updates on strategies that have been successfully employed in targeting MDA-9/Syntenin resulting in profound pre-clinical anti-cancer activity.

Molecular Basis of Class III Ligand Recognition by PDZ3 in Murine Protein Tyrosine Phosphatase PTPN13

Protein tyrosine phosphatase PTPN13, also known as PTP-BL in mice, represents a large multi-domain non-transmembrane scaffolding protein that contains five consecutive PDZ domains. Here, we report the solution structures of the extended murine PTPN13 PDZ3 domain in its apo form and in complex with its physiological ligand, the carboxy-terminus of protein kinase C-related kinase-2 (PRK2), determined by multidimensional NMR spectroscopy. Both in its ligand-free state and when complexed to PRK2, PDZ3 of PTPN13 adopts the classical compact, globular D/E fold. PDZ3 of PTPN13 binds five carboxy-terminal amino acids of PRK2 via a groove located between the EB-strand and the DB-helix. The PRK2 peptide resides in the canonical PDZ3 binding cleft in an elongated manner and the amino acid side chains in position P0 and P-2, cysteine and aspartate, of the ligand face the groove between EB-strand and DB-helix, whereas the PRK2 side chains of tryptophan and alanine located in position P-1 and P-3 point away from the binding cleft. These structures are rare examples of selective class III ligand recognition by a PDZ domain and now provide a basis for the detailed structural investigation of the promiscuous interaction between the PDZ domains of PTPN13 and their ligands. They will also lead to a better understanding of the proposed scaffolding function of these domains in multi-protein complexes assembled by PTPN13 and could ultimately contribute to low molecular weight antagonists that might even act on the PRK2 signaling pathway to modulate rearrangements of the actin cytoskeleton.

Lock and chop: A novel method for the generation of a PICK1 PDZ domain and piperidine-based inhibitor co-crystal structure

The membrane protein interacting with kinase C1 (PICK1) plays a trafficking role in the internalization of neuron receptors such as the amino-3-hydroxyl-5-methyl-4-isoxazole-propionate (AMPA) receptor. Reduction of surface AMPA type receptors on neurons reduces synaptic communication leading to cognitive impairment in progressive neurodegenerative diseases such as Alzheimer disease. The internalization of AMPA receptors is mediated by the PDZ domain of PICK1 which binds to the GluA2 subunit of AMPA receptors and targets the receptor for internalization through endocytosis, reducing synaptic communication. We planned to block the PICK1-GluA2 protein-protein interaction with a small molecule inhibitor to stabilize surface AMPA receptors as a therapeutic possibility for neurodegenerative diseases. Using a fluorescence polarization assay, we identified compound BIO124 as a modest inhibitor of the PICK1-GluA2 interaction. We further tried to improve the binding affinity of BIO124 using structure-aided drug design but were unsuccessful in producing a co-crystal structure using previously reported crystallography methods for PICK1. Here, we present a novel method through which we generated a co-crystal structure of the PDZ domain of PICK1 bound to BIO124.

ZO-2, a tight junction protein involved in gene expression, proliferation, apoptosis, and cell size regulation

ZO-2 is a peripheral tight junction protein that belongs to the membrane-associated guanylate kinase protein family. Here, we explain the modular and supramodular organization of ZO-2 that allows it to interact with a wide variety of molecules, including cell-cell adhesion proteins, cytoskeletal components, and nuclear factors. We also describe how ZO proteins evolved through metazoan evolution and analyze the intracellular traffic of ZO-2, as well as the roles played by ZO-2 at the plasma membrane and nucleus that translate into the regulation of proliferation, cell size, and apoptosis. In addition, we focus on the impact of ZO-2 expression on male fertility and on maladies like cancer, cholestasis, and hearing loss.

Protein Engineering by Combined Computational and In Vitro Evolution Approaches

Two alternative strategies are commonly used to study protein-protein interactions (PPIs) and to engineer protein-based inhibitors. In one approach, binders are selected experimentally from combinatorial libraries of protein mutants that are displayed on a cell surface. In the other approach, computational modeling is used to explore an astronomically large number of protein sequences to select a small number of sequences for experimental testing. While both approaches have some limitations, their combination produces superior results in various protein engineering applications. Such applications include the design of novel binders and inhibitors, the enhancement of affinity and specificity, and the mapping of binding epitopes. The combination of these approaches also aids in the understanding of the specificity profiles of various PPIs.

Application of wavelet transform for PDZ domain classification

PDZ domains have been identified as part of an array of signaling proteins that are often unrelated, except for the well-conserved structural PDZ domain they contain. These domains have been linked to many disease processes including common Avian influenza, as well as very rare conditions such as Fraser and Usher syndromes. Historically, based on the interactions and the nature of bonds they form, PDZ domains have most often been classified into one of three classes (class I, class II and others - class III), that is directly dependent on their binding partner. In this study, we report on three unique feature extraction approaches based on the bigram and trigram occurrence and existence rearrangements within the domain's primary amino acid sequences in assisting PDZ domain classification. Wavelet packet transform (WPT) and Shannon entropy denoted by wavelet entropy (WE) feature extraction methods were proposed. Using 115 unique human and mouse PDZ domains, the existence rearrangement approach yielded a high recognition rate (78.34%), which outperformed our occurrence rearrangements based method. The recognition rate was (81.41%) with validation technique. The method reported for PDZ domain classification from primary sequences proved to be an encouraging approach for obtaining consistent classification results. We anticipate that by increasing the database size, we can further improve feature extraction and correct classification.

The PDZ domain protein Mcc is a novel effector of non-canonical Wnt signaling during convergence and extension in zebrafish

During vertebrate gastrulation, a complex set of mass cellular rearrangements shapes the embryonic body plan and appropriately positions the organ primordia. In zebrafish and Xenopus, convergence and extension (CE) movements simultaneously narrow the body axis mediolaterally and elongate it from head to tail. This process is governed by polarized cell behaviors that are coordinated by components of the non-canonical, β-catenin-independent Wnt signaling pathway, including Wnt5b and the transmembrane planar cell polarity (PCP) protein Vangl2. However, the intracellular events downstream of Wnt/PCP signals are not fully understood. Here, we show that zebrafish mutated in colorectal cancer (mcc), which encodes an evolutionarily conserved PDZ domain-containing putative tumor suppressor, is required for Wnt5b/Vangl2 signaling during gastrulation. Knockdown of mcc results in CE phenotypes similar to loss of vangl2 and wnt5b, whereas overexpression of mcc robustly rescues the depletion of wnt5b, vangl2 and the Wnt5b tyrosine kinase receptor ror2. Biochemical experiments establish a direct physical interaction between Mcc and the Vangl2 cytoplasmic tail. Lastly, CE defects in mcc morphants are suppressed by downstream activation of RhoA and JNK. Taken together, our results identify Mcc as a novel intracellular effector of non-canonical Wnt5b/Vangl2/Ror2 signaling during vertebrate gastrulation.

Periaxin and AHNAK nucleoprotein 2 form intertwined homodimers through domain swapping

Periaxin (PRX) is an abundant protein in the peripheral nervous system, with an important role in myelination. PRX participates in large molecular complexes, most likely through the interactions of its N-terminal PSD-95/Discs-large/ZO-1 (PDZ)-like domain. We present the crystal structures of the PDZ-like domains from PRX and its homologue AHNAK nucleoprotein 2 (AHNAK2). The unique intertwined, domain-swapped dimers provide a structural basis for the homodimerization of both proteins. The core of the homodimer is formed by a 6-stranded antiparallel β sheet, with every other strand from a different chain. The AHNAK2 PDZ domain structure contains a bound class III ligand peptide. The binding pocket is preformed, and the peptide-PDZ interactions have unique aspects, including two salt bridges and weak recognition of the peptide C terminus. Tight homodimerization may be central to the scaffolding functions of PRX and AHNAK2 in molecular complexes linking the extracellular matrix to the cytoskeletal network.

Discovery and confirmation of ligand binding specificities of the Schistosoma japonicum polarity protein Scribble

Background

Schistosomiasis is a chronic debilitating parasitic disease that afflicts more than 200 million individuals worldwide. Long-term administration of chemotherapy with the single available drug, praziquantel, has led to growing concerns about drug resistance. The PSD-95/Dlg/ZO-1 (PDZ) domain is an important module found in many scaffolding proteins, which has been recognized as promising targets for the development of novel drugs. However, the parasite-derived PDZ domains and their associated functions are still largely unknown.

Methodology/Principal Findings

The gene encoding the Schistosoma japonicum Scribble protein (SjScrib) was identified by homologous search with the S. mansoni Scrib sequence. By screening an arbitrary peptide library in yeast two-hybrid (Y2H) assays, we identified and confirmed the ligand binding specificity for each of the four PDZ domains of SjScrib. Both SjScrib-PDZ1 and SjScrib-PDZ3 recognize type I C-terminal PDZ-domain binding motifs (PBMs), which can be deduced as consensus sequences of -[Φ][x][E][TS][x][ILF] and -[x][RKx][E_TS][T][WΦ][ILV], respectively. SjScrib-PDZ2 prefers stringent type II C-terminal PBMs, which significantly differs from that of its human ortholog. SjScrib-PDZ4 binds to typical II C-terminal PBMs with a consensus sequence -[x][F_W][x][LI][x][LIV], in which the aromatic residue Phe is predominantly selected at position -4. The irregular and unconventional internal ligand binding specificities for the PDZ domains of SjScrib were confirmed by point mutations of the key amino acids within the ligand binding motifs. We also compared the differences in ligand specificities between SjScrib-PDZs and hScrib-PDZs, and explored the structural basis for the ligand binding properties of SjScrib-PDZs.

Conclusions/Significance

In this study, we characterized and confirmed the ligand binding specificities of all four PDZ domains of SjScrib for the first time. We denoted the differential ligand binding specificities between SjScrib-PDZs and hScrib-PDZs as well as the structural basis for these properties. This work may provide a fundamental basis for the rational design of novel anti-schistosomal drugs.

Schistosomiasis japonica remains a major public health problem in China and Southeast Asia. The long-term of treatments with the only available drug, praziquantel, has raised the concerns about drug resistance. Protein-protein interactions (PPIs), for highly discriminating specificities, are thought to be the innovative targets for a generation of new drugs. The PDZ domain is one of the most important modules for PPIs. A number of compounds screened based on binding specificities of PDZ domains have shown their potential therapeutic power in several disease models with less side effects. Although domain loss events are widespread in S. japonicum, a panel of PDZ domains is conserved in this species. So far, however, little is known about ligand binding specificities and the molecular functions of parasite-derived PDZ domain-containing proteins. In this study, by yeast two-hybrid screening of a random library, we confirmed the ligand binding properties of a multiple PDZ domain-containing protein Scribble of S. japonicum for the first time. Divergent ligand specificities between the homologous PDZ domains of S. japonicum and human Scribble orthologs were revealed. Internal motif recognition and irregular ligand interaction models for the SjScrib-PDZ domains were identified. These results provide an important basis for the rational discovery of anti-schistosomal drugs.

The periplasmic PDZ domain-containing protein Prc modulates full virulence, envelops stress responses, and directly interacts with dipeptidyl peptidase of Xanthomonas oryzae pv. oryzae

PDZ domain-containing proteases, also known as HtrA family proteases, play important roles in bacterial cells by modulating disease pathogenesis and cell-envelope stress responses. These proteases have diverse functions through proteolysis- and nonproteolysis-dependent modes. Here, we report that the genome of the causative agent of rice bacterial blight, Xanthomonas oryzae pv. oryzae, encodes seven PDZ domain-containing proteins. Systematic inactivation of their encoding genes revealed that PXO_01122 and PXO_04290 (prc) are involved in virulence. prc encodes a putative HtrA family protease that localizes in the bacterial periplasm. Mutation of prc also resulted in susceptibility to multiple environmental stresses, including H2O2, sodium dodecylsulfate, and osmolarity stresses. Comparative subproteomic analyses showed that the amounts of 34 periplasmic proteins were lower in the prc mutant than in wild-type. These proteins were associated with proteolysis, biosynthesis of macromolecules, carbohydrate or energy metabolism, signal transduction, and protein translocation or folding. We provide in vivo and in vitro evidence demonstrating that Prc stabilizes and directly binds to one of these proteins, DppP, a dipeptidyl peptidase contributing to full virulence. Taken together, our results suggest that Prc contributes to bacterial virulence by acting as a periplasmic modulator of cell-envelope stress responses.

Cluster based prediction of PDZ-peptide interactions

Background

PDZ domains are one of the most promiscuous protein recognition modules that bind with short linear peptides and play an important role in cellular signaling. Recently, few high-throughput techniques (e.g. protein microarray screen, phage display) have been applied to determine in-vitro binding specificity of PDZ domains. Currently, many computational methods are available to predict PDZ-peptide interactions but they often provide domain specific models and/or have a limited domain coverage.

Results

Here, we composed the largest set of PDZ domains derived from human, mouse, fly and worm proteomes and defined binding models for PDZ domain families to improve the domain coverage and prediction specificity. For that purpose, we first identified a novel set of 138 PDZ families, comprising of 548 PDZ domains from aforementioned organisms, based on efficient clustering according to their sequence identity. For 43 PDZ families, covering 226 PDZ domains with available interaction data, we built specialized models using a support vector machine approach. The advantage of family-wise models is that they can also be used to determine the binding specificity of a newly characterized PDZ domain with sufficient sequence identity to the known families. Since most current experimental approaches provide only positive data, we have to cope with the class imbalance problem. Thus, to enrich the negative class, we introduced a powerful semi-supervised technique to generate high confidence non-interaction data. We report competitive predictive performance with respect to state-of-the-art approaches.

Conclusions

Our approach has several contributions. First, we show that domain coverage can be increased by applying accurate clustering technique. Second, we developed an approach based on a semi-supervised strategy to get high confidence negative data. Third, we allowed high order correlations between the amino acid positions in the binding peptides. Fourth, our method is general enough and will easily be applicable to other peptide recognition modules such as SH2 domains and finally, we performed a genome-wide prediction for 101 human and 102 mouse PDZ domains and uncovered novel interactions with biological relevance. We make all the predictive models and genome-wide predictions freely available to the scientific community.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-S1-S5) contains supplementary material, which is available to authorized users.

Vaccinia virus F11 promotes viral spread by acting as a PDZ-containing scaffolding protein to bind myosin-9A and inhibit RhoA signaling

The vaccinia F11 protein promotes viral spread by modulating the cortical actin cytoskeleton by inhibiting RhoA signaling via an unknown mechanism. PDZ domains are widely conserved protein interaction modules whose occurrence in viral proteins is unprecedented. We found that F11 contains a central PDZ-like domain that is required to downregulate RhoA signaling and enhance viral spread. The PDZ-like domain interacts with the PDZ binding motif of the Rho GTPase-activating protein (GAP) Myosin-9A. In the absence of Myosin-9A, RhoA signaling is not inhibited, resulting in fewer actin tails and reduced virus release concomitant with less viral spread. The loss of Myosin-9A GAP activity or its ability to bind F11 also reduces actin tail formation. Furthermore, the ability of Myosin-9A to promote viral spread depends on F11 binding RhoA. Thus, F11 acts as a functional PDZ-containing scaffolding protein to inhibit RhoA signaling by binding Myosin-9A.

Identification of methyllysine peptides binding to chromobox protein homolog 6 chromodomain in the human proteome

Methylation is one of the important post-translational modifications that play critical roles in regulating protein functions. Proteomic identification of this post-translational modification and understanding how it affects protein activity remain great challenges. We tackled this problem from the aspect of methylation mediating protein-protein interaction. Using the chromodomain of human chromobox protein homolog 6 as a model system, we developed a systematic approach that integrates structure modeling, bioinformatics analysis, and peptide microarray experiments to identify lysine residues that are methylated and recognized by the chromodomain in the human proteome. Given the important role of chromobox protein homolog 6 as a reader of histone modifications, it was interesting to find that the majority of its interacting partners identified via this approach function in chromatin remodeling and transcriptional regulation. Our study not only illustrates a novel angle for identifying methyllysines on a proteome-wide scale and elucidating their potential roles in regulating protein function, but also suggests possible strategies for engineering the chromodomain-peptide interface to enhance the recognition of and manipulate the signal transduction mediated by such interactions.

Post-synaptic density-95 (PSD-95) binding capacity of G-protein-coupled receptor 30 (GPR30), an estrogen receptor that can be identified in hippocampal dendritic spines

The estrogen 17β-estradiol (E2) modulates dendritic spine plasticity in the cornu ammonis 1 (CA1) region of the hippocampus, and GPR30 (G-protein coupled estrogen receptor 1 (GPER1)) is an estrogen-sensitive G-protein-coupled receptor (GPCR) that is expressed in the mammalian brain and in specific subregions that are responsive to E2, including the hippocampus. The subcellular localization of hippocampal GPR30, however, remains unclear. Here, we demonstrate that GPR30 immunoreactivity is detected in dendritic spines of rat CA1 hippocampal neurons in vivo and that GPR30 protein can be found in rat brain synaptosomes. GPR30 immunoreactivity is identified at the post-synaptic density (PSD) and in the adjacent peri-synaptic zone, and GPR30 can associate with the spine scaffolding protein PSD-95 both in vitro and in vivo. This PSD-95 binding capacity of GPR30 is specific and determined by the receptor C-terminal tail that is both necessary and sufficient for PSD-95 interaction. The interaction with PSD-95 functions to increase GPR30 protein levels residing at the plasma membrane surface. GPR30 associates with the N-terminal tandem pair of PDZ domains in PSD-95, suggesting that PSD-95 may be involved in clustering GPR30 with other receptors in the hippocampus. We demonstrate that GPR30 has the potential to associate with additional post-synaptic GPCRs, including the membrane progestin receptor, the corticotropin releasing hormone receptor, and the 5HT1a serotonin receptor. These data demonstrate that GPR30 is well positioned in the dendritic spine compartment to integrate E2 sensitivity directly onto multiple inputs on synaptic activity and might begin to provide a molecular explanation as to how E2 modulates dendritic spine plasticity.

Computational structural analysis of protein interactions and networks

Protein interactions have been at the focus of computational biology in recent years. In particular, interest has come from two different communities--structural and systems biology. Here, we will discuss key systems and structural biology methods that have been used for analysis and prediction of protein-protein interactions and the insight these approaches have provided on the nature and organization of protein-protein interactions inside cells.

Molecular characterization and ligand binding specificity of the PDZ domain-containing protein GIPC3 from Schistosoma japonicum

Background

Schistosomiasis is a serious global health problem that afflicts more than 230 million people in 77 countries. Long-term mass treatments with the only available drug, praziquantel, have caused growing concerns about drug resistance. PSD-95/Dlg/ZO-1 (PDZ) domain-containing proteins are recognized as potential targets for the next generation of drug development. However, the PDZ domain-containing protein family in parasites has largely been unexplored.

Methods

We present the molecular characteristics of a PDZ domain-containing protein, GIPC3, from Schistosoma japonicum (SjGIPC3) according to bioinformatics analysis and experimental approaches. The ligand binding specificity of the PDZ domain of SjGIPC3 was confirmed by screening an arbitrary peptide library in yeast two-hybrid (Y2H) assays. The native ligand candidates were predicted by Tailfit software based on the C-terminal binding specificity, and further validated by Y2H assays.

Results

SjGIPC3 is a single PDZ domain-containing protein comprised of 328 amino acid residues. Structural prediction revealed that a conserved PDZ domain was presented in the middle region of the protein. Phylogenetic analysis revealed that SjGIPC3 and other trematode orthologues clustered into a well-defined cluster but were distinguishable from those of other phyla. Transcriptional analysis by quantitative RT-PCR revealed that the SjGIPC3 gene was relatively highly expressed in the stages within the host, especially in male adult worms. By using Y2H assays to screen an arbitrary peptide library, we confirmed the C-terminal binding specificity of the SjGIPC3-PDZ domain, which could be deduced as a consensus sequence, -[SDEC]-[STIL]-[HSNQDE]-[VIL]*. Furthermore, six proteins were predicted to be native ligand candidates of SjGIPC3 based on the C-terminal binding properties and other biological information; four of these were confirmed to be potential ligands using the Y2H system.

Conclusions

In this study, we first characterized a PDZ domain-containing protein GIPC3 in S. japonicum. The SjGIPC3-PDZ domain is able to bind both type I and II ligand C-terminal motifs. The identification of native ligand will help reveal the potential biological function of SjGIPC3. These data will facilitate the identification of novel drug targets against S. japonicum infections.

Reduced false positives in PDZ binding prediction using sequence and structural descriptors

Abstract—Identifying the binding partners of proteins is a problem of fundamental importance in computational biology. The PDZ is one of the most common and well-studied protein binding domains, hence it is a perfect model system for designing protein binding predictors. The standard approach to identifying the binding partners of PDZ domains uses multiple sequence alignments to infer the set of contact residues that are used in a predictive model. We expand on the sequence alignment approach by incorporating structural information to generate descriptors of the binding site geometry. Furthermore, we generate a real-value score for binary predictions by applying a filter based on models that predict the probability distributions of contact residues at each of the canonical PDZ ligand binding positions. Under training cross validation, our model produced an order of magnitude more predictions at a false positive proportion (FPP) of 10 percent than our benchmark model chosen from the literature. Evaluated using an independent cross validation, with computationally predicted structures, our model was able to make five times as many predictions as the benchmark model, with a Matthews' correlation coefficient (MCC) of 0.33. In addition, our model achieved a false positive proportion of 0.14, while the benchmark model had a 0.25 false positive proportion.

In silico and in vitro elucidation of BH3 binding specificity toward Bcl-2

Interactions between Bcl-2-like proteins and BH3 domains play a key role in the regulation of apoptosis. Despite the overall structural similarity of their interaction with helical BH3 domains, Bcl-2-like proteins exhibit an intricate spectrum of binding specificities whose underlying basis is not well understood. Here, we characterize these interactions using Rosetta FlexPepBind, a protocol for the prediction of peptide binding specificity that evaluates the binding potential of different peptides based on structural models of the corresponding peptide-receptor complexes. For two prominent players, Bcl-xL and Mcl-1, we obtain good agreement with a large set of experimental SPOT array measurements and recapitulate the binding specificity of peptides derived by yeast display in a previous study. We extend our approach to a third member of this family, Bcl-2: we test our blind prediction of the binding of 180 BIM-derived peptides with a corresponding experimental SPOT array. Both prediction and experiment reveal a Bcl-2 binding specificity pattern that resembles that of Bcl-xL. Finally, we extend this application to accurately predict the specificity pattern of additional human BH3-only derived peptides. This study characterizes the distinct patterns of binding specificity of BH3-only derived peptides for the Bcl-2 like proteins Bcl-xL, Mcl-1, and Bcl-2 and provides insight into the structural basis of determinants of specificity.

Identification of new binding partners of the chemosensory signaling protein Gγ13 expressed in taste and olfactory sensory cells

Tastant detection in the oral cavity involves selective receptors localized at the apical extremity of a subset of specialized taste bud cells called taste receptor cells (TRCs). The identification of the genes coding for the taste receptors involved in this process have greatly improved our understanding of the molecular mechanisms underlying detection. However, how these receptors signal in TRCs, and whether the components of the signaling cascades interact with each other or are organized in complexes is mostly unexplored. Here we report on the identification of three new binding partners for the mouse G protein gamma 13 subunit (Gγ13), a component of the bitter taste receptors signaling cascade. For two of these Gγ13 associated proteins, namely GOPC and MPDZ, we describe the expression in taste bud cells for the first time. Furthermore, we demonstrate by means of a yeast two-hybrid interaction assay that the C terminal PDZ binding motif of Gγ13 interacts with selected PDZ domains in these proteins. In the case of the PDZ domain-containing protein zona occludens-1 (ZO-1), a major component of the tight junction defining the boundary between the apical and baso-lateral region of TRCs, we identified the first PDZ domain as the site of strong interaction with Gγ13. This association was further confirmed by co-immunoprecipitation experiments in HEK 293 cells. In addition, we present immunohistological data supporting partial co-localization of GOPC, MPDZ, or ZO-1, and Gγ13 in taste buds cells. Finally, we extend this observation to olfactory sensory neurons (OSNs), another type of chemosensory cells known to express both ZO-1 and Gγ13. Taken together our results implicate these new interaction partners in the sub-cellular distribution of Gγ13 in olfactory and gustatory primary sensory cells.

Elucidation of the binding preferences of peptide recognition modules: SH3 and PDZ domains

Peptide-binding domains play a critical role in regulation of cellular processes by mediating protein interactions involved in signalling. In recent years, the development of large-scale technologies has enabled exhaustive studies on the peptide recognition preferences for a number of peptide-binding domain families. These efforts have provided significant insights into the binding specificities of these modular domains. Many research groups have taken advantage of this unprecedented volume of specificity data and have developed a variety of new algorithms for the prediction of binding specificities of peptide-binding domains and for the prediction of their natural binding targets. This knowledge has also been applied to the design of synthetic peptide-binding domains in order to rewire protein-protein interaction networks. Here, we describe how these experimental technologies have impacted on our understanding of peptide-binding domain specificities and on the elucidation of their natural ligands. We discuss SH3 and PDZ domains as well characterized examples, and we explore the feasibility of expanding high-throughput experiments to other peptide-binding domains.

Putting into practice domain-linear motif interaction predictions for exploration of protein networks

PDZ domains recognise short sequence motifs at the extreme C-termini of proteins. A model based on microarray data has been recently published for predicting the binding preferences of PDZ domains to five residue long C-terminal sequences. Here we investigated the potential of this predictor for discovering novel protein interactions that involve PDZ domains. When tested on real negative data assembled from published literature, the predictor displayed a high false positive rate (FPR). We predicted and experimentally validated interactions between four PDZ domains derived from the human proteins MAGI1 and SCRIB and 19 peptides derived from human and viral C-termini of proteins. Measured binding intensities did not correlate with prediction scores, and the high FPR of the predictor was confirmed. Results indicate that limitations of the predictor may arise from an incomplete model definition and improper training of the model. Taking into account these limitations, we identified several novel putative interactions between PDZ domains of MAGI1 and SCRIB and the C-termini of the proteins FZD4, ARHGAP6, NET1, TANC1, GLUT7, MARCH3, MAS, ABC1, DLL1, TMEM215 and CYSLTR2. These proteins are localised to the membrane or suggested to act close to it and are often involved in G protein signalling. Furthermore, we showed that, while extension of minimal interacting domains or peptides toward tandem constructs or longer peptides never suppressed their ability to interact, the measured affinities and inferred specificity patterns often changed significantly. This suggests that if protein fragments interact, the full length proteins are also likely to interact, albeit possibly with altered affinities and specificities. Therefore, predictors dealing with protein fragments are promising tools for discovering protein interaction networks but their application to predict binding preferences within networks may be limited.

Characterization of PDZ domain-peptide interaction interface based on energetic patterns

PDZ domain is one of the abundant modular domains that recognize short peptide sequences to mediate protein-protein interactions. To decipher the binding specificity of PDZ domain, we analyzed the interactions between 11 mouse PDZ domains and 217 [corrected] peptides using a method called MIEC-SVM, which energetically characterizes the domain-peptide interaction using molecular interaction energy components (MIECs) and predicts binding specificity using support vector machine (SVM). Cross-validation and leave-one-domain-out test showed that the MIEC-SVM using all 44 PDZ-peptide residue pairs at the interaction interface outperformed the sequence-based methods in the literature. A further feature (residue pair) selection procedure illustrated that 16 residue pairs were uninformative to the binding specificity, even though they contributed significantly (~50%) to the binding energy. If only using the 28 informative residue pairs, the performance of the MIEC-SVM on predicting the PDZ binding specificity was significantly improved. This analysis suggests that the informative and uninformative residue interactions between the PDZ domain and the peptide may represent those contributing to binding specificity and affinity, respectively. We performed additional structural and energetic analyses to shed light on understanding how the PDZ-peptide recognition is established. The success of the MIEC-SVM method on PDZ domains in this study and SH3 domains in our previous studies illustrates its generality on characterizing protein-peptide interactions and understanding protein recognition from a structural and energetic viewpoint.

Genome-wide analysis of PDZ domain binding reveals inherent functional overlap within the PDZ interaction network

Binding selectivity and cross-reactivity within one of the largest and most abundant interaction domain families, the PDZ family, has long been enigmatic. The complete human PDZ domain complement (the PDZome) consists of 267 domains and we applied here a Bayesian selectivity model to predict hundreds of human PDZ domain interactions, using target sequences of 22,997 non-redundant proteins. Subsequent analysis of these binding scores shows that PDZs can be divided into two genome-wide clusters that coincide well with the division between canonical class 1 and 2 PDZs. Within the class 1 PDZs we observed binding overlap at unprecedented levels, mediated by two residues at positions 1 and 5 of the second α-helix of the binding pocket. Eight PDZ domains were subsequently selected for experimental binding studies and to verify the basics of our predictions. Overall, the PDZ domain class 1 cross-reactivity identified here implies that auxiliary mechanisms must be in place to overcome this inherent functional overlap and to minimize cross-selectivity within the living cell. Indeed, when we superimpose PDZ domain binding affinities with gene ontologies, network topology data and the domain position within a PDZ superfamily protein, functional overlap is minimized and PDZ domains position optimally in the binding space. We therefore propose that PDZ domain selectivity is achieved through cellular context rather than inherent binding specificity.