Calorimetric Dissection of Colicin DNase−Immunity Protein Complex Specificity

Evolutionary paths that link orthogonal pairs of binding proteins

Some protein binding pairs exhibit extreme specificities that functionally insulate them from homologs. Such pairs evolve mostly by accumulating single-point mutations, and mutants are selected if their affinity exceeds the threshold required for function1-4. Thus, homologous and high-specificity binding pairs bring to light an evolutionary conundrum: how does a new specificity evolve while maintaining the required affinity in each intermediate5,6? Until now, a fully functional single-mutation path that connects two orthogonal pairs has only been described where the pairs were mutationally close thus enabling experimental enumeration of all intermediates2. We present an atomistic and graph-theoretical framework for discovering low molecular strain single-mutation paths that connect two extant pairs, enabling enumeration beyond experimental capability. We apply it to two orthogonal bacterial colicin endonuclease-immunity pairs separated by 17 interface mutations7. We were not able to find a strain-free and functional path in the sequence space defined by the two extant pairs. But including mutations that bridge amino acids that cannot be exchanged through single-nucleotide mutations led us to a strain-free 19-mutation trajectory that is completely viable in vivo. Our experiments show that the specificity switch is remarkably abrupt, resulting from only one radical mutation on each partner. Furthermore, each of the critical specificity-switch mutations increases fitness, demonstrating that functional divergence could be driven by positive Darwinian selection. These results reveal how even radical functional changes in an epistatic fitness landscape may evolve.

The Thermodynamic Fingerprints of Ultra-Tight Nanobody-Antigen Binding Probed via Two-Color Single-Molecule Coincidence Detection

Life on the molecular scale is based on a versatile interplay of biomolecules, a feature that is relevant for the formation of macromolecular complexes. Fluorescence-based two-color coincidence detection is widely used to characterize molecular binding and was recently improved by a brightness-gated version which gives more accurate results. We developed and established protocols which make use of coincidence detection to quantify binding fractions between interaction partners labeled with fluorescence dyes of different colors. Since the applied technique is intrinsically related to single-molecule detection, the concentration of diffusing molecules for confocal detection is typically in the low picomolar regime. This makes the approach a powerful tool for determining bi-molecular binding affinities, in terms of KD values, in this regime. We demonstrated the reliability of our approach by analyzing very strong nanobody-EGFP binding. By measuring the affinity at different temperatures, we were able to determine the thermodynamic parameters of the binding interaction. The results show that the ultra-tight binding is dominated by entropic contributions.

Affinity of disordered protein complexes is modulated by entropy-energy reinforcement

Intrinsically disordered proteins (IDPs), which are very common and essential to many biological activities, sometimes function via interaction with another IDP and form a fuzzy complex, which can be highly stable. It is unclear what the biophysical forces are that govern their thermodynamics and specificity, which are essential for de novo fuzzy complex design. Here, we explored the fuzzy complex formed between ProTα and H1, which are oppositely charged IDPs, by swapping the charges between them, generating variants that have either greater polyampholytic or polyelectrolytic nature as well as different charge patterns. Charge swapping and shuffling dramatically change the affinity of the fuzzy complex, which is contributed to by both enthalpy and entropy, where the latter is dominated by counterion release.

The association between two intrinsically disordered proteins (IDPs) may produce a fuzzy complex characterized by a high binding affinity, similar to that found in the ultrastable complexes formed between two well-structured proteins. Here, using coarse-grained simulations, we quantified the biophysical forces driving the formation of such fuzzy complexes. We found that the high-affinity complex formed between the highly and oppositely charged H1 and ProTα proteins is sensitive to electrostatic interactions. We investigated 52 variants of the complex by swapping charges between the two oppositely charged proteins to produce sequences whose negatively or positively charged residue content was more homogeneous or heterogenous (i.e., polyelectrolytic or polyampholytic, having higher or lower absolute net charges, respectively) than the wild type. We also changed the distributions of oppositely charged residues within each participating sequence to produce variants in which the charges were segregated or well mixed. Both types of changes significantly affect binding affinity in fuzzy complexes, which is governed by both enthalpy and entropy. The formation of H1–ProTa is supported by an increase in configurational entropy and by entropy due to counterion release. The latter can be twice as large as the former, illustrating the dominance of counterion entropy in modulating the binding thermodynamics. Complexes formed between proteins with greater absolute net charges are more stable, both enthalpically and entropically, indicating that enthalpy and entropy have a mutually reinforcing effect. The sensitivity of the thermodynamics of the complex to net charge and the charge pattern within each of the binding constituents may provide a means to achieve binding specificity between IDPs.

Unraveling the Thermodynamics of Ultra-tight Binding of Intrinsically Disordered Proteins

