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.     

Truncation,modification, and optimization of MIG6segment 2 peptide to target lung cancer-related EGFR

Human epidermal growth factor receptor (EGFR) plays a central role in the pathological progression and metastasis of lung cancer; the development and clinical application of therapeutic agents that target the receptor provide important insights for new lung cancer therapies. The tumor-suppressor protein MIG6 is a negative regulator of EGFR, which can bind at the activation interface of asymmetric dimer of EGFR kinase domains to disrupt dimerization and then inactivate the kinase (Zhang X. et al. Nature 2007, 450: 741–744). The protein adopts two separated segments, i.e. MIG6^{segment 1} and MIG6^{segment 2}, to directly interact with EGFR. Here, computational modeling and analysis of the intermolecular interaction between EGFR kinase domain and MIG6^{segment 2} peptide revealed that the peptide is folded into a two-stranded β-sheet composed of β-strand 1 and β-strand 2; only the β-strand 2 can directly interact with EGFR activation loop, while leaving β-strand 1 apart from the kinase. A C-terminal island within the β-strand 2 is primarily responsible for peptide binding, which was truncated from the MIG6^{segment 2} and exhibited weak affinity to EGFR kinase domain. Structural and energetic analysis suggested that phosphorylation at residues Tyr394 and Tyr395 of truncated peptide can considerably improve EGFR affinity, and mutation of other residues can further optimize the peptide binding capability. Subsequently, three derivative versions of the truncated peptide, including phosphorylated and dephosphorylated peptides as well as a double-point mutant were synthesized and purified, and their affinities to the recombinant protein of human EGFR kinase domain were determined by fluorescence anisotropy titration. As expected theoretically, the dephosphorylated peptide has no observable binding to the kinase, and phosphorylation and mutation can confer low and moderate affinities to the peptide, respectively, suggesting a good consistence between the computational analysis and experimental assay.     

Preparation of orthogonally protected (2S,3R)-2-amino-3-methyl-4-phosphonobutyric acid (Pmab) as a phosphatase-stable phosphothreonine mimetic and its use in the synthesis of polo-box domain-binding peptides

Reported herein is the first stereoselective synthesis of (2S,3R)-4-[bis-(tert-butyloxy)phosphinyl]-2-[(9H-fluoren-9-ylmethoxy)carbonyl]amino-3-methylbutanoic acid [(N-Fmoc, O,O-(bis-(tert-butyl))-Pmab), 4] as a hydrolytically-stable phosphothreonine mimetic bearing orthogonal protection compatible with standard solid-phase protocols. The synthetic approach used employs Evans' oxazolidinone for chiral induction. Also presented is the application of 4 in the solid-phase synthesis of polo-like kinase 1 (Plk1) polo box domain (PBD)-binding peptides. These Pmab-containing peptides retain PBD binding efficacy similar to a parent pThr containing peptide. Reagent 4 should be a highly useful reagent for the preparation of signal transduction-directed peptides.     

The ATCUN domain as a probe of intermolecular interactions: application to calmodulin-peptide complexes

The utility of a three-residue Cu2+ binding motif (ATCUN domain) for studying intermolecular interactions is demonstrated. By comparing a set of 1H-15N correlation spectra recorded on complexes of calmodulin (CaM) and peptides with the ATCUN tag in the presence and absence of Cu2+ the two possible canonical binding orientations of the peptide can be rapidly distinguished. The methodology is confirmed with studies of complexes of CaM and peptides from myosin light chain kinase and CaM kinase kinase, for which high-resolution structures are available, and then applied to a complex with CaM kinase I for which structural data has not been obtained. The orientation of the CaM kinase I and myosin light chain kinase peptides are shown to be identical. In the case of a complex of CaM with a peptide for which structural information is not available, the present methodology, in combination with 1H-15N residual dipolar couplings measured on CaM, and the database of existing CaM-peptide structures, allows a homology model to be built rapidly and with confidence.     

Revisiting the structural basis and energetic landscape of susceptibility difference between HLA isotypes to allergic rhinitis

