Simultaneous prediction of binding free energy and specificity for PDZ domain–peptide interactions

Interactions between protein domains and linear peptides underlie many biological processes. Among these interactions, the recognition of C-terminal peptides by PDZ domains is one of the most ubiquitous. In this work, we present a mathematical model for PDZ domain–peptide interactions capable of predicting both affinity and specificity of binding based on X-ray crystal structures and comparative modeling with Rosetta. We developed our mathematical model using a large phage display dataset describing binding specificity for a wild type PDZ domain and 91 single mutants, as well as binding affinity data for a wild type PDZ domain binding to 28 different peptides. Structural refinement was carried out through several Rosetta protocols, the most accurate of which included flexible peptide docking and several iterations of side chain repacking and backbone minimization. Our findings emphasize the importance of backbone flexibility and the energetic contributions of side chain-side chain hydrogen bonds in accurately predicting interactions. We also determined that predicting PDZ domain–peptide interactions became increasingly challenging as the length of the peptide increased in the N-terminal direction. In the training dataset, predicted binding energies correlated with those derived through calorimetry and specificity switches introduced through single mutations at interface positions were recapitulated. In independent tests, our best performing protocol was capable of predicting dissociation constants well within one order of magnitude of the experimental values and specificity profiles at the level of accuracy of previous studies. To our knowledge, this approach represents the first integrated protocol for predicting both affinity and specificity for PDZ domain–peptide interactions.     

A flexible docking procedure for the exploration of peptide binding selectivity to known structures and homology models of PDZ domains

PDZ domains are important scaffolding modules that typically bind to the C-termini of their interaction partners. Several structures of such complexes have been solved, revealing a conserved binding site in the PDZ domain and an extended conformation of the bound peptide. A compendium of information regarding PDZ complexes demonstrates that dissimilar C-terminal peptides bind to the same PDZ domain, and different PDZ domains can bind the same peptides. A detailed understanding of the PDZ-peptide recognition is needed to elucidate this complexity. To this end, we have designed a family of docking protocols for PDZ domains (termed PDZ-DocScheme) that is based on simulated annealing molecular dynamics and rotamer optimization, and is applicable to the docking of long peptides (20-40 rotatable bonds) to both known PDZ structures and to the more complicated problem of homology models of these domains. The resulting protocol reproduces the structures of PDZ complexes with peptides 4-8 amino acids long within 1-2 A from the experimental structure when the docking is performed to the original structure. If the structure of the target PDZ domain is an apo structure or a homology model, the docking protocol yields structures within 3 A in 9 out of 12 test cases. The automated docking procedure PDZ-DocScheme can serve in the generation of a structural context for validation of PDZ domain specificity from mutagenesis and ligand binding data.     

A fluorescence polarization based screening assay for identification of small molecule inhibitors of the PICK1 PDZ domain

