Human PDZ Domains Research Articles

Abstract 2584: Investigating the effects of non-motif residues on the binding affinity of CFTR and HPV16 to human PDZ domains

Abstract 2584: Investigating the effects of non-motif residues on the binding affinity of CFTR and HPV16 to human PDZ domains

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

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

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

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

Structural basis of coronavirus E protein interactions with human PALS1 PDZ domain

SARS-CoV-2 infection leads to coronavirus disease 2019 (COVID-19), which is associated with severe and life-threatening pneumonia and respiratory failure. However, the molecular basis of these symptoms remains unclear. SARS-CoV-1 E protein interferes with control of cell polarity and cell-cell junction integrity in human epithelial cells by binding to the PALS1 PDZ domain, a key component of the Crumbs polarity complex. We show that C-terminal PDZ binding motifs of SARS-CoV-1 and SARS-CoV-2 E proteins bind the PALS1 PDZ domain with 29.6 and 22.8 μM affinity, whereas the related sequence from MERS-CoV did not bind. We then determined crystal structures of PALS1 PDZ domain bound to both SARS-CoV-1 and SARS-CoV-2 E protein PDZ binding motifs. Our findings establish the structural basis for SARS-CoV-1/2 mediated subversion of Crumbs polarity signalling and serve as a platform for the development of small molecule inhibitors to suppress SARS-CoV-1/2 mediated disruption of polarity signalling in epithelial cells.

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.

Interactomic affinity profiling by holdup assay: Acetylation and distal residues impact the PDZome-binding specificity of PTEN phosphatase.

Protein domains often recognize short linear protein motifs composed of a core conserved consensus sequence surrounded by less critical, modulatory positions. PTEN, a lipid phosphatase involved in phosphatidylinositol 3-kinase (PI3K) pathway, contains such a short motif located at the extreme C-terminus capable to recognize PDZ domains. It has been shown that the acetylation of this motif could modulate the interaction with several PDZ domains. Here we used an accurate experimental approach combining high-throughput holdup chromatographic assay and competitive fluorescence polarization technique to measure quantitative binding affinity profiles of the PDZ domain-binding motif (PBM) of PTEN. We substantially extended the previous knowledge towards the 266 known human PDZ domains, generating the full PDZome-binding profile of the PTEN PBM. We confirmed that inclusion of N-terminal flanking residues, acetylation or mutation of a lysine at a modulatory position significantly altered the PDZome-binding profile. A numerical specificity index is also introduced as an attempt to quantify the specificity of a given PBM over the complete PDZome. Our results highlight the impact of modulatory residues and post-translational modifications on PBM interactomes and their specificity.

Structural characterization and computational analysis of PDZ domains in Monosiga brevicollis.

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

Dual Specificity PDZ- and 14-3-3-Binding Motifs: A Structural and Interactomics Study

Protein-protein interaction motifs are often alterable by post-translational modifications. For example, 19% of predicted human PDZ domain-binding motifs (PBMs) have been experimentally proven to be phosphorylated, and up to 82% are theoretically phosphorylatable. Phosphorylation of PBMs may drastically rewire their interactomes, by altering their affinities for PDZ domains and 14-3-3 proteins. The effect of phosphorylation is often analyzed by performing "phosphomimetic" mutations. Here, we focused on the PBMs of HPV16-E6 viral oncoprotein and human RSK1 kinase. We measured the binding affinities of native, phosphorylated, and phosphomimetic variants of both PBMs toward the 266 human PDZ domains. We co-crystallized all the motif variants with a selected PDZ domain to characterize the structural consequence of the different modifications. Finally, we elucidated the structural basis of PBM capture by 14-3-3 proteins. This study provides novel atomic and interactomic insights into phosphorylatable dual specificity motifs and the differential effects of phosphorylation and phosphomimetic approaches.

EF-X in Silico - Modeling Protein Dynamics in an Electric Field