The human leukocyte antigen class II (HLA II) molecules are implicated in the immunopathogenesis of allergic rhinitis (AR). The HLA II contains three allelic isotypes HLA-DR, -DQ, and -QP that exhibit considerably different susceptibility to AR. Here, we investigated the structural basis and energetic landscape of the susceptibility difference between the three HLA II isotypes to AR by combining computational analysis and experimental assay. Multiple sequence alignment revealed a low conservation among the three subtypes with sequence identity of ∼10% between them, suggesting that the peptide repertoires presented by HLA-DR, -DP and -DQ are not overlapped to each other, and they may be involved in different immune functions and dysfunctions. Structural analysis imparted that the antigenic peptides are rooted on the peptide-binding groove of HLA molecules and hold in a PPII-like helical conformation. Subsequently, the interaction behavior of 17 AR allergen-derived peptides with HLA-DR, −DP and −DQ was investigated using a statistics-based quantitative structure-activity relationship (QSAR) predictor. It was found a significant difference between the binding capabilities of these antigenic peptides to HLA-DR and to HLA-DP/-DQ; the former showed a generally higher affinity than the latter with p-value of 0.02 obtained from 2-tailed Student’s t-test. The computational findings were then confirmed by HLA II–peptide stability assay, which demonstrated that the AR allergen-derived peptides have a high in vitro selectivity for HLA-DR over HLA-DP/-DQ. Thus, the HLA-DR isotype, rather than HLA-DP and −DQ, is expected to associate with the pathological process of AR.     

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.     

Design and synthesis of peptides that bind alpha-bungarotoxin with high affinity

BACKGROUND: Alpha-bungarotoxin (alpha-BTX) is a highly toxic snake venom alpha-neurotoxin that binds to acetylcholine receptor (AChR) at the neuromuscular junction, and is a potent inhibitor of this receptor. We describe the design and synthesis of peptides that bind alpha-BTX with high affinity, and inhibit its interaction with AChR with an IC(50) of 2 nM. The design of these peptides was based on a lead peptide with an IC(50) of 3x10(-7) M, previously identified by us [M. Balass et al., Proc. Natl. Acad. Sci. USA 94 (1997) 6054] using a phage-display peptide library. RESULTS: Employing nuclear magnetic resonance-derived structural information [T. Scherf et al., Proc. Natl. Acad. Sci. USA 94 (1997) 6059] of the complex of alpha-BTX with the lead peptide, as well as structure-function analysis of the ligand-binding site of AChR, a systematic residue replacement of the lead peptide, one position at a time, yielded 45 different 13-mer peptides. Of these, two peptides exhibited a one order of magnitude increase in inhibitory potency in comparison to the lead peptide. The design of additional peptides, with two or three replacements, resulted in peptides that exhibited a further increase in inhibitory potency (IC(50) values of 2 nM), that is more than two orders of magnitude better than that of the original lead peptide, and better than that of any known peptide derived from AChR sequence. The high affinity peptides had a protective effect on mice against alpha-BTX lethality. CONCLUSIONS: Synthetic peptides with high affinity to alpha-BTX may be used as potential lead compounds for developing effective antidotes against alpha-BTX poisoning. Moreover, the procedure employed in this study may serve as a general approach for the design and synthesis of peptides that interact with high affinity with any desired biological target.     

Discovery of a small-molecule inhibitor of Dvl–CXXC5 interaction by computational approaches

The Wnt/β-catenin signaling pathway plays a significant role in the control of osteoblastogenesis and bone formation. CXXC finger protein 5 (CXXC5) has been recently identified as a negative feedback regulator of osteoblast differentiation through a specific interaction with Dishevelled (Dvl) protein. It was reported that targeting the Dvl–CXXC5 interaction could be a novel anabolic therapeutic target for osteoporosis. In this study, complex structure of Dvl PDZ domain and CXXC5 peptide was simulated with molecular dynamics (MD). Based on the structural analysis of binding modes of MD-simulated Dvl PDZ domain with CXXC5 peptide and crystal Dvl PDZ domain with synthetic peptide–ligands, we generated two different pharmacophore models and applied pharmacophore-based virtual screening to discover potent inhibitors of the Dvl–CXXC5 interaction for the anabolic therapy of osteoporosis. Analysis of 16 compounds selected by means of a virtual screening protocol yielded four compounds that effectively disrupted the Dvl–CXXC5 interaction in the fluorescence polarization assay. Potential compounds were validated by fluorescence spectroscopy and nuclear magnetic resonance. We successfully identified a highly potent inhibitor, BMD4722, which directly binds to the Dvl PDZ domain and disrupts the Dvl–CXXC5 interaction. Overall, CXXC5–Dvl PDZ domain complex based pharmacophore combined with various traditional and simple computational methods is a promising approach for the development of modulators targeting the Dvl–CXXC5 interaction, and the potent inhibitor BMD4722 could serve as a starting point to discover or design more potent and specific the Dvl–CXXC5 interaction disruptors.     

High performance liquid chromatography coupled on-line to capillary electrophoresis with laser-induced fluorescence detection for detecting inhibitors of Src homology 2 domain-phosphopeptide binding in mixtures