Protein interactions mediated by the intrinsically disordered proteins (IDPs) are generally associated with lower affinities compared to those between globular proteins. Here, we characterize the association between the intrinsically disordered HigA2 antitoxin and its globular target HigB2 toxin from Vibrio cholerae using competition ITC experiments. We demonstrate that this interaction reaches one of the highest affinities reported for IDP-target systems (K_D = 3 pM) and can be entirely attributed to a short, 20-residue-long interaction motif that folds into α-helix upon binding. We perform an experimentally based decomposition of the IDP-target association parameters into folding and binding contributions, which allows a direct comparison of the binding contribution with those from globular ultra-high affinity binders. We find that the HigA2-HigB2 interface is energy optimized to a similar extent as the interfaces of globular ultra-high affinity complexes, such as barnase-barstar. Evaluation of other ultra-high affinity IDP-target systems shows that a strategy based on entropy optimization can also achieve comparably high, picomolar affinities. Taken together, these examples show how IDP-target interactions achieve picomolar affinities either through enthalpy optimization (HigA2-HigB2), resembling the ultra-high affinity binding of globular proteins, or via bound-state fuzziness and entropy optimization (CcdA-CcdB, histone H1-prothymosin α).

Structural design principles for specific ultra-high affinity interactions between colicins/pyocins and immunity proteins

The interactions of the antibiotic proteins colicins/pyocins with immunity proteins is a seminal model system for studying protein–protein interactions and specificity. Yet, a precise and quantitative determination of which structural elements and residues determine their binding affinity and specificity is still lacking. Here, we used comparative structure-based energy calculations to map residues that substantially contribute to interactions across native and engineered complexes of colicins/pyocins and immunity proteins. We show that the immunity protein α1–α2 motif is a unique structurally-dissimilar element that restricts interaction specificity towards all colicins/pyocins, in both engineered and native complexes. This motif combines with a diverse and extensive array of electrostatic/polar interactions that enable the exquisite specificity that characterizes these interactions while achieving ultra-high affinity. Surprisingly, the divergence of these contributing colicin residues is reciprocal to residue conservation in immunity proteins. The structurally-dissimilar immunity protein α1–α2 motif is recognized by divergent colicins similarly, while the conserved immunity protein α3 helix interacts with diverse colicin residues. Electrostatics thus plays a key role in setting interaction specificity across all colicins and immunity proteins. Our analysis and resulting residue-level maps illuminate the molecular basis for these protein–protein interactions, with implications for drug development and rational engineering of these interfaces.

Development of a Force Field for the Simulation of Single-Chain Proteins and Protein-Protein Complexes

The accuracy of atomistic physics-based force fields for the simulation of biological macromolecules has typically been benchmarked experimentally using biophysical data from simple, often single-chain systems. In the case of proteins, the careful refinement of force field parameters associated with torsion-angle potentials and the use of improved water models have enabled a great deal of progress toward the highly accurate simulation of such monomeric systems in both folded and, more recently, disordered states. In living organisms, however, proteins constantly interact with other macromolecules, such as proteins and nucleic acids, and these interactions are often essential for proper biological function. Here, we show that state-of-the-art force fields tuned to provide an accurate description of both ordered and disordered proteins can be limited in their ability to accurately describe protein-protein complexes. This observation prompted us to perform an extensive reparameterization of one variant of the Amber protein force field. Our objective involved refitting not only the parameters associated with torsion-angle potentials but also the parameters used to model nonbonded interactions, the specification of which is expected to be central to the accurate description of multicomponent systems. The resulting force field, which we call DES-Amber, allows for more accurate simulations of protein-protein complexes, while still providing a state-of-the-art description of both ordered and disordered single-chain proteins. Despite the improvements, calculated protein-protein association free energies still appear to deviate substantially from experiment, a result suggesting that more fundamental changes to the force field, such as the explicit treatment of polarization effects, may simultaneously further improve the modeling of single-chain proteins and protein-protein complexes.

Computational studies of protein-protein dissociation by statistical potential and coarse-grained simulations: a case study on interactions between colicin E9 endonuclease and immunity proteins

Proteins carry out their diverse functions in cells by forming interactions with each other. The dynamics of these interactions are quantified by the measurement of association and dissociation rate constants. Relative to the efforts made to model the association of biomolecules, little has been studied to understand the principles of protein complex dissociation. Using the interaction between colicin E9 endonucleases and immunity proteins as a test system, here we develop a coarse-grained simulation method to explore the dissociation mechanisms of protein complexes. The interactions between proteins in the complex are described by the knowledge-based potential that was constructed by the statistics from available protein complexes in the structural database. Our study provides the supportive evidences to the dual recognition mechanism for the specificity of binding between E9 DNase and immunity proteins, in which the conserved residues of helix III of Im2 and Im9 proteins act as the anchor for binding, while the sequence variations in helix II make positive or negative contributions to specificity. Beyond that, we further suggest that this binding specificity is rooted in the process of complex dissociation instead of association. While we increased the flexibility of protein complexes, we further found that they are less prone to dissociation, suggesting that conformational fluctuations of protein complexes play important functional roles in regulating their binding and dissociation. Our studies therefore bring new insights to the molecule mechanisms of protein-protein interactions, while the method can serve as a new addition to a suite of existing computational tools for the simulations of protein complexes.

