The functional binding epitope of a high affinity variant of human growth hormone mapped by shotgun alanine-scanning mutagenesis: insights into the mechanisms responsible for improved affinity

A high-affinity variant of human growth hormone (hGH(v)) contains 15 mutations within site 1 and binds to the hGH receptor (hGHR) approximately 400-fold tighter than does wild-type (wt) hGH (hGH(wt)). We used shotgun scanning combinatorial mutagenesis to dissect the energetic contributions of individual residues within the hGH(v) binding epitope and placed them in context with previously determined structural information. In all, the effects of alanine substitutions were determined for 35 hGH(v) residues that are directly contained in or closely border the binding interface. We found that the distribution of binding energy in the functional epitope of hGH(v) differs significantly from that of hGH(wt). The residues that contributed the majority of the binding energy in the wt interaction (the so-called binding "hot spot") remain important, but their contributions are attenuated in the hGH(v) interaction, and additional binding energy is acquired from residues on the periphery of the original hotspot. Many interactions that inhibited the binding of hGH(wt) are replaced by interactions that make positive contributions to the binding of hGH(v). These changes produce an expanded and diffused hot spot in which improved affinity results from numerous small contributions distributed broadly over the interface. The mutagenesis results are consistent with previous structural studies, which revealed widespread structural differences between the wt and variant hormone-receptor interfaces. Thus, it appears that the improved binding affinity of hGH(v) site 1 was not achieved through minor adjustments to the wt interface, but rather, results from a wholesale reconfiguration of many of the original binding elements.     

Dominant Saccharomyces cerevisiae msh6 mutations cause increased mispair binding and decreased dissociation from mispairs by Msh2-Msh6 in the presence of ATP

A previous study described four dominant msh6 mutations that interfere with both the Msh2-Msh6 and Msh2-Msh3 mismatch recognition complexes (Das Gupta, R., and Kolodner, R. D. (2000) Nat. Genet. 24, 53-56). Modeling predicted that two of the amino acid substitutions (G1067D and G1142D) interfere with protein-protein interactions at the ATP-binding site-associated dimer interface, one (S1036P) similarly interferes with protein-protein interactions and affects the Msh2 ATP-binding site, and one (H1096A) affects the Msh6 ATP-binding site. The ATPase activity of the Msh2-Msh6-G1067D and Msh2-Msh6-G1142D complexes was inhibited by GT, +A, and +AT mispairs, and these complexes showed increased binding to GT and +A mispairs in the presence of ATP. The ATPase activity of the Msh2-Msh6-S1036P complex was inhibited by a GT mispair, and it bound the GT mispair in the presence of ATP, whereas its interaction with insertion mispairs was unchanged compared with the wild-type complex. The ATPase activity of the Msh2-Msh6-H1096A complex was generally attenuated, and its mispair-binding behavior was unaffected. These results are in contrast to those obtained with the wild-type Msh2-Msh6 complex, which showed mispair-stimulated ATPase activity and ATP inhibition of mispair binding. These results indicate that the dominant msh6 mutations cause more stable binding to mispairs and suggest that there may be differences in how base base and insertion mispairs are recognized.     

The role of protein dynamics in increasing binding affinity for an engineered protein-protein interaction established by H/D exchange mass spectrometry

It is generally accepted that protein and solvation dynamics play fundamental roles in the mechanisms of protein-protein binding; however, assessing their contribution meaningfully has not been straightforward. Here, hydrogen/deuterium exchange mass spectrometry (H/D-Ex) was employed to assess the role of dynamics for a high-affinity human growth hormone variant (hGHv) and the wild-type growth hormone (wt-hGH) each binding to the extracellular domain of their receptor (hGHbp). Comparative analysis of the transient fluctuations in the bound and unbound states revealed that helix-1 of hGHv undergoes significant transient unfolding in its unbound state, a characteristic that was not found in wt-hGH or apparent in the temperature factor data from the X-ray analysis of the unbound hGHv structure. In addition, upon hormone binding, an overall increase in stability was observed for the beta-sheet structure of hGHbp which included sites distant from the binding interface. On the basis of the stability, binding kinetics, and thermodynamic data presented, the increase in the binding free energy of hGHv is primarily generated by factors that appear to increase the energy of the unbound state relative to the free energy of the bound complex. This implies that an alternate route to engineer new interactions aiming to increase protein-protein association energies may be achieved by introducing certain mutations that destabilize one of the interacting molecules without destabilizing the resulting bound complex. Importantly, although the hGHv molecule is less stable than its wt-hGH counterpart, its resulting active ternary complex with two copies of hGHbp has comparable stability to the wt complex.     

