Previous studies of the N-terminal PDZ tandem from PSD-95 produced divergent models and failed to identify interdomain contacts stabilizing the structure. We used ensemble and single-molecule FRET along with replica-exchange molecular dynamics to fully characterize the energy landscape. Simulations and experiments identified two conformations: an open-like conformation with a small contact interface stabilized by salt bridges, and a closed-like conformation with a larger contact interface stabilized by surface-exposed hydrophobic residues. Both interfaces were confirmed experimentally. Proximity of interdomain contacts to the binding pockets may explain the observed coupling between conformation and binding. The low-energy barrier between conformations allows submillisecond dynamics, which were time-averaged in previous NMR and FRET studies. Moreover, the small contact interfaces were likely overridden by lattice contacts as crystal structures were rarely sampled in simulations. Our hybrid approach can identify transient interdomain interactions, which are abundant in multidomain proteins yet often obscured by dynamic averaging.
- Award ID(s):
- 2137558
- NSF-PAR ID:
- 10329882
- Date Published:
- Journal Name:
- International Journal of Molecular Sciences
- Volume:
- 23
- Issue:
- 3
- ISSN:
- 1422-0067
- Page Range / eLocation ID:
- 1454
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
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
-
Protein tyrosine phosphatases (PTPs) are emerging drug targets for many diseases, including cancer, autoimmunity, and neurological disorders. A high degree of structural similarity between their catalytic domains, however, has hindered the development of selective pharmacological agents. Our previous research uncovered two unfunctionalized terpenoid inhibitors that selectively inhibit PTP1B over T-cell PTP (TCPTP), two PTPs with high sequence conservation. Here, we use molecular modeling, with supporting experimental validation, to study the molecular basis of this unusual selectivity. Molecular dynamics (MD) simulations suggest that PTP1B and TCPTP share a h-bond network that connects the active site to a distal allosteric pocket; this network stabilizes the closed conformation of the catalytically essential WPD loop, which it links to the L–11 loop and neighboring α3 and α7 helices on the other side of the catalytic domain. Terpenoid binding to either of two proximal C-terminal sites─an α site and a β site─can disrupt the allosteric network; however, binding to the α site forms a stable complex only in PTP1B. In TCPTP, two charged residues disfavor binding at the α site in favor of binding at the β site, which is conserved between the two proteins. Our findings thus indicate that minor amino acid differences at the poorly conserved α site enable selective binding, a property that might be enhanced with chemical elaboration, and illustrate more broadly how minor differences in the conservation of neighboring─yet functionally similar─allosteric sites can affect the selectivity of inhibitory scaffolds (e.g., fragments).more » « less
-
more » « less
We develop integrated co-evolution and dynamic coupling (ICDC) approach to identify, mutate, and assess distal sites to modulate function. We validate the approach first by analyzing the existing mutational fitness data of TEM-1 β-lactamase and show that allosteric positions co-evolved and dynamically coupled with the active site significantly modulate function. We further apply ICDC approach to identify positions and their mutations that can modulate binding affinity in a lectin, cyanovirin-N (CV-N), that selectively binds to dimannose, and predict binding energies of its variants through Adaptive BP-Dock. Computational and experimental analyses reveal that binding enhancing mutants identified by ICDC impact the dynamics of the binding pocket, and show that rigidification of the binding residues compensates for the entropic cost of binding. This work suggests a mechanism by which distal mutations modulate function through dynamic allostery and provides a blueprint to identify candidates for mutagenesis in order to optimize protein function.
-
more » « less
Abstract 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 representativeM. 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 fromM. 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.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.3390/ijms23031454