Ultrahigh specificity in a network of computationally designed protein-interaction pairs

Protein networks in all organisms comprise homologous interacting pairs. In these networks, some proteins are specific, interacting with one or a few binding partners, whereas others are multispecific and bind a range of targets. We describe an algorithm that starts from an interacting pair and designs dozens of new pairs with diverse backbone conformations at the binding site as well as new binding orientations and sequences. Applied to a high-affinity bacterial pair, the algorithm results in 18 new ones, with cognate affinities from pico- to micromolar. Three pairs exhibit 3-5 orders of magnitude switch in specificity relative to the wild type, whereas others are multispecific, collectively forming a protein-interaction network. Crystallographic analysis confirms design accuracy, including in new backbones and polar interactions. Preorganized polar interaction networks are responsible for high specificity, thus defining design principles that can be applied to program synthetic cellular interaction networks of desired affinity and specificity.

The molecular basis of ultrahigh specificity in protein-protein interactions remains obscure. The authors present a computational method to design atomically accurate new pairs exhibiting >100,000-fold specificity switches, generating a large and complex interaction network.

Structures of the Ultra-High-Affinity Protein-Protein Complexes of Pyocins S2 and AP41 and Their Cognate Immunity Proteins from Pseudomonas aeruginosa

How ultra-high-affinity protein-protein interactions retain high specificity is still poorly understood. The interaction between colicin DNase domains and their inhibitory immunity (Im) proteins is an ultra-high-affinity interaction that is essential for the neutralisation of endogenous DNase catalytic activity and for protection against exogenous DNase bacteriocins. The colicin DNase-Im interaction is a model system for the study of high-affinity protein-protein interactions. However, despite the fact that closely related colicin-like bacteriocins are widely produced by Gram-negative bacteria, this interaction has only been studied using colicins from Escherichia coli. In this work, we present the first crystal structures of two pyocin DNase-Im complexes from Pseudomonas aeruginosa, pyocin S2 DNase-ImS2 and pyocin AP41 DNase-ImAP41. These structures represent divergent DNase-Im subfamilies and are important in extending our understanding of protein-protein interactions for this important class of high-affinity protein complex. A key finding of this work is that mutations within the immunity protein binding energy hotspot, helix III, are tolerated by complementary substitutions at the DNase-Immunity protein binding interface. Im helix III is strictly conserved in colicins where an Asp forms polar interactions with the DNase backbone. ImAP41 contains an Asp-to-Gly substitution in helix III and our structures show the role of a co-evolved substitution where Pro in DNase loop 4 occupies the volume vacated and removes the unfulfilled hydrogen bond. We observe the co-evolved mutations in other DNase-Immunity pairs that appear to underpin the split of this family into two distinct groups.

Globular and disordered-the non-identical twins in protein-protein interactions

In biology proteins from different structural classes interact across and within classes in ways that are optimized to achieve balanced functional outputs. The interactions between intrinsically disordered proteins (IDPs) and other proteins rely on changes in flexibility and this is seen as a strong determinant for their function. This has fostered the notion that IDP's bind with low affinity but high specificity. Here we have analyzed available detailed thermodynamic data for protein-protein interactions to put to the test if the thermodynamic profiles of IDP interactions differ from those of other protein-protein interactions. We find that ordered proteins and the disordered ones act as non-identical twins operating by similar principles but where the disordered proteins complexes are on average less stable by 2.5 kcal mol(-1).

A "fuzzy"-logic language for encoding multiple physical traits in biomolecules

To carry out their activities, biological macromolecules balance different physical traits, such as stability, interaction affinity, and selectivity. How such often opposing traits are encoded in a macromolecular system is critical to our understanding of evolutionary processes and ability to design new molecules with desired functions. We present a framework for constraining design simulations to balance different physical characteristics. Each trait is represented by the equilibrium fractional occupancy of the desired state relative to its alternatives, ranging from none to full occupancy, and the different traits are combined using Boolean operators to effect a "fuzzy"-logic language for encoding any combination of traits. In another paper, we presented a new combinatorial backbone design algorithm AbDesign where the fuzzy-logic framework was used to optimize protein backbones and sequences for both stability and binding affinity in antibody-design simulation. We now extend this framework and find that fuzzy-logic design simulations reproduce sequence and structure design principles seen in nature to underlie exquisite specificity on the one hand and multispecificity on the other hand. The fuzzy-logic language is broadly applicable and could help define the space of tolerated and beneficial mutations in natural biomolecular systems and design artificial molecules that encode complex characteristics.

