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

The conformation of the nSrc specificity-determining loop in the Src SH3 domain is modulated by a WX conserved sequence motif found in SH3 domains

Cellular signaling networks are modulated by multiple protein-protein interaction domains that coordinate extracellular inputs and processes to regulate cellular processes. Several of these domains recognize short linear motifs, or SLiMs, which are often highly conserved and are closely regulated. One such domain, the Src homology 3 (SH3) domain, typically recognizes proline-rich SLiMs and is one of the most abundant SLiM-binding domains in the human proteome. These domains are often described as quite versatile, and indeed, SH3 domains can bind ligands in opposite orientations dependent on target sequence. Furthermore, recent work has identified diverse modes of binding for SH3 domains and a wide variety of sequence motifs that are recognized by various domains. Specificity is often attributed to the RT and nSrc loops near the peptide-binding cleft in this domain family, particularly for Class I binding, which is defined as RT and nSrc loop interactions with the N-terminus of the ligand. Here, we used the Src and Abl SH3 domains as a model to further investigate the role of the RT and nSrc loops in SH3 specificity. We created chimeric domains with both the RT and nSrc loop sequences swapped between these SH3 domains, and used fluorescence anisotropy assays to test how relative binding affinities were affected for Src SH3- and Abl SH3-specific ligands. We also used Alphafold-Multimer to model our SH3:peptide complexes in combination with molecular dynamics simulations. We identified a position that contributes to the nSrc loop conformation in Src SH3, the amino acid immediately following a highly conserved Trp that creates a hydrophobic pocket critical for SH3 ligand recognition. We defined this as the WX motif, where X = Trp for Src and Cys for Abl. A broad importance of this position for modulating nSrc loop conformation in SH3 domains is suggested by analyses of previously deposited SH3 structures, multiple sequence alignment of SH3 domains in the human proteome, and our biochemical and computational data of mutant Src and Abl SH3 domains. Overall, our work uses experimental approaches and structural modeling to better understand specificity determinants in SH3 domains.

Receptor Determinants for β-Arrestin Functional Specificity at C-X-C Chemokine Receptor 5

β-arrestins are multifaceted adaptor proteins that mediate G protein-coupled receptor (GPCR) desensitization, internalization, and signaling. It is emerging that receptor-specific determinants specify these divergent functions at GPCRs, yet this remains poorly understood. Here, we set out to identify the receptor determinants responsible for β-arrestin-mediated regulation of the chemokine receptor C-X-C motif chemokine receptor 5 (CXCR5). Using bioluminescence resonance energy transfer, we show that β-arrestin1 and β-arrestin2 are dose-dependently recruited to CXCR5 by its cognate ligand C-X-C motif chemokine ligand 13 (CXCL13). The carboxy-terminal tail of CXCR5 contains several serine/threonine residues that can be divided into three discrete phospho-site clusters based on their position relative to transmembrane domain 7. Mutagenesis experiments revealed that the distal and medial phospho-site clusters, but not the proximal, are required for agonist-stimulated β-arrestin1 or β-arrestin2 recruitment to CXCR5. Consistent with this, we provide evidence that the distal and medial, but not proximal, phospho-site clusters are required for receptor desensitization. Surprisingly, the individual phospho-site clusters are not required for agonist-stimulated internalization of CXCR5. Further, we show that CXCL13-stimulated CXCR5 internalization and ERK1/2 phosphorylation, but not desensitization, remain intact in human embryonic kidney 293 cells lacking β-arrestin1 and β-arrestin2. Our study provides evidence that β-arrestins are recruited to CXCR5 and are required for desensitization but are dispensable for internalization or signaling, suggesting that discrete receptor determinants specify the divergent functions of β-arrestins. SIGNIFICANCE STATEMENT: C-X-C motif ligand 13 (CXCL13) and C-X-C motif chemokine receptor 5 (CXCR5) are important in the immune system and are linked to diseases, yet regulation of CXCR5 signaling remains poorly understood. We provide evidence that a phospho-site cluster located at the extreme distal carboxyl-terminal tail of the receptor is responsible for β-arrestin recruitment and receptor desensitization. β-arrestins are not required for CXCL13-stimulated internalization or signaling, indicating that β-arrestins perform only one of their functions at CXCR5 and that discrete receptor determinants specify the divergent functions of β-arrestins.

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.

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.

