skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.
Attention:The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 7:00 AM ET to 7:30 AM ET on Friday, April 24 due to maintenance. We apologize for the inconvenience.


Title: 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.  more » « less
Award ID(s):
1904711
PAR ID:
10350863
Author(s) / Creator(s):
; ; ; ; ; ;
Date Published:
Journal Name:
Molecules
Volume:
26
Issue:
19
ISSN:
1420-3049
Page Range / eLocation ID:
6034
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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 quiteversatile, 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. 
    more » « less
  2. Dynamic allosterism allows the propagation of signal throughout a protein. The PDZ (PSD-95/Dlg1/ZO-1) family has been named as a classic example of dynamic allostery in small modular domains. While the PDZ family consists of more than 200 domains, previous efforts have primarily focused on a few well-studied PDZ domains, including PTP-BL PDZ2, PSD-95 PDZ3, and Par6 PDZ. Taken together, experimental and computational studies have identified regions of these domains that are dynamically coupled to ligand binding. These regions include the αA helix, the αB lower-loop, and the αC helix. In this review, we summarize the specific residues on the αA helix, the αB lower-loop, and the αC helix of PTP-BL PDZ2, PSD-95 PDZ3, and Par6 PDZ that have been identified as participants in dynamic allostery by either experimental or computational approaches. This review can serve as an index for researchers to look back on the previously identified allostery in the PDZ family. Interestingly, our summary of previous work reveals clear consistencies between the domains. While the PDZ family has a low sequence identity, we show that some of the most consistently identified allosteric residues within PTP-BL PDZ2 and PSD-95 PDZ3 domains are evolutionarily conserved. These residues include A46/A347, V61/V362, and L66/L367 on PTP-BL PDZ2 and PSD-95 PDZ3, respectively. Finally, we expose a need for future work to explore dynamic allostery within (1) PDZ domains with multiple binding partners and (2) multidomain constructs containing a PDZ domain. 
    more » « less
  3. The PDZ family has drawn attention as possible drug targets because of the domains’ wide ranges of function and highly conserved binding pockets. The PICK1 PDZ domain has been proposed as a possible drug target because the interactions between the PICK1 PDZ domain and the GluA2 subunit of the AMPA receptor have been shown to progress neurodegenerative diseases. BIO124 has been identified as a sub µM inhibitor of the PICK1–GluA2 interaction. Here, we use all-atom molecular dynamics simulations to reveal the atomic-level interaction pattern between the PICK1 PDZ domain and BIO124. Our simulations reveal three unique binding conformations of BIO124 in the PICK1 PDZ binding pocket, referred to here as state 0, state 1, and state 2. Each conformation is defined by a unique hydrogen bonding network and a unique pattern of hydrophobic interactions between BIO124 and the PICK1 PDZ domain. Interestingly, each conformation of BIO124 results in different dynamic changes to the PICK1 PDZ domain. Unlike states 1 and 2, state 0 induces dynamic coupling between BIO124 and the αA helix. Notably, this dynamic coupling with the αA helix is similar to what has been observed in other PDZ–ligand complexes. Our analysis indicates that the interactions formed between BIO124 and I35 may be the key to inducing dynamic coupling with the αA helix. Lastly, we suspect that the conformational shifts observed in our simulations may affect the stability and thus the overall effectiveness of BIO124. We propose that a physically larger inhibitor may be necessary to ensure sufficient interactions that permit stable binding between a drug and the PICK1 PDZ domain. 
    more » « less
  4. Introduction: The PICK1 PDZ domain has been identified as a potential drug target forneurological disorders. After many years of effort, a few inhibitors, such as TAT-C5 and mPD5,have been discovered experimentally to bind to the PDZ domain with a relatively high bindingaffinity. With the rapid growth of computational research, there is an urgent need for more efficientcomputational methods to design viable ligands that target proteins.Method: Recently, a newly developed program called AfDesign (part of ColabDesign) at https://github.com/sokrypton/ColabDesign), an open-source software built on AlphaFold, has beensuggested to be capable of generating ligands that bind to targeted proteins, thus potentially facilitatingthe ligand development process. To evaluate the performance of this program, we exploredits ability to target the PICK1 PDZ domain, given our current understanding of it. We found thatthe designated length of the ligand and the number of recycles play vital roles in generating ligandswith optimal properties.Results: Utilizing AfDesign with a sequence length of 5 for the ligand produced the highest comparableligands to that of prior identified ligands. Moreover, these designed ligands displayed significantlylower binding energy compared to manually created sequences.Conclusion: This work demonstrated that AfDesign can potentially be a powerful tool to facilitatethe exploration of the ligand space for the purpose of targeting PDZ domains. 
    more » « less
  5. null (Ed.)
    Abstract Proteins are tiny players involved in the activation and deactivation of multiple signaling cascades through interactions in cells. The TNFR1 and MADD interact with each other and mediate downstream protein signaling pathways which cause neuronal cell death and Alzheimer’s disease. In the current study, a molecular docking approach was employed to explore the interactive behavior of TNFR1 and MADD proteins and their role in the activation of downstream signaling pathways. The computational sequential and structural conformational results revealed that Asp400, Arg58, Arg59 were common residues of TNFR1 and MADD which are involved in the activation of downstream signaling pathways. Aspartic acid in negatively charged residues is involved in the biosynthesis of protein. However, arginine is a positively charged residue with the potential to interact with oppositely charged amino acids. Furthermore, our molecular dynamic simulation results also ensured the stability of the backbone of TNFR1 and MADD death domains (DDs) in binding interactions. This DDs interaction mediates some conformational changes in TNFR1 which leads to the activation of mediators proteins in the cellular signaling pathways. Taken together, a better understanding of TNFR1 and MADD receptors and their activated signaling cascade may help treat Alzheimer’s disease. The death domains of TNFR1 and MADD could be used as a novel pharmacological target for the treatment of Alzheimer’s disease by inhibiting the MAPK pathway. 
    more » « less