A force-activated trip switch triggers rapid dissociation of a colicin from its immunity protein

A single-molecule force study shows that rapid dissociation of a high-affinity protein interaction can be triggered by site-specific remodelling of one protein partner, and that prevention of remodelling maintains avidity.

Colicins are protein antibiotics synthesised by Escherichia coli strains to target and kill related bacteria. To prevent host suicide, colicins are inactivated by binding to immunity proteins. Despite their high avidity (K_d≈fM, lifetime ≈4 days), immunity protein release is a pre-requisite of colicin intoxication, which occurs on a timescale of minutes. Here, by measuring the dynamic force spectrum of the dissociation of the DNase domain of colicin E9 (E9) and immunity protein 9 (Im9) complex using an atomic force microscope we show that application of low forces (<20 pN) increases the rate of complex dissociation 10⁶-fold, to a timescale (lifetime ≈10 ms) compatible with intoxication. We term this catastrophic force-triggered increase in off-rate a trip bond. Using mutational analysis, we elucidate the mechanism of this switch in affinity. We show that the N-terminal region of E9, which has sparse contacts with the hydrophobic core, is linked to an allosteric activator region in E9 (residues 21–30) whose remodelling triggers immunity protein release. Diversion of the force transduction pathway by the introduction of appropriately positioned disulfide bridges yields a force resistant complex with a lifetime identical to that measured by ensemble techniques. A trip switch within E9 is ideal for its function as it allows bipartite complex affinity, whereby the stable colicin:immunity protein complex required for host protection can be readily converted to a kinetically unstable complex whose dissociation is necessary for cellular invasion and competitor death. More generally, the observation of two force phenotypes for the E9:Im9 complex demonstrates that force can re-sculpt the underlying energy landscape, providing new opportunities to modulate biological reactions in vivo; this rationalises the commonly observed discrepancy between off-rates measured by dynamic force spectroscopy and ensemble methods.

Many proteins interact with other proteins as part of their function. One method of modulating the activity of protein complexes is to break them apart. Some complexes, however, are extremely kinetically stable and it is unclear how these can dissociate on a biologically relevant timescale. In this study we address this question using protein complexes between colicin E9 (a bacterial toxin) and its immunity protein Im9. These highly avid complexes (with a lifetime of days) must be broken apart for colicin to be activated. By using single-molecule force methods we show that pulling on one end of colicin E9 drastically destabilises the complex so that it dissociates a million-fold faster than its intrinsic rate. We then show that preventing this destabilisation (by the insertion of cross-links that pin the N-terminus of E9 in place) yields a kinetically stable complex. It has previously been postulated that force can destabilise a protein complex by partially unfolding one or more binding partners. Our work provides new experimental evidence that shows this is the case and provides a mechanism for this phenomenon, which we term a trip bond. For the E9:Im9 complex, trip bond behaviour allows a stable complex to be rapidly dissociated by application of a surprisingly small force.

Coarse-grained simulations of protein-protein association: an energy landscape perspective

Understanding protein-protein association is crucial in revealing the molecular basis of many biological processes. Here, we describe a theoretical simulation pipeline to study protein-protein association from an energy landscape perspective. First, a coarse-grained model is implemented and its applications are demonstrated via molecular dynamics simulations for several protein complexes. Second, an enhanced search method is used to efficiently sample a broad range of protein conformations. Third, multiple conformations are identified and clustered from simulation data and further projected on a three-dimensional globe specifying protein orientations and interacting energies. Results from several complexes indicate that the crystal-like conformation is favorable on the energy landscape even if the landscape is relatively rugged with metastable conformations. A closer examination on molecular forces shows that the formation of associated protein complexes can be primarily electrostatics-driven, hydrophobics-driven, or a combination of both in stabilizing specific binding interfaces. Taken together, these results suggest that the coarse-grained simulations and analyses provide an alternative toolset to study protein-protein association occurring in functional biomolecular complexes.

Conformational properties of the unfolded state of Im7 in nondenaturing conditions