Crystal structure and site 1 binding energetics of human placental lactogen

In primates, placental lactogen (PL) is a pituitary hormone with fundamental roles during pregnancy involving fetal growth, metabolism, and stimulating lactation in the mother. Human placental lactogen (hPL) is highly conserved with human growth hormone (hGH) and both hormones bind to the hPRLR extracellular domain (ECD), the first step in receptor homodimerization, in a Zn2+-dependent manner. A modified surface plasmon resonance method was developed to measure the kinetics for hPL and hGH binding to the hPRLR ECD, with and without Zn2+ and showed that hPL has about a tenfold higher affinity for the hPRLR ECD1 than hGH. The crystal structure of the free state of hPL has been determined to 2.0 A resolution showing the molecule possesses an overall structure similar to other long chain four-helix bundle cytokines. Comparison of the free hPL structure with the 1:1 complex structure of hGH bound to the hPRLR ECD1 suggests that two surface loops undergo conformational changes >10 A upon binding. An 18 residue Ala-scan was used to characterize the binding energy epitope for the site 1 interface of hPL. Individual alanine substitutions at five positions reduced binding affinity by a DeltaDeltaG > or = 3 kcal mol(-1). A comparison of the hPL site 1 epitope with that previously determined for hGH indicates contributions of individual residues track reasonably well between hPL and hGH. In particular, residues involved in the zinc-binding site and Lys172 constitute the principal binding determinants for both hormones. However, several residues that are identical between hPL and hGH contribute quite differently to the binding of the hPRLR ECD1. Additionally, the overall magnitudes of the DeltaDeltaG changes observed from the Ala-scan of hPL were markedly larger than those determined in the comparative scan of hGH to the hPRLR ECD1. The structural and biophysical data presented here show that subtle changes in the structural context of an interaction can lead to significantly different effects at the individual residue level.     

Intramolecular cooperativity in a protein binding site assessed by combinatorial shotgun scanning mutagenesis

Combinatorial shotgun alanine-scanning was used to assess intramolecular cooperativity in the high affinity site (site 1) of human growth hormone (hGH) for binding to its receptor. A total of 19 side-chains were analyzed and statistically significant data were obtained for 145 of the 171 side-chain pairs. The analysis revealed that 90% of the side-chain pairs exhibited no statistically significant pair interactions, and the remaining 10% of side-chain pairs exhibited only small interactions corresponding to cooperative interaction energies with magnitudes less than 0.4 kcal/mol. The statistical predictions were tested by measuring affinities for purified mutant proteins and were found to be accurate for five of six side-chain pairs tested. The results reveal that hGH site 1 behaves in a highly additive manner and suggest that shotgun scanning should be useful for assessing cooperative effects in other protein-protein interactions.     

Hot-spot residues at the E9/Im9 interface help binding via different mechanisms

Protein-protein association involves many interface interactions, but they do not contribute equally. Ala scanning experiments reveal that only a few mutations significantly lower binding affinity. These key residues, which appear to drive protein-protein association, are called hot-spot residues. Molecular dynamics simulations of the Colicin E9/Im9 complex show Im9 Glu41 and Im9 Ser50, both hot-spots, bind via different mechanisms. The results suggest that Im9 Ser50 restricts Glu41 in a conformation auspicious for salt-bridge formation across the interface. This type of model may be helpful in engineering hot-spot clusters at protein-protein interfaces and, consequently, the design of specificity. (c) 2008 Wiley Periodicals, Inc. Biopolymers 89: 916-920, 2008.This article was originally published online as an accepted preprint. The "Published Online" date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com.     

