Abstract: Multiple type I interferons (IFN-alpha/beta) elicit Jak/Stat activation, rapid gene induction, and pleiotropic effects, such as differentiation, antiviral protection, and blocks in proliferation, which are dependent on the IFN subtype and the cellular context. To date, ligand- and receptor-specific molecular determinants underlying IFN-alpha/beta differential activities or potencies have been well characterized. To analyze cellular determinants that impact subtype-specific potency, human fibrosarcoma U5A-derived clones, exhibiting a gradient of IFN sensitivity by virtue of increasing receptor levels, were monitored for Jak/Stat signaling, gene induction, cell cycle lengthening, and apoptosis. In cells with scarce receptors, IFN-beta was more potent than IFN-alpha2 in antiproliferative activity, while the two subtypes were equipotent in all other readouts. Conversely, in cells with abundant receptors, IFN-alpha2 matched or even surpassed IFN-beta in all readouts tested. Our results suggest that the differential activities of the IFN subtypes are dictated not only by the intrinsic ligand/receptor binding kinetics but also by the density of cell surface receptor components.

Abstract: The three-dimensional structures of proteins are stabilized by the interactions between amino acid residues. Here we report a method where four distances are calculated between any two side chains to provide an exact spatial definition of their bonds. The data were binned into a four-dimensional grid and compared to a random model, from which the preference for specific four-distances was calculated. A clear relation between the quality of the experimental data and the tightness of the distance distribution was observed, with crystal structure data providing far tighter distance distributions than NMR data. Since the four-distance data have higher information content than classical bond descriptions, we were able to identify many unique inter-residue features not found previously in proteins. For example, we found that the side chains of Arg, Glu, Val and Leu are not symmetrical in respect to the interactions of their head groups. The described method may be developed into a function, which computationally models accurately protein structures.

Abstract:

Abstract: The association reaction between pairs of proteins proceeds through an encounter complex that develops into the final complex. Here, we combined Brownian dynamics simulations with experimental studies to analyze the structures of the encounter complexes along the association reaction between TEM1-beta-lactamase and its inhibitor, beta-lactamase-inhibitor protein. The encounter complex can be considered as an ensemble of short-lived low free-energy states that are stabilized primarily by electrostatic forces and desolvation. For the wild-type, the simulation showed two main encounter regions located outside the physical binding site. One of these regions was located near the experimentally determined transition state. To validate whether these encounters are fruitful or futile, we examined three groups of mutations that altered the encounter. The first group consisted of mutations that increased the experimental rate of association through electrostatic optimization. This resulted in an increase in the size of the encounter region located near the experimentally determined transition state, as well as a decrease in the energy of this region and an increase in the number of successful trajectories (i.e., encounters that develop into complex). A second group of mutations was specifically designed to either increase or decrease the size and energy of the second encounter complex, but either way it did not affect k(on). A third group of mutations consisted of residues that increased k(on) without significantly affecting the encounter complexes. These results indicate that the size and energy of the encounter regions are only two of several parameters that lead to fruitful association, and that electrostatic optimization is a major driving force in fast association.

Abstract: Methods for protein modeling and design advanced rapidly in recent years. At the heart of these computational methods is an energy function that calculates the free energy of the system. Many of these functions were also developed to estimate the consequence of mutation on protein stability or binding affinity. In the current study, we chose six different methods that were previously reported as being able to predict the change in protein stability (DeltaDeltaG) upon mutation: CC/PBSA, EGAD, FoldX, I-Mutant2.0, Rosetta and Hunter. We evaluated their performance on a large set of 2156 single mutations, avoiding for each program the mutations used for training. The correlation coefficients between experimental and predicted DeltaDeltaG values were in the range of 0.59 for the best and 0.26 for the worst performing method. All the tested computational methods showed a correct trend in their predictions, but failed in providing the precise values. This is not due to lack in precision of the experimental data, which showed a correlation coefficient of 0.86 between different measurements. Combining the methods did not significantly improve prediction accuracy compared to a single method. These results suggest that there is still room for improvement, which is crucial if we want forcefields to perform better in their various tasks.

Abstract: Antizyme (Az) is a highly conserved key regulatory protein bearing a major role in regulating polyamine levels in the cell. It has the ability to bind and inhibit ornithine decarboxylase (ODC), targeting it for degradation. Az inhibitor (AzI) impairs the activity of Az. In this study, we mapped the binding sites of ODC and AzI on Az using Ala scan mutagenesis and generated models of the two complexes by constrained computational docking. In order to scan a large number of mutants in a short time, we developed a workflow combining high-throughput mutagenesis, small-scale parallel partial purification of His-tagged proteins and their immobilization on a tris-nitrilotriacetic-acid-coated surface plasmon resonance chip. This combination of techniques resulted in a significant reduction in time for production and measurement of large numbers of mutant proteins. The data-driven docking results suggest that both proteins occupy the same binding site on Az, with Az binding within a large groove in AzI and ODC. However, single-mutant data provide information concerning the location of the binding sites only, not on their relative orientations. Therefore, we generated a large number of double-mutant cycles between residues on Az and ODC and used the resulting interaction energies to restrict docking. The model of the complex is well defined and accounts for the mutant data generated here, and previously determined biochemical data for this system. Insights on the structure and function of the complexes, as well as general aspects of the method, are discussed.

Abstract: Studies of protein-protein interactions, carried out in polymer solutions, are designed to mimic the crowded environment inside living cells. It was shown that crowding enhances oligomerization and polymerization of macromolecules. Conversely, we have shown that crowding has only a small effect on the rate of association of protein complexes. Here, we investigated the equilibrium effects of crowding on protein heterodimerization of TEM1-beta-lactamase with beta-lactamase inhibitor protein (BLIP) and barnase with barstar. We also contrasted these with the effect of crowding on the weak binding pair CyPet-YPet. We measured the association and dissociation rates as well as the affinities and thermodynamic parameters of these interactions in polyethylene glycol and dextran solutions. For TEM1-BLIP and for barnase-barstar, only a minor reduction in association rate constants compared to that expected based on solution viscosity was found. Dissociation rate constants showed similar levels of reduction. Overall, this resulted in a binding affinity that is quite similar to that in aqueous solutions. On the other hand, for the CyPet-YPet pair, aggregation, and not enhanced dimerization, was detected in polyethylene glycol solutions. The results suggest that typical crowding agents have only a small effect on specific protein-protein dimerization reactions. Although crowding in the cell results from proteins and other macromolecules, one may still speculate that binding in vivo is not very different from that measured in dilute solutions.

Abstract:

Abstract: Type I interferons (IFNs) are multifunctional cytokines that activate cellular responses by binding a common receptor consisting of two subunits, IFNAR-1 and IFNAR-2. Although the binding of IFNs to IFNAR-2 is well characterized, the binding to the lower affinity IFNAR-1 remains less well understood. Previous reports identified a region of human IFN-alpha2 on the B and C helices ("site 1A": N65, L80, Y85, Y89) that plays a key role in binding IFNAR-1 and contributes strongly to differential activation by various type I IFNs. The current studies demonstrate that residues on the D helix are also involved in IFNAR-1 binding. In particular, residue 120 (Arg in IFN-alpha2; Lys in IFN-alpha2/alpha1) appears to be a "hot-spot" residue: substitution by alanine significantly decreased biological activity, and the charge-reversal mutation of residue 120 to Glu caused drastic loss of antiviral and antiproliferative activity for both IFN-alpha2 and IFN-alpha2/alpha1. Mutations in residues of helix D maintained their affinity for IFNAR-2 but had decreased affinity for IFNAR-1. Single-site or multiple-site mutants in the IFNAR-1 binding site that had little or no detectable in vitro biological activity were capable of blocking in vitro antiviral and antiproliferative activity of native IFN-alpha2; i.e., they are type I IFN antagonists. These prototype IFN antagonists can be developed further for possible therapeutic use in systemic lupus erythematosus, and analogous molecules can be designed for use in animal models.

Abstract: A new method is presented for the redesign of protein-protein interfaces, resulting in specificity of the designed pair while maintaining high affinity. The design is based on modular interface architecture and was carried out on the interaction between TEM1 beta-lactamase and its inhibitor protein, beta-lactamase inhibitor protein. The interface between these two proteins is composed of several mostly independent modules. We previously showed that it is possible to delete a complete module without affecting the overall structure of the interface. Here, we replace a complete module with structure fragments taken from nonrelated proteins. Nature-optimized fragments were chosen from 10(7) starting templates found in the Protein Data Bank. A procedure was then developed to identify sets of interacting template residues with a backbone arrangement mimicking the original module. This generated a final list of 361 putative replacement modules that were ranked using a novel scoring function based on grouped atom-atom contact surface areas. The top-ranked designed complex exhibited an affinity of at least the wild-type level and a mode of binding that was remarkably specific despite the absence of negative design in the procedure. In retrospect, the combined application of three factors led to the success of the design approach: utilizing the modular construction of the interface, capitalizing on native rather than artificial templates, and ranking with an accurate atom-atom contact surface scoring function.

