Achieving Peptide Binding Specificity and Promiscuity by Loops: Case of the Forkhead-Associated Domain

The regulation of a series of cellular events requires specific protein–protein interactions, which are usually mediated by modular domains to precisely select a particular sequence from diverse partners. However, most signaling domains can bind to more than one peptide sequence. How do proteins create promiscuity from precision? Moreover, these complex interactions typically occur at the interface of a well-defined secondary structure, α helix and β sheet. However, the molecular recognition primarily controlled by loop architecture is not fully understood. To gain a deep understanding of binding selectivity and promiscuity by the conformation of loops, we chose the forkhead-associated (FHA) domain as our model system. The domain can bind to diverse peptides via various loops but only interact with sequences containing phosphothreonine (pThr). We applied molecular dynamics (MD) simulations for multiple free and bound FHA domains to study the changes in conformations and dynamics. Generally, FHA domains share a similar folding structure whereby the backbone holds the overall geometry and the variety of sidechain atoms of multiple loops creates a binding surface to target a specific partner. FHA domains determine the specificity of pThr by well-organized binding loops, which are rigid to define a phospho recognition site. The broad range of peptide recognition can be attributed to different arrangements of the loop interaction network. The moderate flexibility of the loop conformation can help access or exclude binding partners. Our work provides insights into molecular recognition in terms of binding specificity and promiscuity and helpful clues for further peptide design.     

The anchorless adhesin Eap (extracellular adherence protein) from Staphylococcus aureus selectively recognizes extracellular matrix aggregates but binds promiscuously to monomeric matrix macromolecules.

Besides a number of cell wall-anchored adhesins, the majority of Staphylococcus aureus strains produce anchorless, cell wall-associated proteins, such as Eap (extracellular adherence protein). Eap contains four to six tandem repeat (EAP)-domains. Eap mediates diverse biological functions, including adherence and immunomodulation, thus contributing to S. aureus pathogenesis. Eap binding to host macromolecules is unusually promiscuous and includes matrix or matricellular proteins as well as plasma proteins. The structural basis of this promiscuity is poorly understood. Here, we show that in spite of the preferential location of the binding epitopes within triple helical regions in some collagens there is a striking specificity of Eap binding to different collagen types. Collagen I, but not collagen II, is a binding substrate in monomolecular form. However, collagen I is virtually unrecognized by Eap when incorporated into banded fibrils. By contrast, microfibrils containing collagen VI as well as basement membrane-associated networks containing collagen IV, or aggregates containing fibronectin bound Eap as effectively as the monomeric proteins. Therefore, Eap-binding to extracellular matrix ligands is promiscuous at the molecular level but not indiscriminate with respect to supramolecular structures containing the same macromolecules. In addition, Eap bound to banded fibrils after their partial disintegration by matrix-degrading proteinases, including matrix metalloproteinase 1. Therefore, adherence to matrix suprastructures by S. aureus can be supported by inflammatory reactions.     

Drug Promiscuity in PDB: Protein Binding Site Similarity Is Key

Drug repositioning applies established drugs to new disease indications with increasing success. A pre-requisite for drug repurposing is drug promiscuity (polypharmacology) – a drug’s ability to bind to several targets. There is a long standing debate on the reasons for drug promiscuity. Based on large compound screens, hydrophobicity and molecular weight have been suggested as key reasons. However, the results are sometimes contradictory and leave space for further analysis. Protein structures offer a structural dimension to explain promiscuity: Can a drug bind multiple targets because the drug is flexible or because the targets are structurally similar or even share similar binding sites? We present a systematic study of drug promiscuity based on structural data of PDB target proteins with a set of 164 promiscuous drugs. We show that there is no correlation between the degree of promiscuity and ligand properties such as hydrophobicity or molecular weight but a weak correlation to conformational flexibility. However, we do find a correlation between promiscuity and structural similarity as well as binding site similarity of protein targets. In particular, 71% of the drugs have at least two targets with similar binding sites. In order to overcome issues in detection of remotely similar binding sites, we employed a score for binding site similarity: LigandRMSD measures the similarity of the aligned ligands and uncovers remote local similarities in proteins. It can be applied to arbitrary structural binding site alignments. Three representative examples, namely the anti-cancer drug methotrexate, the natural product quercetin and the anti-diabetic drug acarbose are discussed in detail. Our findings suggest that global structural and binding site similarity play a more important role to explain the observed drug promiscuity in the PDB than physicochemical drug properties like hydrophobicity or molecular weight. Additionally, we find ligand flexibility to have a minor influence.     

Assisted peptide folding by surface pattern recognition