Structure of a phage display-derived variant of human growth hormone complexed to two copies of the extracellular domain of its receptor: evidence for strong structural coupling between receptor binding sites

The structure of the ternary complex between the phage display- optimized, high-affinity Site 1 variant of human growth hormone (hGH) and two copies of the extracellular domain (ECD) of the hGH receptor (hGHR) has been determined at 2.6 A resolution. There are widespread and significant structural differences compared to the wild-type ternary hGH hGHR complex. The hGH variant (hGH(v)) contains 15 Site 1 mutations and binds>10(2) tighter to the hGHR ECD (hGH(R1)) at Site 1. It is biologically active and specific to hGHR. The hGH(v) Site 1 interface is somewhat smaller and 20% more hydrophobic compared to the wild-type (wt) counterpart. Of the ten hormone-receptor H-bonds in the site, only one is the same as in the wt complex. Additionally, several regions of hGH(v) structure move up to 9A in forming the interface. The contacts between the C-terminal domains of two receptor ECDs (hGH(R1)- hGH(R2)) are conserved; however, the large changes in Site 1 appear to cause global changes in the domains of hGH(R1) that affect the hGH(v)-hGH(R2) interface indirectly. This coupling is manifested by large changes in the conformation of groups participating in the Site 2 interaction and results in a structure for the site that is reorganized extensively. The hGH(v)- hGH(R2) interface contains seven H-bonds, only one of which is found in the wt complex. Several groups on hGH(v) and hGH(R2) undergo conformational changes of up to 8 A. Asp116 of hGH(v) plays a central role in the reorganization of Site 2 by forming two new H-bonds to the side-chains of Trp104(R2) and Trp169(R2), which are the key binding determinants of the receptor. The fact that a different binding solution is possible for Site 2, where there were no mutations or binding selection pressures, indicates that the structural elements found in these molecules possess an inherent functional plasticity that enables them to bind to a wide variety of binding surfaces.     

Shotgun alanine scanning shows that growth hormone can bind productively to its receptor through a drastically minimized interface

The high affinity binding site (Site1) of the human growth hormone (hGH) binds to its cognate receptor (hGHR) via a concave surface patch containing about 35 residues. Using 167 sequences from a shotgun alanine scanning analysis of Site1, we have determined that over half of these residues can be simultaneously changed to an alanine or a non-isosteric amino acid while still retaining a high affinity interaction. Among these hGH variants the distribution of the mutation is highly variable throughout the interface, although helix 4 is more conserved than the other binding elements. Kinetic and thermodynamic analyses were performed on 11 representative hGH Site1 variants that contained 14-20 mutations. Generally, the tightest binding variants showed similar associated rate constants (k(on)) as the wild-type (wt) hormone, indicating that their binding proceeds through a similar transition state intermediate. However, calorimetric analyses indicate very different thermodynamic partitioning: wt-hGH binding exhibits favorable enthalpy and entropy contributions, whereas the variants display highly favorable enthalpy and highly unfavorable entropy contributions. The heat capacities (DeltaCp) on binding measured for wt-hGH and its variants are significantly larger than normally seen for typical protein-protein interactions, suggesting large conformational or solvation effects. The multiple Site1 mutations are shown to indirectly affect binding of the second receptor at Site2 through an allosteric mechanism. We show that the stability of the ternary hormone-receptor complex reflects the affinity of the Site2 binding and is surprisingly exempt from changes in Site1 affinity, directly demonstrating that dissociation of the active signaling complex is a stepwise process.     

Optimization of peptidic HIV‐1 fusion inhibitor T20 by phage display