PDZ (PSD-95/Discs-large/ZO-1 homology) domains represent putative targets in several diseases including cancer, stroke, addiction and neuropathic pain. Here we describe the application of a simple and fast screening assay based on fluorescence polarization (FP) to identify inhibitors of the PDZ domain in PICK1 (protein interacting with C kinase 1). We screened 43,380 compounds for their ability to inhibit binding of an Oregon Green labeled C-terminal dopamine transporter peptide (OrG-DAT C13) to purified PICK1 in solution. The assay was highly reliable with excellent screening assay parameters (Z'≈0.7 and Z≈0.6). Out of ~200 compounds that reduced FP to less than 80% of the control wells, six compounds were further characterized. The apparent affinities of the compounds were determined in FP competition binding experiments and ranged from ~5.0 μM to ~193 μM. Binding to the PICK1 PDZ domain was confirmed for five of the compounds (CSC-03, CSC-04, CSC-43, FSC-231 and FSC-240) in a non-fluorescence based assay by their ability to inhibit pull-down of PICK1 by a C-terminal DAT GST fusion protein. CSC-03 displayed the highest apparent affinity (5.0 μM) in the FP assay, and was according to fluorescence resonance energy transfer (FRET) experiments capable of inhibiting the interaction between the C-terminus of the GluR2 subunit of the AMPA-type glutamate receptor and PICK1 in live cells. Additional experiments suggested that CSC-03 most likely is an irreversible inhibitor but with specificity for PICK1 since it did not bind three different PDZ domains of PSD-95. Summarized, our data suggest that FP based screening assays might be a widely applicable tool in the search for small molecule inhibitors of PDZ domain interactions.     

Novel rearranged ions observed for

C-terminal rearrangement ions [b(n-1) + H2O] (where n refers to the total number of residues of peptides) are frequently observed for peptides which contain basic amino acid(s), especially arginine, at or near their N termini in low- and high-energy collision-induced dissociation or post-source decay (PSD) spectra. Here we report a novel rearrangement, associated with PSD for serine- or threonine-containing peptides that are susceptible to C-terminal rearrangement. Based on PSD analyses of serine- or threonine-containing bradykinin and its analogs, which have been ethyl-esterified or 18O labeled at their C termini, the [b(k) + H2O] (where k denotes the position adjacent to the left of the Ser/Thr residue) ion is generally thought to be formed by the transfer of the hydroxyl moiety of a serine or threonine residue to the carbonyl group of the residue to its left accompanied by the loss of the remaining C-terminal portion of the peptide. When the Ser/Thr is at or near the C terminus, the present [b(k) + H2O] ion could be formed via two pathways, i.e., the Ser/Thr-related rearrangement and the conventional C-terminal rearrangement, which has been clearly verified by 18O labeling at the C terminus. In addition, the ions which are formally designated as [y(m)b(l) + H2O], where y(m)b(l) denotes a b-type internal ion, are also briefly described.     

A systematic family-wide investigation reveals that ~30% of mammalian PDZ domains engage in PDZ-PDZ interactions

PDZ domains are independently folded modules that typically mediate protein-protein interactions by binding to the C termini of their target proteins. However, in a few instances, PDZ domains have been reported to dimerize with other PDZ domains. To investigate this noncanonical-binding mode further, we used protein microarrays comprising virtually every mouse PDZ domain to systematically query all possible PDZ-PDZ pairs. We then used fluorescence polarization to retest and quantify interactions and coaffinity purification to test biophysically validated interactions in the context of their full-length proteins. Overall, we discovered 37 PDZ-PDZ interactions involving 46 PDZ domains (~30% of all PDZ domains tested), revealing that dimerization is a more frequently used binding mode than was previously appreciated. This suggests that many PDZ domains evolved to form multiprotein complexes by simultaneously interacting with more than one ligand.     

Predicting protein-protein interactions using graph invariants and a neural network

The PDZ domain of proteins mediates a protein-protein interaction by recognizing the hydrophobic C-terminal tail of the target protein. One of the challenges put forth by the DREAM (Discussions on Reverse Engineering Assessment and Methods) 2009 Challenge consists of predicting a position weight matrix (PWM) that describes the specificity profile of five PDZ domains to their target peptides. We consider the primary structures of each of the five PDZ domains as a numerical sequence derived from graph-theoretic models of each of the individual amino acids in the protein sequence. Using available PDZ domain databases to obtain known targets, the graph-theoretic based numerical sequences are then used to train a neural network to recognize their targets. Given the challenge sequences, the target probabilities are computed and a corresponding position weight matrix is derived. In this work we present our method. The results of our method placed second in the DREAM 2009 challenge.     

Design and verification of halogen-bonding system at the complex interface of human fertilization-related MUP PDZ5 domain with CAMK’s C-terminal peptide

Calmodulin-dependent protein kinase (CAMK) is physiologically activated in fertilized human oocytes and is involved in the Ca²⁺ response pathways that link the fertilization calmodulin signal to meiosis resumption and cortical granule exocytosis. The kinase has an unstructured C-terminal tail that can be recognized and bound by the PDZ5 domain of its cognate partner, the multi-PDZ domain protein (MUP). In the current study, we reported a rational biomolecular design of halogen-bonding system at the complex interface of CAMK’s C-terminal peptide with MUP PDZ5 domain by using high-level computational approaches. Four organic halogens were employed as atom probes to explore the structural geometry and energetic property of designed halogen bonds in the PDZ5–peptide complex. It was found that the heavier halogen elements such as bromine Br and iodine I can confer stronger halogen bond but would cause bad atomic contacts and overlaps at the complex interface, while fluorine F cannot form effective halogen bond in the complex. In addition, the halogen substitution at different positions of peptide’s aromatic ring would result in distinct effects on the halogen-bonding system. The computational findings were then verified by using fluorescence analysis; it is indicated that the halogen type and substitution position play critical role in the interaction strength of halogen bonds, and thus the PDZ5–peptide binding affinity can be improved considerably by optimizing their combination.     

Secondary Amino Alcohols: Traceless Cleavable Linkers for Use in Affinity Capture and Release

Capture and release of peptides is often a critical operation in the pathway to discovering materials with novel functions. However, the best methods for efficient capture impede facile release. To overcome this challenge, we report linkers based on secondary amino alcohols for the release of peptides after capture. These amino alcohols are based on serine (seramox) or isoserine (isoseramox) and can be incorporated into peptides during solid‐phase peptide synthesis through reductive amination. Both linkers are quantitatively cleaved within minutes under NaIO₄ treatment. Cleavage of isoseramox produced a native peptide N‐terminus. This linker also showed broad substrate compatibility; incorporation into a synthetic peptide library resulted in the identification of all sequences by nanoLC‐MS/MS. The linkers are cell compatible; a cell‐penetrating peptide that contained this linker was efficiently captured and identified after uptake into cells. These findings suggest that such secondary amino alcohol based linkers might be suitable tools for peptide‐discovery platforms.     

Protein assembly by orthogonal chemical ligation methods

Chemical synthesis harbors the potential to provide ready access to natural proteins as well as to create nonnatural ones. The Staudinger ligation of a peptide containing a C-terminal phosphinothioester with a peptide containing an N-terminal azide gives an amide with no residual atoms. This method for amide bond formation is orthogonal and complementary to other ligation methods. Herein, we describe the first use of the Staudinger ligation to couple peptides on a solid support. The fragment thus produced is used to assemble functional ribonuclease A via native chemical ligation. The synthesis of a protein by this route expands the versatility of chemical approaches to protein production.     

A One‐Pot Chemically Cleavable Bis‐Linker Tether Strategy for the Synthesis of Heterodimeric Peptides

Heterodimeric peptides linked by disulfide bonds are attractive drug targets. However, their chemical assembly can be tedious, time‐consuming, and low yielding. Inspired by the cellular synthesis of pro‐insulin in which the two constituent peptide chains are expressed as a single‐chain precursor separated by a connecting C‐peptide, we have developed a novel chemically cleavable bis‐linker tether which allows the convenient assembly of two peptide chains as a single “pro”‐peptide on the same solid support. Following the peptide cleavage and post‐synthetic modifications, this bis‐linker tether can be removed in one‐step by chemical means. This method was used to synthesize a drug delivery‐cargo conjugate, TAT‐PKCi peptide, and a two‐disulfide bridged heterodimeric peptide, thionin (7‐19)‐(24‐32R), a thionin analogue. To our knowledge, this is the first report of a one‐pot chemically cleavable bis‐linker strategy for the facile synthesis of cross‐bridged two‐chain peptides.     

A general method for designing combinatorial peptide libraries decodable by amino acid analysis

Herein we describe an algorithm for designing combinatorial peptide libraries for split-and-mix synthesis on solid support that are decodable by amino acid analysis (AAA) of the beads. AAA is a standard service analysis available in most biochemical laboratories, and it allows one to control the quality of the peptide on each bead, an important feature that is missing from most library decoding protocols. In the algorithm, each AA is assigned to two variable positions in the sequence grouped in a "unique pair". This arrangement limits sequence design because both the number of unique pairs U (setting the maximum number of variable AA) and the maximum number S of different AA per variable position depend on the peptide length N (U=N(N-1)/2), S=N-1). The method is therefore only suitable for focused libraries. An application example is shown for the selection of peptides with N-terminal proline or hydroxyproline catalyzing an aldol reaction from a combinatorial library of 65536 octapeptides. A simple enumeration program is available to help design combinatorial libraries decodable by amino acid analysis. The method applies to linear and cyclic peptides, can be used for nonnatural building blocks, including beta-amino acids, and should help to explore the vast chemistry of linear and cyclic peptide for catalysis and bioactivity.     

