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
This content will become publicly available on August 26, 2024
Direct observation of negative cooperativity in a detoxification enzyme at the atomic level by Electron Paramagnetic Resonance spectroscopy and simulation
The catalytic activity of human glutathione S‐transferase A1‐1 (hGSTA1‐1), a homodimeric detoxification enzyme, is dependent on the conformational dynamics of a key C‐terminal helix α9 in each monomer. However, the structural details of how the two monomers interact upon binding of substrates is not well understood and the structure of the ligand‐free state of the hGSTA1‐1 homodimer has not been resolved. Here, we used a combination of electron paramagnetic resonance (EPR) distance measurements and weighted ensemble (WE) simulations to characterize the conformational ensemble of the ligand‐free state at the atomic level. EPR measurements reveal a broad distance distribution between a pair of Cu(II) labels in the ligand‐free state that gradually shifts and narrows as a function of increasing ligand concentration. These shifts suggest changes in the relative positioning of the two α9 helices upon ligand binding. WE simulations generated unbiased pathways for the seconds‐timescale transition between alternate states of the enzyme, leading to the generation of atomically detailed structures of the ligand‐free state. Notably, the simulations provide direct observations of negative cooperativity between the monomers of hGSTA1‐1, which involve the mutually exclusive docking of α9 in each monomer as a lid over the active site. We identify key interactions between residues that lead to this negative cooperativity. Negative cooperativity may be essential for interaction of hGSTA1‐1 with a wide variety of toxic substrates and their subsequent neutralization. More broadly, this work demonstrates the power of integrating EPR distances with WE rare‐events sampling strategy to gain mechanistic information on protein function at the atomic level.
more » « less- Award ID(s):
- 2112871
- NSF-PAR ID:
- 10466021
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- Protein Science
- Volume:
- 32
- Issue:
- 10
- ISSN:
- 0961-8368
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
Abstract Structure and functions of S100 proteins are regulated by two distinct calcium binding EF hand motifs. In this work, we used solution‐state NMR spectroscopy to investigate the cooperativity between the two calcium binding sites and map the allosteric changes at the target binding site. To parse the contribution of the individual calcium binding events, variants of S100A12 were designed to selectively bind calcium to either the EF‐I (N63A) or EF‐II (E31A) loop, respectively. Detailed analysis of the backbone chemical shifts for wildtype protein and its mutants indicates that calcium binding to the canonical EF‐II loop is the principal trigger for the conformational switch between ‘closed’ apo to the ‘open’ Ca2+‐bound conformation of the protein. Elimination of binding in S100‐specific EF‐I loop has limited impact on the calcium binding affinity of the EF‐II loop and the concomitant structural rearrangement. In contrast, deletion of binding in the EF‐II loop significantly attenuates calcium affinity in the EF‐I loop and the structure adopts a ‘closed’ apo‐like conformation. Analysis of experimental amide nitrogen (15N) relaxation rates (
R 1,R 2, and15N–{1H} NOE) and molecular dynamics (MD) simulations demonstrate that the calcium bound state is relatively floppy with pico–nanosecond motions induced in functionally relevant domains responsible for target recognition such as the hinge domain and the C‐terminal residues. Experimental relaxation studies combined with MD simulations show that while calcium binding in the EF‐I loop alone does not induce significant motions in the polypeptide chain, EF‐I regulates fluctuations in the polypeptide in the presence of bound calcium in the EF‐II loop. These results offer novel insights into the dynamic regulation of target recognition by calcium binding and unravels the role of cooperativity between the two calcium binding events in S100A12. -
more » « less
Abstract For decades, researchers have elucidated essential enzymatic functions on the atomic length scale by tracing atomic positions in real-time. Our work builds on possibilities unleashed by mix-and-inject serial crystallography (MISC) at X-ray free electron laser facilities. In this approach, enzymatic reactions are triggered by mixing substrate or ligand solutions with enzyme microcrystals. Here, we report in atomic detail (between 2.2 and 2.7 Å resolution) by room-temperature, time-resolved crystallography with millisecond time-resolution (with timepoints between 3 ms and 700 ms) how the
Mycobacterium tuberculosis enzyme BlaC is inhibited by sulbactam (SUB). Our results reveal ligand binding heterogeneity, ligand gating, cooperativity, induced fit, and conformational selection all from the same set of MISC data, detailing how SUB approaches the catalytic clefts and binds to the enzyme noncovalently before reacting to atrans- enamine. This was made possible in part by the application of singular value decomposition to the MISC data using a program that remains functional even if unit cell parameters change up to 3 Å during the reaction. -
more » « less
Abstract The multidrug resistance (MDR) system effectively expels antibiotics out of bacteria causing serious issues during bacterial infection. In addition to drug, indole, a common metabolic waste of bacteria, is expelled by MDR system of gram‐negative bacteria for their survival. Experimental results suggest that AcrB, one of the key components of MDR system, undergoes large scale conformation changes during the pumping due to proton‐motive process. However, due to extremely short time scale, it is difficult to observe (experimentally) those changes in the AcrB, which might facilitate the pumping process. Molecular simulations can shed light to understand the conformational changes for transport of indole in AcrB. Examination of conformational changes using all‐atom simulation is, however, impractical. Here, we develop a hybrid coarse‐grained force field to study the conformational changes of AcrB in presence of indole in the porter domain of monomer II. Using the coarse‐grained force field, we investigated the conformational changes of AcrB for a number of model systems considering the effect of protonation in aspartic acid (Asp) residues Asp407 and Asp408 in the transmembrane domain of monomer II. Our results show that in the presence of indole, protonation of Asp408 or Asp407 residue causes conformational changes from binding state to extrusion state in monomer II, while remaining two monomers (I and III) approach access state in AcrB protein. We also observed that all three AcrB monomers prefer to go back to access state in the absence of indole. Steered molecular dynamics simulations were performed to demonstrate the feasibility of indole transport mechanism for protonated systems. Identification of indole transport pathway through AcrB can be very helpful in understanding the drug efflux mechanism used by the MDR bacteria.
-
null (Ed.)We present a quantitative study comparing the binding of 4-methoxypyridine, MeOPy, ligand to Co( ii )octaethylporphyrin, CoOEP, at the phenyloctane/HOPG interface and in toluene solution. Scanning tunneling microscopy (STM) was used to study the ligand binding to the porphyrin receptors adsorbed on graphite. Electronic spectroscopy was employed for examining this process in fluid solution. The on surface coordination reaction was completely reversible and followed a simple Langmuir adsorption isotherm. Ligand affinities (or Δ G ) for the binding processes in the two different chemical environments were determined from the respective equilibrium constants. The free energy value of −13.0 ± 0.3 kJ mol −1 for the ligation reaction of MeOPy to CoOEP at the solution/HOPG interface is less negative than the Δ G for cobalt porphyrin complexed to the ligand in solution, −16.8 ± 0.2 kJ mol −1 . This result indicates that the MeOPy–CoOEP complex is more stable in solution than on the surface. Additional thermodynamic values for the formation of the surface ligated species (Δ H c = −50 kJ mol −1 and Δ S c = −120 J mol −1 ) were extracted from temperature dependent STM measurements. Density functional computational methods were also employed to explore the energetics of both the solution and surface reactions. At high concentrations of MeOPy the monolayer was observed to be stripped from the surface. Computational results indicate that this is not because of a reduction in adsorption energy of the MeOPy–CoOEP complex. Nearest neighbor analysis of the MeOPy–CoOEP in the STM images revealed positive cooperative ligand binding behavior. Our studies bring new insights to the general principles of affinity and cooperativity in the ligand–receptor interactions at the solution/solid interface. Future applications of STM will pave the way for new strategies designing highly functional multisite receptor systems for sensing, catalysis, and pharmacological applications.more » « less