The unfolded ensemble in aqueous solution represents the starting point of protein folding. Characterisation of this species is often difficult since the native state is usually predominantly populated at equilibrium. Previous work has shown that the four-helix protein, Im7 (immunity protein 7), folds via an on-pathway intermediate. While the transition states and folding intermediate have been characterised in atomistic detail, knowledge of the unfolded ensemble under the same ambient conditions remained sparse. Here, we introduce destabilising amino acid substitutions into the sequence of Im7, such that the unfolded state becomes predominantly populated at equilibrium in the absence of denaturant. Using far- and near-UV CD, fluorescence, urea titration and heteronuclear NMR experiments, we show that three amino acid substitutions (L18A-L19A-L37A) are sufficient to prevent Im7 folding, such that the unfolded state is predominantly populated at equilibrium. Using measurement of chemical shifts, (15)N transverse relaxation rates and sedimentation coefficients, we show that the unfolded species of L18A-L19A-L37A deviates significantly from random-coil behaviour. Specifically, we demonstrate that this unfolded species is compact (R(h)=25 Å) relative to the urea-denatured state (R(h)≥30 Å) and contains local clusters of hydrophobic residues in regions that correspond to the four helices in the native state. Despite these interactions, there is no evidence for long-range stabilising tertiary interactions or persistent helical structure. The results reveal an unfolded ensemble that is conformationally restricted in regions of the polypeptide chain that ultimately form helices I, II and IV in the native state.

Structure of the ultra-high-affinity colicin E2 DNase--Im2 complex

How proteins achieve high-affinity binding to a specific protein partner while simultaneously excluding all others is a major biological problem that has important implications for protein design. We report the crystal structure of the ultra-high-affinity protein-protein complex between the endonuclease domain of colicin E2 and its cognate immunity (Im) protein, Im2 (K(d)∼10(-)(15) M), which, by comparison to previous structural and biophysical data, provides unprecedented insight into how high affinity and selectivity are achieved in this model family of protein complexes. Our study pinpoints the role of structured water molecules in conjoining hotspot residues that govern stability with residues that control selectivity. A key finding is that a single residue, which in a noncognate context massively destabilizes the complex through frustration, does not participate in specificity directly but rather acts as an organizing center for a multitude of specificity interactions across the interface, many of which are water mediated.

Nuclease colicins and their immunity proteins

It is more than 80 years since Gratia first described 'a remarkable antagonism between two strains of Escherichia coli'. Shown subsequently to be due to the action of proteins (or peptides) produced by one bacterium to kill closely related species with which it might be cohabiting, such bacteriocins have since been shown to be commonplace in the internecine warfare between bacteria. Bacteriocins have been studied primarily from the twin perspectives of how they shape microbial communities and how they penetrate bacteria to kill them. Here, we review the modes of action of a family of bacteriocins that cleave nucleic acid substrates in E. coli, known collectively as nuclease colicins, and the specific immunity (inhibitor) proteins that colicin-producing organisms make in order to avoid committing suicide. In a process akin to targeting in mitochondria, nuclease colicins engage in a variety of cellular associations in order to translocate their cytotoxic domains through the cell envelope to the cytoplasm. As well as informing on the process itself, the study of nuclease colicin import has also illuminated functional aspects of the host proteins they parasitize. We also review recent studies where nuclease colicins and their immunity proteins have been used as model systems for addressing fundamental problems in protein folding and protein-protein interactions, areas of biophysics that are intimately linked to the role of colicins in bacterial competition and to the import process itself.

A structure-based benchmark for protein-protein binding affinity

We have assembled a nonredundant set of 144 protein-protein complexes that have high-resolution structures available for both the complexes and their unbound components, and for which dissociation constants have been measured by biophysical methods. The set is diverse in terms of the biological functions it represents, with complexes that involve G-proteins and receptor extracellular domains, as well as antigen/antibody, enzyme/inhibitor, and enzyme/substrate complexes. It is also diverse in terms of the partners' affinity for each other, with K(d) ranging between 10(-5) and 10(-14) M. Nine pairs of entries represent closely related complexes that have a similar structure, but a very different affinity, each pair comprising a cognate and a noncognate assembly. The unbound structures of the component proteins being available, conformation changes can be assessed. They are significant in most of the complexes, and large movements or disorder-to-order transitions are frequently observed. The set may be used to benchmark biophysical models aiming to relate affinity to structure in protein-protein interactions, taking into account the reactants and the conformation changes that accompany the association reaction, instead of just the final product.

CAPRI targets T29-T42: proving ground for new docking procedures

The critical assessment of protein interactions (CAPRI) experiment provides a unique opportunity for unbiased assessment of docking procedures. The recent CAPRI targets T29-T42 entailed docking of bound, unbound, and modeled structures, presenting a wide range of prediction difficulty. We submitted accurate predictions for targets T40, T41, and T42, a good prediction for T32 and acceptable predictions for T29 and T34. The accuracy of our docking results generally matched the prediction difficulty; hence, docking of modeled proteins produced less accurate results. However, there were interesting exceptions: an accurate prediction was submitted for the dimer of modeled tetratricopeptide repeat (T42) and only an acceptable prediction for the bound/unbound case T29. The ensembles of docking models produced in the scans included an acceptable or better prediction for every target. We show here that our recently developed postscan reevaluation procedure, which tests propensity and solvation measures of the whole interface and the interface core, successfully distinguished these predictions from false docking models. For enzyme-inhibitor targets, we show that the distance of the interface from the enzyme's centroid ranked high native like docking models. Also, for one case we demonstrate that docking of an ensemble of conformers produced by normal modes analysis can improve the accuracy of the prediction.