The HIV fusion inhibitor T20 has been approved to treat those living with HIV/AIDS, but treatment gives rise to resistant viruses. Using combinatorial phage‐displayed libraries, we applied a saturation scan approach to dissect the entire T20 sequence for binding to a prefusogenic five‐helix bundle (5HB) mimetic of HIV‐1 gp41. Our data set compares all possible amino acid substitutions at all positions, and affords a complete view of the complex molecular interactions governing the binding of T20 to 5HB. The scan of T20 revealed that 12 of its 36 positions were conserved for 5HB binding, which cluster into three epitopes: hydrophobic epitopes at the ends and a central dyad of hydrophilic residues. The scan also revealed that the T20 sequence was highly adaptable to mutations at most positions, demonstrating a striking structural plasticity that allows multiple amino acid substitutions at contact points to adapt to conformational changes, and also at noncontact points to fine‐tune the interface. Based on the scan result and structural knowledge of the gp41 fusion intermediate, a library was designed with tailored diversity at particular positions of T20 and was used to derive a variant (T20v1) that was found to be a highly effective inhibitor of infection by multiple HIV‐1 variants, including a common T20‐escape mutant. These findings show that the plasticity of the T20 functional sequence space can be exploited to develop variants that overcome resistance of HIV‐1 variants to T20 itself, and demonstrate the utility of saturation scanning for rapid epitope mapping and protein engineering.     

Modeling Effects of Human Single Nucleotide Polymorphisms on Protein-Protein Interactions

A large set of three-dimensional structures of 264 protein-protein complexes with known nonsynonymous single nucleotide polymorphisms (nsSNPs) at the interface was built using homology-based methods. The nsSNPs were mapped on the proteins' structures and their effect on the binding energy was investigated with CHARMM force field and continuum electrostatic calculations. Two sets of nsSNPs were studied: disease annotated Online Mendelian Inheritance in Man (OMIM) and nonannotated (non-OMIM). It was demonstrated that OMIM nsSNPs tend to destabilize the electrostatic component of the binding energy, in contrast with the effect of non-OMIM nsSNPs. In addition, it was shown that the change of the binding energy upon amino acid substitutions is not related to the conservation of the net charge, hydrophobicity, or hydrogen bond network at the interface. The results indicate that, generally, the effect of nsSNPs on protein-protein interactions cannot be predicted from amino acids' physico-chemical properties alone, since in many cases a substitution of a particular residue with another amino acid having completely different polarity or hydrophobicity had little effect on the binding energy. Analysis of sequence conservation showed that nsSNP at highly conserved positions resulted in a large variance of the binding energy changes. In contrast, amino acid substitutions corresponding to nsSNPs at nonconserved positions, on average, were not found to have a large effect on binding affinity. pKa calculations were performed and showed that amino acid substitutions could change the wild-type proton uptake/release and thus resulting in different pH-dependence of the binding energy.     

Dual requirement for flexibility and specificity for binding of the coiled-coil tropomyosin to its target, actin

The coiled coil is a widespread motif involved in oligomerization and protein-protein interactions, but the structural requirements for binding to target proteins are poorly understood. To address this question, we measured binding of tropomyosin, the prototype coiled coil, to actin as a model system. Tropomyosin binds to the actin filament and cooperatively regulates its function. Our results support the hypothesis that coiled-coil domains that bind to other proteins are flexible. We made mutations that alter interface packing and stability as well as mutations in surface residues in a postulated actin binding site. Actin affinity, measured by cosedimentation, was correlated with coiled-coil stability and local instability and side chain flexibility, analyzed with circular dichroism and fluorescence spectroscopy. The flexibility from interruptions in the stable coiled-coil interface is essential for actin binding. The surface residues in a postulated actin binding site participate in actin binding when the coiled coil within it is poorly packed.     

Probing the effect of point mutations at protein-protein interfaces with free energy calculations