Reversed-phase HPLC was coupled on-line to a rapid, competitive affinity probe capillary electrophoresis (APCE) assay to screen mixtures for compounds that inhibit protein-ligand interactions. The Fyn Src homology 2 (SH2) domain and its phosphopeptide binding partner were used as a model interaction for demonstration of this technique. In the method, mixtures containing possible inhibitors of binding were separated by HPLC at a flow rate of 0.3 mL/min. A small portion of effluent was directed to a fluidic tee where it was mixed on-line with Fyn SH2 domain and a fluorescent phosphopeptide ("affinity probe") known to bind selectively to Fyn SH2 domain. Electropherograms of the reaction mixture were collected on-line at approximately 6s intervals using a flow-gated interface to control injections onto the capillary electrophoresis with laser-induced fluorescence system. The resulting electropherograms contained two peaks, one corresponding to the free affinity probe and the other a complex of the affinity probe and Fyn SH2 domain. Compounds that bound the protein were detected as a decrease in the peak height of the complex and an increase in the peak height of affinity probe with relative standard deviations of <5%. The assay was shown to resolve multiple peptidergic inhibitors and selectively detect them within a complex mixture of peptides. Signals were dependent upon both concentration of active peptide and its potency in binding inhibition. Detection limits were in the range of 2-11 microM depending upon the peptide. Common organic solvents used in HPLC were shown to have minimal effect in the on-line measurement up to approximately 60% in the mobile phase.     

Site selection in tandem arrays of metal-binding domains

The metal binding properties of peptides corresponding to metal-binding sites spanning regions that normally function as linkers in tandem arrays of metal-binding domain-containing proteins were examined. For a peptide with two His residues from one TFIIIA-like zinc finger domain, a canonical TFIIIA-like linker, and two Cys residues from an adjacent zinc domain, the dissociation constant for the 1:1 peptide to cobalt(II) was found to be 15 +/- 10 microM, compared with 60 nM for the corresponding zinc finger domains themselves. Peptides overlapping two sets of metal-binding domains from human TRAF (tumor necrosis factor receptor-associated factor) proteins were examined. In one case, the affinity of the presumed metal-binding domain and that for the linker region were comparable, while in the second case, the affinity of the linker peptide was higher than that for the corresponding presumed metal-binding domain peptide. These studies revealed that cobalt(II) affinities in the micromolar range can occur even for peptides that do not correspond to natural zinc-binding domains and that the degree of distinction between authentic metal-binding domains and the corresponding linker-spanning peptides may be modest, at least for single domain peptide models.     

Calorimetric and structural studies of 1,2,3-trisubstituted cyclopropanes as conformationally constrained peptide inhibitors of Src SH2 domain binding.

Isothermal titration calorimetry and X-ray crystallography have been used to determine the structural and thermodynamic consequences associated with constraining the pTyr residue of the pYEEI ligand for the Src Homology 2 domain of the Src kinase (Src SH2 domain). The conformationally constrained peptide mimics that were used are cyclopropane-derived isosteres whereby a cyclopropane ring substitutes to the N-Calpha-Cbeta atoms of the phosphotyrosine. Comparison of the thermodynamic data for the binding of the conformationally constrained peptide mimics relative to their equivalent flexible analogues as well as a native tetrapeptide revealed an entropic advantage of 5-9 cal mol(-1) K(-1) for the binding of the conformationally constrained ligands. However, an unexpected drop in enthalpy for the binding of the conformationally constrained ligands relative to their flexible analogues was also observed. To evaluate whether these differences reflected conformational variations in peptide binding modes, we have determined the crystal structure of a complex of the Src SH2 domain bound to one of the conformationally constrained peptide mimics. Comparison of this new structure with that of the Src SH2 domain bound to a natural 11-mer peptide (Waksman et al. Cell 1993, 72, 779-790) revealed only very small differences. Hence, cyclopropane-derived peptides are excellent mimics of the bound state of their flexible analogues. However, a rigorous analysis of the structures and of the surface areas at the binding interface, and subsequent computational derivation of the energetic binding parameters, failed to predict the observed differences between the binding thermodynamics of the rigidified and flexible ligands, suggesting that the drop in enthalpy observed with the conformationally constrained peptide mimic arises from sources other than changes in buried surface areas, though the exact origin of the differences remains unclear.     

Design of Decorin‐Based Peptides That Bind to Collagen I and their Potential as Adhesion Moieties in Biomaterials