The structural and energetic basis for high selectivity in a high-affinity protein-protein interaction

High-affinity, high-selectivity protein-protein interactions that are critical for cell survival present an evolutionary paradox: How does selectivity evolve when acquired mutations risk a lethal loss of high-affinity binding? A detailed understanding of selectivity in such complexes requires structural information on weak, noncognate complexes which can be difficult to obtain due to their transient and dynamic nature. Using NMR-based docking as a guide, we deployed a disulfide-trapping strategy on a noncognate complex between the colicin E9 endonuclease (E9 DNase) and immunity protein 2 (Im2), which is seven orders of magnitude weaker binding than the cognate femtomolar E9 DNase-Im9 interaction. The 1.77 A crystal structure of the E9 DNase-Im2 complex reveals an entirely noncovalent interface where the intersubunit disulfide merely supports the crystal lattice. In combination with computational alanine scanning of interfacial residues, the structure reveals that the driving force for binding is so strong that a severely unfavorable specificity contact is tolerated at the interface and as a result the complex becomes weakened through "frustration." As well as rationalizing past mutational and thermodynamic data, comparing our noncognate structure with previous cognate complexes highlights the importance of loop regions in developing selectivity and accentuates the multiple roles of buried water molecules that stabilize, ameliorate, or aggravate interfacial contacts. The study provides direct support for dual-recognition in colicin DNase-Im protein complexes and shows that weakened noncognate complexes are primed for high-affinity binding, which can be achieved by economical mutation of a limited number of residues at the interface.

Energetic determinants of protein binding specificity: insights into protein interaction networks

One of the challenges of the postgenomic era is to provide a more realistic representation of cellular processes by combining a systems biology description of functional networks with information on their interacting components. Here we carried out a systematic large-scale computational study on a structural protein-protein interaction network dataset in order to dissect thermodynamic characteristics of binding determining the interplay between protein affinity and specificity. As expected, interactions involving specific binding sites display higher affinities than those of promiscuous binding sites. Next, in order to investigate a possible role of modular distribution of hot spots in binding specificity, we divided binding sites into modules previously shown to be energetically independent. In general, hot spots that interact with different partners are located in different modules. We further observed that common hot spots tend to interact with partners exhibiting common binding motifs, whereas different hot spots tend to interact with partners with different motifs. Thus, energetic properties of binding sites provide insights into the way proteins modulate interactions with different partners. Knowledge of those factors playing a role in protein specificity is important for understanding how proteins acquire additional partners during evolution. It should also be useful in drug design.

E9-Im9 colicin DNase-immunity protein biomolecular association in water: a multiple-copy and accelerated molecular dynamics simulation study

Protein-protein transient and dynamic interactions underlie all biological processes. The molecular dynamics (MD) of the E9 colicin DNase protein, its Im9 inhibitor protein, and their E9-Im9 recognition complex are investigated by combining multiple-copy (MC) MD and accelerated MD (aMD) explicit-solvent simulation approaches, after validation with crystalline-phase and solution experiments. Im9 shows higher flexibility than its E9 counterpart. Im9 displays a significant reduction of backbone flexibility and a remarkable increase in motional correlation upon E9 association. Im9 loops 23-31 and 54-64 open with respect to the E9-Im9 X-ray structure and show high conformational diversity. Upon association a large fraction (approximately 20 nm2) of E9 and Im9 protein surfaces become inaccessible to water. Numerous salt bridges transiently occurring throughout our six 50 ns long MC-MD simulations are not present in the X-ray model. Among these Im9 Glu31-E9 Arg96 and Im9 Glu41-Lys89 involve interface interactions. Through the use of 10 ns of Im9 aMD simulation, we reconcile the largest thermodynamic impact measured for Asp51Ala mutation with Im9 structure and dynamics. Lys57 acts as an essential molecular switch to shift Im9 surface loop towards an ideal configuration for E9 inhibition. This is achieved by switching Asp60-Lys57 and Asp62-Lys57 hydrogen bonds to Asp51-Lys57 salt bridge. E9-Im9 recognition involves shifts of conformational distributions, reorganization of intramolecular hydrogen bond patterns, and formation of new inter- and intramolecular interactions. The description of key transient biological interactions can be significantly enriched by the dynamic and atomic-level information provided by computer simulations.

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

