Binding free energy landscape of domain-peptide interactions

Peptide recognition domains (PRDs) are ubiquitous protein domains which mediate large numbers of protein interactions in the cell. How these PRDs are able to recognize peptide sequences in a rapid and specific manner is incompletely understood. We explore the peptide binding process of PDZ domains, a large PRD family, from an equilibrium perspective using an all-atom Monte Carlo (MC) approach. Our focus is two different PDZ domains representing two major PDZ classes, I and II. For both domains, a binding free energy surface with a strong bias toward the native bound state is found. Moreover, both domains exhibit a binding process in which the peptides are mostly either bound at the PDZ binding pocket or else interact little with the domain surface. Consistent with this, various binding observables show a temperature dependence well described by a simple two-state model. We also find important differences in the details between the two domains. While both domains exhibit well-defined binding free energy barriers, the class I barrier is significantly weaker than the one for class II. To probe this issue further, we apply our method to a PDZ domain with dual specificity for class I and II peptides, and find an analogous difference in their binding free energy barriers. Lastly, we perform a large number of fixed-temperature MC kinetics trajectories under binding conditions. These trajectories reveal significantly slower binding dynamics for the class II domain relative to class I. Our combined results are consistent with a binding mechanism in which the peptide C terminal residue binds in an initial, rate-limiting step.     

Canonical and Noncanonical Sites Determine NPT2A Binding Selectivity to NHERF1 PDZ1

Na⁺/H⁺ Exchanger Regulatory Factor-1 (NHERF1) is a scaffolding protein containing 2 PDZ domains that coordinates the assembly and trafficking of transmembrane receptors and ion channels. Most target proteins harboring a C-terminus recognition motif bind more-or-less equivalently to the either PDZ domain, which contain identical core-binding motifs. However some substrates such as the type II sodium-dependent phosphate co-transporter (NPT2A), uniquely bind only one PDZ domain. We sought to define the structural determinants responsible for the specificity of interaction between NHERF1 PDZ domains and NPT2A. By performing all-atom/explicit-solvent molecular dynamics (MD) simulations in combination with biological mutagenesis, fluorescent polarization (FP) binding assays, and isothermal titration calorimetry (ITC), we found that in addition to canonical interactions of residues at 0 and -2 positions, Arg at the -1 position of NPT2A plays a critical role in association with Glu43 and His27 of PDZ1 that are absent in PDZ2. Experimentally introduced mutation in PDZ1 (Glu43Asp and His27Asn) decreased binding to NPT2A. Conversely, introduction of Asp183Glu and Asn167His mutations in PDZ2 promoted the formation of favorable interactions yielding micromolar K _Ds. The results describe novel determinants within both the PDZ domain and outside the canonical PDZ-recognition motif that are responsible for discrimination of NPT2A between two PDZ domains. The results challenge general paradigms for PDZ recognition and suggest new targets for drug development.     

Comprehensive Assessment of the Relationship Between Site−2 Specificity and Helix α2 in the Erbin PDZ Domain

PDZ domains are key players in signalling pathways. These modular domains generally recognize short linear C-terminal stretches of sequences in proteins that organize the formation of complex multi-component assemblies. The development of new methodologies for the characterization of the molecular principles governing these interactions is critical to fully understand the functional diversity of the family and to elucidate biological functions for family members. Here, we applied an in vitro evolution strategy to explore comprehensively the capacity of PDZ domains for specific recognition of different amino acids at a key position in C-terminal peptide ligands. We constructed a phage-displayed library of the Erbin PDZ domain by randomizing the binding site⁻² and adjacent residues, which are all contained in helix α2, and we selected for variants binding to a panel of peptides representing all possible position⁻² residues. This approach generated insights into the basis for the common natural class I and II specificities, demonstrated an alternative basis for a rare natural class III specificity for Asp⁻², and revealed a novel specificity for Arg⁻² that has not been reported in natural PDZ domains. A structure of a PDZ-peptide complex explained the minimum requirement for switching specificity from class I ligands containing Thr/Ser⁻² to class II ligands containing hydrophobic residues at position⁻². A second structure explained the molecular basis for the specificity for ligands containing Arg⁻². Overall, the evolved PDZ variants greatly expand our understanding of site⁻² specificities and the variants themselves may prove useful as building blocks for synthetic biology.     

