DexDesign: an OSPREY-based algorithm for designing de novo D-peptide inhibitors.

With over 270 unique occurrences in the human genome, peptide-recognizing PDZ domains play a central role in modulating polarization, signaling, and trafficking pathways. Mutations in PDZ domains lead to diseases such as cancer and cystic fibrosis, making PDZ domains attractive targets for therapeutic intervention. D-peptide inhibitors offer unique advantages as therapeutics, including increased metabolic stability and low immunogenicity. Here, we introduce DexDesign, a novel OSPREY-based algorithm for computationally designing de novo D-peptide inhibitors. DexDesign leverages three novel techniques that are broadly applicable to computational protein design: the Minimum Flexible Set, K*-based Mutational Scan, and Inverse Alanine Scan. We apply these techniques and DexDesign to generate novel D-peptide inhibitors of two biomedically important PDZ domain targets: CAL and MAST2. We introduce a framework for analyzing de novo peptides-evaluation along a replication/restitution axis-and apply it to the DexDesign-generated D-peptides. Notably, the peptides we generated are predicted to bind their targets tighter than their targets' endogenous ligands, validating the peptides' potential as lead inhibitors. We also provide an implementation of DexDesign in the free and open source computational protein design software OSPREY.

Dishevelled2 activates WGEF via its interaction with a unique internal peptide motif of the GEF.

The Wnt-planar cell polarity (Wnt-PCP) pathway is crucial in establishing cell polarity during development and tissue homoeostasis. This pathway is found to be dysregulated in many pathological conditions, including cancer and autoimmune disorders. The central event in Wnt-PCP pathway is the activation of Weak-similarity guanine nucleotide exchange factor (WGEF) by the adapter protein Dishevelled (Dvl). The PDZ domain of Dishevelled2 (Dvl2PDZ) binds and activates WGEF by releasing it from its autoinhibitory state. However, the actual Dvl2PDZ binding site of WGEF and the consequent activation mechanism of the GEF have remained elusive. Using biochemical and molecular dynamics studies, we show that a unique "internal-PDZ binding motif" (IPM) of WGEF mediates the WGEF-Dvl2PDZ interaction to activate the GEF. The residues at P2, P0, P-2 and P-3 positions of IPM play an important role in stabilizing the WGEFpep-Dvl2PDZ interaction. Furthermore, MD simulations of modelled Dvl2PDZ-WGEFIPM peptide complexes suggest that WGEF-Dvl2PDZ interaction may differ from the reported Dvl2PDZ-IPM interactions. Additionally, the apo structure of human Dvl2PDZ shows conformational dynamics different from its IPM peptide bound state, suggesting an induced fit mechanism for the Dvl2PDZ-peptide interaction. The current study provides a model for Dvl2 induced activation of WGEF.

Unconventional PDZ Recognition Revealed in α7 nAChR-PICK1 Complexes.

PDZ domains are modular domains that conventionally bind to C terminal or internal motifs of target proteins to control cellular functions through the regulation of protein complex assemblies. Almost all reported structures of PDZ-target protein complexes rely on fragments or peptides as target proteins. No intact target protein complexed with PDZ was structurally characterized. In this study, we used NMR spectroscopy and other biochemistry and biophysics tools to uncover insights into structural coupling between the PDZ domain of protein interacting with C-kinase 1 (PICK1) and α7 nicotinic acetylcholine receptors (α7 nAChR). Notably, the intracellular domains of both α7 nAChR and PICK1 PDZ exhibit a high degree of plasticity in their coupling. Specifically, the MA helix of α7 nAChR interacts with residues lining the canonical binding site of the PICK1 PDZ, while flexible loops also engage in protein-protein interactions. Both hydrophobic and electrostatic interactions mediate the coupling. Overall, the resulting structure of the α7 nAChR-PICK1 complex reveals an unconventional PDZ binding mode, significantly expanding the repertoire of functionally important PDZ interactions.