Protein-protein association involves many interface interactions, but they do not contribute equally. Ala scanning experiments reveal that only a few mutations significantly lower binding affinity. These key residues, which appear to drive protein-protein association, are called hot-spot residues. Molecular dynamics simulations of the Colicin E9/Im9 complex show Im9 Glu41 and Im9 Ser50, both hot-spots, bind via different mechanisms. The results suggest that Im9 Ser50 restricts Glu41 in a conformation auspicious for salt-bridge formation across the interface. This type of model may be helpful in engineering hot-spot clusters at protein-protein interfaces and, consequently, the design of specificity.

Experimental and computational analyses of the energetic basis for dual recognition of immunity proteins by colicin endonucleases

Colicin endonucleases (DNases) are bound and inactivated by immunity (Im) proteins. Im proteins are broadly cross-reactive yet specific inhibitors binding cognate and non-cognate DNases with K(d) values that vary between 10(-4) and 10(-14) M, characteristics that are explained by a 'dual-recognition' mechanism. In this work, we addressed for the first time the energetics of Im protein recognition by colicin DNases through a combination of E9 DNase alanine scanning and double-mutant cycles (DMCs) coupled with kinetic and calorimetric analyses of cognate Im9 and non-cognate Im2 binding, as well as computational analysis of alanine scanning and DMC data. We show that differential DeltaDeltaGs observed for four E9 DNase residues cumulatively distinguish cognate Im9 association from non-cognate Im2 association. E9 DNase Phe86 is the primary specificity hotspot residue in the centre of the interface, which is coordinated by conserved and variable hotspot residues of the cognate Im protein. Experimental DMC analysis reveals that only modest coupling energies to Im9 residues are observed, in agreement with calculated DMCs using the program ROSETTA and consistent with the largely hydrophobic nature of E9 DNase-Im9 specificity contacts. Computed values for the 12 E9 DNase alanine mutants showed reasonable agreement with experimental DeltaDeltaG data, particularly for interactions not mediated by interfacial water molecules. DeltaDeltaG predictions for residues that contact buried water molecules calculated using solvated rotamer models met with mixed success; however, we were able to predict with a high degree of accuracy the location and energetic contribution of one such contact. Our study highlights how colicin DNases are able to utilise both conserved and variable amino acids to distinguish cognate from non-cognate Im proteins, with the energetic contributions of the conserved residues modulated by neighbouring specificity sites.

Spatial chemical conservation of hot spot interactions in protein-protein complexes

Background

Conservation of the spatial binding organizations at the level of physico-chemical interactions is important for the formation and stability of protein-protein complexes as well as protein and drug design. Due to the lack of computational tools for recognition of spatial patterns of interactions shared by a set of protein-protein complexes, the conservation of such interactions has not been addressed previously.

Results

We performed extensive spatial comparisons of physico-chemical interactions common to different types of protein-protein complexes. We observed that 80% of these interactions correspond to known hot spots. Moreover, we show that spatially conserved interactions allow prediction of hot spots with a success rate higher than obtained by methods based on sequence or backbone similarity. Detection of spatially conserved interaction patterns was performed by our novel MAPPIS algorithm. MAPPIS performs multiple alignments of the physico-chemical interactions and the binding properties in three dimensional space. It is independent of the overall similarity in the protein sequences, folds or amino acid identities. We present examples of interactions shared between complexes of colicins with immunity proteins, serine proteases with inhibitors and T-cell receptors with superantigens. We unravel previously overlooked similarities, such as the interactions shared by the structurally different RNase-inhibitor families.

Conclusion

The key contribution of MAPPIS is in discovering the 3D patterns of physico-chemical interactions. The detected patterns describe the conserved binding organizations that involve energetically important hot spot residues and are crucial for the protein-protein associations.

Energetics of protein homodimerization: Effects of water sequestering on the formation of β-lactoglobulin dimer

Transient protein-protein interactions are functionally relevant as a control mechanism in a variety of biological processes. Analysis of the 3D structure of protein-protein complexes indicates that water molecules trapped at the interface are very common; however, their role in the stability and specificity of protein homodimer interactions has been not addressed yet. To provide new insights into the energetic bases that govern the formation of highly hydrated interfaces, the dissociation process of bovine beta lg variant A at a neutral pH was characterized here thermodynamically by conducting dilution experiments with an isothermal titration calorimeter. Association was enthalpically driven throughout the temperature range spanned. DeltaH and deltaC(p) were significantly more negative than estimates based on surface area changes, suggesting the occurrence of effects additional to the dehydration of the contact surfaces between subunits. Near-UV CD spectra proved to be independent of protein concentration, indicating a rigid body-like association. Furthermore, the process proved not to be coupled to significant changes in the protonation state of ionizable groups or counterion exchange. In contrast, both osmotic stress experiments and a computational analysis of the dimer's 3D structure indicated that a large number of water molecules are incorporated into the interface upon association. Numerical estimates considering the contributions of interface area desolvation and water immobilization accounted satisfactorily for the experimental deltaC(p). Thus, our study highlights the importance of explicitly considering the effects of water sequestering to perform a proper quantitative analysis of the formation of homodimers with highly hydrated interfaces.