Peptide‐Stabilized,Fluorescent Silver Nanoclusters: Solid‐Phase Synthesis and Screening

Few‐atom silver nanoclusters (AgNCs) can exhibit strong fluorescence; however, they require ligands to prevent aggregation into larger nanoparticles. Fluorescent AgNCs in biopolymer scaffolds have so far mainly been synthesized in solution, and peptides have only found limited use compared to DNA. Herein, we demonstrate how solid‐phase methods can increase throughput dramatically in peptide ligand screening and in initial evaluation of fluorescence intensity and chemical stability of peptide‐stabilized AgNCs (P‐AgNCs). 9‐Fluorenylmethyloxycarbonyl (Fmoc) solid‐phase peptide synthesis on a hydroxymethyl‐benzoic acid (HMBA) polyethylene glycol polyacrylamide copolymer (PEGA) resin enabled on‐resin screening and evaluation of a peptide library, leading to identification of novel peptide‐stabilized, fluorescent AgNCs. Using systematic amino acid substitutions, we synthesized and screened a 144‐member library. This allowed us to evaluate the effect of length, charge, and Cys content in peptides used as ligands for AgNC stabilization. The results of this study will enable future spectroscopic studies of these peptide‐stabilized AgNCs for bioimaging and other applications.     

Fmoc-based synthesis of peptide alpha-thioesters using an aryl hydrazine support