Mimicking the binding epitopes of protein–protein interactions by using small peptides is important for generating modular biomimetic systems. A strategy is described for the design of such bioactive peptides without accessible structural data for the targeted interaction, and the effect of incorporating such adhesion peptides in complex biomaterial systems is demonstrated. The highly repetitive structure of decorin was analyzed to identify peptides that are representative of the inner and outer surface, and it was shown that only peptides based on the inner surface of decorin bind to collagen. The peptide with the highest binding affinity for collagen I, LHERHLNNN, served to slow down the diffusion of a conjugated dye in a collagen gel, while its dimer could physically crosslink collagen, thereby enhancing the elastic modulus of the gel by one order of magnitude. These results show the potential of the identified peptides for the design of biomaterials for applications in regenerative medicine.     

Insight into herbicide resistance of W574L mutant Arabidopsis thaliana acetohydroxyacid synthase: molecular dynamics simulations and binding free energy calculations

Acetohydroxyacid synthase (AHAS) is the target enzyme of several classes of herbicides, such as sulfonylureas and imidazolinones. Now many mutant AHASs with herbicide resistance have emerged along with extensive use of herbicides, therefore it is imperative to understand the detailed interaction mechanism and resistance mechanism so as to develop new potent inhibitors for wild-type or resistant AHAS. With the aid of available crystal structures of the Arabidopsis thaliana (At) AHAS-inhibitor complex, molecular dynamics (MD) simulations were used to investigate the interaction and resistance mechanism directly and dynamically at the atomic level. Nanosecond-level MD simulations were performed on six systems consisting of wild-type or W574L mutant AtAHAS in the complex with three sulfonylurea inhibitors, separately, and binding free energy was calculated for each system using the MM-GBSA method. Comprehensive analyses from structural and energetic aspects confirmed the importance of residue W574, and also indicated that W574L mutation might alert the structural charactersistic of the substrate access channel and decrease the binding affinity of inhibitors, which cooperatively weaken the effective channel-blocked effect and finally result in weaker inhibitory effect of inhibitor and corresponding herbicide resistance of W574L mutant. To our knowledge, it is the first report about MD simulations study on the AHAS-related system, which will pave the way to study the interactions between herbicides and wild-type or mutant AHAS dynamically, and decipher the resistance mechanism at the atomic level for better designing new potent anti-resistance herbicides.     

Interactions and uses of antisense peptides in affinity technology.

Antisense peptides, amino acid sequences encoded in the antisense strand of DNA, can interact with significant affinity and selectivity with their corresponding sensepeptides. Experimentally, sense-antisense peptide recognition has been observed repeatedly. However, skepticism about the biological relevance of this phenomenon has persisted. This is due in part to the unexpected and somewhat couterintutive nature of the interaction as well as to its non-universality as an empirical observation. Nonetheless, antisense peptides in several cases investigated so far have been used as immobilized ligands for the successful affinity chromatographic separation of native (sense) peptides and proteins. For example, immobilized antisense peptides corresponding to Arg8-vasopressin (AVP) have been used to separate vasopressin from oxytocin chromatographically as well as to affinity capture AVP-receptor complex. These results, together with improved understanding of the general features of amino acid sequence which drive antisense-sense peptide interactions as well as new ideas for making antisense peptides chimeras, are beginning to suggest improved ways to make antisense-related peptides as affinity agents for separation as well as for other biotechnology applications.     

Structure‐Based Approach To Improve a Small‐Molecule Inhibitor by the Use of a Competitive Peptide Ligand

Structural information about the target–compound complex is invaluable in the early stage of drug discovery. In particular, it is important to know into which part of the initial compound additional interaction sites could be introduced to improve its affinity. Herein, we demonstrate that the affinity of a small‐molecule inhibitor for its target protein could be successfully improved by the constructive introduction of the interaction mode of a competitive peptide. The strategy involved the discrimination of overlapping and non‐overlapping peptide–compound pharmacophores by the use of a ligand‐based NMR spectroscopic approach, INPHARMA. The obtained results enabled the design of a new compound with improved affinity for the platelet receptor glycoprotein VI (GPVI). The approach proposed herein efficiently combines the advantages of compounds and peptides for the development of higher‐affinity druglike ligands.     

Molecular design and structural optimization of potent peptide hydroxamate inhibitors to selectively target human ADAM metallopeptidase domain 17

