skip to main content

This content will become publicly available on August 1, 2023

Title: 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 more » 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. « less
; ;
Award ID(s):
Publication Date:
Journal Name:
Page Range or eLocation-ID:
Sponsoring Org:
National Science Foundation
More Like this
  1. 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, andmore »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.« less
  2. Abstract

    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.

  3. Abstract New drug production, from target identification to marketing approval, takes over 12 years and can cost around $2.6 billion. Furthermore, the COVID-19 pandemic has unveiled the urgent need for more powerful computational methods for drug discovery. Here, we review the computational approaches to predicting protein–ligand interactions in the context of drug discovery, focusing on methods using artificial intelligence (AI). We begin with a brief introduction to proteins (targets), ligands (e.g. drugs) and their interactions for nonexperts. Next, we review databases that are commonly used in the domain of protein–ligand interactions. Finally, we survey and analyze the machine learning (ML) approaches implemented to predict protein–ligand binding sites, ligand-binding affinity and binding pose (conformation) including both classical ML algorithms and recent deep learning methods. After exploring the correlation between these three aspects of protein–ligand interaction, it has been proposed that they should be studied in unison. We anticipate that our review will aid exploration and development of more accurate ML-based prediction strategies for studying protein–ligand interactions.
  4. Mechanistic understanding of DNA recombination in the Cre-loxPsystem has largely been guided by crystallographic structures of tetrameric synaptic complexes. Those studies have suggested a role for protein conformational dynamics that has not been well characterized at the atomic level. We used solution nuclear magnetic resonance (NMR) spectroscopy to discover the link between intrinsic flexibility and function in Cre recombinase. Transverse relaxation-optimized spectroscopy (TROSY) NMR spectra show the N-terminal and C-terminal catalytic domains (CreNTDand CreCat) to be structurally independent. Amide15N relaxation measurements of the CreCatdomain reveal fast-timescale dynamics in most regions that exhibit conformational differences in active and inactive Cre protomers in crystallographic tetramers. However, the C-terminal helix αN, implicated in assembly of synaptic complexes and regulation of DNA cleavage activity viatransprotein–protein interactions, is unexpectedly rigid in free Cre. Chemical shift perturbations and intra- and intermolecular paramagnetic relaxation enhancement (PRE) NMR data reveal an alternative autoinhibitory conformation for the αN region of free Cre, wherein it packsin cisover the protein DNA binding surface and active site. Moreover, binding toloxPDNA induces a conformational change that dislodges the C terminus, resulting in acis-to-transswitch that is likely to enable protein–protein interactions required for assembly of recombinogenic Cre intasomes. These findings necessitate a reexamination ofmore »the mechanisms by which this widely utilized gene-editing tool selects target sites, avoids spurious DNA cleavage activity, and controls DNA recombination efficiency.

    « less
  5. Thiamine pyrophosphate (TPP) riboswitches regulate thiamine metabolism by inhibiting the translation of enzymes essential to thiamine synthesis pathways upon binding to thiamine pyrophosphate in cells across all domains of life. Recent work on the Arabidopsis thaliana TPP riboswitch suggests a multistep TPP binding process involving multiple riboswitch configurational ensembles and Mg 2+ dependence underlies the mechanism of TPP recognition and subsequent transition to the expression-inhibiting state of the aptamer domain followed by changes in the expression platform. However, details of the relationship between TPP riboswitch conformational changes and interactions with TPP and Mg 2+ in the aptamer domain constituting this mechanism are unknown. Therefore, we integrated single-molecule multiparameter fluorescence and force spectroscopy with atomistic molecular dynamics simulations and found that conformational transitions within the aptamer domain's sensor helices associated with TPP and Mg 2+ ligand binding occurred between at least five different ensembles on timescales ranging from µs to ms. These dynamics are orders of magnitude faster than the 10 sec-timescale folding kinetics associated with expression-state switching in the switch helix. Together, our results show that a TPP and Mg 2+ dependent mechanism determines dynamic configurational state ensemble switching of the aptamer domain's sensor helices that regulate the switch helix'smore »stability, which ultimately may lead to the expression-inhibiting state of the riboswitch. Additionally, we propose that two pathways exist for ligand recognition and that this mechanism underlies a kinetic rheostat-like behavior of the Arabidopsis thaliana TPP riboswitch.« less