Change in protein-ligand specificity through binding pocket grafting

Recognition and discrimination of small molecules are crucial for biological processes in living systems. Understanding the mechanisms that underlie binding specificity is of particular interest to synthetic biology, e.g. the engineering of biosensors with de novo ligand affinities. Promising scaffolds for such biosensors are the periplasmic binding proteins (PBPs) due to their ligand-mediated structural change that can be translated into a physically measurable signal.In this study we focused on the two homologous polyamine binding proteins PotF and PotD. Despite their structural similarity, PotF and PotD have different binding specificities for the polyamines putrescine and spermidine. To elucidate how specificity is determined, we grafted the binding site of PotD onto PotF. The introduction of 7 mutations in the first shell of the binding pocket leads to a swap in the binding profile as confirmed by isothermal titration calorimetry. Furthermore, the 1.7 Å crystal structure of the new variant complexed with spermidine reveals the interactions of the specificity determining residues including a defined water network. Altogether our study shows that specificity is encoded in the first shell residues of the PotF binding pocket and that transplantation of these residues allows the swap of the binding specificity.     

The role of disorder in RNA binding affinity and specificity

Most RNA-binding modules are small and bind few nucleotides. RNA-binding proteins typically attain the physiological specificity and affinity for their RNA targets by combining several RNA-binding modules. Here, we review how disordered linkers connecting RNA-binding modules govern the specificity and affinity of RNA–protein interactions by regulating the effective concentration of these modules and their relative orientation. RNA-binding proteins also often contain extended intrinsically disordered regions that mediate protein–protein and RNA–protein interactions with multiple partners. We discuss how these regions can connect proteins and RNA resulting in heterogeneous higher-order assemblies such as membrane-less compartments and amyloid-like structures that have the characteristics of multi-modular entities. The assembled state generates additional RNA-binding specificity and affinity properties that contribute to further the function of RNA-binding proteins within the cellular environment.     

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.     

Origins of specificity and promiscuity in metabolic networks

How enzymes have evolved to their present form is linked to the question of how pathways emerged and evolved into extant metabolic networks. To investigate this mechanism, we have explored the chemical diversity present in a largely unbiased data set of catalytic reactions processed by modern enzymes across the tree of life. In order to get a quantitative estimate of enzyme chemical diversity, we measure enzyme multispecificity or promiscuity using the reaction molecular signatures. Our main finding is that reactions that are catalyzed by a highly specific enzyme are shared by poorly divergent species, suggesting a later emergence of this function during evolution. In contrast, reactions that are catalyzed by highly promiscuous enzymes are more likely to appear uniformly distributed across species in the tree of life. From a functional point of view, promiscuous enzymes are mainly involved in amino acid and lipid metabolisms, which might be associated with the earliest form of biochemical reactions. In this way, results presented in this paper might assist us with the identification of primeval promiscuous catalytic functions contributing to life's minimal metabolism.     

Between promiscuity and specificity: novel roles of EF-hand calcium sensors in neuronal Ca2+ signalling

In recent years, substantial progress has been made towards an understanding of the physiological function of EF-hand calcium sensor proteins of the Calmodulin (CaM) superfamily in neurons. This deeper appreciation is based on the identification of novel target interactions, structural studies and the discovery of novel signalling mechanisms in protein trafficking and synaptic plasticity, in which CaM-like sensor proteins appear to play a role. However, not all interactions are of plausible physiological relevance and in many cases it is not yet clear how the CaM signaling network relates to the proposed function of other EF-hand sensors. In this review, we will summarize these findings and address some of the open questions on the functional role of EF-hand calcium binding proteins in neurons.     

A consensus motif in the RFX DNA binding domain and binding domain mutants with altered specificity.

The RFX DNA binding domain is a novel motif that has been conserved in a growing number of dimeric DNA-binding proteins, having diverse regulatory functions, in eukaryotic organisms ranging from yeasts to humans. To characterize this novel motif, we have performed a detailed dissection of the site-specific DNA binding activity of RFX1, a prototypical member of the RFX family. First, we have performed a site selection procedure to define the consensus binding site of RFX1. Second, we have developed a new mutagenesis-selection procedure to derive a precise consensus motif, and to test the accuracy of a secondary structure prediction, for the RFX domain. Third, a modification of this procedure has allowed us to isolate altered-specificity RFX1 mutants. These results should facilitate the identification both of additional candidate genes controlled by RFX1 and of new members of the RFX family. Moreover, the altered-specificity RFX1 mutants represent valuable tools that will permit the function of RFX1 to be analyzed in vivo without interference from the ubiquitously expressed endogenous protein. Finally, the simplicity, efficiency, and versatility of the selection procedure we have developed make it of general value for the determination of consensus motifs, and for the isolation of mutants exhibiting altered functional properties, for large protein domains involved in protein-DNA as well as protein-protein interactions.     