Solo regulates the localization and activity of PDZ-RhoGEF for actin cytoskeletal remodeling in response to substrate stiffness.

Recent findings indicate that Solo, a RhoGEF, is involved in cellular mechanical stress responses. However, the mechanism of actin cytoskeletal remodeling via Solo remains unclear. Therefore, this study aimed to identify Solo-interacting proteins using the BioID, a proximal-dependent labeling method, and elucidate the molecular mechanisms of function of Solo. We identified PDZ-RhoGEF (PRG) as a Solo-interacting protein. PRG colocalized with Solo in the basal area of cells, depending on Solo localization, and enhanced actin polymerization at the Solo accumulation sites. Additionally, Solo and PRG interaction was necessary for actin cytoskeletal remodeling. Furthermore, the purified Solo itself had little or negligible GEF activity, even its GEF-inactive mutant directly activated the GEF activity of PRG through interaction. Moreover, overexpression of the Solo and PRG binding domains, respectively, had a dominant-negative effect on actin polymerization and actin stress fiber formation in response to substrate stiffness. Therefore, Solo restricts the localization of PRG and regulates actin cytoskeletal remodeling in synergy with PRG in response to the surrounding mechanical environment.

Pseudomonas environmental strain produces a DegQ-derived and PDZ domain containing peptide with protease activity.

In the search of new enzymatic activities with a possible industrial application, we focused on those microorganisms and their molecular mechanisms that allow them to succeed in the environment, particularly in the proteolytic activity and its central role in the microorganisms' successful permanence. The use of highly active serine proteases for industrial applications is a modern need, especially for the formulation of detergents, protein processing, and hair removal from animal skins. This report provides the isolation and identification of a highly proteolytic fragment derived from DegQ produced by a Pseudomonas fluorescens environmental strain isolated from a frog carcass. Zymograms demonstrate that a 10 kDa protein mainly generates the total proteolytic activity of this strain, which is enhanced by the detergent SDS. Mass spectroscopy analysis revealed that the protein derived a couple of peptides, the ones showing the highest coverage belonging to DegQ. Interestingly, this small protein fragment contains a PDZ domain but no obvious residues indicating that it is a protease. Protein model analysis shows that this fragment corresponds to the main PDZ domain from DegQ, and its unique sequence and structure render a proteolytic peptide. The results presented here indicate that a novel DegQ fragment is sufficient for obtaining high protease activity highlighting that the analysis of environmental microorganisms can render new strains or enzymes with helpful biotechnological characteristics.

Structure-based design of peptidomimetic inhibitors of PSD-95 with improved affinity for the PDZ3 domain.

Aberrant brain-derived neurotrophic factor (BDNF) signaling has been proposed to contribute to the pathophysiology of depression and other neurological disorders such as Angelman syndrome. We have previously shown that targeting the tropomyosin receptor kinase B/postsynaptic density protein-95 (PSD-95) nexus in the BDNF signaling pathway by peptidomimetic inhibitors is a promising approach for therapeutic intervention. Here, we used structure-based knowledge to develop a new Syn3 peptidomimetic compound series that fuses peptides derived from the PSD-95-binding protein SynGAP to our prototype compound CN2097. The new compounds target the PSD-95 PDZ3 domain and adjoining αC helix to achieve bivalent binding that results in up to 7-fold stronger affinity compared to CN2097. These compounds were designed to improve CN2097 specificity for the PSD-95 PDZ3 domain, and structure-activity relationship studies were performed to improve their resistance to proteolysis.

A short sequence in the tail of SARS-CoV-2 envelope protein controls accessibility of its PDZ-binding motif to the cytoplasm.