C-Terminal peptide thioesters are key intermediates in the synthesis/semisynthesis of proteins and of cyclic peptides by native chemical ligation. They are prepared by solid-phase peptide synthesis (SPPS) or biosynthetically by protein splicing techniques. Until recently, the chemical synthesis of C-terminal alpha-thioester peptides by SPPS was largely restricted to the use of Boc/Benzyl chemistry due to the poor stability of the thioester bond to the basic conditions required for the deprotection of the N(alpha)-Fmoc group. In the present work, we describe a new method for the SPPS of C-terminal thioesters using Fmoc/t-Bu chemistry. This method is based on the use of an aryl hydrazine linker, which is totally stable to conditions required for Fmoc-SPPS. When the peptide synthesis has been completed, activation of the linker is achieved by mild oxidation. This step converts the acyl hydrazine group into a highly reactive acyl diazene intermediate which reacts with an alpha-amino acid alkyl thioester (H-AA-SR) to yield the corresponding peptide alpha-thioester in good yield. This method has been successfully used to prepare a variety of peptide thioesters, cyclic peptides, and a fully functional Src homology 3 (SH3) protein domain.     

Solid-phase synthesis of peptide and glycopeptide thioesters through side-chain-anchoring strategies

An efficient new strategy for the synthesis of peptide and glycopeptide thioesters is described. The method relies on the side-chain immobilization of a variety of Fmoc-amino acids, protected at their C-termini, on solid supports. Once anchored, peptides were constructed using solid-phase peptide synthesis according to the Fmoc protocol. After unmasking the C-terminal carboxylate, either thiols or amino acid thioesters were coupled to afford, after cleavage, peptide and glycopeptide thioesters in high yields. Using this method a significant proportion of the proteinogenic amino acids could be incorporated as C-terminal amino acid residues, therefore providing access to a large number of potential targets that can serve as acyl donors in subsequent ligation reactions. The utility of this methodology was exemplified in the synthesis of a 28 amino acid glycopeptide thioester, which was further elaborated to an N-terminal fragment of the glycoprotein erythropoietin (EPO) by native chemical ligation.     

Thermodynamic basis for promiscuity and selectivity in protein-protein interactions: PDZ domains, a case study

Like other protein-protein interaction domains, PDZ domains are involved in many key cellular processes. These processes often require that specific multiprotein complexes be assembled, a task that PDZ domains accomplish by binding to specific peptide motifs in target proteins. However, a growing number of experimental studies show that PDZ domains (like other protein-protein interaction domains) can engage in a variety of interactions and bind distinct peptide motifs. Such promiscuity in ligand recognition raises intriguing questions about the molecular and thermodynamic mechanisms that can sustain it. To identify possible sources of promiscuity and selectivity underlying PDZ domain interactions, we performed molecular dynamics simulations of 20 to 25 ns on a set of 12 different PDZ domain complexes (for the proteins PSD-95, Syntenin, Erbin, GRIP, NHERF, Inad, Dishevelled, and Shank). The electrostatic, nonpolar, and configurational entropy binding contributions were evaluated using the MM/PBSA method combined with a quasi-harmonic analysis. The results revealed that PDZ domain interactions are characterized by overwhelmingly favorable nonpolar contributions and almost negligible electrostatic components, a mix that may readily sustain promiscuity. In addition, despite the structural similarity in fold and in recognition modes, the entropic and other dynamical aspects of binding were remarkably variable not only across PDZ domains but also for the same PDZ domain bound to distinct ligands. This variability suggests that entropic and dynamical components can play a role in determining selectivity either of PDZ domain interactions with peptide ligands or of PDZ domain complexes with downstream effectors.     

Selective Alkaline Protease Catalyzed Hydrolysis of Peptide Esters

Procedures for preparing C-terminal free peptides from hydrolysis of its corresponding methyl or benzyl esters catalyzed by alkaline protease has been developed. N-protected peptides having side-chain ester protecting groups or successive hydrophobic amino acid residues in its sequence are hydrolyzed selectively at C-terminal only and leave other bonds (β and γ- ester or peptide bonds) intact. Compounds which cause a side reaction in base mediated saponification could be hydrolyzed safely by this procedure. Products of this hydrolysis are useful intermediates for fragments coupling in the solid phase peptide synthesis.     

A New and Efficient Method for Synthesis of 5′‐Conjugates of Oligonucleotides through Amide‐Bond Formation on Solid Phase