Natively disordered proteins belong to a unique class of biomolecules whose function is related to their flexibility and their ability to adopt desired conformations upon binding to substrates. In some cases these proteins can bind multiple partners, which can lead to distinct structures and promiscuity in functions. In other words, the capacity to recognize molecular patterns on the substrate is often essential for the folding and function of intrinsically disordered proteins. Biomolecular pattern recognition is extremely relevant both in vivo (e.g., for oligomerization, immune response, induced folding, substrate binding, and molecular switches) and in vitro (e.g., for biosensing, catalysis, chromatography, and implantation). Here, we use a minimalist computational model system to investigate how polar/nonpolar patterns on a surface can induce the folding of an otherwise unstructured peptide. We show that a model peptide that exists in the bulk as a molten globular state consisting of many interconverting structures can fold into either a helix-coil-helix or an extended helix structure in the presence of a complementary designed patterned surface at low hydrophobicity (3.7%) or a uniform surface at high hydrophobicity (50%). However, we find that a carefully chosen surface pattern can bind to and catalyze the folding of a natively unfolded protein much more readily or effectively than a surface with a noncomplementary or uniform distribution of hydrophobic residues.     

Identification of specificity and promiscuity of PDZ domain interactions through their dynamic behavior

PDZ domains (PDZs), the most common interaction domain proteins, play critical roles in many cellular processes. PDZs perform their job by binding specific protein partners. However, they are very promiscuous, binding to more than one protein, yet selective at the same time. We examined the binding related dynamics of various PDZs to have insight about their specificity and promiscuity. We used full atomic normal mode analysis and a modified coarse‐grained elastic network model to compute the binding related dynamics. In the latter model, we introduced specificity for each single parameter constant and included the solvation effect implicitly. The modified model, referred to as specific‐Gaussian Network Model (s‐GNM), highlights some interesting differences in the conformational changes of PDZs upon binding to Class I or Class II type peptides. By clustering the residue fluctuation profiles of PDZs, we have shown: (i) binding selectivities can be discriminated from their dynamics, and (ii) the dynamics of different structural regions play critical roles for Class I and Class II specificity. s‐GNM is further tested on a dual‐specific PDZ which showed only Class I specificity when a point mutation exists on the βA‐βB loop. We observe that the binding dynamics change consistently in the mutated and wild type structures. In addition, we found that the binding induced fluctuation profiles can be used to discriminate the binding selectivity of homolog structures. These results indicate that s‐GNM can be a powerful method to study the changes in binding selectivities for mutant or homolog PDZs. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Predicting Protein Function and Binding Profile via Matching of Local Evolutionary and Geometric Surface Patterns

Inferring protein functions from structures is a challenging task, as a large number of orphan protein structures from structural genomics project are now solved without their biochemical functions characterized. For proteins binding to similar substrates or ligands and carrying out similar functions, their binding surfaces are under similar physicochemical constraints, and hence the sets of allowed and forbidden residue substitutions are similar. However, it is difficult to isolate such selection pressure due to protein function from selection pressure due to protein folding, and evolutionary relationship reflected by global sequence and structure similarities between proteins is often unreliable for inferring protein function. We have developed a method, called pevoSOAR (pocket-based evolutionary search of amino acid residues), for predicting protein functions by solving the problem of uncovering amino acids residue substitution pattern due to protein function and separating it from amino acids substitution pattern due to protein folding. We incorporate evolutionary information specific to an individual binding region and match local surfaces on a large scale with millions of precomputed protein surfaces to identify those with similar functions. Our pevoSOAR method also generates a probablistic model called the computed binding a profile that characterizes protein-binding activities that may involve multiple substrates or ligands. We show that our method can be used to predict enzyme functions with accuracy. Our method can also assess enzyme binding specificity and promiscuity. In an objective large-scale test of 100 enzyme families with thousands of structures, our predictions are found to be sensitive and specific: At the stringent specificity level of 99.98%, we can correctly predict enzyme functions for 80.55% of the proteins. The overall area under the receiver operating characteristic curve measuring the performance of our prediction is 0.955, close to the perfect value of 1.00. The best Matthews coefficient is 86.6%. Our method also works well in predicting the biochemical functions of orphan proteins from structural genomics projects.     

Selectivity and promiscuity in the interaction network mediated by protein recognition modules

A substantial fraction of protein interactions in the cell is mediated by families of protein modules binding to relatively short linear peptides. Many of these interactions have a high dissociation constant and are therefore suitable for supporting the formation of dynamic complexes that are assembled and disassembled during signal transduction. Extensive work in the past decade has shown that, although member domains within a family have some degree of intrinsic peptide recognition specificity, the derived interaction networks display substantial promiscuity. We review here recent advances in the methods for deriving the portion of the protein network mediated by these domain families and discuss how specific biological outputs could emerge in vivo despite the observed promiscuity in peptide recognition in vitro.     

