Allosterism in the PDZ Family

SGDAcn is a suppressor for silk gland endoreplication and development

Silkworm silk gland cells undergo multiple rounds of endoreplication, a process in which the genome is duplicated without cell division, leading to cellular polyploidization. This results in the accumulation of genomic DNA, serving as the foundation for rapid silk proteins synthesis. For the first time, we report a previously uncharacterized gene, SGDAcn, in the silkworm silk gland that clustered regularly interspaced short palindromic repeats (CRISPR) / CRISPR-associated nuclease 9-mediated SGDAcn knockout in the posterior silk gland increased cell size and enhanced silk production. SGDAcn knockout facilitated the progression of endoreplication by upregulating the expression of various cyclin genes and promoting energy metabolism, leading to a substantial increase in fibroin gene expression and its transcription factor Dimm, as well as the stimulation of ribosome biogenesis for messenger RNA translation and enhancement of eukaryotic translation initiation factors for protein synthesis. Our findings demonstrate that SGDAcn influences endoreplication, cell growth, and nucleolus size through SGDAcn-EGFR/PI3K/AKT and SGDAcn-NF-κB signaling pathways. Overall, SGDAcn acts as a negative regulator of silk gland development, affecting cell size and protein synthesis, thus modulating silk production. These mechanisms might be conserved in cell growth and progression, making an attractive target for genetic editing to improve silk yield in silkworms and potentially in mammalian cell growth regulation.

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.

Functional protein dynamics in a crystal

Proteins are molecular machines and to understand how they work, we need to understand how they move. New pump-probe time-resolved X-ray diffraction methods open up ways to initiate and observe protein motions with atomistic detail in crystals on biologically relevant timescales. However, practical limitations of these experiments demands parallel development of effective molecular dynamics approaches to accelerate progress and extract meaning. Here, we establish robust and accurate methods for simulating dynamics in protein crystals, a nontrivial process requiring careful attention to equilibration, environmental composition, and choice of force fields. With more than seven milliseconds of sampling of a single chain, we identify critical factors controlling agreement between simulation and experiments and show that simulated motions recapitulate ligand-induced conformational changes. This work enables a virtuous cycle between simulation and experiments for visualizing and understanding the basic functional motions of proteins.

Functional Protein Dynamics in a Crystal

Proteins are molecular machines and to understand how they work, we need to understand how they move. New pump-probe time-resolved X-ray diffraction methods open up ways to initiate and observe protein motions with atomistic detail in crystals on biologically relevant timescales. However, practical limitations of these experiments demands parallel development of effective molecular dynamics approaches to accelerate progress and extract meaning. Here, we establish robust and accurate methods for simulating dynamics in protein crystals, a nontrivial process requiring careful attention to equilibration, environmental composition, and choice of force fields. With more than seven milliseconds of sampling of a single chain, we identify critical factors controlling agreement between simulation and experiments and show that simulated motions recapitulate ligand-induced conformational changes. This work enables a virtuous cycle between simulation and experiments for visualizing and understanding the basic functional motions of proteins.

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.

Investigating the allosteric response of the PICK1 PDZ domain to different ligands with all-atom simulations

The PDZ family is comprised of small modular domains that play critical roles in the allosteric modulation of many cellular signaling processes by binding to the C-terminal tail of different proteins. As dominant modular proteins that interact with a diverse set of peptides, it is of particular interest to explore how different binding partners induce different allosteric effects on the same PDZ domain. Because the PICK1 PDZ domain can bind different types of ligands, it is an ideal test case to answer this question and explore the network of interactions that give rise to dynamic allostery. Here, we use all-atom molecular dynamics simulations to explore dynamic allostery in the PICK1 PDZ domain by modeling two PICK1 PDZ systems: PICK1 PDZ-DAT and PICK1 PDZ-GluR2. Our results suggest that ligand binding to the PICK1 PDZ domain induces dynamic allostery at the αA helix that is similar to what has been observed in other PDZ domains. We found that the PICK1 PDZ-ligand distance is directly correlated with both dynamic changes of the αA helix and the distance between the αA helix and βB strand. Furthermore, our work identifies a hydrophobic core between DAT/GluR2 and I35 as a key interaction in inducing such dynamic allostery. Finally, the unique interaction patterns between different binding partners and the PICK1 PDZ domain can induce unique dynamic changes to the PICK1 PDZ domain. We suspect that unique allosteric coupling patterns with different ligands may play a critical role in how PICK1 performs its biological functions in various signaling networks.

Three Binding Conformations of BIO124 in the Pocket of the PICK1 PDZ Domain

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.

Human Serine Racemase Weakly Binds the Third PDZ Domain of PSD-95

Human serine racemase (hSR) is a pyridoxal-5'-phosphate (PLP)-dependent dimer that catalyzes the formation of D-serine from L-serine, as well as the dehydration of both L- and D-serine to pyruvate and ammonia. As D-serine is a co-agonist of N-methyl-D-aspartate receptors (NMDARs), hSR is a key enzyme in glutamatergic neurotransmission. hSR activity is finely regulated by Mg²⁺, ATP, post-translational modifications, and the interaction with protein partners. In particular, the C-terminus of murine SR binds the third PDZ domain (PDZ3) of postsynaptic density protein 95 (PSD-95), a member of the membrane-associated guanylate kinase (MAGUK) family involved in the trafficking and localization of glutamate receptors. The structural details of the interaction and the stability of the complex have not been elucidated yet. We evaluated the binding of recombinant human PSD-95 PDZ3 to hSR by glutaraldehyde cross-linking, pull-down assays, isothermal titration calorimetry, nuclear magnetic resonance, and enzymatic assays. Overall, a weak interaction was observed, confirming the binding for the human orthologs but supporting the hypothesis that a third protein partner (i.e., stargazin) is required for the regulation of hSR activity by PSD-95 and to stabilize their interaction.