Protein functions result from local and collective atomic motions that span a wide range of time scales. An integrated analysis of experimental and simulation data can shed light on the detailed mechanism of these motions. Applying a high electric field to protein crystals enables conformational changes that can be captured by time-resolved X-ray crystallography. Such an experiment (referred to as EF-X) carried out on a human PDZ domain obtained a series of atomistic “snapshots” of ensemble averages protein dynamics at 50 to 100 ns time intervals. Here, we present a molecular dynamics (MD) study of the same system and provide a detailed picture of the protein dynamics in between the experimental “snapshots”. We replicated the experimental conditions and system geometry in the presence and absence of an electric field. By constructing a model of the protein crystal as a 3x3x3 supercell with a total of 108 individual proteins, we achieved extensive sampling of the protein conformational ensemble at the sub-millisecond time scale. A number of techniques, including principal component analysis and strain analysis, were utilized to quantify the effects of the electric field on the protein structure and crystal symmetry. This study demonstrates how MD simulations can complement information obtained in EF-X experiments by providing the higher spatial and temporal resolution of underlying dynamical processes.

MotifAnalyzer-PDZ: A computational program to investigate the evolution of PDZ-binding target specificity.

Recognition of short linear motifs (SLiMs) or peptides by proteins is an important component of many cellular processes. However, due to limited and degenerate binding motifs, prediction of cellular targets is challenging. In addition, many of these interactions are transient and of relatively low affinity. Here, we focus on one of the largest families of SLiM-binding domains in the human proteome, the PDZ domain. These domains bind the extreme C-terminus of target proteins, and are involved in many signaling and trafficking pathways. To predict endogenous targets of PDZ domains, we developed MotifAnalyzer-PDZ, a program that filters and compares all motif-satisfying sequences in any publicly available proteome. This approach enables us to determine possible PDZ binding targets in humans and other organisms. Using this program, we predicted and biochemically tested novel human PDZ targets by looking for strong sequence conservation in evolution. We also identified three C-terminal sequences in choanoflagellates that bind a choanoflagellate PDZ domain, the Monsiga brevicollis SHANK1 PDZ domain (mbSHANK1), with endogenously-relevant affinities, despite a lack of conservation with the targets of a homologous human PDZ domain, SHANK1. All three are predicted to be signaling proteins, with strong sequence homology to cytosolic and receptor tyrosine kinases. Finally, we analyzed and compared the positional amino acid enrichments in PDZ motif-satisfying sequences from over a dozen organisms. Overall, MotifAnalyzer-PDZ is a versatile program to investigate potential PDZ interactions. This proof-of-concept work is poised to enable similar types of analyses for other SLiM-binding domains (e.g., MotifAnalyzer-Kinase). MotifAnalyzer-PDZ is available at http://motifAnalyzerPDZ.cs.wwu.edu.

Electric-field-stimulated protein mechanics.

The internal mechanics of proteins-the coordinated motions of amino acids and the pattern of forces constraining these motions-connects protein structure to function. Here we describe a new method combining the application of strong electric field pulses to protein crystals with time-resolved X-ray crystallography to observe conformational changes in spatial and temporal detail. Using a human PDZ domain (LNX2PDZ2) as a model system, we show that protein crystals tolerate electric field pulses strong enough to drive concerted motions on the sub-microsecond timescale. The induced motions are subtle, involve diverse physical mechanisms, and occur throughout the protein structure. The global pattern of electric-field-induced motions is consistent with both local and allosteric conformational changes naturally induced by ligand binding, including at conserved functional sites in the PDZ domain family. This work lays the foundation for comprehensive experimental study of the mechanical basis of protein function.

ModPDZpep: a web resource for structure based analysis of human PDZ-mediated interaction networks