We have studied the effect of point mutations of the primary binding residue (P1) at the protein-protein interface in complexes of chymotrypsin and elastase with the third domain of the turkey ovomucoid inhibitor and in trypsin with the bovine pancreatic trypsin inhibitor, using molecular dynamics simulations combined with the linear interaction energy (LIE) approach. A total of 56 mutants have been constructed and docked into their host proteins. The free energy of binding could be reliably calculated for 52 of these mutants that could unambiguously be fitted into the binding sites. We find that the predicted binding free energies are in very good agreement with experimental data with mean unsigned errors between 0.50 and 1.03 kcal/mol. It is also evident that the standard LIE model used to study small drug-like ligand binding to proteins is not suitable for protein-protein interactions. Three different LIE models were therefore tested for each of the series of protein-protein complexes included, and the best models for each system turn out to be very similar. The difference in parameterization between small drug-like compounds and protein point mutations is attributed to the preorganization of the binding surface. Our results clearly demonstrate the potential of free energy calculations for probing the effect of point mutations at protein-protein interfaces and for exploring the principles of specificity of hot spots at the interface.     

The six conserved helicase motifs of the UL5 gene product, a component of the herpes simplex virus type 1 helicase-primase, are essential for its function.

The UL5 protein of herpes simplex virus type 1, one component of the viral helicase-primase complex, contains six sequence motifs found in all members of a superfamily of DNA and RNA helicases. Although this superfamily contains more than 20 members ranging from bacteria to mammalian cells and their viruses, the importance of these motifs has not been addressed experimentally for any one of them. In this study, we have examined the functional significance of these six motifs for the UL5 protein through the introduction of site-specific mutations resulting in single amino acid substitutions of the most highly conserved residues within each motif. A transient replication complementation assay was used to test the effect of each mutation on the function of the UL5 protein in viral DNA replication. In this assay, a mutant UL5 protein expressed from an expression clone is used to complement a replication-deficient null mutant with a mutation in the UL5 gene for the amplification of herpes simplex virus origin-containing plasmids. Eight mutations in conserved regions and three similar mutations in nonconserved regions of the UL5 gene were analyzed, and the results indicate that all six conserved motifs are essential to the function of UL5 protein in viral DNA replication; on the other hand, mutations in nonconserved regions are tolerated. These data provide the first direct evidence for the importance of these conserved regions in any member of the superfamily of DNA and RNA helicases. In addition, three motif mutations were introduced into the viral genome, and the phenotypic analyses of these mutants are consistent with results from the transient replication complementation assay. The ability of these three mutant UL5 proteins to form specific interactions with other members of the helicase-primase complex, UL8 and UL52, indicates that the functional domains required for replication activity of UL5 are separable from domains responsible for protein-protein interactions. It is anticipated that this type of structure-function analysis will lead to the identification of protein domains that contribute not only to the enzymatic activities of helicase or primase but also to protein-protein interactions within members of the complex.     

Cluster conservation as a novel tool for studying protein-protein interactions evolution

Protein-protein interactions networks has come to be a buzzword associated with nets containing edges that represent a pair of interacting proteins (e.g. hormone-receptor, enzyme-inhibitor, antigen-antibody, and a subset of multichain biological machines). Yet, each such interaction composes its own unique network, in which vertices represent amino acid residues, and edges represent atomic contacts. Recent studies have shown that analyses of the data encapsulated in these detailed networks may impact predictions of structure-function correlation. Here, we study homologous families of protein-protein interfaces, which share the same fold but vary in sequence. In this context, we address what properties of the network are shared among relatives with different sequences (and hence different atomic interactions) and which are not. Herein, we develop the general mathematical framework needed to compare the modularity of homologous networks. We then apply this analysis to the structural data of a few interface families, including hemoglobin alpha-beta, growth hormone-receptor, and Serine protease-inhibitor. Our results suggest that interface modularity is an evolutionarily conserved property. Hence, protein-protein interfaces can be clustered down to a few modules, with the boundaries being evolutionarily conserved along homologous complexes. This suggests that protein engineering of protein-protein binding sites may be simplified by varying each module, but retaining the overall modularity of the interface.     

Evolutionary potential of a duplicated repressor-operator pair: simulating pathways using mutation data