The carboxy-terminal tail of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) envelope protein (E) contains a PDZ-binding motif (PBM) which is crucial for coronavirus pathogenicity. During SARS-CoV-2 infection, the viral E protein is expressed within the Golgi apparatus membrane of host cells with its PBM facing the cytoplasm. In this work, we study the molecular mechanisms controlling the presentation of the PBM to host PDZ (PSD-95/Dlg/ZO-1) domain-containing proteins. We show that at the level of the Golgi apparatus, the PDZ-binding motif of the E protein is not detected by E C-terminal specific antibodies nor by the PDZ domain-containing protein-binding partner. Four alanine substitutions upstream of the PBM in the central region of the E protein tail is sufficient to generate immunodetection by anti-E antibodies and trigger robust recruitment of the PDZ domain-containing protein into the Golgi organelle. Overall, this work suggests that the presentation of the PBM to the cytoplasm is under conformational regulation mediated by the central region of the E protein tail and that PBM presentation probably does not occur at the surface of Golgi cisternae but likely at post-Golgi stages of the viral cycle.

The PDLIM family of actin-associated proteins and their emerging role in membrane trafficking.

The PDZ and LIM domain (PDLIM) proteins are associated with the actin cytoskeleton and have conserved in roles in metazoan actin organisation and function. They primarily function as scaffolds linking various proteins to actin and its binding partner α-actinin via two conserved domains; an N-terminal postsynaptic density 95, discs large and zonula occludens-1 (PDZ) domain, and either single or multiple C-terminal LIN-11, Isl-1 and MEC-3 (LIM) domains in the actinin-associated LIM protein (ALP)- and Enigma-related proteins, respectively. While their role in actin organisation, such as in stress fibres or in the Z-disc of muscle fibres is well known, emerging evidence also suggests a role in actin-dependent membrane trafficking in the endosomal system. This is mediated by a recently identified interaction with the sorting nexin 17 (SNX17) protein, an adaptor for the trafficking complex Commander which is itself intimately linked to actin-directed formation of endosomal recycling domains. In this review we focus on the currently understood structural basis for PDLIM function. The PDZ domains mediate direct binding to distinct classes of PDZ-binding motifs (PDZbms), including α-actinin and other actin-associated proteins, and a highly specific interaction with the type III PDZbm such as the one found in the C-terminus of SNX17. The structures of the LIM domains are less well characterised and how they engage with their ligands is completely unknown. Despite the lack of experimental structural data, we find that recently developed machine learning-based structure prediction methods provide insights into their potential interactions and provide a template for further studies of their molecular functions.

Autoregulation of the LIM kinases by their PDZ domain.

LIM domain kinases (LIMK) are important regulators of actin cytoskeletal remodeling. These protein kinases phosphorylate the actin depolymerizing factor cofilin to suppress filament severing, and are key nodes between Rho GTPase cascades and actin. The two mammalian LIMKs, LIMK1 and LIMK2, contain consecutive LIM domains and a PDZ domain upstream of the C-terminal kinase domain. The roles of the N-terminal regions are not fully understood, and the function of the PDZ domain remains elusive. Here, we determine the 2.0 Å crystal structure of the PDZ domain of LIMK2 and reveal features not previously observed in PDZ domains including a core-facing arginine residue located at the second position of the 'x-Φ-G-Φ' motif, and that the expected peptide binding cleft is shallow and poorly conserved. We find a distal extended surface to be highly conserved, and when LIMK1 was ectopically expressed in yeast we find targeted mutagenesis of this surface decreases growth, implying increased LIMK activity. PDZ domain LIMK1 mutants expressed in yeast are hyperphosphorylated and show elevated activity in vitro. This surface in both LIMK1 and LIMK2 is critical for autoregulation independent of activation loop phosphorylation. Overall, our study demonstrates the functional importance of the PDZ domain to autoregulation of LIMKs.

HCMV UL8 interaction with ß-catenin and DVL2 regulates viral reactivation in CD34+ hematopoietic progenitor cells.