Characterization of PDZ domain-peptide interaction interface based on energetic patterns

PDZ domain is one of the abundant modular domains that recognize short peptide sequences to mediate protein-protein interactions. To decipher the binding specificity of PDZ domain, we analyzed the interactions between 11 mouse PDZ domains and 2387 peptides using a method called MIEC-SVM, which energetically characterizes the domain-peptide interaction using molecular interaction energy components (MIECs) and predicts binding specificity using support vector machine (SVM). Cross-validation and leave-one-domain-out test showed that the MIEC-SVM using all 44 PDZ-peptide residue pairs at the interaction interface outperformed the sequence-based methods in the literature. A further feature (residue pair) selection procedure illustrated that 16 residue pairs were uninformative to the binding specificity, even though they contributed significantly (~50%) to the binding energy. If only using the 28 informative residue pairs, the performance of the MIEC-SVM on predicting the PDZ binding specificity was significantly improved. This analysis suggests that the informative and uninformative residue interactions between the PDZ domain and the peptide may represent those contributing to binding specificity and affinity, respectively. We performed additional structural and energetic analyses to shed light on understanding how the PDZ-peptide recognition is established. The success of the MIEC-SVM method on PDZ domains in this study and SH3 domains in our previous studies illustrates its generality on characterizing protein-peptide interactions and understanding protein recognition from a structural and energetic viewpoint.     

HIV-1跨膜蛋白gp41与N-取代吡咯衍生物的结合模式研究

采用分子对接,分子动力学(MD)模拟和分子力学/泊松-波尔兹曼溶剂可有面积方法与分子力学/广义伯恩溶剂可及面积方法(MM-PBSA/MM-GBSA),预测两种N-取代吡咯衍生物与HIV-1 跨膜蛋白gp41疏水口袋的结合模式与作用机理．分子对接采用多种受体构象,并从结果中选取几种可能的结合模式进行MD 模拟,然后通过MM-PBSA计算结合能的方法识别最优的结合模式. MM-PBSA计算结果表明,范德华相互作用是结合的主要驱动力,而极性相互作用决定了配体在结合过程中的取向．进一步的结合能分解显示,配体的羧基与gp41残基Arg579的静电相互作用对结合有重要贡献．上述工作为进一步优化N-取代吡咯衍生物类的HIV-1融合抑制剂建立了良好的理论基础．     

Targeting Protein-Protein Interactions with Trimeric Ligands: High Affinity Inhibitors of the MAGUK Protein Family

PDZ domains in general, and those of PSD-95 in particular, are emerging as promising drug targets for diseases such as ischemic stroke. We have previously shown that dimeric ligands that simultaneously target PDZ1 and PDZ2 of PSD-95 are highly potent inhibitors of PSD-95. However, PSD-95 and the related MAGUK proteins contain three consecutive PDZ domains, hence we envisioned that targeting all three PDZ domains simultaneously would lead to more potent and potentially more specific interactions with the MAGUK proteins. Here we describe the design, synthesis and characterization of a series of trimeric ligands targeting all three PDZ domains of PSD-95 and the related MAGUK proteins, PSD-93, SAP-97 and SAP-102. Using our dimeric ligands targeting the PDZ1-2 tandem as starting point, we designed novel trimeric ligands by introducing a PDZ3-binding peptide moiety via a cysteine-derivatized NPEG linker. The trimeric ligands generally displayed increased affinities compared to the dimeric ligands in fluorescence polarization binding experiments and optimized trimeric ligands showed low nanomolar inhibition towards the four MAGUK proteins, thus being the most potent inhibitors described. Kinetic experiments using stopped-flow spectrometry showed that the increase in affinity is caused by a decrease in the dissociation rate of the trimeric ligand as compared to the dimeric ligands, likely reflecting the lower probability of simultaneous dissociation of all three PDZ ligands. Thus, we have provided novel inhibitors of the MAGUK proteins with exceptionally high affinity, which can be used to further elucidate the therapeutic potential of these proteins.     