BackgroundPDZ domains recognize short sequence stretches usually present in C-terminal of their interaction partners. Because of the involvement of PDZ domains in many important biological processes, several attempts have been made for developing bioinformatics tools for genome-wide identification of PDZ interaction networks. Currently available tools for prediction of interaction partners of PDZ domains utilize machine learning approach. Since, they have been trained using experimental substrate specificity data for specific PDZ families, their applicability is limited to PDZ families closely related to the training set. These tools also do not allow analysis of PDZ-peptide interaction interfaces.ResultsWe have used a structure based approach to develop modPDZpep, a program to predict the interaction partners of human PDZ domains and analyze structural details of PDZ interaction interfaces. modPDZpep predicts interaction partners by using structural models of PDZ-peptide complexes and evaluating binding energy scores using residue based statistical pair potentials. Since, it does not require training using experimental data on peptide binding affinity, it can predict substrates for diverse PDZ families. Because of the use of simple scoring function for binding energy, it is also fast enough for genome scale structure based analysis of PDZ interaction networks. Benchmarking using artificial as well as real negative datasets indicates good predictive power with ROC-AUC values in the range of 0.7 to 0.9 for a large number of human PDZ domains. Another novel feature of modPDZpep is its ability to map novel PDZ mediated interactions in human protein-protein interaction networks, either by utilizing available experimental phage display data or by structure based predictions.ConclusionsIn summary, we have developed modPDZpep, a web-server for structure based analysis of human PDZ domains. It is freely available at http://www.nii.ac.in/modPDZpep.html or http://202.54.226.235/modPDZpep.html.ReviewersThis article was reviewed by Michael Gromiha and Zoltán Gáspári.Electronic supplementary materialThe online version of this article (doi:10.1186/s13062-016-0151-4) contains supplementary material, which is available to authorized users.

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.

Structures of the human Pals1 PDZ domain with and without ligand suggest gated access of Crb to the PDZ peptide-binding groove.

Many components of epithelial polarity protein complexes possess PDZ domains that are required for protein interaction and recruitment to the apical plasma membrane. Apical localization of the Crumbs (Crb) transmembrane protein requires a PDZ-mediated interaction with Pals1 (protein-associated with Lin7, Stardust, MPP5), a member of the p55 family of membrane-associated guanylate kinases (MAGUKs). This study describes the molecular interaction between the Crb carboxy-terminal motif (ERLI), which is required for Drosophila cell polarity, and the Pals1 PDZ domain using crystallography and fluorescence polarization. Only the last four Crb residues contribute to Pals1 PDZ-domain binding affinity, with specificity contributed by conserved charged interactions. Comparison of the Crb-bound Pals1 PDZ structure with an apo Pals1 structure reveals a key Phe side chain that gates access to the PDZ peptide-binding groove. Removal of this side chain enhances the binding affinity by more than fivefold, suggesting that access of Crb to Pals1 may be regulated by intradomain contacts or by protein-protein interaction.

Discovery of novel interacting partners of PSMD9, a proteasomal chaperone: Role of an Atypical and versatile PDZ-domain motif interaction and identification of putative functional modules

PSMD9 (Proteasome Macropain non-ATPase subunit 9), a proteasomal assembly chaperone, harbors an uncharacterized PDZ-like domain. Here we report the identification of five novel interacting partners of PSMD9 and provide the first glimpse at the structure of the PDZ-domain, including the molecular details of the interaction. We based our strategy on two propositions: (a) proteins with conserved C-termini may share common functions and (b) PDZ domains interact with C-terminal residues of proteins. Screening of C-terminal peptides followed by interactions using full-length recombinant proteins, we discovered hnRNPA1 (an RNA binding protein), S14 (a ribosomal protein), CSH1 (a growth hormone), E12 (a transcription factor) and IL6 receptor as novel PSMD9-interacting partners. Through multiple techniques and structural insights, we clearly demonstrate for the first time that human PDZ domain interacts with the predicted Short Linear Sequence Motif (SLIM) at the C-termini of the client proteins. These interactions are also recapitulated in mammalian cells. Together, these results are suggestive of the role of PSMD9 in transcriptional regulation, mRNA processing and editing, hormone and receptor activity and protein translation. Our proof-of-principle experiments endorse a novel and quick method for the identification of putative interacting partners of similar PDZ-domain proteins from the proteome and for discovering novel functions.