Human ADAMs (a disintegrin and metalloproteinases) have been established as an attractive therapeutic target of inflammatory disorders such as inflammatory bowel disease (IBD). The ADAM metallopeptidase domain 17 (ADAM17 or TACE) and its close relative ADAM10 are two of the most important ADAM members that share high conservation in sequence, structure and function, but exhibit subtle difference in regulation of downstream cell signaling events. Here, we described a systematic protocol that combined computational modeling and experimental assay to discover novel peptide hydroxamate derivatives as potent and selective inhibitors for ADAM17 over ADAM10. In the procedure, a virtual combinatorial library of peptide hydroxamate compounds was generated by exploiting intermolecular interactions involved in crystal and modeled structures. The library was examined in detail to identify few promising candidates with both high affinity to ADAM17 and low affinity to ADAM10, which were then tested in vitro with enzyme inhibition assay. Consequently, two peptide hydroxamates Hxm-Phe-Ser-Asn and Hxm-Phe-Arg-Gln were found to exhibit potent inhibition against ADAM17 (K_i = 92 and 47 nM, respectively) and strong selectivity for ADAM17 over ADAM10 (∼7-fold and ∼5-fold, S = 0.86 and 0.71, respectively). The structural basis and energetic property of ADAM17 and ADAM10 interactions with the designed inhibitors were also investigated systematically. It is found that the exquisite network of nonbonded interactions involving the side chains of peptide hydroxamates is primarily responsible for inhibitor selectivity, while the coordination interactions and hydrogen bonds formed by the hydroxamate moiety and backbone of peptide hydroxamates confer high affinity to inhibitor binding.     

Molecular dynamics and thermodynamics of protein-RNA interactions: mutation of a conserved aromatic residue modifies stacking interactions and structural adaptation in the U1A-stem loop 2 RNA complex.

Molecular dynamics (MD) simulations and free energy component analysis have been performed to evaluate the molecular origins of the 5.5 kcal/mol destabilization of the complex formed between the N-terminal RNP domain of U1A and stem loop 2 of U1 snRNA upon mutation of a conserved aromatic residue, Phe56, to Ala. MD simulations, including counterions and water, have been carried out on the wild type and Phe56Ala peptide-stem loop 2 RNA complexes, the free wild type and Phe56Ala peptides, and the free stem loop 2 RNA. The MD structure of the Phe56Ala-stem loop 2 complex is similar to that of the wild type complex except the stacking interaction between Phe56 and A6 of stem loop 2 is absent and loop 3 of the peptide is more dynamic. However, the MD simulations predict large changes in the structure and dynamics of helix C and increased dynamic range of loop 3 for the free Phe56Ala peptide compared to the wild type peptide. Since helix C and loop 3 are highly variable regions of RNP domains, this indicates that a significant contribution to the reduced affinity of the Phe56Ala peptide for RNA results from cooperation between highly conserved and highly variable regions of the RNP domain of U1A. Surprisingly, these structural effects, which are manifested as cooperative free energy changes, occur in the free peptide, rather than in the complex, and are revealed only by study of both the initial and final states of the complexation process. Free energy component analysis correctly accounts for the destabilization of the Phe56Ala-stem loop 2 complex, and indicates that approximately 80% of the destabilization is due to the loss of the stacking interaction and approximately 20% is due to differences in U1A adaptation.     

Structure-based design and optimization of antihypertensive peptides to obtain high inhibitory potency against both renin and angiotensin I-converting enzyme

The human renin–angiotensin system (RAS) plays an essential role in regulating blood pressure and systemic vascular resistance. Renin and angiotensin I-converting enzyme (ACE) are two key enzymes in RAS and have long been recognized as attractive antihypertensive targets. Here, a synthetic strategy was proposed integrating quantitative structure–activity relationship (QSAR), molecular dynamics (MD) simulation and binding free energy analysis to discover novel dual renin and ACE peptidic inhibitors. With the strategy a number of candidates were generated virtually, from which eight promising peptides were selected and synthesized for biological assay. Consequently, three peptides (RYLP, YTAWVP and YRAWVL) were successfully identified to have satisfactory inhibitory profile against both renin and ACE with IC₅₀ values of <1 mM and <10 μM, respectively. Structural analysis and energetic dissection revealed different binding modes of peptide to renin and ACE; a peptide only inserts its C-terminus into the active site of ACE, whereas the whole peptide packs tightly against renin. In addition, when limited to structural diversity it is hard to reconcile the renin and ACE inhibitory activities of short peptides such as dipeptides. These findings can be used to guide peptide optimization with improved biological activity.     

A mini-review and perspective on multicyclic peptide mimics of antibodies

This review gave a brief introduction on recent development in monocyclic and multicyclic peptide mimics of antibodies and provides a perspective on screening and design of multicyclic peptide mimics of antibodies in the future.