Energetics of peptide recognition by the second PDZ domain of human protein tyrosine phosphatase 1E

Formation of protein-protein assemblies is essential in maintaining cell structure and function. Conservation of structural motifs and binding sites is the result of evolutionary pressure for solutions compatible with both molecular economy and regulation. PDZ domains are a typical example: A conserved fold governs specificity toward recognition of C-terminal protein sequences by small sequential and/or structural deviations within a canonical binding mode. The energetic principles underlying the strength and specificity of PDZ-protein interactions are practically unknown. We use the second PDZ domain (PDZ2) of the human protein tyrosine phosphatase (hPTP1E) as a model to study the energetics of peptide binding to a class I PDZ domain. Calorimetric experiments reveal the enthalpy, entropy, and heat capacity changes accompanying PDZ2 binding to the C-terminal pentadecapeptide derived from the guanine nucleotide exchange factor RA-GEF2. Association is driven by favorable enthalpy and entropy changes below 18 degrees C. Above that temperature the entropy change opposes complex formation. Structure-based predictions poorly reproduce the observed thermodynamic profile of the PDZ-peptide complex. On the basis of MD simulations and experimental findings by others we suggest that changes in the dynamics of the PDZ domain upon peptide binding make a large contribution to the observed thermodynamic parameters. Possible impacts of subtle, ligand-induced structural "stiffening" of PDZ domains are discussed. In our hands, the C-terminal segment of the tumor suppressor APC binds much less tightly to PDZ2 than what has been proposed earlier from surface plasmon resonance experiments.     

Analysis of PDZ domain-ligand interactions using carboxyl-terminal phage display

PDZ domains mediate protein-protein interactions at specialized subcellular sites, such as epithelial cell tight junctions and neuronal post-synaptic densities. Because most PDZ domains bind extreme carboxyl-terminal sequences, the phage display method has not been amenable to the study of PDZ domain binding specificities. For the first time, we demonstrate the functional display of a peptide library fused to the carboxyl terminus of the M13 major coat protein. We used this library to analyze carboxyl-terminal peptide recognition by two PDZ domains. For each PDZ domain, the library provided specific ligands with sub-micromolar binding affinities. Synthetic peptides and homology modeling were used to dissect and rationalize the binding interactions. Our results establish carboxyl-terminal phage display as a powerful new method for mapping PDZ domain binding specificity.     

Structure, dynamics and binding characteristics of the second PDZ domain of PTP-BL

The PDZ domains of the protein tyrosine phosphatase PTP-BL mediate interactions by binding to specific amino acid sequences in target proteins. The solution structure of the second PDZ domain of PTP-BL, PDZ2, displays a compact fold with six β strands and two α-helices. A unique feature of this domain compared to the canonical PDZ fold is an extended flexible loop at the base of the binding pocket, termed L1, that folds back onto the protein backbone, a feature that is shared by both the murine and human orthologues. The structure of PDZ2 differs significantly from the orthologous human structure. A comparison of structural quality indicators clearly demonstrates that the PDZ2 ensemble is statistically more reasonable than that of the human orthologue. The analysis of ¹⁵N relaxation data for PDZ2 shows a normal pattern, with more rigid secondary structures and more flexible loop structures. Close to the binding pocket, Leu85 and Thr88 display greater mobility when compared to surrounding residues. Peptide binding studies demonstrated a lack of interaction between murine PDZ2 and the C terminus of the murine Fas/CD95 receptor, suggesting that the Fas/CD95 receptor is not an in vivo target for PDZ2. In addition, PDZ2 specifically binds the C termini of both human Fas/CD95 receptor and the RIL protein, despite RIL containing a non-canonical PDZ-interacting sequence of E-x-V. A model of PDZ2 with the RIL peptide reveals that the PDZ2 binding pocket is able to accommodate the bulkier side-chain of glutamic acid while maintaining crucial protein to peptide hydrogen bond interactions.     

Cooperative phosphoinositide and peptide binding by PSD-95/discs large/ZO-1 (PDZ) domain of polychaetoid, Drosophila zonulin