The Human PDZome: A Gateway to PSD95-Disc Large-Zonula Occludens (PDZ)-mediated Functions

Protein–protein interactions organize the localization, clustering, signal transduction, and degradation of cellular proteins and are therefore implicated in numerous biological functions. These interactions are mediated by specialized domains able to bind to modified or unmodified peptides present in binding partners. Among the most broadly distributed protein interaction domains, PSD95-disc large-zonula occludens (PDZ) domains are usually able to bind carboxy-terminal sequences of their partners. In an effort to accelerate the discovery of PDZ domain interactions, we have constructed an array displaying 96% of the human PDZ domains that is amenable to rapid two-hybrid screens in yeast. We have demonstrated that this array can efficiently identify interactions using carboxy-terminal sequences of PDZ domain binders such as the E6 oncoviral protein and protein kinases (PDGFRβ, BRSK2, PCTK1, ACVR2B, and HER4); this has been validated via mass spectrometry analysis. Taking advantage of this array, we show that PDZ domains of Scrib and SNX27 bind to the carboxy-terminal region of the planar cell polarity receptor Vangl2. We also have demonstrated the requirement of Scrib for the promigratory function of Vangl2 and described the morphogenetic function of SNX27 in the early Xenopus embryo. The resource presented here is thus adapted for the screen of PDZ interactors and, furthermore, should facilitate the understanding of PDZ-mediated functions.

Predicting PDZ domain mediated protein interactions from structure

BackgroundPDZ domains are structural protein domains that recognize simple linear amino acid motifs, often at protein C-termini, and mediate protein-protein interactions (PPIs) in important biological processes, such as ion channel regulation, cell polarity and neural development. PDZ domain-peptide interaction predictors have been developed based on domain and peptide sequence information. Since domain structure is known to influence binding specificity, we hypothesized that structural information could be used to predict new interactions compared to sequence-based predictors.ResultsWe developed a novel computational predictor of PDZ domain and C-terminal peptide interactions using a support vector machine trained with PDZ domain structure and peptide sequence information. Performance was estimated using extensive cross validation testing. We used the structure-based predictor to scan the human proteome for ligands of 218 PDZ domains and show that the predictions correspond to known PDZ domain-peptide interactions and PPIs in curated databases. The structure-based predictor is complementary to the sequence-based predictor, finding unique known and novel PPIs, and is less dependent on training–testing domain sequence similarity. We used a functional enrichment analysis of our hits to create a predicted map of PDZ domain biology. This map highlights PDZ domain involvement in diverse biological processes, some only found by the structure-based predictor. Based on this analysis, we predict novel PDZ domain involvement in xenobiotic metabolism and suggest new interactions for other processes including wound healing and Wnt signalling.ConclusionsWe built a structure-based predictor of PDZ domain-peptide interactions, which can be used to scan C-terminal proteomes for PDZ interactions. We also show that the structure-based predictor finds many known PDZ mediated PPIs in human that were not found by our previous sequence-based predictor and is less dependent on training–testing domain sequence similarity. Using both predictors, we defined a functional map of human PDZ domain biology and predict novel PDZ domain function. Users may access our structure-based and previous sequence-based predictors at http://webservice.baderlab.org/domains/POW.

The Role of Flexibility and Conformational Selection in the Binding Promiscuity of PDZ Domains

