Synthesis and evaluation of bidentate ligands designed to interact with PDZ domains

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.

Emerging Themes in PDZ Domain Signaling: Structure, Function, and Inhibition

Post-synaptic density-95, disks-large and zonula occludens-1 (PDZ) domains are small globular protein-protein interaction domains widely conserved from yeast to humans. They are composed of ∼90 amino acids and form a classical two α-helical/six β-strand structure. The prototypical ligand is the C-terminus of partner proteins; however, they also bind internal peptide sequences. Recent findings indicate that PDZ domains also bind phosphatidylinositides and cholesterol. Through their ligand interactions, PDZ domain proteins are critical for cellular trafficking and the surface retention of various ion channels. In addition, PDZ proteins are essential for neuronal signaling, memory, and learning. PDZ proteins also contribute to cytoskeletal dynamics by mediating interactions critical for maintaining cell-cell junctions, cell polarity, and cell migration. Given their important biological roles, it is not surprising that their dysfunction can lead to multiple disease states. As such, PDZ domain-containing proteins have emerged as potential targets for the development of small molecular inhibitors as therapeutic agents. Recent data suggest that the critical binding function of PDZ domains in cell signaling is more than just glue, and their binding function can be regulated by phosphorylation or allosterically by other binding partners. These studies also provide a wealth of structural and biophysical data that are beginning to reveal the physical features that endow this small modular domain with a central role in cell signaling.

Structure-based design of PDZ ligands as inhibitors of 5-HT(2A) receptor/PSD-95 PDZ1 domain interaction possessing anti-hyperalgesic activity

Disrupting the interaction between the PDZ protein PSD-95 and the C-terminal domain of the 5-HT2A serotonin receptor has been shown to reduce hyperalgesia in a rodent model of neuropathic pain. Here, we designed and synthesized PDZ ligands capable of binding to the first PDZ domain (PDZ1) of the PSD-95 protein and evaluated their biological activity in vitro and in vivo. A series of substituted indoles was identified by docking simulations, and six novel analogues were synthesized. Three analogues displayed strong interactions with the first PDZ domain (PDZ1) of PDZ-95 in (1)H-(15)N heteronuclear single-quantum coherence (HSQC) experiments and two of them were able to inhibit the interaction between PSD-95 and the 5-HT2A receptor in vitro. We identified compound 8b as the analogue able to significantly suppress mechanical hyperalgesia in an experimental model of traumatic neuropathic pain in the rat. This effect was suppressed by the coadministration of the 5-HT2A receptor antagonist M100907, consistent with an inhibitory effect upon 5-HT2A receptor/PSD-95 interaction. Finally, we determined an NMR-restraint driven model structure for the PSD95 PDZ1/8b complex, which confirms that indole 8b binds to the putative PDZ-ligand binding site.

Accidental interaction between PDZ domains and diclofenac revealed by NMR-assisted virtual screening

In silico approaches have become indispensable for drug discovery as well as drug repositioning and adverse effect prediction. We have developed the eF-seek program to predict protein–ligand interactions based on the surface structure of proteins using a clique search algorithm. We have also developed a special protein structure prediction pipeline and accumulated predicted 3D models in the Structural Atlas of the Human Genome (SAHG) database. Using this database, genome-wide prediction of non-peptide ligands for proteins in the human genome was performed, and a subset of predicted interactions including 14 PDZ domains was then confirmed by NMR titration. Surprisingly, diclofenac, a non-steroidal anti-inflammatory drug, was found to be a non-peptide PDZ domain ligand, which bound to 5 of 15 tested PDZ domains. The critical residues for the PDZ–diclofenac interaction were also determined. Pharmacological implications of the accidental PDZ–diclofenac interaction are further discussed.

An in silico analysis of the binding modes and binding affinities of small molecule modulators of PDZ-peptide interactions

Inhibitors of PDZ-peptide interactions have important implications in a variety of biological processes including treatment of cancer and Parkinson’s disease. Even though experimental studies have reported characterization of peptidomimetic inhibitors of PDZ-peptide interactions, the binding modes for most of them have not been characterized by structural studies. In this study we have attempted to understand the structural basis of the small molecule-PDZ interactions by in silico analysis of the binding modes and binding affinities of a set of 38 small molecules with known K_i or K_d values for PDZ2 and PDZ3 domains of PSD-95 protein. These two PDZ domains show differential selectivity for these compounds despite having a high degree of sequence similarity and almost identical peptide binding pockets. Optimum binding modes for these ligands for PDZ2 and PDZ3 domains were identified by using a novel combination of semi-flexible docking and explicit solvent molecular dynamics (MD) simulations. Analysis of the binding modes revealed most of the peptidomimectic ligands which had high K_i or K_d moved away from the peptide binding pocket, while ligands with high binding affinities remained in the peptide binding pocket. The differential specificities of the PDZ2 and PDZ3 domains primarily arise from differences in the conformation of the loop connecting βB and βC strands, because this loop interacts with the N-terminal chemical moieties of the ligands. We have also computed the MM/PBSA binding free energy values for these 38 compounds with both the PDZ domains from multiple 5 ns MD trajectories on each complex i.e. a total of 228 MD trajectories of 5 ns length each. Interestingly, computational binding free energies show good agreement with experimental binding free energies with a correlation coefficient of approximately 0.6. Thus our study demonstrates that combined use of docking and MD simulations can help in identification of potent inhibitors of PDZ-peptide complexes.

Identification of PDZ ligands by docking-based virtual screening for the development of novel analgesic agents

Disrupting the interaction between the PDZ protein, PSD-95, and its target ligands (such as the glutamate NMDA receptor or the serotonin 5-HT2A receptor) was found to reduce hyperalgesia in various models of neuropathic pain. Here, we set out to identify lead molecules which would interact with PSD-95, and hence, would potentially display analgesic activity. We describe the virtual screening of the Asinex and Cambridge databases which together contain almost one million molecules. Using three successive docking filters and visual inspection, we identified three structural classes of molecules and synthesized a potential lead compound from each class. The binding of the molecules with the PDZ domains of PSD-95 was assessed by (1)H-(15)N HSQC NMR experiments. The analgesic activity of the best ligand, quinoline 2, was evaluated in vivo in a model of neuropathic pain and showed promising results.