PDZ domains are well known protein-protein interaction modules that, as part of multidomain proteins, assemble molecular complexes. Some PDZ domains have been reported to interact with membrane lipids, in particular phosphatidylinositol phosphates, but few studies have been aimed at elucidating the prevalence or the molecular details of such interactions. We screened 46 Drosophila PDZ domains for phosphoinositide-dependent cellular localization and discovered that the second PDZ domain of polychaetoid (Pyd PDZ2) interacts with phosphatidylinositol 4,5-bisphosphate (PtdIns(4,5)P(2)) at the plasma membrane. Surface plasmon resonance binding experiments with recombinant protein established that Pyd PDZ2 interacts with phosphatidylinositol phosphates with apparent affinities in the micromolar range. Electrostatic interactions involving an extended positively charged surface of Pyd PDZ2 are crucial for the PtdIns(4,5)P(2)-dependent membrane interactions as shown by a combination of three-dimensional modeling, mutagenesis, binding, and localization studies. In vivo localization studies further suggested that both lipid and peptide binding contribute to membrane localization. We identified the transmembrane protein Crumbs as a Pyd PDZ2 ligand and probed the relation between peptide and PtdIns(4,5)P(2) binding. Contrary to the prevalent view on PDZ/peptide/lipid binding, we did not find competition between peptide and lipid ligands. Instead, preloading the protein with the 10-mer Crb3 peptide increased the apparent affinity of Pyd PDZ2 for PtdIns(4,5)P(2) 6-fold. Our results suggest that membrane localization of Pyd PDZ2 may be driven by a combination of peptide and PtdIns(4,5)P(2) binding, which raises the intriguing possibility that the domain may coordinate protein- and phospholipid-mediated signals.     

The kinetics of PDZ domain-ligand interactions and implications for the binding mechanism

PDZ domains are protein adapter modules present in a few hundred human proteins. They play important roles in scaffolding and signal transduction. PDZ domains usually bind to the C termini of their target proteins. To assess the binding mechanism of this interaction we have performed the first in-solution kinetic study for PDZ domains and peptides corresponding to target ligands. Both PDZ3 from postsynaptic density protein 95 and PDZ2 from protein tyrosine phosphatase L1 bind their respective target peptides through an apparent A + B --> A.B mechanism without rate-limiting conformational changes. But a mutant with a fluorescent probe (Trp) outside of the binding pocket suggests that slight changes in the structure take place upon binding in protein tyrosine phosphatase-L1 PDZ2. For PDZ3 from postsynaptic density protein 95 the pH dependence of the binding reaction is consistent with a one-step mechanism with one titratable group. The salt dependence of the interaction shows that the formation of electrostatic interactions is rate-limiting for the association reaction but not for dissociation of the complex.     

Probing Binding Sites and Mechanisms of Action of an IKs Activator by Computations and Experiments

The slow delayed rectifier (I_Ks) channel is composed of the KCNQ1 channel and KCNE1 auxiliary subunit, and functions to repolarize action potentials in the human heart. I_Ks activators may provide therapeutic efficacy for treating long QT syndromes. Here, we show that a new KCNQ1 activator, ML277, can enhance I_Ks amplitude in adult guinea pig and canine ventricular myocytes. We probe its binding site and mechanism of action by computational analysis based on our recently reported KCNQ1 and KCNQ1/KCNE1 3D models, followed by experimental validation. Results from a pocket analysis and docking exercise suggest that ML277 binds to a side pocket in KCNQ1 and the KCNE1-free side pocket of KCNQ1/KCNE1. Molecular-dynamics (MD) simulations based on the most favorable channel/ML277 docking configurations reveal a well-defined ML277 binding space surrounded by the S2-S3 loop and S4-S5 helix on the intracellular side, and by S4–S6 transmembrane helices on the lateral sides. A detailed analysis of MD trajectories suggests two mechanisms of ML277 action. First, ML277 restricts the conformational dynamics of the KCNQ1 pore, optimizing K⁺ ion coordination in the selectivity filter and increasing current amplitudes. Second, ML277 binding induces global motions in the channel, including regions critical for KCNQ1 gating transitions. We conclude that ML277 activates I_Ks by binding to an intersubunit space and allosterically influencing pore conductance and gating transitions. KCNE1 association protects KCNQ1 from an arrhythmogenic (constitutive current-inducing) effect of ML277, but does not preclude its current-enhancing effect.     