Structures of designed armadillo repeat proteins binding to peptides fused to globular domains

Designed armadillo repeat proteins (dArmRP) are α‐helical solenoid repeat proteins with an extended peptide binding groove that were engineered to develop a generic modular technology for peptide recognition. In this context, the term “peptide” not only denotes a short unstructured chain of amino acids, but also an unstructured region of a protein, as they occur in termini, loops, or linkers between folded domains. Here we report two crystal structures of dArmRPs, in complex with peptides fused either to the N‐terminus of Green Fluorescent Protein or to the C‐terminus of a phage lambda protein D. These structures demonstrate that dArmRPs bind unfolded peptides in the intended conformation also when they constitute unstructured parts of folded proteins, which greatly expands possible applications of the dArmRP technology. Nonetheless, the structures do not fully reflect the binding behavior in solution, that is, some binding sites remain unoccupied in the crystal and even unexpected peptide residues appear to be bound. We show how these differences can be explained by restrictions of the crystal lattice or the composition of the crystallization solution. This illustrates that crystal structures have to be interpreted with caution when protein–peptide interactions are characterized, and should always be correlated with measurements in solution.     

Efficiency, selectivity, and robustness of nucleocytoplasmic transport

All materials enter or exit the cell nucleus through nuclear pore complexes (NPCs), efficient transport devices that combine high selectivity and throughput. NPC-associated proteins containing phenylalanine–glycine repeats (FG nups) have large, flexible, unstructured proteinaceous regions, and line the NPC. A central feature of NPC-mediated transport is the binding of cargo-carrying soluble transport factors to the unstructured regions of FG nups. Here, we model the dynamics of nucleocytoplasmic transport as diffusion in an effective potential resulting from the interaction of the transport factors with the flexible FG nups, using a minimal number of assumptions consistent with the most well-established structural and functional properties of NPC transport. We discuss how specific binding of transport factors to the FG nups facilitates transport, and how this binding and competition between transport factors and other macromolecules for binding sites and space inside the NPC accounts for the high selectivity of transport. We also account for why transport is relatively insensitive to changes in the number and distribution of FG nups in the NPC, providing an explanation for recent experiments where up to half the total mass of the FG nups has been deleted without abolishing transport. Our results suggest strategies for the creation of artificial nanomolecular sorting devices.     

Nucleocytoplasmic transport: a role for nonspecific competition in karyopherin-nucleoporin interactions

Nucleocytoplasmic transport occurs through the nuclear pore complex (NPC), which in yeast is a ~50 MDa complex consisting of ~30 different proteins. Small molecules can freely exchange through the NPC, but macromolecules larger than ~40 kDa must be aided across by transport factors, most of which belong to a related family of proteins termed karyopherins (Kaps). These transport factors bind to the disordered phenylalanine-glycine (FG) repeat domains in a family of NPC proteins termed FG nups, and this specific binding allows the transport factors to cross the NPC. However, we still know little in terms of the molecular and kinetic details regarding how this binding translates to selective passage of transport factors across the NPC. Here we show that the specific interactions between Kaps and FG nups are strongly modulated by the presence of a cellular milieu whose proteins appear to act as very weak competitors that nevertheless collectively can reduce Kap/FG nup affinities by several orders of magnitude. Without such modulation, the avidities between Kaps and FG nups measured in vitro are too tight to be compatible with the rapid transport kinetics observed in vivo. We modeled the multivalent interactions between the disordered repeat binding sites in the FG nups and multiple cognate binding sites on Kap, showing that they should indeed be sensitive to even weakly binding competitors; the introduction of such competition reduces the availability of these binding sites, dramatically lowering the avidity of their specific interactions and allowing rapid nuclear transport.     

VASP-E: Specificity Annotation with a Volumetric Analysis of Electrostatic Isopotentials

Algorithms for comparing protein structure are frequently used for function annotation. By searching for subtle similarities among very different proteins, these algorithms can identify remote homologs with similar biological functions. In contrast, few comparison algorithms focus on specificity annotation, where the identification of subtle differences among very similar proteins can assist in finding small structural variations that create differences in binding specificity. Few specificity annotation methods consider electrostatic fields, which play a critical role in molecular recognition. To fill this gap, this paper describes VASP-E (Volumetric Analysis of Surface Properties with Electrostatics), a novel volumetric comparison tool based on the electrostatic comparison of protein-ligand and protein-protein binding sites. VASP-E exploits the central observation that three dimensional solids can be used to fully represent and compare both electrostatic isopotentials and molecular surfaces. With this integrated representation, VASP-E is able to dissect the electrostatic environments of protein-ligand and protein-protein binding interfaces, identifying individual amino acids that have an electrostatic influence on binding specificity. VASP-E was used to examine a nonredundant subset of the serine and cysteine proteases as well as the barnase-barstar and Rap1a-raf complexes. Based on amino acids established by various experimental studies to have an electrostatic influence on binding specificity, VASP-E identified electrostatically influential amino acids with 100% precision and 83.3% recall. We also show that VASP-E can accurately classify closely related ligand binding cavities into groups with different binding preferences. These results suggest that VASP-E should prove a useful tool for the characterization of specific binding and the engineering of binding preferences in proteins.     