In molecular recognition, it is often the case that ligand binding is coupled to conformational change in one or both of the binding partners. Two hypotheses describe the limiting cases involved; the first is the induced fit and the second is the conformational selection model. The conformational selection model requires that the protein adopts conformations that are similar to the ligand-bound conformation in the absence of ligand, whilst the induced-fit model predicts that the ligand-bound conformation of the protein is only accessible when the ligand is actually bound. The flexibility of the apo protein clearly plays a major role in these interpretations. For many proteins involved in signaling pathways there is the added complication that they are often promiscuous in that they are capable of binding to different ligand partners. The relationship between protein flexibility and promiscuity is an area of active research and is perhaps best exemplified by the PDZ domain family of proteins. In this study we use molecular dynamics simulations to examine the relationship between flexibility and promiscuity in five PDZ domains: the human Dvl2 (Dishevelled-2) PDZ domain, the human Erbin PDZ domain, the PDZ1 domain of InaD (inactivation no after-potential D protein) from fruit fly, the PDZ7 domain of GRIP1 (glutamate receptor interacting protein 1) from rat and the PDZ2 domain of PTP-BL (protein tyrosine phosphatase) from mouse. We show that despite their high structural similarity, the PDZ binding sites have significantly different dynamics. Importantly, the degree of binding pocket flexibility was found to be closely related to the various characteristics of peptide binding specificity and promiscuity of the five PDZ domains. Our findings suggest that the intrinsic motions of the apo structures play a key role in distinguishing functional properties of different PDZ domains and allow us to make predictions that can be experimentally tested.

Calcyon Forms a Novel Ternary Complex with Dopamine D1 Receptor through PSD-95 Protein and Plays a Role in Dopamine Receptor Internalization

Calcyon, once known for interacting directly with the dopamine D(1) receptor (D(1)DR), is implicated in various neuropsychiatric disorders including schizophrenia, bipolar disorder, and attention deficit hyperactivity disorder. Although its direct interaction with D(1)DR has been shown to be misinterpreted, it still plays important roles in D(1)DR signaling. Here, we found that calcyon interacts with the PSD-95 and subsequently forms a ternary complex with D(1)DR through PSD-95. Calcyon is phosphorylated on Ser-169 by the PKC activator phorbol 12-myristate 13-acetate or by the D(1)DR agonist SKF-81297, and its phosphorylation increases its association with PSD-95 and recruitment to the cell surface. Interestingly, the internalization of D(1)DR at the cell surface was enhanced by phorbol 12-myristate 13-acetate and SKF-81297 in the presence of calcyon, but not in the presence of its S169A phospho-deficient mutant, suggesting that the phosphorylation of calcyon and the internalization of the surface D(1)DR are tightly correlated. Our results suggest that calcyon regulates D(1)DR trafficking by forming a ternary complex with D(1)DR through PSD-95 and thus possibly linking glutamatergic and dopamine receptor signalings. This also raises the possibility that a novel ternary complex could represent a potential therapeutic target for the modulation of related neuropsychiatric disorders.

Two PDZ binding motifs within NS5 have roles in Tick-borne encephalitis virus replication

The flavivirus genus includes important human neurotropic pathogens like Tick-borne encephalitis virus (TBEV) and West-Nile virus (WNV). Flavivirus replication occurs at replication complexes, where the NS5 protein provides both RNA cap methyltransferase and RNA-dependent RNA polymerase activities. TBEVNS5 contains two PDZ binding motifs (PBMs) important for specific targeting of human PDZ proteins including Scribble, an association important for viral down regulation of cellular defense systems and neurite outgrowth.To determine whether the PBMs of TBEVNS5 affects virus replication we constructed a DNA based sub-genomic TBEV replicon expressing firefly luciferase. The PBMs within NS5 were mutated individually and in concert and the replicons were assayed in cell culture. Our results show that the replication rate was impaired in all mutants, which indicates that PDZ dependent host interactions influence TBEV replication. We also find that the C-terminal PBMs present in TBEVNS5 and WNVNS5 are targeting various human PDZ domain proteins. TBEVNS5 has affinity to Zonula occludens-2 (ZO-2), GIAP C-terminus interacting protein (GIPC), calcium/calmodulin-dependent serine protein kinase (CASK), glutamate receptor interacting protein 2, (GRIP2) and Interleukin 16 (IL-16). A different pattern was observed for WNVNS5 as it associate with a broader repertoire of putative host PDZ proteins.