IMPORTANCE: CD34+ hematopoietic progenitor cells (HPCs) are an important cellular reservoir for latent human cytomegalovirus (HCMV). Several HCMV genes are expressed during latency that are involved with the maintenance of the viral genome in CD34+ HPC. However, little is known about the process of viral reactivation in these cells. Here, we describe a viral protein, pUL8, and its interaction and stabilization with members of the Wnt/ß-catenin pathway as an important component of viral reactivation. We further define that pUL8 and ß-catenin interact with DVL2 via a PDZ-binding domain, and loss of UL8 interaction with ß-catenin-DVL2 restricts viral reactivation. Our findings will be instrumental in understanding the molecular processes involved in HCMV reactivation in order to design new antiviral therapeutics.

Identification of small molecule antivirals against HTLV-1 by targeting the hDLG1-Tax-1 protein-protein interaction.

Human T-cell leukemia virus type-1 (HTLV-1) is the first pathogenic retrovirus discovered in human. Although HTLV-1-induced diseases are well-characterized and linked to the encoded Tax-1 oncoprotein, there is currently no strategy to target Tax-1 functions with small molecules. Here, we analyzed the binding of Tax-1 to the human homolog of the drosophila discs large tumor suppressor (hDLG1/SAP97), a multi-domain scaffolding protein involved in Tax-1-transformation ability. We have solved the structures of the PDZ binding motif (PBM) of Tax-1 in complex with the PDZ1 and PDZ2 domains of hDLG1 and assessed the binding of 10 million molecules by virtual screening. Among the 19 experimentally confirmed compounds, one systematically inhibited the Tax-1-hDLG1 interaction in different biophysical and cellular assays, as well as HTLV-1 cell-to-cell transmission in a T-cell model. Thus, our work demonstrates that interactions involving Tax-1 PDZ-domains are amenable to small-molecule inhibition, which provides a framework for the design of targeted therapies for HTLV-1-induced diseases.

Suppression of the gene encoding PDZ domain-containing protein decreases cold tolerance and overwintering survival of the mosquito, Culex pipiens (Culicidae: Diptera).

In diapausing mosquitoes, cold tolerance and prolonged lifespan are important features that are crucial for overwintering success. In the mosquito Culex pipiens, we suggest that PDZ domain-containing protein (PDZ) (post synaptic density protein [PSD95], drosophila disc large tumor suppressor [Dlg1], and zonula occludens-1 protein [zo-1]) domain-containing protein is involved with these diapause features for overwintering survival in Culex mosquitoes. The expression level of pdz was significantly higher in diapausing adult females in the early stage in comparison to their nondiapausing counterparts. Suppression of the gene that encodes PDZ by RNA interference significantly decreased actin accumulation in the midgut of early-stage adult diapausing females. Inhibition of pdz also significantly reduced the survivability of diapausing females which indicates that this protein could play a key role in preserving the midgut tissues during early diapause.

GIPC3 couples to MYO6 and PDZ domain proteins, and shapes the hair cell apical region.

GIPC3 has been implicated in auditory function. Here, we establish that GIPC3 is initially localized to the cytoplasm of inner and outer hair cells of the cochlea and then is increasingly concentrated in cuticular plates and at cell junctions during postnatal development. Early postnatal Gipc3KO/KO mice had mostly normal mechanotransduction currents, but had no auditory brainstem response at 1 month of age. Cuticular plates of Gipc3KO/KO hair cells did not flatten during development as did those of controls; moreover, hair bundles were squeezed along the cochlear axis in mutant hair cells. Junctions between inner hair cells and adjacent inner phalangeal cells were also severely disrupted in Gipc3KO/KO cochleas. GIPC3 bound directly to MYO6, and the loss of MYO6 led to altered distribution of GIPC3. Immunoaffinity purification of GIPC3 from chicken inner ear extracts identified co-precipitating proteins associated with adherens junctions, intermediate filament networks and the cuticular plate. Several of immunoprecipitated proteins contained GIPC family consensus PDZ-binding motifs (PBMs), including MYO18A, which bound directly to the PDZ domain of GIPC3. We propose that GIPC3 and MYO6 couple to PBMs of cytoskeletal and cell junction proteins to shape the cuticular plate.