The relationship between zeros and factors of binding polynomials and cooperativity in protein-ligand binding

Cooperativity in the protein-ligand binding process is discussed in terms of the zeros of the binding polynomial and the corresponding possible factorizations of the binding polynomial into polynomials having non-negative coefficients. Particular attention is paid to the case in which the real parts of all zeros are negative (Hurwitz polynomial) and the case in which the binding polynomial admits no positive factorization (positive irreducible polynomial). Such factorizations are then interpreted as site linkage patterns and related to cooperativity. The possible combinations of zeros of the binding polynomials for the MWC and KNF tetrahedral, square and linear models are determined and the corresponding factorization and linkage patterns analyzed. An application and interpretation are then made for data obtained from Trout I hemoglobin.     

Multiple roles of a conserved GAF domain tyrosine residue in cyanobacterial and plant phytochromes

The phytochrome family of red/far-red photoreceptors has been optimized to support photochemical isomerization of a bound bilin chromophore, a process that triggers a conformational change and modulates biochemical output from the surrounding protein scaffold. Recent studies have established that the efficiency of this photochemical process is profoundly altered by mutation of a conserved tyrosine residue (Tyr176) within the bilin-binding GAF domain of the cyanobacterial phytochrome Cph1 [Fischer, A. J., and Lagarias, J. C. (2004) Harnessing phytochrome's glowing potential, Proc. Natl. Acad. Sci. U.S.A. 101, 17334-17339]. Here, we show that the equivalent mutation in plant phytochromes behaves similarly, indicating that the function of this tyrosine in the primary photochemical mechanism is conserved. Saturation mutagenesis of Tyr176 in Cph1 establishes that no other residue can support comparably efficient photoisomerization. The spectroscopic consequences of Tyr176 mutations also reveal that Tyr176 regulates the conversion of the porphyrin-like conformation of the bilin precursor to a more extended conformation. The porphyrin-binding ability of the Tyr176Arg mutant protein indicates that Tyr176 also regulates the ligand-binding specificity of apophytochrome. On the basis of the hydrogen-bonding ability of Tyr176 substitutions that support the nonphotochemical C15-Z,syn to C15-Z,anti interconversion, we propose that Tyr176 orients the carboxyl side chain of a conserved acidic residue to stabilize protonation of the bilin chromophore. A homology model of the GAF domain of Cph1 predicts a C5-Z,syn, C10-Z,syn, C15-Z,anti configuration for the chromophore and implicates Glu189 as the proposed acidic residue stabilizing the extended conformation, an interpretation consistent with site-directed mutagenesis of this conserved acidic residue.     

The N-terminal domain determines the affinity and specificity of H1 binding to chromatin

Linker histone H1, one of the most abundant nuclear proteins in multicellular eukaryotes, is a key component of the chromatin structure mainly due to its role in the formation and maintenance of the 30nm chromatin fiber. It has a three-domain structure; a central globular domain flanked by a short N-terminal domain and a long, highly basic C-terminal domain. Previous studies have shown that the binding abilities of H1 are at large determined by the properties of the C-terminal domain; much less attention has been paid to role of the N-terminal domain. We have previously shown that H1 can be reconstituted via cytoplasmic mRNA injection in Xenopus oocytes, cells that lack somatic H1. The heterologously expressed H1 proteins are incorporated into in vivo assembled chromatin at specific sites and the binding event is monitored as an increase in nucleosomal repeat length (NRL). Using this setup we have here compared the binding properties of wt-H1.4 and hH1.4 devoid of its N-terminal domain (ΔN-hH1.4). The ΔN-hH1.4 displays a drastically lower affinity for chromatin binding as compared to the wild type hH1.4. Our data also indicates that ΔN-hH1.4 is more prone to unspecific chromatin binding than the wild type. We conclude that the N-terminal domain of H1 is an important determinant of affinity and specificity of H1-chromatin interactions.     

A Conserved residue in the yeast Bem1p SH3 domain maintains the high level of binding specificity required for function

The yeast Bem1p SH3b and Nbp2p SH3 domains are unusual because they bind to peptides containing the same consensus sequence, yet they perform different functions and display low sequence similarity. In this work, by analyzing the interactions of these domains with six biologically relevant peptides containing the consensus sequence, they are shown to possess finely tuned and distinct binding specificities. We also identify a residue in the Bem1p SH3b domain that inhibits binding, yet is highly conserved for the purpose of preventing nonspecific interactions. Substitution of this residue results in a marked reduction of in vivo function that is caused by titration of the domain away from its proper targets through nonspecific interactions with other proteins. This work provides a clear illustration of the importance of intrinsic binding specificity for the function of protein-protein interaction modules, and the key role of "negative" interactions in determining the specificity of a domain.     

The transmembrane domain provides nucleotide binding specificity to the bacterial conjugation protein TrwB

In order to understand the functional significance of the transmembrane domain of TrwB, an integral membrane protein involved in bacterial conjugation, the protein was purified in the native, and also as a truncated soluble form (TrwBΔN70). The intact protein (TrwB) binds preferentially purine over pyrimidine nucleotides, NTPs over NDPs, and ribo- over deoxyribonucleotides. In contrast, TrwBΔN70 binds uniformly all tested nucleotides. The transmembrane domain has the general effect of making the nucleotide binding site(s) less accessible, but more selective. This is in contrast to other membrane proteins in which most of the protein mass, including the catalytic domain, is outside the membrane, but whose activity is not modified by the presence or absence of the transmembrane segment.     

The intermolecular interaction between the PH domain and the C-terminal domain of Arabidopsis dynamin-like 6 determines lipid binding specificity

Dynamin and its related proteins are a group of mechanochemical proteins involved in the modulation of lipid membranes in various biological processes. Here we investigate the nature of membrane binding of the Arabidopsis dynamin-like 6 (ADL6) involved in vesicle trafficking from the trans-Golgi network to the central vacuole. Fractionation experiments by continuous sucrose gradients and gel filtration revealed that the majority of ADL6 is associated with membranes in vivo. Amino acid sequence analysis revealed that ADL6 has a putative pleckstrin homology (PH) domain. In vitro lipid binding assays demonstrated that ADL6 showed high affinity binding to phosphatidylinositol 3-phosphate (PtdIns-3-P) and that the PH domain was responsible for this interaction. However, the PH domain alone binds equally well to both PtdIns-3-P and phosphatidylinositol 4-phosphate (PtdIns-4-P). Interestingly, the high affinity binding of the PH domain to PtdIns-3-P was restored by a protein-protein interaction between the PH domain and the C-terminal region. In addition, deletion of the inserted regions within the PH domain results in high affinity binding of the PH domain to PtdIns-3-P. These results suggest that ADL6 binds specifically to PtdIns-3-P and that the lipid binding specificity is determined by the interaction between the PH domain and the C-terminal domain of ADL6.     

The salt-dependence of a protein-ligand interaction: ion-protein binding energetics

Using the binding of a nucleotide inhibitor (guanosine-3'-monophosphate) to a ribonuclease (ribonuclease Sa) as a model system, we show that the salt-dependence of the interaction arises due to specific ion binding at the site of nucleotide binding. The presence of specific ion-protein binding is concluded from a combination of differential scanning calorimetry and NMR data. Isothermal titration calorimetry data are then fit to determine the energetic profile (enthalpy, entropy, and heat capacity) for both the ion-protein and nucleotide-protein interactions. The results provide insight into the energetics of charge-charge interactions, and have implications for the interpretation of an observed salt-dependence. Further, the presence of specific ion-binding leads to a system behavior as a function of temperature that is drastically different from that predicted from Poisson-Boltzmann calculations.     

A new measure of cooperativity in protein-ligand binding

An allosteric binding system consisting of a single ligand and a nondissociating macromolecule having multiple binding sites can be represented by a binding polynomial. Various properties of the binding process can be obtained by analyzing the coefficients of the binding polynomial and such functions as the binding curve and the Hill plot. The Hill plot has an asymptote of unit slope at each end and the departure of the slope from unity at any point can be used to measure the effective interaction free energy at that point. Of particular interest in detecting and measuring cooperativity are extrema of the Hill slope and its value at the half-saturation point. If the binding polynomial is symmetric, then there is an extremum of the Hill slope at the half-saturation point. This value, the Hill coefficient, is a convenient measure of cooperativity. The purpose of this paper is to express the Hill coefficient for symmetric binding polynomials in terms of the roots of the polynomial and to give an interpretation of cooperativity in terms of the geometric pattern of the roots in the complex plane. This interpretation is then applied to the binding polynomials for the MWC (Monod-Wyman-Changeux) and KNF (Koshland-Nemethy-Filmer) models.     

The ATPase domain of hsp70 possesses a unique binding specificity for 3'-sulfogalactolipids

The region(s) of hsp70 critical for sulfogalactolipid (SGL) recognition has been defined through deletion analysis and site-directed mutagenesis. Truncated polymerase chain reaction products of hsp70 generated N-terminal fragments of 43, 35, 29, and 22 kDa. The C terminus substrate-binding domain (28 kDa) was also expressed. The N-terminal ATPase domain (rP43) shared the binding specificity of hsp70, because only sulfogalactosyl ceramide and sulfogalactosyl glycerolipid were recognized by both TLC overlay and RELISA. The C-terminal domain showed no binding. SGL binding of rP29 and rP22 was severely reduced. The loss of SGL binding for rP35 by RELISA but not TLC overlay was considered as a function of receptor presentation. The truncation of rP43 to rP35 demonstrates that residues 318-387 (the base of the ATP binding cleft) are critical for high affinity SGL binding. Mutagenesis showed that Arg(342) and Phe(198) are crucial for this process. SGL binding, mediated by these conserved residues within the ATPase domain of hsp70, implies that this binding specificity is evolutionarily conserved.     

The roles of conserved amino acids on substrate binding and conformational integrity of ClpB N-terminal domain

Escherichia coli heat shock protein ClpB disaggregates denatured protein in cooperation with the DnaK chaperone system. Several studies showed that the N-terminal domain is essential for the chaperone activity, but its role is still largely unknown. The N-terminal domain contains two structurally similar subdomains, and conserved amino acids Thr7 and Ser84 share the same position in two apparent sequence repeats. T7A and S84A substitutions affected chaperone activity of ClpB without significantly changing the native conformation [Liu, Z. et al. (2002) J. Mol. Biol. 321, 111-120]. In this study, we aimed to better understand the roles of several conserved amino acid residues, including Thr7 and Ser84, in the N-terminal domain. We investigated the effects of mutagenesis on substrate binding and conformational states of ClpB N-terminal domain fragment (ClpBN). Fluorescence polarization analysis showed that the T7A and S84A substitutions enhanced the interaction between ClpBN and protein aggregates. Interestingly, further analyses suggested that the mechanisms by which they do so are quite different. For T7A substitution, the increased substrate affinity could be due to a conformational change in the hydrophobic core as revealed by NMR spectroscopy. In contrast, for S84A, increased substrate binding would be explained by a unique conformational state of this mutant as revealed by pressure perturbation analysis. The thermal transition temperature of the S84A mutant, monitored by DSC, was 6.1 degrees C lower than that of wild-type. Our results revealed that conserved amino acids Thr7 and Ser84 both participated in maintaining the conformational integrity of the ClpB N-terminal domain.     

A single residue of GDF-5 defines binding specificity to BMP receptor IB

Growth and differentiation factor 5 (GDF-5), a member of the TGF-beta superfamily, is involved in many developmental processes, like chondrogenesis and joint formation. Mutations in GDF-5 lead to diseases, e.g. chondrodysplasias like Hunter-Thompson, Grebe and DuPan syndromes and brachydactyly. Similar to other TGF-beta superfamily members, GDF-5 transmits signals through binding to two different types of membrane-bound serine-/threonine-kinase receptors termed type I and type II. In contrast to the large number of ligands, only seven type I and five type II receptors have been identified to date, implicating a limited promiscuity in ligand-receptor interaction. However, in contrast to other members of the TGF-beta superfamily, GDF-5 shows a pronounced specificity in type I receptor interaction in cross-link experiments binding only to BMP receptor IB (BMPR-IB). In mice, deletion of either GDF-5 or BMPR-IB results in a similar phenotype, indicating that GDF-5 signaling is highly dependent on BMPR-IB. Here, we demonstrate by biosensor analysis that GDF-5 also binds to BMP receptor IA (BMPR-IA) but with approximately 12-fold lower affinity. Structural and mutational analyses revealed a single residue of GDF-5, Arg57 located in the pre-helix loop, being solely responsible for the high binding specificity to BMPR-IB. In contrast to wild-type GDF-5, variant GDF-5R57A interacts with BMPR-IA and BMPR-IB with a comparable high binding affinity. These results provide important insights into how receptor-binding specificity is generated at the molecular level and might be useful for the generation of receptor subtype specific activators or inhibitors.