Ample evidence has accumulated for the evolutionary importance of duplication events. However, little is known about the ensuing step-by-step divergence process and the selective conditions that allow it to progress. Here we present a computational study on the divergence of two repressors after duplication. A central feature of our approach is that intermediate phenotypes can be quantified through the use of in vivo measured repression strengths of Escherichia coli lac mutants. Evolutionary pathways are constructed by multiple rounds of single base pair substitutions and selection for tight and independent binding. Our analysis indicates that when a duplicated repressor co-diverges together with its binding site, the fitness landscape allows funneling to a new regulatory interaction with early increases in fitness. We find that neutral mutations do not play an essential role, which is important for substantial divergence probabilities. By varying the selective pressure we can pinpoint the necessary ingredients for the observed divergence. Our findings underscore the importance of coevolutionary mechanisms in regulatory networks, and should be relevant for the evolution of protein-DNA as well as protein-protein interactions.     

Non-Redundant Unique Interface Structures as Templates for Modeling Protein Interactions

Improvements in experimental techniques increasingly provide structural data relating to protein-protein interactions. Classification of structural details of protein-protein interactions can provide valuable insights for modeling and abstracting design principles. Here, we aim to cluster protein-protein interactions by their interface structures, and to exploit these clusters to obtain and study shared and distinct protein binding sites. We find that there are 22604 unique interface structures in the PDB. These unique interfaces, which provide a rich resource of structural data of protein-protein interactions, can be used for template-based docking. We test the specificity of these non-redundant unique interface structures by finding protein pairs which have multiple binding sites. We suggest that residues with more than 40% relative accessible surface area should be considered as surface residues in template-based docking studies. This comprehensive study of protein interface structures can serve as a resource for the community. The dataset can be accessed at http://prism.ccbb.ku.edu.tr/piface.     

Computational design of a new hydrogen bond network and at least a 300-fold specificity switch at a protein-protein interface

The redesign of protein-protein interactions is a stringent test of our understanding of molecular recognition and specificity. Previously we engineered a modest specificity switch into the colicin E7 DNase-Im7 immunity protein complex by identifying mutations that are disruptive in the native complex, but can be compensated by mutations on the interacting partner. Here we extend the approach by systematically sampling alternate rigid body orientations to optimize the interactions in a binding mode specific manner. Using this protocol we designed a de novo hydrogen bond network at the DNase-immunity protein interface and confirmed the design with X-ray crystallographic analysis. Subsequent design of the second shell of interactions guided by insights from the crystal structure on tightly bound water molecules, conformational strain, and packing defects yielded new binding partners that exhibited specificities of at least 300-fold between the cognate and the non-cognate complexes. This multi-step approach should be applicable to the design of polar protein-protein interactions and contribute to the re-engineering of regulatory networks mediated by protein-protein interactions.     

A novel method for protein-protein interaction site prediction using phylogenetic substitution models

Protein-protein binding events mediate many critical biological functions in the cell. Typically, functionally important sites in proteins can be well identified by considering sequence conservation. However, protein-protein interaction sites exhibit higher sequence variation than other functional regions, such as catalytic sites of enzymes. Consequently, the mutational behavior leading to weak sequence conservation poses significant challenges to the protein-protein interaction site prediction. Here, we present a phylogenetic framework to capture critical sequence variations that favor the selection of residues essential for protein-protein binding. Through the comprehensive analysis of diverse protein families, we show that protein binding interfaces exhibit distinct amino acid substitution as compared with other surface residues. On the basis of this analysis, we have developed a novel method, BindML, which utilizes the substitution models to predict protein-protein binding sites of protein with unknown interacting partners. BindML estimates the likelihood that a phylogenetic tree of a local surface region in a query protein structure follows the substitution patterns of protein binding interface and nonbinding surfaces. BindML is shown to perform well compared to alternative methods for protein binding interface prediction. The methodology developed in this study is very versatile in the sense that it can be generally applied for predicting other types of functional sites, such as DNA, RNA, and membrane binding sites in proteins.