Discovery of a PDZ Domain Inhibitor Targeting the Syndecan/Syntenin Protein-Protein Interaction: A Semi-Automated "Hit Identification-to-Optimization" Approach.

The rapid identification of early hits by fragment-based approaches and subsequent hit-to-lead optimization represents a challenge for drug discovery. To address this challenge, we created a strategy called "DOTS" that combines molecular dynamic simulations, computer-based library design (chemoDOTS) with encoded medicinal chemistry reactions, constrained docking, and automated compound evaluation. To validate its utility, we applied our DOTS strategy to the challenging target syntenin, a PDZ domain containing protein and oncology target. Herein, we describe the creation of a "best-in-class" sub-micromolar small molecule inhibitor for the second PDZ domain of syntenin validated in cancer cell assays. Key to the success of our DOTS approach was the integration of protein conformational sampling during hit identification stage and the synthetic feasibility ranking of the designed compounds throughout the optimization process. This approach can be broadly applied to other protein targets with known 3D structures to rapidly identify and optimize compounds as chemical probes and therapeutic candidates.

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.

Peptide-binding specificity prediction using fine-tuned protein structure prediction networks.

Peptide-binding proteins play key roles in biology, and predicting their binding specificity is a long-standing challenge. While considerable protein structural information is available, the most successful current methods use sequence information alone, in part because it has been a challenge to model the subtle structural changes accompanying sequence substitutions. Protein structure prediction networks such as AlphaFold model sequence-structure relationships very accurately, and we reasoned that if it were possible to specifically train such networks on binding data, more generalizable models could be created. We show that placing a classifier on top of the AlphaFold network and fine-tuning the combined network parameters for both classification and structure prediction accuracy leads to a model with strong generalizable performance on a wide range of Class I and Class II peptide-MHC interactions that approaches the overall performance of the state-of-the-art NetMHCpan sequence-based method. The peptide-MHC optimized model shows excellent performance in distinguishing binding and non-binding peptides to SH3 and PDZ domains. This ability to generalize well beyond the training set far exceeds that of sequence-only models and should be particularly powerful for systems where less experimental data are available.

Inhibitors against Two PDZ Domains of MDA-9 Suppressed Migration of Breast Cancer Cells.

Melanoma differentiation-associated gene 9 (MDA-9) is a small adaptor protein with tandem PDZ domains that promotes tumor progression and metastasis in various human cancers. However, it is difficult to develop drug-like small molecules with high affinity due to the narrow groove of the PDZ domains of MDA-9. Herein, we identified four novel hits targeting the PDZ1 and PDZ2 domains of MDA-9, namely PI1A, PI1B, PI2A, and PI2B, using a protein-observed nuclear magnetic resonance (NMR) fragment screening method. We also solved the crystal structure of the MDA-9 PDZ1 domain in complex with PI1B and characterized the binding poses of PDZ1-PI1A and PDZ2-PI2A, guided by transferred paramagnetic relaxation enhancement. The protein-ligand interaction modes were then cross-validated by the mutagenesis of the MDA-9 PDZ domains. Competitive fluorescence polarization experiments demonstrated that PI1A and PI2A blocked the binding of natural substrates to the PDZ1 and PDZ2 domains, respectively. Furthermore, these inhibitors exhibited low cellular toxicity, but suppressed the migration of MDA-MB-231 breast carcinoma cells, which recapitulated the phenotype of MDA-9 knockdown. Our work has paved the way for the development of potent inhibitors using structure-guided fragment ligation in the future.

Allosteric Signaling in PDZ Energetic Networks: Embedding Error Analysis.