A systematic study of traditional Chinese medicine treating hepatitis B virus-related hepatocellular carcinoma based on target-driven reverse network pharmacology

Hepatocellular carcinoma (HCC) is a serious global health problem, and hepatitis B virus (HBV) infection remains the leading cause of HCC. It is standard care to administer antiviral treatment for HBV-related HCC patients with concurrent anti-cancer therapy. However, a drug with repressive effects on both HBV infection and HCC has not been discovered yet. In addition, drug resistance and side effects have made existing therapeutic regimens suboptimal. Traditional Chinese medicine (TCM) has multi-ingredient and multi-target advantages in dealing with multifactorial HBV infection and HCC. TCM has long been served as a valuable source and inspiration for discovering new drugs. In present study, a target-driven reverse network pharmacology was applied for the first time to systematically study the therapeutic potential of TCM in treating HBV-related HCC. Firstly, 47 shared targets between HBV and HCC were screened as HBV-related HCC targets. Next, starting from 47 targets, the relevant chemical components and herbs were matched. A network containing 47 targets, 913 chemical components and 469 herbs was established. Then, the validated results showed that almost 80% of the herbs listed in chronic hepatitis B guidelines and primary liver cancer guidelines were included in the 469 herbs. Furthermore, functional analysis was conducted to understand the biological processes and pathways regulated by these 47 targets. The docking results indicated that the top 50 chemical components bound well to targets. Finally, the frequency statistical analysis results showed the 469 herbs against HBV-related HCC were mainly warm in property, bitter in taste, and distributed to the liver meridians. Taken together, a small library of 913 chemical components and 469 herbs against HBV-related HCC were obtained with a target-driven approach, thus paving the way for the development of therapeutic modalities to treat HBV-related HCC.

Peptide Targeting of PDZ-Dependent Interactions as Pharmacological Intervention in Immune-Related Diseases

PDZ (postsynaptic density (PSD95), discs large (Dlg), and zonula occludens (ZO-1)-dependent interactions are widely distributed within different cell types and regulate a variety of cellular processes. To date, some of these interactions have been identified as targets of small molecules or peptides, mainly related to central nervous system disorders and cancer. Recently, the knowledge of PDZ proteins and their interactions has been extended to various cell types of the immune system, suggesting that their targeting by viral pathogens may constitute an immune evasion mechanism that favors viral replication and dissemination. Thus, the pharmacological modulation of these interactions, either with small molecules or peptides, could help in the control of some immune-related diseases. Deeper structural and functional knowledge of this kind of protein–protein interactions, especially in immune cells, will uncover novel pharmacological targets for a diversity of clinical conditions.

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.

Sequence variation in the β7-β8 loop of bacterial class A sortase enzymes alters substrate selectivity

Gram-positive bacteria contain sortase enzymes on their cell surfaces that catalyze transpeptidation reactions critical for proper cellular function. In vitro, sortases are used in sortase-mediated ligation (SML) reactions for a variety of protein engineering applications. Historically, sortase A from Staphylococcus aureus (saSrtA) has been the enzyme of choice to catalyze SML reactions. However, the stringent specificity of saSrtA for the LPXTG sequence motif limits its uses. Here, we describe the impact on substrate selectivity of a structurally conserved loop with a high degree of sequence variability in all classes of sortases. We investigate the contribution of this β7–β8 loop by designing and testing chimeric sortase enzymes. Our chimeras utilize natural sequence variation of class A sortases from eight species engineered into the SrtA sequence from Streptococcus pneumoniae. While some of these chimeric enzymes mimic the activity and selectivity of the WT protein from which the loop sequence was derived (e.g., that of saSrtA), others results in chimeric Streptococcus pneumoniae SrtA enzymes that are able to accommodate a range of residues in the final position of the substrate motif (LPXTX). Using mutagenesis, structural comparisons, and sequence analyses, we identify three interactions facilitated by β7–β8 loop residues that appear to be broadly conserved or converged upon in class A sortase enzymes. These studies provide the foundation for a deeper understanding of sortase target selectivity and can expand the sortase toolbox for future SML applications.

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.

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

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

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.