Dynamics Govern Specificity of a Protein-Protein Interface: Substrate Recognition by Thrombin

Biomolecular recognition is crucial in cellular signal transduction. Signaling is mediated through molecular interactions at protein-protein interfaces. Still, specificity and promiscuity of protein-protein interfaces cannot be explained using simplistic static binding models. Our study rationalizes specificity of the prototypic protein-protein interface between thrombin and its peptide substrates relying solely on binding site dynamics derived from molecular dynamics simulations. We find conformational selection and thus dynamic contributions to be a key player in biomolecular recognition. Arising entropic contributions complement chemical intuition primarily reflecting enthalpic interaction patterns. The paradigm “dynamics govern specificity” might provide direct guidance for the identification of specific anchor points in biomolecular recognition processes and structure-based drug design.     

Structure and Mode of Peptide Binding of Pheromone Receptor PrgZ

We present the crystal structure of the pheromone receptor protein PrgZ from Enterococcus faecalis in complex with the heptapeptide cCF10 (LVTLVFV), which is used in signaling between conjugative recipient and donor cells. Comparison of PrgZ with homologous oligopeptide-binding proteins (AppA and OppA) explains the high specificity of PrgZ for hydrophobic heptapeptides versus the promiscuity of peptide binding in the homologous proteins.     

Protein kinase C binding partners

Members of the protein kinase C family respond to second messengers and are involved in controlling a broad array of cellular functions. The overlapping specificity and promiscuity of these proteins has promoted the view that specific binding proteins constrain individual family members to create the appropriate specificity of action. It is speculated that such protein kinase C-regulator protein interactions affect substrate availability as well as exposure to allosteric activator(s) and that consequent interactions specify cellular location and impose integration with other signaling systems. These predicted features have been realized in the identification of many protein kinase C interacting proteins and examples of these are discussed.     

Ubiquitin-binding domains

Ubiquitin-binding domains (UBDs) are a collection of modular protein domains that non-covalently bind to ubiquitin. These recently discovered motifs interpret and transmit information conferred by protein ubiquitylation to control various cellular events. Detailed molecular structures are known for a number of UBDs, but to understand their mechanism of action, we also need to know how binding specificity is determined, how ubiquitin binding is regulated, and the function of UBDs in the context of full-length proteins. Such knowledge will be key to our understanding of how ubiquitin regulates cellular proteins and processes.     

Bim, Bad and Bmf: intrinsically unstructured BH3-only proteins that undergo a localized conformational change upon binding to prosurvival Bcl-2 targets

All BH3-only proteins, key initiators of programmed cell death, interact tightly with multiple binding partners and have sequences of low complexity, properties that are the hallmark of intrinsically unstructured proteins (IUPs). We show, using spectroscopic methods, that the BH3-only proteins Bim, Bad and Bmf are unstructured in the absence of binding partners. Detailed sequence analyses are consistent with this observation and suggest that most BH3-only proteins are unstructured. When Bim binds and inactivates prosurvival proteins, most residues remain disordered, only the BH3 element becomes structured, and the short alpha-helical molecular recognition element can be considered to behave as a 'bead on a string'. Coupled folding and binding is typical of many IUPs that have important signaling roles, such as BH3-only proteins, as the inherent structural plasticity favors interaction with multiple targets. This understanding offers promise for the development of BH3 mimetics, as multiple modes of binding are tolerated.     

Role of Denatured-State Properties in Chaperonin Action Probed by Single-Molecule Spectroscopy

The bacterial chaperonin GroEL/GroES assists folding of a broad spectrum of denatured and misfolded proteins. Here, we explore the limits of this remarkable promiscuity by mapping two denatured proteins with very different conformational properties, rhodanese and cyclophilin A, during binding and encapsulation by GroEL/GroES with single-molecule spectroscopy, microfluidic mixing, and ensemble kinetics. We find that both proteins bind to GroEL with high affinity in a reaction involving substantial conformational adaptation. However, whereas the compact denatured state of rhodanese is encapsulated efficiently upon addition of GroES and ATP, the more expanded and unstructured denatured cyclophilin A is not encapsulated but is expelled into solution. The origin of this surprising disparity is the weaker interactions of cyclophilin A with a transiently formed GroEL-GroES complex, which may serve as a crucial checkpoint for substrate discrimination.