All-Atom Monte Carlo Approach to Protein-Peptide Binding

We develop a procedure for exploring the free energy landscape of protein-peptide binding at atomic detail and apply it to PDZ domain-peptide interactions. The procedure involves soft constraints on receptor proteins providing limited chain flexibility, including backbone motions. Peptide chains are left fully flexible and kept in spatial proximity of the protein through periodic boundary conditions. By extensive Monte Carlo simulations, full representative conformational ensembles at temperatures where bound and unbound states coexist are obtained. To make this approach computationally feasible, we develop an effective all-atom energy function centering on hydrophobicity, hydrogen bonding, and electrostatic interactions. Our initial focus is a set of 11 PDZ domain-peptide pairs with experimentally determined complex structures. Minimum-energy conformations are found to be highly similar to the respective native structures in eight of the cases (all-atom peptide RMSDs < 6 Å). Having achieved that, we turn to a more complete characterization of the bound peptide state through a clustering scheme applied on the full ensembles of peptide structures. We find a significant diversity among bound peptide conformations for several PDZ domains, in particular involving the N terminal side of the peptide chains. Our computational model is then tested further on a set of nine PDZ domain-peptide pairs where the peptides are not originally present in the experimentally determined structures. We find a similar success rate in terms of the nativeness of minimum-energy conformations. Finally, we investigate the ability of our approach to capture variations in binding affinities for different peptide sequences. This is done in particular for a set of related sequences binding to the third PDZ domain of PSD-95 with encouraging results.     

Degenerate specificity of PDZ domains from RhoA-specific nucleotide exchange factors PDZRhoGEF and LARG

PDZ domains are ubiquitous protein-protein interaction modules which bind short, usually carboxyterminal fragments of receptors, other integral or membrane-associated proteins, and occasionally cytosolic proteins. Their role in organizing multiprotein complexes at the cellular membrane is crucial for many signaling pathways, but the rules defining their binding specificity are still poorly understood and do not readily explain the observed diversity of their known binding partners. Two homologous RhoA-specific, multidomain nucleotide exchange factors PDZRhoGEF and LARG contain PDZ domains which show a particularly broad recognition profile, as suggested by the identification of five diverse biological targets. To investigate the molecular roots of this phenomenon, we constructed a phage display library of random carboxyterminal hexapeptides. Peptide variants corresponding to the sequences identified in library selection were synthesized and their affinities for both PDZ domains were measured and compared with those of peptides derived from sequences of natural partners. Based on the analysis of the binding sequences identified for PDZRhoGEF, we propose a sequence for an 'optimal' binding partner. Our results support the hypothesis that PDZ-peptide interactions may be best understood when one considers the sum of entropic and dynamic effects for each peptide as a whole entity, rather than preferences for specific residues at a given position.     

A Method for Dissecting Two Site Receptor Interactions in Ligand Binding Studies

Abstract