Allosteric signaling in proteins has been known for some half a century, yet how the signal traverses the protein remains an active area of research. Recently, the importance of electrostatics to achieve long-range signaling has become increasingly appreciated. Our laboratory has been working on developing network approaches to capture such interactions. In this study, we turn our attention to the well-studied allosteric model protein, PDZ. We study the allosteric dynamics on a per-residue basis in key constructs involving the PDZ domain, its allosteric effector, and its peptide ligand. We utilize molecular dynamics trajectories to create the networks for the constructs to explore the allosteric effect by plotting the heat kernel results onto axes defined by principal components. We introduce a new metric to quantitate the volume sampled by a residue in the latent space. We relate our findings to PDZ and the greater field of allostery.

Structural analysis of human papillomavirus E6 interactions with Scribble PDZ domains.

The cell polarity regulator Scribble has been shown to be a critical regulator of the establishment and development of tissue architecture, and its dysregulation promotes or suppresses tumour development in a context-dependent manner. Scribble activity is subverted by numerous viruses. This includes human papillomaviruses (HPVs), who target Scribble via the E6 protein. Binding of E6 from high-risk HPV strains to Scribble via a C-terminal PDZ-binding motif leads to Scribble degradation in vivo. However, the precise molecular basis for Scribble-E6 interactions remains to be defined. We now show that Scribble PDZ1 and PDZ3 are the major interactors of HPV E6 from multiple high-risk strains, with each E6 protein displaying a unique interaction profile. We then determined crystal structures of Scribble PDZ1 and PDZ3 domains in complex with the PDZ-binding motif (PBM) motifs of E6 from HPV strains 16, 18 and 66. Our findings reveal distinct interaction patterns for each E6 PBM motif from a given HPV strain, suggesting that a complex molecular interplay exists that underpins the overt Scribble-HPV E6 interaction and controls E6 carcinogenic potential.

Detrimental Role of PDZ-RhoGEF in Pathological Cardiac Hypertrophy.

BACKGROUND: Postsynaptic density 95/disk-large/ZO-1 Rho guanine nucleotide exchange factor (PDZ-RhoGEF, PRG) functions as a RhoGEF for activated Gα13 and transmits activation signals to downstream signaling pathways in various pathological processes. Although the prohypertrophic effect of activated Gα13 (guanine nucleotide binding protein alpha 13; a heterotrimeric G protein) is well-established, the role of PDZ-RhoGEF in pathological cardiac hypertrophy is still obscure. METHODS: Genetically engineered mice and neonatal rat ventricular myocytes were generated to investigate the function of PRG in pathological myocardial hypertrophy. The prohypertrophic stimuli-induced alternations in the morphology and intracellular signaling were measured in myocardium and neonatal rat ventricular myocytes. Furthermore, multiple molecular methodologies were used to identify the precise molecular mechanisms underlying PDZ-RhoGEF function. RESULTS: Increased PDZ-RhoGEF expression was documented in both hypertrophied hearts and neonatal rat ventricular myocytes. Upon prohypertrophic stimuli, the PDZ-RhoGEF-deficient hearts displayed alleviated cardiomyocyte enlargement and attenuated collagen deposition with improved cardiac function, whereas the adverse hypertrophic responses in hearts and neonatal rat ventricular myocytes were markedly exaggerated by PDZ-RhoGEF overexpression. Mechanistically, RhoA (ras homolog family member A)-dependent signaling pathways may function as the downstream effectors of PDZ-RhoGEF in hypertrophic remodeling, as confirmed by rescue experiments using a RhoA inhibitor and dominant-negative RhoA. Furthermore, PDZ-RhoGEF is associated with activated Gα13 and contributes to Gα13-mediated activation of RhoA-dependent signaling. CONCLUSIONS: Our data provide the first evidence that PDZ-RhoGEF promotes pathological cardiac hypertrophy by linking activated Gα13 to RhoA-dependent signaling pathways. Therefore, PDZ-RhoGEF has the potential to be a diagnostic marker or therapeutic target for pathological cardiac hypertrophy.