An efficient method for synthesis of oligonucleotide 5′‐conjugates through amide‐bond formation on solid phase is described. Protected oligonucleotides containing a 5′‐carboxylic acid function were obtained by use of a novel non‐nucleosidic phosphoramidite building block, where the carboxylic acid moiety was protected by a 2‐chlorotrityl group. The protecting group is stable to the phosphoramidite coupling conditions used in solid‐phase oligonucleotide assembly, but is easily deprotected by mild acidic treatment. The protecting group may be removed also by ammonolysis. 5′‐Carboxylate‐modified oligonucleotides were efficiently conjugated on solid support under normal peptide‐coupling conditions to various amines or to the N‐termini of small peptides to yield products of high purity. The method is well‐suited in principle for the synthesis of peptide‐oligonucleotide conjugates containing an amide linkage between the 5′‐end of an oligonucleotide and the N‐terminus of a peptide.     

Cysteine/Penicillamine Ligation Independent of Terminal Steric Demands for Chemical Protein Synthesis

The chemical ligation of two unprotected peptides to generate a natural peptidic linkage specifically at the C‐ and N‐termini is a desirable goal in chemical protein synthesis but is challenging because it demands high reactivity and selectivity (chemo‐, regio‐, and stereoselectivity). We report an operationally simple and highly effective chemical peptide ligation involving the ligation of peptides with C‐terminal salicylaldehyde esters to peptides with N‐terminal cysteine/penicillamine. The notable features of this method include its tolerance of steric hinderance from the side groups on either ligating terminus, thereby allowing flexible disconnection at sites that are otherwise difficult to functionalize. In addition, this method can be expanded to selective desulfurization and one‐pot ligation‐desulfurization reactions. The effectiveness of this method was demonstrated by the synthesis of VISTA (216‐311), PD‐1 (192‐288) and Eglin C.     

Scaffold proteins and the regeneration of visual pigments

CRALBP, cellular retinaldehyde-binding protein, is a retinoid-binding protein necessary for efficient regeneration of rod and cone visual pigments. The C terminus of CRALBP binds to the PDZ domains of EBP50/NHERF-1, which in turn bind to ezrin and actin, proteins localized to the apical processes of the retinal pigment epithelium. In this study, we examined structural features associated with the interaction of the two proteins. The C-terminal amino-acid sequence of 11 orthologous CRALBPs is either ENTAL, ENTAF or EDTAL. Peptides ending in each of these sequences inhibited the interaction of CRALBP and EBP50/NHERF-1 with the use of an overlay assay. Molecular modeling showed that both NTAL and NTAF formed similar networks of H bonds with PDZ1 of EBP50/ NHERF-1, and the side chains of both C-terminal Leu and Phe fit into the peptide-binding groove of PDZ1x CRALBP.11-cis-retinal and EBP50/NHERF-1 migrated as single components when analyzed individually by gel filtration and as a complex when mixed together before gel filtration. Complex formation was abolished by preincubation of EBP50/NHERF-1 with peptide EVENTAL. The ligand absorption spectrum of the complex was identical with that of CRALBP x 11-cis-retinal, demonstrating that complex formation did not perturb the ligand-binding domain of CRALBP.     

Uncovering quantitative protein interaction networks for mouse PDZ domains using protein microarrays

One of the principal challenges in systems biology is to uncover the networks of protein-protein interactions that underlie most biological processes. To date, experimental efforts directed at this problem have largely produced only qualitative networks that are replete with false positives and false negatives. Here, we describe a domain-centered approach--compatible with genome-wide investigations--that enables us to measure the equilibrium dissociation constant (K(D)) of recombinant PDZ domains for fluorescently labeled peptides that represent physiologically relevant binding partners. Using a pilot set of 22 PDZ domains, 4 PDZ domain clusters, and 20 peptides, we define a gold standard dataset by determining the K(D) for all 520 PDZ-peptide combinations using fluorescence polarization. We then show that microarrays of PDZ domains identify interactions of moderate to high affinity (K(D) < or = 10 microM) in a high-throughput format with a false positive rate of 14% and a false negative rate of 14%. By combining the throughput of protein microarrays with the fidelity of fluorescence polarization, our domain/peptide-based strategy yields a quantitative network that faithfully recapitulates 85% of previously reported interactions and uncovers new biophysical interactions, many of which occur between proteins that are co-expressed. From a broader perspective, the selectivity data produced by this effort reveal a strong concordance between protein sequence and protein function, supporting a model in which interaction networks evolve through small steps that do not involve dramatic rewiring of the network.