A general model has been developed describing the relationship between the measured (IC₅₀) and absolute affinities (K_I), observed in radioligand binding studies when two ligands, one radioactive, interact with two receptors or binding sites. The model shows the dependence of the IC₅₀'s upon the concentration of radioligand for any combinations of the absolute affinities of the radioligand (K_d's) and the displacing ligand (K_I's). By constraining the affinities of the two ligands for the sites, five special cases of the general model can be described that model all possible 'selectivities' the ligands may have for the sites. The properties of these five cases can be exploited experimentally to probe the nature of the ligand/site interactions by the simple expedient of constructing a number of displacement curves at different radioligand concentrations. The method has been tested experimentally in three situations where two ligand/two site interactions occur, and is shown to be a useful technique to qualitatively examine the underlying binding reactions.     

A Structural Portrait of the PDZ Domain Family

PDZ (PSD-95/Discs-large/ZO1) domains are interaction modules that typically bind to specific C-terminal sequences of partner proteins and assemble signaling complexes in multicellular organisms. We have analyzed the existing database of PDZ domain structures in the context of a specificity tree based on binding specificities defined by peptide-phage binding selections. We have identified 16 structures of PDZ domains in complex with high-affinity ligands and have elucidated four additional structures to assemble a structural database that covers most of the branches of the PDZ specificity tree. A detailed comparison of the structures reveals features that are responsible for the diverse specificities across the PDZ domain family. Specificity differences can be explained by differences in PDZ residues that are in contact with the peptide ligands, but these contacts involve both side-chain and main-chain interactions. Most PDZ domains bind peptides in a canonical conformation in which the ligand main chain adopts an extended β-strand conformation by interacting in an antiparallel fashion with a PDZ β-strand. However, a subset of PDZ domains bind peptides with a bent main-chain conformation and the specificities of these non-canonical domains could not be explained based on canonical structures. Our analysis provides a structural portrait of the PDZ domain family, which serves as a guide in understanding the structural basis for the diverse specificities across the family.     

High‐energy water sites determine peptide binding affinity and specificity of PDZ domains

PDZ domains have well known binding preferences for distinct C‐terminal peptide motifs. For most PDZ domains, these motifs are of the form [S/T]‐W‐[I/L/V]. Although the preference for S/T has been explained by a specific hydrogen bond interaction with a histidine in the PDZ domain and the (I/L/V) is buried in a hydrophobic pocket, the mechanism for Trp specificity at the second to last position has thus far remained unknown. Here, we apply a method to compute the free energies of explicit water molecules and predict that potency gained by Trp binding is due to a favorable release of high‐energy water molecules into bulk. The affinities of a series of peptides for both wild‐type and mutant forms of the PDZ domain of Erbin correlate very well with the computed free energy of binding of displaced waters, suggesting a direct relationship between water displacement and peptide affinity. Finally, we show a correlation between the magnitude of the displaced water free energy and the degree of Trp‐sensitivity among subtypes of the HTRA PDZ family, indicating a water‐mediated mechanism for specificity of peptide binding.     

Peptide binding to the PDZ3 domain by conformational selection

The PDZ domains, a large family of peptide recognition proteins, bind to the C‐terminal segment of membrane ion channels and receptors thereby mediating their localization. The peptide binding process is not known in detail and seems to differ among different PDZ domains. For the third PDZ domain of the synaptic protein PSD‐95 (PDZ3), a lock‐and‐key mechanism was postulated on the basis of the almost perfect overlap of the crystal structures in the presence and absence of its peptide ligand. Here, peptide binding to PDZ3 is investigated by explicit solvent molecular dynamics (MD) simulations (for a total of 1.3 μs) and the cut‐based free energy profile method for determining free energy barriers and basins. The free energy landscape of apo PDZ3 indicates that there are multiple basins within the native state. These basins differ by the relative orientation of the α2 helix and β2 strand, the two secondary structure elements that make up the peptide binding site. Only the structure with the smallest aperture of the binding site is populated in the MD simulations of the complex whose analysis reveals that the peptide ligand binds to PDZ3 by selecting one of three conformations. Thus, the dynamical information obtained by the atomistic simulations increment the static, that is, partial, picture of the PDZ3 binding mechanism based on the X‐ray crystallography data. Importantly, the simulation results show for the first time that conformational selection is a possible mechanism of peptide binding by PDZ domains in general. Proteins 2012. © 2012 Wiley Periodicals, Inc.     

Phage display can select over-hydrophobic sequences that may impair prediction of natural domain-peptide interactions

MOTIVATION: The phage display peptide selection approach is widely used for defining binding specificities of globular domains. PDZ domains recognize partner proteins via C-terminal motifs and are often used as a model for interaction predictions. Here, we investigated to which extent phage display data that were recently published for 54 human PDZ domains can be applied to the prediction of human PDZ-peptide interactions. RESULTS: Promising predictions were obtained for one-third of the 54 PDZ domains. For the other two-thirds, we detected in the phage display peptides an important bias for hydrophobic amino acids that seemed to impair correct predictions. Therefore, phage display-selected peptides may be over-hydrophobic and of high affinity, while natural interaction motifs are rather hydrophilic and mostly combine low affinity with high specificity. We suggest that potential amino acid composition bias should systematically be investigated when applying phage display data to the prediction of specific natural domain-linear motif interactions.