Biomimetic Interactions of Proteins with Functionalized Nanoparticles: A Thermodynamic Study

Gold nanoparticles (NPs) functionalized with L-amino acid-terminated monolayers provide an effective platform for the recognition of protein surfaces. Isothermal titration calorimetry (ITC) was used to quantify the binding thermodynamics of these functional NPs with alpha-chymotrypsin (ChT), histone, and cytochrome c (CytC). The enthalpy and entropy changes for the complex formation depend upon the nanoparticle structure and the surface characteristics of the proteins, e.g., distributions of charged and hydrophobic residues on the surface. Enthalpy-entropy compensation studies on these NP-protein systems indicate an excellent linear correlation between DeltaH and TDeltaS with a slope (alpha) of 1.07 and an intercept (TDeltaS0) of 35.2 kJ mol(-1). This behavior is closer to those of native protein-protein systems (alpha = 0.92 and TDeltaS0 = 41.1 kJ mol(-1)) than other protein-ligand and synthetic host-guest systems.

Mutations designed to destabilize the receptor-bound conformation increase MICA-NKG2D association rate and affinity

MICA is a major histocompatibility complex-like protein that undergoes a structural transition from disorder to order upon binding its immunoreceptor, NKG2D. We redesigned the disordered region of MICA with RosettaDesign to increase NKG2D binding. Mutations that stabilize this region were expected to increase association kinetics without changing dissociation kinetics, increase affinity of interaction, and reduce entropy loss upon binding. MICA mutants were stable in solution, and they were amenable to surface plasmon resonance evaluation of NKG2D binding kinetics and thermodynamics. Several MICA mutants bound NKG2D with enhanced affinity, kinetic changes were primarily observed during association, and thermodynamic changes in entropy were as expected. However, none of the 15 combinations of mutations predicted to stabilize the receptor-bound MICA conformation enhanced NKG2D affinity, whereas all 10 mutants predicted to be destabilized bound NKG2D with increased on-rates. Five of these had affinities enhanced by 0.9-1.8 kcal/mol over wild type by one to three non-contacting substitutions. Therefore, in this case, mutations designed to mildly destabilize a protein enhanced association and affinity.

Colicin biology

Colicins are proteins produced by and toxic for some strains of Escherichia coli. They are produced by strains of E. coli carrying a colicinogenic plasmid that bears the genetic determinants for colicin synthesis, immunity, and release. Insights gained into each fundamental aspect of their biology are presented: their synthesis, which is under SOS regulation; their release into the extracellular medium, which involves the colicin lysis protein; and their uptake mechanisms and modes of action. Colicins are organized into three domains, each one involved in a different step of the process of killing sensitive bacteria. The structures of some colicins are known at the atomic level and are discussed. Colicins exert their lethal action by first binding to specific receptors, which are outer membrane proteins used for the entry of specific nutrients. They are then translocated through the outer membrane and transit through the periplasm by either the Tol or the TonB system. The components of each system are known, and their implication in the functioning of the system is described. Colicins then reach their lethal target and act either by forming a voltage-dependent channel into the inner membrane or by using their endonuclease activity on DNA, rRNA, or tRNA. The mechanisms of inhibition by specific and cognate immunity proteins are presented. Finally, the use of colicins as laboratory or biotechnological tools and their mode of evolution are discussed.

The molecular architecture of protein-protein binding sites

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?

CH···π Interaction for Rhenium-Based Rectangles: An Interaction That Is Rarely Designed into a Host−Guest Pair

Alkoxy- and thiolato-bridged Re(I) molecular rectangles [{(CO)3Re(mu-ER)2Re(CO)3}2(mu-bpy)2] (ER = SC4H9, 1a; SC8H17, 1b; OC4H9, 2a; OC12H25, 2b; bpy = 4,4'-bipyridine) exhibit strong interactions with several planar aromatic molecules. The nature of their binding was studied by spectral techniques and verified by X-ray diffraction analysis. Standard absorption and fluorescence titrations showed that a relatively strong 1:1 interaction occurs between aromatic guests such as pyrene and these rectangles. The results of a single-crystal X-ray diffraction analysis show that the recognition of 1 with a pyrene molecule is mainly due to CH...pi interactions and the face of the guest pyrene is located over the edges of the bpy linkers of 1. This is a fairly novel example of an interaction that is rarely designed into a host-guest pair. Furthermore, the interaction of 1 with Ag+ results in the self-organization of supramolecular arrays, as revealed by solid-state data.