Abstract: Protein-water interactions have long been recognized as a major determinant of chain folding, conformational stability, binding specificity and catalysis. However, the detailed effects of water on stabilizing protein-protein interactions remain elusive. A way to test experimentally the contribution of water-mediated interactions is by applying double mutant cycle analysis on pairs of residues that do not form direct interactions, but are bridged by water. Seven such interactions within the interface between TEM1 and BLIP proteins were evaluated. No significant interaction free energy was found between either of them. Water can bridge interactions, but also stabilize the structure of the monomer. To distinguish between these, we performed a bioinformatic analysis using AQUAPROT (http://bioinfo.weizmann.ac.il/aquaprot) to determine the degree of water conservation between the bound and unbound states. 29 structures of twelve complexes and 20 related monomers were analyzed. Of the 262 water molecules located within the interfaces, 145 were conserved between the unbound and bound structures. Strikingly, all 50 buried or partially buried waters in the monomer structures were conserved at the same location in the bound structures. Thus, buried waters have an important role in stabilizing the monomer fold rather than contributing to protein-protein binding, and are not replaced by residues from the incoming protein. Taking together the experimental and bioinformatics evidence suggests that exposed waters within the interface may be good sites for protein engineering, while buried or mostly buried waters should be left unchanged.

Abstract: Type I interferons (IFNs) signal for their diverse biological effects by binding a common receptor on target cells, composed of the two transmembrane IFNAR1 and IFNAR2 proteins. We have previously differentially enhanced the antiproliferative activity of IFN by increasing the weak binding affinity of IFN to IFNAR1. In this study, we further explored the affinity interdependencies between the two receptor subunits and the role of IFNAR1 in differential IFN activity. For this purpose, we generated a panel of mutations targeting the IFNAR2 binding site on the background of the IFNalpha2 YNS mutant, which increases the affinity to IFNAR1 by 60-fold, resulting in IFNAR2-to-IFNAR1 binding affinity ratios ranging from 1000:1 to 1:1000. Both the antiproliferative and antiviral potencies of the interferon mutants clearly correlated to the in situ binding IC(50) values, independently of the relative contributions of the individual receptors, thus relating to the integral lifetime of the complex. However, the antiproliferative potency correlated throughout the entire range of affinities, as well as with prolonged IFNAR1 receptor down-regulation, whereas the antiviral potency reached a maximum at binding affinities equivalent to that of wild-type IFNalpha2. Our data suggest that (i) the specific activity of interferon is related to the ternary complex binding affinity and not to affinity toward individual receptor components and (ii) although the antiviral pathway is strongly dependent on pSTAT1 activity, the cytostatic effect requires additional mechanisms that may involve IFNAR1 down-regulation. This differential interferon response is ultimately mediated through distinct gene expression profiling.

Abstract: 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.

Abstract: Proteins fold into a well-defined structure as a result of the collapse of the polypeptide chain, while transient protein-complex formation mainly is a result of binding of two folded individual monomers. Therefore, a protein-protein interface does not resemble the core of monomeric proteins, but has a more polar nature. Here, we address the question of whether the physico-chemical characteristics of intraprotein versus interprotein bonds differ, or whether interfaces are different from folded monomers only in the preference for certain types of interactions. To address this question we assembled a high resolution, nonredundant, protein-protein interaction database consisting of 1374 homodimer and 572 heterodimer complexes, and compared the physico-chemical properties of these interactions between protein interfaces and monomers. We performed extensive statistical analysis of geometrical properties of interatomic interactions of different types: hydrogen bonds, electrostatic interactions, and aromatic interactions. Our study clearly shows that there is no significant difference in the chemistry, geometry, or packing density of individual interactions between interfaces and monomeric structures. However, the distribution of different bonds differs. For example, side-chain-side-chain interactions constitute over 62% of all interprotein interactions, while they make up only 36% of the bonds stabilizing a protein structure. As on average, properties of backbone interactions are different from those of side chains, a quantitative difference is observed. Our findings clearly show that the same knowledge-based potential can be used for protein-binding sites as for protein structures. However, one has to keep in mind the different architecture of the interfaces and their unique bond preference.

Abstract: Proteins bind one another in aqua's solution to form tight and specific complexes. Previously we have shown that this is achieved through the modular architecture of the interaction network formed by the interface residues, where tight cooperative interactions are found within modules but not between them. Here we extend this study to cover the entire interface of TEM1 beta-lactamase and its protein inhibitor BLIP using an improved method for deriving interaction maps based on REDUCE to add hydrogen atoms and then by evaluating the interactions using modifications of the programs PROBE, NCI and PARE. An extensive mutagenesis study of the interface residues indeed showed that each module is energetically independent on other modules, and that cooperativity is found only within a module. By solving the X-ray structure of two interface mutations affecting two different modules, we demonstrated that protein-protein binding occur via the structural reorganization of the binding modules, either by a "lock and key" or an induced fit mechanism. To explain the cooperativity within a module, we performed multiple-mutant cycle analysis of cluster 2 resulting in a high-resolution energy map of this module. Mutant studies are usually done in reference to alanine, which can be regarded as a deletion of a side-chain. However, from a biological perspective, there is a major interest to understand non-Ala substitutions, as they are most common. Using X-ray crystallography and multiple-mutant cycle analysis we demonstrated the added complexity in understanding non-Ala mutations. Here, a double mutation replacing the wild-type Glu,Tyr to Tyr,Asn on TEM1 (res id 104,105) caused a major backbone structural rearrangement of BLIP, changing the composition of two modules but not of other modules within the interface. This shows the robustness of the modular approach, yet demonstrates the complexity of in silico protein design.

Abstract: All type I IFNs act through a single cell surface receptor composed of the IFNAR1 and IFNAR2 subunits and two associated cytoplasmic tyrosine kinases of the Janus family, Tyk2 and Jak1. A central issue in type I IFN biology is to understand how a multitude of subtypes can generate similar signaling outputs but also govern specific cellular responses. This review summarizes results from the last decade that contributed to our current state of knowledge of IFN-receptor complex structure and assembly.

Abstract: Site-specific conjugation of proteins to surfaces, spectroscopic probes, or other functional units is a key task for implementing biochemical assays. The streptavidin-biotin interaction has proven a highly versatile tool for detection, quantification, and functional analysis of proteins. We have developed an approach for site-specific reversible biotinylation of recombinant proteins through their histidine tag using biotin conjugated to the multivalent chelator trisnitrilotriacetic acid (BTtris-NTA). Stable binding of BTtris-NTA to His-tagged proteins was demonstrated, which is readily reversed by addition of imidazole, enabling versatile conjugation schemes in solution as well as at interfaces. Gel filtration experiments revealed that His-tagged proteins bind to streptavidin doped with BTtris-NTA in a 2:1 stoichiometry. Furthermore, an increased binding affinity toward His-tagged proteins was observed for BTtris-NTA linked to streptavidin compared to tris-NTA in solution and on surfaces. These results indicate an efficient cooperative interaction of two adjacent tris-NTA moieties with a single His-tag, yielding an extremely tight complex with a lifetime of several days. We demonstrate several applications of BTtris-NTA including multiplexed capturing of proteins to biosensor surfaces, cell surface labeling, and Western blot detection. The remarkable selectivity of the His-tag-specific biotinylation, as well as the highly stable, yet reversible complex provides the basis for numerous further applications for functional protein analysis.

Abstract: In a typical cell, proteins function in the crowded cytoplasmic environment where 30% of the space is occupied by macromolecules of varying size and nature. This environment may be simulated in vitro using synthetic polymers. Here, we followed the association and diffusion rates of TEM1-beta-lactamase (TEM) and the beta-lactamase inhibitor protein (BLIP) in the presence of crowding agents of varying molecular mass, from monomers (ethylene glycol, glycerol, or sucrose) to polymeric agents such as different polyethylene glycols (PEGs, 0.2-8 kDa) and Ficoll. An inverse linear relation was found between translational diffusion of the proteins and viscosity in all solutions tested, in accordance with the Stokes-Einstein (SE) relation. Conversely, no simple relation was found between either rotational diffusion rates or association rates (k(on)) and viscosity. To assess the translational diffusion-independent steps along the association pathway, we introduced a new factor, alpha, which corrects the relative change in k(on) by the relative change in solution viscosity, thus measuring the deviations of the association rates from SE behavior. We found that these deviations were related to the three regimes of polymer solutions: dilute, semidilute, and concentrated. In the dilute regime PEGs interfere with TEM-BLIP association by introducing a repulsive force due to solvophobic preferential hydration, which results in slower association than predicted by the SE relation. Crossing over from the dilute to the semidilute regime results in positive deviations from SE behavior, i.e., relatively faster association rates. These can be attributed to the depletion interaction, which results in an effective attraction between the two proteins, winning over the repulsive force. In the concentrated regime, PEGs again dramatically slow down the association between TEM and BLIP, an effect that does not depend on the physical dimensions of PEGs, but rather on their mass concentration. This is probably a manifestation of the monomer-like repulsive depletion effect known to occur in concentrated polymer solutions. As a transition from moderate to high crowding agent concentration can occur in the cellular milieu, this behavior may modulate protein association in vivo, thereby modulating biological function.

Abstract: The development of bioinformatic tools by individual labs results in the abundance of parallel programs for the same task. For example, identification of binding site regions between interacting proteins is done using: ProMate, WHISCY, PPI-Pred, PINUP and others. All servers first identify unique properties of binding sites and then incorporate them into a predictor. Obviously, the resulting prediction would improve if the most suitable parameters from each of those predictors would be incorporated into one server. However, because of the variation in methods and databases, this is currently not feasible. Here, the protein-binding site prediction server is extended into a general protein-binding sites research tool, ProMateus. This web tool, based on ProMate's infrastructure enables the easy exploration and incorporation of new features and databases by the user, providing an evaluation of the benefit of individual features and their combination within a set framework. This transforms the individual research into a community exercise, bringing out the best from all users for optimized predictions. The analysis is demonstrated on a database of protein protein and protein-DNA interactions. This approach is basically different from that used in generating meta-servers. The implications of the open-research approach are discussed. ProMateus is available at http://bip.weizmann.ac.il/promate.

Abstract: Interferons (IFNs) stand in the frontline of defense against viral infections. In this study, we aimed at characterizing the gene expression profile specific to the antiviral effect out of the hundreds of genes involved also in other IFN activities. We found that the IFN-induced antiviral state is maintained for a prolonged time even after IFN occlusion. This was achieved through the active expression of a small set of <40 genes long after IFN was occluded, from which two groups are distinguished: one includes genes participating in direct inhibition of viral replication, such as Mx and OAS; the second group is related to antigen presentation, including all genes involved in the proteasome-to-immunoproteasome switch and class I MHC genes. Transcription of these genes continued after IFN removal and was Stat1 independent, suggesting the involvement of other signaling elements in addition to the canonical signal transduction pathway. Not less important were genes whose upregulation, in cases by many fold, is terminated once IFN is removed. Among these are viral sensing genes, such as retinoic acid-inducible gene-I protein (RIG-I), melanoma differentiation-associated gene 5 (MDA5) and toll-like receptor (TLR), cytokines, and apoptotic-related genes. Our findings provide a systemwide depiction of prolonged intracellular antiviral protection without the need for ongoing IFN stimulation.

Abstract: EmrE is an Escherichia coli H(+)-coupled multidrug transporter that provides a unique experimental paradigm because of its small size and stability, and because its activity can be studied in detergent solution. In this work, we report a study of the transient kinetics of substrate binding and substrate-induced proton release in EmrE. For this purpose, we measured transient changes in the tryptophan fluorescence upon substrate binding and the rates of substrate-induced proton release. The fluorescence of the essential and fully conserved Trp residue at position 63 is sensitive to the occupancy of the binding site with either protons or substrate. The maximal rate of binding to detergent-solubilized EmrE of TPP(+), a high-affinity substrate, is 2 x 10(7) M(-1).s(-1), a rate typical of diffusion-limited reactions. Rate measurements with medium- and low-affinity substrates imply that the affinity is determined mainly by the k(off) of the substrate. The rates of substrate binding and substrate-induced release of protons are faster at basic pHs and slower at lower pHs. These findings imply that the substrate-binding rates are determined by the generation of the species capable of binding; this is controlled by the high affinity to protons of the glutamate at position 14, because an Asp replacement with a lower pK is faster at the same pHs.

Abstract: The process of protein-protein association starts with their random collision, which may develop into an encounter complex followed by a transition state and final complex formation. Here we aim to experimentally characterize the nature of the transition state of protein-protein association for three different protein-protein interactions; Barnase-Barstar, TEM1-BLIP and IFNalpha2-IFNAR2, and use the data to model the transition state structures. To model the transition state, we determined inter-protein distance-constraints of the activated complex by using double mutant cycles (DMC) assuming that interacting residues are spatially close. Significant DeltaDeltaG(double dagger)(int) values were obtained only between residues on Barnase and Barstar. However, introducing specific mutations that optimize the charge complementarity between BLIP and TEM1 resulted in the introduction of significant DeltaDeltaG(double dagger)(int) values also between residues of these two proteins. While electrostatic interactions make major contributions towards stabilizing the transition state, we show two examples where steric hindrance exerts an effect on the transition state as well. To model the transition-state structures from the experimental DeltaDeltaG(double dagger)(int) values, we introduced a method for structure perturbation, searching for those inter-protein orientations that best support the experimental DeltaDeltaG(double dagger)(int) values. Two types of transition states were found, specific as observed for Barnase-Barstar and the electrostatically optimized TEM1-BLIP mutants, and diffusive as shown for wild-type TEM1-BLIP and IFNalpha2-IFNAR2. The specific transition states are characterized by defined inter-protein orientations, which cannot be modeled for the diffusive transition states. Mutations introduced through rational design can change the transition state from diffusive to specific. Together, these data provide a structural view of the mechanism allowing rates of association to differ by five orders of magnitude between different protein complexes.

Abstract: The 3D structures of macromolecules are difficult to grasp and also to communicate. By their nature, movies or animations are particularly useful for highlighting key features by offering a 'guided tour' of structures and conformation changes. However, high-quality movies are rarely seen because they are currently difficult and time consuming to make. By adopting the traditional movie 'storyboard' concept, which gives guidance and direction to filming, eMovie makes the creation of lengthy molecular animations much easier. This tool is a plug-in for the open-source molecular graphics program PyMOL, and enables experts and novices alike to produce informative and high-quality molecular animations.

Abstract: The formation of specific protein interactions plays a crucial role in most, if not all, biological processes, including signal transduction, cell regulation, the immune response and others. Recent advances in our understanding of the molecular architecture of protein-protein binding sites, which facilitates such diversity in binding affinity and specificity, are enabling us to address key questions. What is the amino acid composition of binding sites? What are interface hotspots? How are binding sites organized? What are the differences between tight and weak interacting complexes? How does water contribute to binding? Can the knowledge gained be translated into protein design? And does a universal code for binding exist, or is it the architecture and chemistry of the interface that enable diverse but specific binding solutions?

Abstract: All alpha-interferons (IFNalpha) bind the IFNAR1 receptor subunit with low affinity. Increasing the binding affinity was shown to specifically increase the antiproliferative potency of IFNalpha2. Here, we constructed a phage display library by randomizing three positions on IFNalpha2 previously shown to confer weak binding to IFNAR1. The tightest binding variant selected, comprised of mutations H57Y, E58N, and Q61S (YNS), was shown to bind IFNAR1 60-fold tighter compared with wild-type IFNalpha2, and 3-fold tighter compared with IFNbeta. Binding of YNS to IFNAR2 was comparable with wild-type IFNalpha2. The YNS mutant conferred a 150-fold higher antiproliferative potency in WISH cells compared with wild-type IFNalpha2, whereas its antiviral activity was increased by only 3.5-fold. The high antiproliferative activity was related to an induction of apoptosis, as demonstrated by annexin V binding assays, and to specific gene induction, particularly TRAIL. To determine the potency of the YNS mutant in a xenograft cancer model, we injected it twice a week to nude mice carrying transplanted MDA231 human breast cancer cells. After 5 weeks, no tumors remained in mice treated with YNS, whereas most mice treated with wild-type IFNalpha2 showed visible tumors. Histological analysis of these tumors showed a significant anti-angiogenic effect of YNS, compared with wild-type IFNalpha2. This work demonstrates the application of detailed biophysical understanding in the process of protein engineering, yielding an interferon variant with highly increased biological potency.

Abstract: MOTIVATION: Secondary structures are key descriptors of a protein fold and its topology. In recent years, they facilitated intensive computational tasks for finding structural homologues, fold prediction and protein design. Their popularity stems from an appealing regularity in patterns of geometry and chemistry. However, the definition of secondary structures is of subjective nature. An unsupervised de-novo discovery of these structures would shed light on their nature, and improve the way we use these structures in algorithms of structural bioinformatics. METHODS: We developed a new method for unsupervised partitioning of undirected graphs, based on patterns of small recurring network motifs. Our input was the network of all H-bonds and covalent interactions of protein backbones. This method can be also used for other biological and non-biological networks. RESULTS: In a fully unsupervised manner, and without assuming any explicit prior knowledge, we were able to rediscover the existence of conventional alpha-helices, parallel beta-sheets, anti-parallel sheets and loops, as well as various non-conventional hybrid structures. The relation between connectivity and crystallographic temperature factors establishes the existence of novel secondary structures.

Abstract: The expression, purification, detection, and assay of recombinant proteins have been made more convenient and rapid by the use of small affinity tags. To facilitate the purification of interferon-alpha2c (IFN-alpha2c) by metal chelate affinity chromatography, N-terminal 6-histidine tag was introduced via genetic manipulation. Two preparations of IFN material were purified; one contained IFN-alpha2c with the 6-histidine tag, and the other contained IFN-alpha2c without the 6-histidine tag. The antigenic properties of the human IFN-alpha2c subvariant with and without the 6-histidine tag, as well as the effects of the N-terminal 6-histidine tag on IFN-alpha2c interaction with the extracellular domain of human IFN-alpha receptor chain 2 (IFNAR2-EC) were examined. For the purposes of this study, IFNs were characterized by Western blots with anti-IFN monoclonal antibodies (mAb) and bioassays. Immunoblot analyses showed differences between IFN-alpha2c-6-histidine tag and IFN-alpha2a, b, c in their interaction with IFNAR2-EC. We also observed differences between IFN-alpha2c-6-histidine tag and IFN-alpha2a, b, c in bioactivities. This study is the first report that shows that an N-terminal 6-histidine tag on IFN-alpha2c can affect its interaction with receptor and cause a different bioactivity.

Abstract: Alpha and beta interferons (IFN-alpha and IFN-beta) are multifunctional cytokines that exhibit differential activities through a common receptor composed of the subunits IFNAR1 and IFNAR2. Here we combined biophysical and functional studies to explore the mechanism that allows the alpha and beta IFNs to act differentially. For this purpose, we have engineered an IFN-alpha2 triple mutant termed the HEQ mutant that mimics the biological properties of IFN-beta. Compared to wild-type (wt) IFN-alpha2, the HEQ mutant confers a 30-fold higher binding affinity towards IFNAR1, comparable to that measured for IFN-beta, resulting in a much higher stability of the ternary complex as measured on model membranes. The HEQ mutant, like IFN-beta, promotes a differentially higher antiproliferative effect than antiviral activity. Both bring on a down-regulation of the IFNAR2 receptor upon induction, confirming an increased ternary complex stability of the plasma membrane. Oligonucleotide microarray experiments showed similar gene transcription profiles induced by the HEQ mutant and IFN-beta and higher levels of gene induction or repression than those for wt IFN-alpha2. Thus, we show that the differential activities of IFN-beta are directly related to the binding affinity for IFNAR1. Conservation of the residues mutated in the HEQ mutant within IFN-alpha subtypes suggests that IFN-alpha has evolved to bind IFNAR1 weakly, apparently to sustain differential levels of biological activities compared to those induced by IFN-beta.

Abstract: De novo design and redesign of proteins and protein complexes have made promising progress in recent years. Here, we give an overview of how to use available computer-based tools to design proteins to bind faster and tighter to their protein-complex partner by electrostatic optimization between the two proteins. Electrostatic optimization is possible because of the simple relation between the Debye-Huckel energy of interaction between a pair of proteins and their rate of association. This can be used for rapid, structure-based calculations of the electrostatic attraction between the two proteins in the complex. Using these principles, we developed two computer programs that predict the change in k(on), and as such the affinity, on introducing charged mutations. The two programs have a web interface that is available at <webref type="url">www.weizmann.ac.il/home/bcges/PARE.html</webref> and <webref type="url">http://bip.weizmann.ac.il/hypare</webref>. When mutations leading to charge optimization are introduced outside the physical binding site, the rate of dissociation is unchanged and therefore the change in k(on) parallels that of the affinity. This design method was evaluated on a number of different protein complexes resulting in binding rates and affinities of hundreds of fold faster and tighter compared to wild type. In this chapter, we demonstrate the procedure and go step by step over the methodology of using these programs for protein-association design. Finally, the way to easily implement the principle of electrostatic design for any protein complex of choice is shown.

Abstract: Type I interferons (IFNs) elicit antiviral, antiproliferative and immunomodulatory properties in cells. All of them bind to the same receptor proteins, IFNAR1 and IFNAR2, with different affinities. While the 13 known IFNalphas are highly conserved, the C-terminal unstructured tail was found to have large variation in its net charge, from neutral to +4. This led us to speculate that the tail may have a role in modulation of the IFN biological activity, through fine-tuning the binding to IFNAR2. To evaluate this hypothesis, we replaced the tail of IFNalpha2 with that of IFNalpha8 and IFNbeta tails, or deleted the last five residues of this segment. Mutations to the more positively charged tail of IFNalpha8 resulted in a 20-fold higher affinity to IFNAR2, which results in a higher antiviral and antiproliferative activity. Double and multiple mutant cycle analysis placed the tail near a negatively charged loop on IFNAR2, comprising of residues Glu 132-134. Deleting the tail resulted in only twofold reduction in binding compared to the wild-type. Next, we modeled the location of the tail using a two-step procedure: first we generated 200 models of the tail docked on IFNAR2 using HADDOCK, second the models were scored according to the fit between experimentally determined rates of association of nine mutant complexes, and their calculated rates using the PARE software. From the results we suggest that the unstructured tail of IFNalpha is gaining a specific structure in the bound state, binding to a groove below the 132-134 loop in IFNAR2.

Abstract: Type I interferons activate cellular responses by forming a ternary complex with two receptor components, IFNAR1 and IFNAR2. While the binding of the IFNAR2 receptor to interferon is of high affinity and well characterized, the binding to IFNAR1 is weak, transient, and poorly understood. Here, we mapped the complete binding region of IFNAR1 on IFNalpha2 by creating a panel of 21 single alanine mutant proteins, and determined their binding affinities. The IFNAR1 binding site on IFNalpha2 maps to the center of the B and C helices, opposite to the binding site for IFNAR2. No hot spots for binding were found in the interface, with individual mutations having an up to fivefold effect on binding. Of the nine residues that affected binding, three adjacent conserved residues, located on the B helix, conferred an increase in the binding affinity to IFNAR1, as well as an increase in the biological activity of the interferon mutant. This suggests that binding of alpha interferons to the IFNAR1 receptor is sub-optimal. A correlation between binding affinity and biological activity was found, albeit not across the whole range of affinities. In WISH cells, but not DAUDI cells, the anti-proliferative activity was markedly affected by fluctuations in the IFNalpha2 affinity towards the IFNAR1 receptor. On the other hand, the antiviral activity of interferons on WISH cells seems to change in accordance to the binding affinity towards IFNAR1 only as long as the binding affinity is not beyond twofold of the wild-type. In accordance, the biological roles of the two interferon-receptor subunits are discussed.

Abstract: The rate of association of a protein complex is a function of an intrinsic basal rate and of the magnitude of electrostatic steering. In the present study we analyze the contribution of electrostatics towards the association rate of proteins in a database of 68 transient hetero-protein-protein complexes. Our calculations are based on an upgraded version of the computer algorithm PARE, which was shown to successfully predict the impact of mutations on k(on) by calculating the difference in Columbic energy of interaction of a pair of proteins. HyPare (http://bip.weizmann.ac.il/HyPare), automatically calculates the impact of mutations on a per-residue basis for all residues of a protein-protein interaction, achieving a precision similar to that of PARE. Our calculations show that electrostatics play a marginal role (<10 fold) in determining the rate of association for about half of the complexes in the database. Strong electrostatic steering, which results in an increase of over 100-fold in k(on), was calculated for about 25% of the complexes. Applying HyPare to all 68 complexes in the database shows that a small number of residues are hotspots for association. About 40% of the hotspots are calculated to increase the rate of association upon mutation, and thus increase binding affinity. This is a much higher ratio than found for hotspots for dissociation, where the large majority cause weaker binding. About 40% of the hotspots are located outside the physical boundary of the binding site, making them ideal candidates for protein engineering. Our data shows that a majority of protein-protein complexes are not optimized for fast association. Hotspots are not evenly distributed between all types of amino acids. About 75% of all hotspots are of charged residues. This is understandable, as a charge-reverse mutant changes the total charge by 2. The small number of hydrophobic residues that are hotspots upon mutation probably relates to their location and surrounding. For 18 out of the 68 complexes in the database, experimental values of k(on) are available. For these, a basal rate of association was calculated to be in the range of 10(4)M(-1)s(-1) to 10(7)M(-1)s(-1). Some of these rates were verified independently from experimental mutant data. The basal rates were correlated with the size of the proteins and the shape of the interface.

Abstract: The principal goal of the Israel Structural Proteomics Center (ISPC) is to determine the structures of proteins related to human health in their functional context. Emphasis is on the solution of structures of proteins complexed with their natural partner proteins and/or with DNA. To date, the ISPC has solved the structures of 14 proteins, including two protein complexes. It has adopted automated high-throughput (HTP) cloning and expression techniques and is now expressing in Escherichia coli, Pichia pastoris and baculovirus, and in a cell-free E. coli system. Protein expression in E. coli is the primary system of choice in which different parameters are tested in parallel. Much effort is being devoted to development of automated refolding of proteins expressed as inclusion bodies in E. coli. The current procedure utilizes tagged proteins from which the tag can subsequently be removed by TEV protease, thus permitting streamlined purification of a large number of samples. Robotic protein crystallization screens and optimization utilize both the batch method under oil and vapour diffusion. In order to record and organize the data accumulated by the ISPC, a laboratory information-management system (LIMS) has been developed which facilitates data monitoring and analysis. This permits optimization of conditions at all stages of protein production and structure determination. A set of bioinformatics tools, which are implemented in our LIMS, is utilized to analyze each target.

Abstract: The association of two proteins is preceded by a mutual diffusional search in solution. The role of translational and rotational diffusion in this process has been studied theoretically for many years. However, systematic experimental verification of theoretical results is still lacking. We report here measurements of association rates of the proteins beta-lactamase (TEM) and beta-lactamase inhibitor protein (BLIP) in solutions of glycerol and poly(ethylene glycol) of increasing viscosity. We also measured translational and rotational diffusion in the same solutions, using fluorescence correlation spectroscopy and fluorescence anisotropy, respectively. It is found that in glycerol both translational and rotational diffusion rates are inversely dependent on viscosity, as predicted by the classical Stokes-Einstein relations, while the association rate depends nonlinearly on viscosity. In contrast, the association rate depends only weakly on the viscosity of the polymer solutions, which results in a similar weak dependence of k(on) on viscosity. The data are modeled using the theory of diffusion-limited association. Deviations from the theory are explained by a short-range solute-induced repulsion between the proteins in glycerol solution and an attractive depletion interaction generated by the polymers. These results open the way to the creation of a unified framework for all nonspecific effects involved in the protein association process, as well as to better theoretical understanding of these effects. Further, they reflect on the complex factors controlling protein association within the crowded environment of cells and suggest that a high concentration of macromolecules does not significantly impede protein association.

Abstract: Protein-protein interactions are essential for life. Yet, our understanding of the general principles governing binding is not complete. In the present study, we show that the interface between proteins is built in a modular fashion; each module is comprised of a number of closely interacting residues, with few interactions between the modules. The boundaries between modules are defined by clustering the contact map of the interface. We show that mutations in one module do not affect residues located in a neighboring module. As a result, the structural and energetic consequences of the deletion of entire modules are surprisingly small. To the contrary, within their module, mutations cause complex energetic and structural consequences. Experimentally, this phenomenon is shown on the interaction between TEM1-beta-lactamase and beta-lactamase inhibitor protein (BLIP) by using multiple-mutant analysis and x-ray crystallography. Replacing an entire module of five interface residues with Ala created a large cavity in the interface, with no effect on the detailed structure of the remaining interface. The modular architecture of binding sites, which resembles human engineering design, greatly simplifies the design of new protein interactions and provides a feasible view of how these interactions evolved.

Abstract:

Abstract: Association of two proteins can be described as a two-step process, with the formation of an encounter complex followed by desolvation and establishment of a tight complex. Here, by using the computer algorithm PARE, we designed a set of mutants of the Ras effector protein Ral guanine nucleotide dissociation stimulator (RalGDS) with optimized electrostatic steering. The fastest binding RalGDS mutant, M26K,D47K,E54K, binds Ras 14-fold faster and 25-fold tighter compared with WT. A linear correlation was found between the calculated and experimental data, with a correlation coefficient of 0.97 and a slope of 0.65 for the 24 mutants produced. The data suggest that increased electrostatic steering specifically stabilizes the encounter complex and transition state. This conclusion is backed up by Phi analysis of the encounter complex and transition state of the RalGDS(M26K,D47K,E54K)/Ras complex, with both values being close to 1. Upon further formation of the final complex, the increased Coulombic interactions are probably counterbalanced by the cost of desolvation of charges, keeping the dissociation rate constant almost unchanged. This mechanism is also reflected by the mutual compensation of enthalpy and entropy changes quantified by isothermal titration calorimetry. The binding constants of the faster binding RalGDS mutants toward Ras are similar to those of Raf, the most prominent Ras effector, suggesting that the design methodology may be used to switch between signal transduction pathways.

Abstract: Prolonging the circulatory half-life of low mass protein drugs can be achieved either by administration of a pro-drug or through co-injection with a carrier protein that will slowly release the active protein. The rate of release is concentration and affinity dependent. The optimal relationship between these two in prolonging the half-life of a pro-drug is the focus of this work. Interferon (IFN) beta is one of the most widely used protein drugs in the clinic. Here, we show that the circulatory half-life of IFNbeta can be significantly extended by co-administration with the extracellular domain of the IFN receptor ifnar2 (ifnar2-EC). To investigate the concentration/affinity relation, a range of tighter binding ifnar2-EC mutants was designed that bind IFNbeta, but not IFNalpha2, up to 50-fold tighter compared with the wild-type ifnar2-EC. This increased affinity is related to a slower dissociation rate, whereas the association of IFNbeta with ifnar2-EC is already near optimum. Using the wild-type and mutant receptors, we investigated their potential in occluding IFNbeta from circulation in a tissue culture assay, as well as in rats. To determine the potential of ifnar2-EC as a carrier protein, we co-administered a mixture of IFNbeta and ifnar2-EC to rats both intravenously and subcutaneously, and followed the blood plasma concentrations of IFNbeta over time. The tighter binding ifnar2-EC mutant had a clear advantage in prolonging the half-life of IFNbeta in circulation, even when lower protein concentrations were administered. A numerical simulation of the in vivo data demonstrates that the optimal binding affinity of a carrier protein is around the concentration needed to obtain optimal activity of the ligand.

Abstract: Many peptide and protein drugs have a short circulatory half-life in vivo. The covalent attachment of polyethylene glycol (PEG) chains (PEGylation) can overcome this deficiency, but pegylated peptides and proteins are often inactive. In this study, we present a novel PEG-IFNalpha2 conjugate, PEG(40)-FMS-IFNalpha2, capable of regenerating native interferon alpha2 (IFNalpha2) at a slow rate under physiological conditions. A 2-sulfo-9-fluorenylmethoxycarbonyl (FMS) containing bifunctional reagent, MAL-FMS-NHS, has been synthesized, enabling the linkage of a 40 kDa PEG-SH to IFNalpha2 through a slowly hydrolyzable bond. By use of a BIAcore binding assay, the in vitro rate of regeneration of native interferon was estimated to have a half-life of 65 h. Following subcutaneous administration to rats and monitoring circulating antiviral activity, active IFNalpha2 levels peaked at 50 h, with substantial levels still being detected 200 h after administration. This value contrasts with a half-life of about 1 h measured for unmodified interferon. The concentration of active IFNalpha2 scaled linearly with the quantity injected. Comparing subcutaneous to intravenous administration of PEG(40)-FMS-IFNalpha2, we found that the long circulatory lifetime of IFNalpha2 was affected both by the slow rate of absorption of the PEGylated protein from the subcutaneous volume and by the slow rate of discharge from the PEG in circulation. A numerical simulation of the results was in good agreement with the results observed in vivo. The pharmacokinetic profile of this novel IFNalpha2 conjugate combines a prolonged maintenance in vivo with the regeneration of active-native IFNalpha2, ensuring ready access to peripheral tissues and thus an overall advantage over currently used formulations.

Abstract: Is the whole protein surface available for interaction with other proteins, or are specific sites pre-assigned according to their biophysical and structural character? And if so, is it possible to predict the location of the binding site from the surface properties? These questions are answered quantitatively by probing the surfaces of proteins using spheres of radius of 10 A on a database (DB) of 57 unique, non-homologous proteins involved in heteromeric, transient protein-protein interactions for which the structures of both the unbound and bound states were determined. In structural terms, we found the binding site to have a preference for beta-sheets and for relatively long non-structured chains, but not for alpha-helices. Chemically, aromatic side-chains show a clear preference for binding sites. While the hydrophobic and polar content of the interface is similar to the rest of the surface, hydrophobic and polar residues tend to cluster in interfaces. In the crystal, the binding site has more bound water molecules surrounding it, and a lower B-factor already in the unbound protein. The same biophysical properties were found to hold for the unbound and bound DBs. All the significant interface properties were combined into ProMate, an interface prediction program. This was followed by an optimization step to choose the best combination of properties, as many of them are correlated. During optimization and prediction, the tested proteins were not used for data collection, to avoid over-fitting. The prediction algorithm is fully automated, and is used to predict the location of potential binding sites on unbound proteins with known structures. The algorithm is able to successfully predict the location of the interface for about 70% of the proteins. The success rate of the predictor was equal whether applied on the unbound DB or on the disjoint bound DB. A prediction is assumed correct if over half of the predicted continuous interface patch is indeed interface. The ability to predict the location of protein-protein interfaces has far reaching implications both towards our understanding of specificity and kinetics of binding, as well as in assisting in the analysis of the proteome.

Abstract: Physiological media constitutes a crowded environment that serves as the field of action for protein-protein interaction in vivo. Measuring protein-protein interaction in crowded solutions can mimic this environment. In this work we follow the process of protein-protein association and its rate constants (k(on)) of the beta-lactamase (TEM)-beta-lactamase inhibitor protein (BLIP) complex in crowded solution using both low and high molecular mass crowding agents. In all crowded solutions (0-40% (w/w) of ethylene glycol (EG), poly(ethylene glycol) (PEG) 200, 1000, 3350, 8000 Da Ficoll-70 and Haemaccel the measured absolute k(on), but not k(off) values, were found to be slower as compared to buffer. However, there is a fundamental difference between low and high mass crowding agents. In the presence of low mass crowding agents and Haemaccel k(on) depends inversely on the solution viscosity. In high mass polymer solutions k(on) changes only slightly, even at viscosities 12-fold higher than water. The border between low and high molecular mass polymers is sharp and is dictated by the ratio between the polymer length (L) and its persistence length (Lp). Polymers that are long enough to form a flexible coil (L/Lp > 2) behave as high molecular mass polymers and those who are unable to do so (L/Lp < 2) behave as low molecular mass polymers. We concluded that although polymers solution are crowded, this property is not uniform; i.e. there are areas in the solution that contain bulk water, and in these areas proteins can diffuse and associate almost as if they were in diluted environment. This porous medium may be taken as mimicking some aspects of the cellular environment, where many of the macromolecules are organized along membranes and the cytoskeleton. To determine the contribution of electrostatic attraction between proteins in crowded milieu, we followed k(on) of wt-TEM and three BLIP analogs with up to 100-fold increased values of k(on) due to electrostatic steering. Faster associating BLIP variants keep their relative advantage in all crowded solutions, including Haemaccel. This result suggests that faster associating protein complexes keep their advantage also in complex environment.

Abstract: Docking algorithms produce many possible structures of a protein-protein complex. In most cases some of them resemble the correct structure within an r.m.s.d. of <3 A. A major challenge in the field of docking is to extract the correct structure out of this pool, the so-called 'scoring'. Here, we introduce a new scoring function, which discriminates between the many wrong and few true conformations. The scoring function is based on measuring the tightness of fit of the two docked proteins at a predicted binding interface. The location of the binding interface is identified using the recently developed computer algorithm ProMate. The new scoring function does not rely on energy considerations. It is therefore tolerant to low-resolution descriptions of the interface. A linear relation between the score and the r.m.s.d. relative to the 'true structure' is found in most of the cases evaluated. The function was tested on the docking results of 21 complexes in their unbound form. It was found to be successful in 77% of the examined cases, defining success as scoring a 'true' result with a p value of better than 0.1.

Abstract: The potent antiviral and antiproliferative activities of human type I interferons (IFNs) are mediated by a single receptor comprising two subunits, IFNAR1 and IFNAR2. The structure of the IFNAR2 IFN binding ectodomain (IFNAR2-EC), the first helical cytokine receptor structure determined in solution, reveals the molecular basis for IFN binding. The atypical perpendicular orientation of its two fibronectin domains explains the lack of C domain involvement in ligand binding. A model of the IFNAR2-EC/IFNalpha2 complex based on double mutant cycle-derived constraints uncovers an extensive and predominantly aliphatic hydrophobic patch on the receptor that interacts with a matching hydrophobic surface of IFNalpha2. An adjacent motif of alternating charged side chains guides the two proteins into a tight complex. The binding interface may account for crossreactivity and ligand specificity of the receptor. This molecular description of IFN binding should be invaluable for study and design of IFN-based biomedical agents.

Abstract: The structure of a protein-protein interaction, its affinity and thermodynamic characteristics depict a 'frozen' state of a complex. This picture ignores the kinetic nature of complex formation and dissociation, which are of major biological and biophysical interest. This review highlights recent advances in deciphering the kinetic pathway of protein-protein complexation, the nature of the encounter complex, transition state and intermediate along the reaction, and the effects of mutation, viscosity, pH and salt on association.

Abstract: The human interferon receptor (IFNAR) mediates the antiviral and antiproliferative activities of type I interferons (IFNs). This receptor is comprised of subunits IFNAR1 and IFNAR2, the latter exhibiting nanomolar affinity for IFNs. Here the extracellular domain of IFNAR2 (IFNAR2-EC), a soluble 25 kDa IFN-binding polypeptide, and its complex with IFN-alpha 2 were studied using multidimensional NMR. IFNAR2-EC is comprised of two fibronectin-III (FN-III) domains connected by a helical hinge region. The deduced global fold was utilized to improve the alignment of IFNAR2-EC against structurally related receptors and to model its structure. A striking feature of IFNAR2-EC is the limited and localized deviations in chemical shifts exhibited upon ligand binding, observed for only 15% of its backbone (1)H and (15)N nuclei. Analysis of these deviations maps the IFN-alpha 2 binding site upon IFNAR2-EC to a contiguous surface on the N-terminal domain, including the S3-S4 loop (residues 44-53), the S5-S6 loop and S6 beta-strand (residues 74-82), and the S7 beta-strand and the hinge region (residues 95-105). The C-terminal domain contributes only marginally to ligand binding, and no change in the hypothesized interdomain interface is observed. The proposed binding domain encompasses all residues implicated by mutagenesis studies in IFN binding, and suggests adjacent residues cooperate in forming the binding surface. D(2)O-exchange experiments indicate that binding of IFN-alpha2 induces tightening of the N-terminal domain of IFNAR2-EC. This increase in receptor rigidity may play an important role in initiating the intracellular stage of the IFN signaling cascade.

Abstract: Association of a protein complex follows a two-step mechanism, with the first step being the formation of an encounter complex that evolves into the final complex. Here, we analyze recent experimental data of the association of TEM1-beta-lactamase with BLIP using theoretical calculations and simulation. We show that the calculated Debye-Hückel energy of interaction for a pair of proteins during association resembles an energy funnel, with the final complex at the minima. All attraction is lost at inter-protein distances of 20 A, or rotation angles of >60 degrees from the orientation of the final complex. For faster-associating protein complexes, the energy funnel deepens and its volume increases. Mutations with the largest impact on association (hotspots for association) have the largest effect on the size and depth of the energy funnel. Analyzing existing evidence, we suggest that the transition state along the association pathway is the formation of the final complex from the encounter complex. Consequently, pairs of proteins forming an encounter complex will tend to dissociate more readily than to evolve into the final complex. Increasing directional diffusion by increasing favorable electrostatic attraction results in a faster forming and slower dissociating encounter complex. The possible applicability of electrostatic calculations for protein-protein docking is discussed.

Abstract: Polypeptide drugs are generally short-lived species in circulation. In this study, we have covalently linked seven moieties of 2-sulfo-9-fluorenylmethoxycarbonyl (FMS) to the amino groups of human interferon-alpha2. The derivative thus obtained (FMS(7)-IFN-alpha2) has approximately 4% the biological potency and 33 +/- 4% the receptor binding capacity of the native cytokine. Upon incubation, FMS(7)-IFN-alpha2 undergoes time-dependent spontaneous hydrolysis, generating active interferon with t(1/2) values of 24 +/- 2 h at pH 8.5 and 98 +/- 10 h at pH 7.4. When native IFN-alpha2 is intravenously administered to mice, circulating antiviral activity is maintained for a short duration and then declines with t(1/2) = 4 +/- 0.5 h, reaching undetectable values at approximately 18 h after administration. With intravenously administered FMS(7)-IFN-alpha2, there is a lag period of 2 h, followed by a progressive elevation in circulating antiviral-active protein, which peaked at 20 h and declined with t(1/2) = 35 +/- 4 h. FMS(7)-IFN-alpha2 is resistant to alpha-chymotrypsin digest and to proteolytic inactivation by human serum proteases in vitro. We have thus introduced here an inactive IFN-alpha2 derivative, which is resistant to in situ inactivation and has the capability of slowly reverting to the native active protein at physiological conditions in vivo and in vitro. Having these attributes, FMS(7)-IFN-alpha2 maintains prolonged circulating antiviral activity in mice, exceeding 7-8 times the activity of intravenously administered native cytokine.

Abstract: The pleiotropic activity of type I interferons has been attributed to the specific interaction of IFN with the cell-surface receptor components ifnar1 and ifnar2. To date, the structure of IFN has been solved, but not that of the receptor or the complex. In this study, the structure of the IFN-alpha 2-ifnar2 complex was generated with a docking procedure, using nuclear Overhauser effect-like distance constraints obtained from double-mutant cycle experiments. The interaction free energy between 13 residues of the ligand and 11 of the receptor was measured by double-mutant cycles. Of the 100 pairwise interactions probed, five pairs of residues were found to interact. These five interactions were incorporated as distance constraints into the flexible docking program prodock by using fixed and movable energy-gradient grids attached to the receptor and ligand, respectively. Multistart minimization and Monte Carlo minimization docking of IFN-alpha 2 onto ifnar2 converged to a well-defined average structure, with the five distance constraints being satisfied. Furthermore, no structural artifacts or intraloop energy strain were observed. The mutual binding sites on IFN-alpha 2 and ifnar2 predicted from the model showed an almost complete superposition with the ones determined from mutagenesis studies. Based on this structure, differences in IFN-alpha 2 versus IFN-beta binding are discussed.

Abstract: Association of a protein complex follows a two step reaction mechanism, with the first step being the formation of an encounter complex which evolves into the final complex. Here we present new experimental data for the association of the bacterial ribonuclease barnase and its polypeptide inhibitor barstar which shed light on the thermodynamics and structure of the transition state and preceding encounter complex of association at diminishing electrostatic attraction. We show that the activation entropy at the transition state is close to zero, with the activation enthalpy being equal to the free energy of binding. This observation was independent of the magnitude of the mutual electrostatic attraction, which were altered by mutagenesis or by addition of salt. The low activation entropy implies that the transition state is mostly solvated at all ionic strengths. The structure of the transition state was probed by measuring pairwise interaction energies using double-mutant-cycles. While at low ionic strength all proximal charge-pairs form contacts, at high salt only a subset of these interactions are maintained. More specifically, charge-charge interactions between partially buried residues are lost, while exposed charged residues maintain their ability to form specific interactions even at the highest salt concentration. Uncharged residues do not interact at any ionic strength. The results presented here suggest that the barnase-barstar binding sites are correctly aligned during the transition state even at diminishing electrostatic attraction, although specific short range interactions of uncharged residues are not yet formed. Furthermore, most of the interface desolvation (which contributes to the entropy of the system) has not yet occurred. This picture seems to be valid at low and high salt. However, at high salt, interactions of the activated complex are limited to a more restricted set of residues which are easier approached during diffusion, prior to final docking. This suggest that the steering region at high salt is more limited, albeit maintaining its specificity.

Abstract:

Abstract: Investigating protein-protein interactions by mutational analysis requires practical techniques for quantifying rate constants and equilibrium constants over several orders of magnitude with reasonably high sample throughput. We have employed spectroscopic interferometry for label-free monitoring of the interaction between the cytokine interferon alpha2 (IFNalpha2) and the extracellular domain of its receptor ifnar2 (ifnar2-EC). We implemented a versatile surface chemistry for the glass substrate of this transducer for covalent immobilization of proteins. Affinity capturing with a monoclonal anti-ifnar2-EC antibody (mAb) followed by crosslinking with a second, noncompetitive mAb provided stable, but still reversible, immobilization of ifnar2-EC. We measured kinetics and affinity of numerous of mutants of IFNalpha2 and ifnar2-EC. Dissociation rate constants up to 0.3 s(-1) and association rate constants up to 3 x 10(6) M(-)1 s(-1) were resolved by the system. Dissociation constants down to 200 microM were measured with protein concentrations up to 50 microM without no background signal or nonspecific binding. The instrument detection limit is approximately 10 pm without the need for temperature stabilization or referencing channels. The system proved effective for large-scale mutational analysis involving alanine scanning mutagenesis and double mutant cycles.

Abstract: An experimental approach to evaluate the net binding free energy of buried hydrogen bonds and salt bridges is presented. The approach, which involves a modified multiple-mutant cycle protocol, was applied to selected interactions between TEM-1-beta-lactamase and its protein inhibitor, BLIP. The selected interactions (two salt bridges and two hydrogen bonds) all involving BLIP-D49, define a distinct binding unit. The penta mutant, where all side-chains constructing the binding unit were mutated to Ala, was used as a reference state to which combinations of side-chains were introduced. At first, pairs of interacting residues were added allowing the determination of interaction energies in the absence of neighbors, using double mutant cycles. Addition of neighboring residues allowed the evaluation of their cooperative effects on the interaction. The two isolated salt bridges were either neutral or repulsive whereas the two hydrogen bonds contribute 0.3 kcal mol(-1 )each. Conversely, a double mutant cycle analysis of these interactions in their native environment showed that they all stabilize the complex by 1-1.5 kcal mol(-1). Examination of the effects of neighboring residues on each of the interactions revealed that the formation of a salt bridge triad, which involves two connected salt bridges, had a strong cooperative effect on stabilizing the complex independent of the presence or absence of additional neighbors. These results demonstrate the importance of forming net-works of buried salt bridges. We present theoretical electrostatic calculations which predict the observed mode of cooperativity, and suggest that the cooperative networking effect results from the favorable contribution of the protein to the interaction. Furthermore, a good correlation between calculated and experimentally determined interaction energies for the two salt bridges, and to a lesser extent for the two hydrogen bonds, is shown. The data analysis was performed on values of DeltaDeltaG(double dagger)K(d) which reflect the strength of short range interactions, while DeltaDeltaG(o)K(D) values which include the effects of long range electrostatic forces that alter specifically DeltaDeltaG(double dagger)k(a) were treated separately.

Abstract: A protein design strategy was developed to specifically enhance the rate of association (k(on)) between a pair of proteins without affecting the rate of dissociation (k(off)). The method is based on increasing the electrostatic attraction between the proteins by incorporating charged residues in the vicinity of the binding interface. The contribution of mutations towards the rate of association was calculated using a newly developed computer algorithm, which predicted accurately the rate of association of mutant protein complexes relative to the wild type. Using this design strategy, the rate of association and the affinity between TEM1 beta-lactamase and its protein inhibitor BLIP was enhanced 250-fold, while the dissociation rate constant was unchanged. The results emphasize that long range electrostatic forces specifically alter k(on), but do not effect k(off). The design strategy presented here is applicable for increasing rates of association and affinities of protein complexes in general.

Abstract: Type I interferons bind to two cell surface receptors, ifnar1 and ifnar2, as the first step in the activation of several signal transduction pathways that elicit an anti-viral state and an anti-proliferative response. Here, we quantitatively mapped the complete binding region of ifnar2 on interferon (IFN)alpha2 by 35 individual mutations to alanine and isosteric residues. Of the six "hot-spot" residues identified (Leu-30, Arg-33, Arg-144, Ala-145, Met-148, and Arg-149), four are located on the E-helix, which is located at the center of the binding site flanked by residues on the A-helix and the AB-loop. The contribution of residues of the D-helix, which have been previously implicated in binding, proved to be marginal for the interaction with the extracellular domain of ifnar2. Interestingly, the ifnar2 binding site overlaps the largest continuous hydrophobic patch on IFNalpha2. Thus, hydrophobic interactions seem to play a significant role stabilizing this interaction, with the charged residues contributing toward the rapid association of the complex. Relating the anti-viral and anti-proliferative activity of the various interferon mutants with their affinity toward ifnar2 results in linear function over the whole range of affinities investigated, suggesting that ifnar2 binding is the rate-determining step in cellular activation. Dose-time analysis of the anti-viral response revealed that shortening the incubation time of low-level activation cannot be compensated by higher IFN doses. Considering the strict dependence of the cellular response on affinity, these results suggest that for maintaining transcription of IFN-responsive genes over a longer time period, low but continuous signaling through the IFN receptor is essential.

Abstract: BLIP is a secreted protein from Streptomyces clavuligerus that inhibits a wide range of beta-lactamases. Here we investigate the tight interaction of BLIP, expressed heterologousely in E. coli, with TEM-1. Kinetic and thermodynamic constants were determined using methods with the proteins either in a homogeneous or in a heterogeneous phase. While values of Delta DeltaG(mut-wt) are similar whether measured by fluorescence quench, enzyme inhibition, or surface plasmon resonance, absolute values of DeltaG and kinetic constants vary. Association and dissociation rate constants of 10(5) M-1 s-1 and 10(-)4 s-1, respectively, and a nanomolar affinity were determined for the wild-type proteins. The highest affinity is measured at pH 7.5, with a decreasing association rate constant at higher pH values, and an increasing dissociation rate constant at lower pH values. The marginal effect of salt on the kinetics of binding, as well as the calculated surface potentials, suggests a limited role for electrostatic forces in guiding this reaction. Still, mutations of interfacial residues affect the rate of association significantly, so that an increase in the net negative charge on either protein reduces the association rate constant. We show that simple electrostatic rules can explain this behavior. BLIP inhibits the catalytic activity of TEM-1 by binding its active site. Yet, mutations of active site residues on TEM-1 only have a moderate though cooperative effect on the binding energy. This can be explained in light of the peripheral location of the active site in the interface between the two proteins.

Abstract: Type I interferons (IFN) exert pleiotropic activities through binding to two cell surface receptors, ifnar1 and ifnar2. We are investigating the biophysical basis of IFN signaling by characterizing the complex of the extra-cellular domain of ifnar2 (ifnar2-EC) with IFNs on the level of purified recombinant proteins in vitro. Here, we present a detailed mutational study on the functional epitopes on both IFN and ifnar2. Kinetic and thermodynamic parameters were determined by label-free heterogeneous phase detection. On IFNalpha2, a relatively small functional epitope comprising ten amino acid residues was localized, which is nearly entirely formed by residues on the AB loop. Two hot-spot residues, L30 and R33, account for two-thirds of the total interaction energy. Comparing the anti-viral potency of the various mutants to the binding affinity towards ifnar2 revealed a proportional correlation between the two, suggesting a rate-limiting role of ifnar2 binding in IFN signaling. On ifnar2, residues T46, I47 and M48 were identified as hot-spots in the interaction with IFNalpha2. For another ten residues on ifnar2, significant contribution of interaction energy was determined. Based on these data, the functional epitope on ifnar2 was defined according to a homology model based on other members of the class II hCR family in good agreement with the complementary functional epitope on IFNalpha2. Although IFNalpha2 and IFNbeta bind competitively to the same functional epitope, mutational analysis revealed distinct centers of binding for these IFNs on ifnar2. This small shift of the binding site may result in different angular orientation, which can be critically coupled to cytoplasmic signaling.

Abstract: The rate of association of proteins is dictated by diffusion, but can be enhanced by favorable electrostatic forces. Here the relationship between the electrostatic energy of interaction, and the kinetics of protein-complex formation was analyzed for the protein pairs of: hirudin-thrombin, acetylcholinesterase-fasciculin and barnase-barstar, and for a panel of point mutants of these proteins. Electrostatic energies of interaction were calculated as the difference between the electrostatic energy of the complex and the sum of the energies of the two individual proteins, using the computer simulation package DelPhi. Calculated electrostatic energies of interaction were compared to experimentally determined rates of association. One kcal/mol of Coulombic interaction energy increased the rate of association by a factor of 2.8, independent of the protein-complex or mutant analyzed. Electrostatic energies of interaction were also determined from the salt dependence of the association rate constant, using the same basic equation as for the theoretical calculation. A Brönsted analysis of the electrostatic energies of interactions plotted versus experimentally determined ln(rate)s of association shows a linear relation between the two, with a beta value close to 1. This is interpreted as the energy of the transition state varies according to the electrostatic interaction energy, fitting a two state model for the association reaction. Calculating electrostatic rate enhancement from the electrostatic interaction energy can be used as a powerful tool to design protein complexes with altered rates of association and affinities.

Abstract: Type I interferons are cytokines which activate an anti-viral response by binding to two specific cell surface receptors, ifnar1 and ifnar2. Here, we report purification and refolding of the extracellular part of human ifnar2 (ifnar2-EC) expressed in Escherichia coli and its characterization with respect to its interaction with interferon alpha2 (IFNalpha2). The 25 kDa, non-glycosylated ifnar2-EC is a stable, fully active protein, which inhibits antiviral activity of IFNalpha2. The stoichiometry of binding IFNalpha2 is 1:1, as determined by gel filtration, chemical cross-linking and solid-phase detection. The affinity of this interaction is 10 nM, which is similar to the affinity measured for the cell surface-bound ifnar2 receptor. No difference in affinity was found throughout various assays using optical detection as BIAcore or reflectometric interference spectorscopy. However, the binding kinetics as measured in homogeneous phase by fluorescence de-quenching was about three times faster than that measured on a sensor surface. The rate of complex formation is relatively high compared to other cytokine-receptor interactions. The salt dependence of the association kinetics suggest a limited but significant contribution of electrostatic forces towards the rate of complex formation. The dissociation constant increases with decreasing pH according to the protonation of a base with a pKa of 6.7. The surface properties of the IFNalpha2 binding surface on ifnar2 were interpreted according to the pH and salt dependence of the interaction.

Abstract: The electrostatic enhancement of the association rate of barnase and barstar is calculated using a transition-state theory like expression and atomic-detail modeling of the protein molecules. This expression predicts that the rate enhancement is simply the average Boltzmann factor in the region of configurational space where association occurs instantaneously in the diffusion-controlled limit. Based on experimental evidence, this "transition state" is defined by configurations in which, relative to the stereospecifically bound complex, the two proteins are shifted apart by approximately 8 A (so a layer of water can be accommodated in the interface) and the two binding surfaces are rotated away by 0 degrees to 3 degrees. The values of the average Boltzmann factor, calculated by solving the Poisson-Boltzmann equation, for the wild-type complex and 16 complexes with single mutations are found to correlate well with experimental results for the electrostatic rate enhancement. The predicted rate enhancement is found to be somewhat insensitive to the precise definition of the transition state, due to the long-range nature of electrostatic interactions. The experimental ionic strength dependence of the rate enhancement is also reasonably reproduced.

Abstract: Barnase, a small extracellular ribonuclease from Bacillus amyloliquefaciens and its intracellular inhibitor barstar have co-evolved to bind tightly and rapidly. Barnase has also evolved to be catalytically active. The active site of barnase and its binding site for barstar use the same subset of amino acids. The exception is Glu73 (the general base in catalysis), which although located at the centre of the binding site, is separated by three ordered water molecules from barstar. We examined in this work the contribution of Glu73 to both catalysis and barstar binding. Truncation mutants of the general base (Glu73 --> Ala or Ser) retain a residual RNase activity of about 0.3% while mutants with larger hydrophobic replacements (Glu 73 --> Trp or Phe) have virtually no catalytic activity. This, and binding data of 3'-GMP with the different barnase mutants suggest that the loss in activity results from the elimination of the general base, which can be substituted to some extent by water or other polar side-chains in truncation mutants. All of the Glu73 mutations lead to a weakening of the free energy of complex formation with barstar by 1.4 to 3.0 kcal/mol (including Gln). This is surprising, since Glu73 does not interact directly with barstar and there is an electrostatic repulsion between Glu73 on barnase and the negatively charged binding surface of barstar. A newly developed method of constructing double mutant cycles between multiple mutations at the same site appears to pinpoint a favourable interaction between Glu73 and one of its nearest neighbours in barstar, Asp39. The coupling energy between those residues is presumably indirect: the carboxylate of Glu73 organizes neighbouring positively charged groups in barnase, Lys27, Arg83, and Arg87 to interact with Asp39 in barstar. This emphasizes that an apparent interaction between a pair of residues as measured with double mutant cycles is the sum of their direct and indirect interactions.

Abstract: We have studied the thermodynamics of the interaction between the ribonuclease barnase and its natural polypeptide inhibitor barstar. The contribution of specific residues and interactions within the barnase-barstar interface to the enthalpy of binding has been examined using isothermal titration calorimetry and protein engineering. The enthalpy of association of the wild-type proteins is -18.9 (+/-0.1) kcal/mol at pH 8 and at 25 degrees C. The enthalpy of binding remains favourable for 31 different combinations of mutations in the interface. The effects on the binding enthalpy upon replacing a side-chain involved in the interaction of barnase and barstar are, however, always unfavourable and in most cases larger than the effects on the free energy of binding. Interaction enthalpies calculated by double mutant cycle analysis are in some cases much larger than the interaction free energies. The interaction enthalpies for complexes between different barnase mutants with amino acid substitutions of the general base residue glutamic acid 73 and a barstar variant (D39A) vary by as much as 8.3 kcal/mol while the coupling free energies differ only by 1 kcal/mol. The use of enthalpies for the analysis of structure-activity relationships appears to be complicated by enthalpy-entropy compensation of weak intermolecular interactions. These tend to cancel out in measurements of free energy, which is thus the preferred quantity for simple analysis of interactions.

Abstract: We have documented the folding pathway of the 10-kDa protein barstar from the first few microseconds at the resolution of individual residues from its well characterized denatured state. The denatured state had been shown from NMR to have flickering native-like structure in the first two of its four alpha-helices. phi-value analysis shows that the first helix becomes substantially consolidated as the intermediate is formed in a few hundred microseconds, as does the second to a lesser extent. A native-like structure then is formed in a few hundred milliseconds as the whole structure consolidates. Peptide fragments corresponding to sequences containing the first two helices separately and together as a helix-loop-helix motif have little helical structure under conditions that favor folding. The early stages of folding fit the nucleation-condensation model that was proposed for the smaller chymotrypsin inhibitor 2, which is a single module of structure and folds by two-state kinetics. The early stages of the multistate folding of the larger, multimodular, barnase have proved experimentally inaccessible. The folding pathway of barstar links those of CI2 and barnase to give a unified scheme for folding.

Abstract: The rapid association of barnase and its intracellular inhibitor barstar has been analysed from the effects of mutagenesis and electrostatic screening. A basal association rate constant of 10(5) M(-1) s(-1) is increased to over 5 x 10(9) M(-1) s(-1) by electrostatic forces. The association between the oppositely charged proteins proceeds through the rate-determining formation of an early, weakly specific complex, which is dominated by long-range electrostatic interactions, followed by precise docking to form the high affinity complex. This mode of binding is likely to be used widely in nature to increase association rate constants between molecules and its principles may be used for protein design.

Abstract: The polypeptide inhibitor of the ribonuclease barnase, barstar, has two cysteine residues in positions 40 and 82. These have been proposed to form a disulfide bridge leading to an increase in stability without changing the inhibitory activity of the protein. Barstar and a mutant (E80A) were oxidized in vitro and the biochemical and physico-chemical properties of the oxidized monomers were analysed. The oxidized proteins show no inhibition of barnase using a plate assay and are significantly destabilized. CD spectra indicate a loss of secondary structure. The amino acid substitution E80 --> A stabilizes the oxidized barstar to about the same extent as it does the reduced protein, indicating, however, that the helical region which it is in is intact.

Abstract: ppGpp serves as an alarmon in prokaryotes, distributing and coordinating different cellular processes according to the nutritional potential of the growth medium. This work is interpreted as favoring the view that, in addition to its previously documented role in regulating the rate of ribosome synthesis, ppGpp participates in coordinating DNA replication and cell division. We studied the effects of ppGpp on the cell division cycle, using cells containing plasmid pSM11 that codes for the 55-kDa truncated RelA protein under the inducible Ptac promoter. In this system it was found that the rate of initiation of new rounds of DNA replication is inversely correlated with the intracellular level of ppGpp. Furthermore, ppGpp levels similar to those found during the activation of stringent control inhibited replication initiation, in a manner comparable to that resulting from inhibition of protein synthesis by amino acid starvation or by chloramphenicol addition. However, in contrast to chloramphenicol treatment, elevated ppGpp levels did not block septum formation, and, in fact, there is some evidence for enhanced septation. As a result, the residual cell division following elevation in ppGpp levels was higher than after chloramphenicol treatment, resulting in cells with a size similar to that of stationary phase cells.

Abstract: The temperature-induced unfolding of single, double, and triple mutants of barstar, the specific intracellular protein inhibitor of barnase from Bacillus amyloliquefaciens, has been studied by high-sensitivity differential scanning calorimetry. The thermal unfolding of barstar mutants, where at least one of the two cysteine residues in the molecule had been replaced by alanine, follows a two-state mechanism at neutral and alkaline pH. The unfolding enthalpy and heat capacity changes are slightly lower than those accepted for highly compact, small, globular proteins. We have found that at pH 2.5, where barstar seems to be in a molten globule state, the protein has a heat capacity between that of the native and the unfolded states and shows some tendency for association. Scanning calorimetry experiments were also extended to the barstar--barnase complex in the neutral and alkaline pH range. The binding constants obtained from DSC studies are similar to those already obtained from other (kinetic) studies. The interaction of barstar and barnase was also investigated by isothermal calorimetry in various buffers within the pH range 6.0-10.0 and a temperature range of 15-35 degrees C. The favorable enthalpy contribution to the binding is about 4 times higher than the entropic one at 25 degrees C. The overall data analysis of the combined calorimetric results has led to the thermodynamic characterization of barstar unfolding and the interaction of barstar and barnase over a wide range of temperatures.

Abstract: The interaction of barnase, an extracellular RNase of Bacillus amylolique-faciens, with its intracellular inhibitor barstar is a suitable paradigm for protein-protein interactions, since the structures of both the free and the complexed proteins are available at high resolution. The contributions of residues from both proteins to the energetics of kinetics and thermodynamics of binding were measured by double mutant cycle analysis. Such cycles reveal whether the contributions from a pair of residues are additive, or the effects of mutations are coupled. The aim of the study was to determine which of the interactions are co-operative. Double mutant cycles were constructed between a subset of five barnase and seven barstar residues, which were shown by structural and mutagenesis studies to be important in stabilising the complex. The coupling energy between two residues was found to decrease with the distance between them. Generally, residues separated by less than 7 A interact co-operatively. At greater separations, the effects of mutation are additive, and the energetics of the interactions are independent of each other. The highest coupling energies are found between pairs of charged residues (1.6 to 7 kcal mol-1). Three of the six most important interactions detected by double mutant cycle analysis (with coupling energies of more than 3.0 kcal mol-1) had not been noted previously from examination of the crystal structure. The effects of mutation on the kinetics of association are all additive, apart from charged residues located at distances of up to 10 A apart, which are co-operative. This can be explained by the fact that the transition state for association occurs before most interactions are formed.

Abstract: BACKGROUND: Barstar is the intracellular inhibitor of barnase, an extracellular RNAse of Bacillus amyloliquefaciens. The dissociation constant of the barnase-barstar complex is 10(-14) M with an association rate constant between barnase and barstar of 3.7 x 10(8) s-1 M-1. The rapid association arises in part from the clustering of four acidic residues (Asp35, Asp39, Glu76 and Glu80) on the barnase-binding surface of barstar. The negatively charged barnase-binding surface of barstar effectively 'steers' the inhibitor towards the positively charged active site of barnase. RESULTS: Mutating any one of the four acidic side chains of barstar to an alanine results in an approximately two-fold decrease in the association rate constant, while the dissociation rate constant increases from five orders of magnitude for Asp39-->Ala, to no significant change for Glu80-->Ala. The stability of barstar is increased by all four mutations, the increase ranging from 0.3 kcal mol-1 for Asp35-->Ala or Asp39-->Ala, to 2.1 kcal mol-1 for Glu80-->Ala. CONCLUSIONS: The evolutionary pressure on barstar for rapid binding of barnase is so strong that glutamate is preferred over alanine at position 80, even though it does not directly interact with barnase in the complex and significantly destabilizes the inhibitor structure. This, and other examples from the literature, suggest that proteins evolve primarily to optimize their function in vivo, with relatively little evolutionary pressure to increase stability above a certain threshold, thus allowing greater latitude in the evolution of enzyme activity.

Abstract: Transcriptional control of the himA and the himD/hip genes coding for the two subunits of the integration host factor (IHF) was investigated. The promoters for the two genes were identified by the use of primer extension and S1 analysis. Expression from both promoters was found to increase as the cells enter stationary phase. Mutation in rpoS, known to be induced upon entry to stationary phase, dramatically reduced the growth-phase response of the himA P4 promoter but had only a small effect on the induction of the himD/hip promoter. The increased activity of both promoters required the presence of the relA and spoT genes, suggesting that ppGpp plays a major role in the response to stationary phase. An artificial increase in ppGpp in exponentially growing cells induced a rapid increase in himA P4 and himD/hip mRNA levels. Experiments with a mutant defective in rpoS showed that the response of the himA P4 promoter to high ppGpp levels was greatly reduced while that of himD/hip was only slightly affected. Therefore, it seems that different mechanisms involving RpoS and ppGpp regulate the growth-phase response of the two promoters. We propose that the effect of ppGpp on himA P4 is mediated via RpoS whereas the himD/hip promoter is affected by ppGpp independently of RpoS. Expression of the himD/hip and himA genes was found to be subject to negative autoregulation. IHF-binding sites, implicated in autoregulation, were found to overlap both the himD/hip and himA P4 promoters. An additional IHF-binding site was found upstream of the himD/hip promoter. All three sites show low binding affinity to IHF suggesting that autoregulation can take place only after sufficiently high levels of IHF accumulate in the cell.

Abstract: We have solved, refined, and analyzed the 2.0-Ã¥ resolution crystal structure of a 1:1 complex between the bacterial ribonuclease, barnase, and a Cys-->Ala(40,82) double mutant of its intracellular polypeptide inhibitor, barstar. Barstar inhibits barnase by sterically blocking the active site with a helix and adjacent loop segment. Almost half of the 14 hydrogen bonds between barnase and barstar involve two charged residues, and a third involve one charged partner. The electrostatic contribution to the overall binding energy is considerably greater than for other protein-protein interactions. Consequently, the very high rate constant for the barnase-barstar association (10(8) s-1 M-1) is most likely due to electrostatic steering effects. The barnase active-site residue His102 is located in a pocket on the surface of barstar, and its hydrogen bonds with Asp39 and Gly31 residues of barstar are directly responsible for the pH dependence of barnase-barstar binding. There is a high degree of complementarity both of the shape and of the charge of the interacting surfaces, but neither is perfect. The surface complementarity is slightly poorer than in protease-inhibitor complexes but a little better than in antibody-antigen interactions. However, since the burial of solvent in the barnase-barstar interface improves the fit significantly by filling in the majority of gaps, as well as stabilizing unfavorable electrostatic interactions, its role seems to be more important than in other protein-protein complexes. The electrostatic interactions between barnase and barstar are very similar to those between barnase and the tetranucleotide d(CGAC). In the barnase-barstar complex, the two phosphate-binding sites in the barnase active site are occupied by Asp39 and Gly43 of barstar. However, barstar has no equivalent for a guanine base of an RNA substrate, resulting in the occupation of the guanine recognition site in the barnase-barstar complex by nine ordered water molecules. Upon barnase-barstar binding, entropy losses resulting from the immobilization of segments of the protein chain and the energetic costs of conformational changes are minimized due to the essentially preformed active site of barnase. However, a certain degree of flexibility within the barnase active site is required to allow for the structural differences between barnase-barstar binding and barnase-RNA binding. A comparison between the bound and the free barstar structure shows that the overall structural response to barnase-binding is significant. This response can be best described as outwardly oriented, rigid-body movements of the four alpha-helices of barstar, resulting in the structure of bound barstar being somewhat expanded.

Abstract: Barnase, an extracellular ribonuclease of Bacillus amyloliquefaciens, forms a very tight complex with its intracellular polypeptide inhibitor barstar. At pH 8, the values for the rate constants k1 (association) and k-1 (dissociation) are 6.0 x 10(8) s-1 M-1 and 8.0 x 10(-6) s-1, respectively. The value of Ki, the dissociation constant of barstar and barnase, calculated from the ratio k-1/k1 is 1.3 x 10(-14) M, which corresponds to a delta G of -18.9 kcal/mol at 25 degrees C. The dissociation constant increases with decreasing pH according to the ionization of an acid in free barnase of pKa 6.4, with very weak, if any, binding to the protonated form. This pH dependence for dissociation of the complex can be attributed almost entirely to residue His102 in barnase, as determined by a His102-->Ala mutation. Analysis of the pH dependence of the kinetic constants indicates that binding is, at least, a two-step process. The first, and rate-determining, step is association at close to the diffusion-controlled rate. There is then the precise docking of the complex. The value of Ki increases to 2.4 x 10(-11) M in the presence of 500 mM NaCl, and to 1.6 x 10(-11) M at pH 5 (100 mM NaCl). The binding site of barstar on barnase was mapped by measuring the values of Ki for a broad range of site-specific mutants of barnase. Mutagenesis of residues Lys27, Arg59, Arg87, and His102 to Ala increases the values of Ki by a factor of 10(4).(ABSTRACT TRUNCATED AT 250 WORDS)

Abstract: Barstar, a small globular protein which undergoes reversible unfolding, is a good candidate for studies on protein folding. It possesses two cysteine residues that complicate folding studies by forming a variable mixture of disulfide-bridged forms. We have constructed and analyzed, therefore, a double mutant Cys40-->Ala,Cys82-->Ala. Equilibrium unfolding with urea follows a simple two-step mechanism. The midpoint for unfolding ([U]1/2) is 3.87 +/- 0.03 M urea, with m(delta delta G/delta [urea]) = 1.25 +/- 0.04 kcal/mol2. The free energy of unfolding, delta GU-FH2O, is 4.84 +/- 0.18 kcal/mol. Identical results were found on monitoring the intrinsic tryptophan fluorescence or the circular dichroism signal at 221 nm, showing that the transition is due to the global denaturation of the protein. Barstar contains two proline residues, one of which (Pro48) has a cis N-aminoacyl bond conformation in the folded state. A transiently generated form of the unfolded protein, which contains the proline residues in their native conformations, has a rate constant for refolding (31 s-1) similar to that for refolding of the equilibrium-unfolded protein, which results in a "misfolded" form of the protein (32 s-1). The two refolded states are different: the free energies of unfolding measured from kinetic constants for the native and misfolded variants are 5.4 +/- 0.3 and 2.85 +/- 0.1 kcal/mol, respectively. The rate constant for the unfolding in water of the misfolded protein is 0.87 s-1, compared with 0.068 s-1 for the unfolding of the native protein. This difference can be explained by a nonnative trans peptidyl-proline bond at position 48 in the misfolded protein.(ABSTRACT TRUNCATED AT 250 WORDS)

Abstract: Intracellular levels of guanosine 3',5'-bispyrophosphate (ppGpp) governed by the relA gene are normally regulated by aminoacyl-tRNA availability for protein synthesis. An experimental system is described in which cellular levels of ppGpp are controlled instead by induction of plasmid pKK223-3 derivatives with the relA structural gene, or portions thereof, under control of the Ptac promoter. In amino acid-rich media, isopropyl-1-thio-beta-D-galactopyranoside induction of transcription of the wild type relA gene in pSM10 yields about a 100-fold overexpression of a metabolically stable, full length (743 amino acid) RelA protein to levels approximating the number of cellular ribosomes. This overexpression is accompanied by a roughly parallel and relC-dependent elevation of ppGpp levels. Induction of a relA gene deletion mutant in pSM11 containing 455 amino-terminal amino acids results in much lower levels of expression of a metabolically unstable 55-kDa protein and elevated ppGpp levels that are almost equivalent to induced pSM10 and are relC-independent. Induction of a larger deletion in pSM12 containing 331 amino-terminal amino acids does not provoke ppGpp accumulation. We are able to elicit high levels of ppGpp without changing nutritional abundance and without massive overexpression of the RelA protein by inducing the metabolically unstable, truncated RelA protein. We find the effects of elevated ppGpp levels to include a slowing of growth, an inhibition of stable RNA accumulation, an inhibition of cellular rrn P1 promoter activities as measured by primer extension, and changes in the pattern of gene expression viewed by two-dimensional electrophoresis of cellular proteins.

Abstract: The most widely studied "relaxed" mutant of the relA locus, the relA1 allele, is shown here to consist of an IS2 insertion between the 85th and 86th codons of the otherwise wild-type relA structural gene, which normally encodes a 743-amino acid (84 kDa) protein. The RelA protein is a ribosome-dependent ATP:GTP (GDP) pyrophosphoryltransferase that is activated during the stringent response to amino acid starvation and thereby occasions the accumulation of guanosine 3',5'-bispyrophosphate (ppGpp). We propose that the IS2 insertion functionally splits the RelA protein into two (alpha and beta) peptide fragments which can complement each other in trans to yield residual ppGpp synthetic activity; neither fragment shows this activity when expressed alone. Cell strains with a single copy relA null allele show physiological behavior that is much the same as relA1 mutant strains. Both relA1 and relA null strains accumulate ppGpp during glucose starvation and do not accumulate ppGpp during the stringent response. The presence of ppGpp in verifiable relA null strains is interpreted as unequivocal evidence for an alternate route of ppGpp synthesis that exists in addition to the relA-dependent reaction.

Abstract: The relA gene product of Escherichia coli is known to be responsible for the synthesis of guanosine 3',5'-bispyrophosphate (ppGpp) during the stringent response to amino acid starvation. This report presents the sequence of the relA gene region and assignment of its 743-codon open reading frame by the following criteria: 1) genetic complementation of ppGpp synthesis in a relaxed (relA1) mutant during the stringent response; 2) changes in 3-aminotriazole resistance during growth to mimic a relA+ phenotype; 3) verification of the presence of an amber codon at the normal carboxyl terminus of the relA gene; and 4) immunological assays of expression of the RelA protein. The apparent molecular mass of the cloned relA gene product is calculated to be 83,856 daltons and as visualized by immunoblotting is identical to that of the previously characterized protein. A promoter has been identified that directs relA gene transcription towards the pyrG gene, in a counterclockwise direction on the E. coli chromosome. Genomic Southern blot analyses verify that the relA regions cloned and subjected to nucleotide sequence analysis correspond to homologous regions on the E. coli chromosome.