Protein-protein and protein-salt interactions in aqueous protein solutions containing concentrated electrolytes

Protein-protein and protein-salt interactions have been obtained for ovalbumin in solutions of ammonium sulfate and for lysozyme in solutions of ammonium sulfate, sodium chloride, potassium isothiocyanate, and potassium chloride. The two-body interactions between ovalbumin molecules in concentrated ammonium-sulfate solutions can be described by the DLVO potentials plus a potential that accounts for the decrease in free volume available to the protein due to the presence of the salt ions. The interaction between ovalbumin and ammonium sulfate is unfavorable, reflecting the kosmotropic nature of sulfate anions. Lysozyme-lysozyme interactions cannot be described by the above potentials because anion binding to lysozyme alters these interactions. Lysozyme-isothiocyanate complexes are strongly attractive due to electrostatic interactions resulting from bridging by the isothiocyanate ion. Lysozyme-lysozyme interactions in sulfate solutions are more repulsive than expected, possibly resulting from a larger excluded volume of a lysozyme-sulfate bound complex or perhaps, hydration forces between the lysozyme-sulfate complexes.     

Potential of mean force for protein-protein interaction studies.

Calculating protein-protein interaction energies is crucial for understanding protein-protein associations. On the basis of the methodology of mean-field potential, we have developed an empirical approach to estimate binding free energy for protein-protein interactions. This knowledge-based approach has been used to derive distance-dependent free energies of protein complexes from a nonredundant training set in the Protein Data Bank (PDB), with a careful treatment of homology. We calculate atom pair potentials for 16 pair interactions, which can reflect the importance of hydrophobic interactions and specific hydrogen-bonding interactions. The derived potentials for hydrogen-bonding interactions show a valley of favorable interactions at a distance of approximately 3 A, corresponding to that of an established hydrogen bond. For the test set of 28 protein complexes, the calculated energies have a correlation coefficient of 0.75 compared with experimental binding free energies. The performance of the method in ranking the binding energies of different protein-protein complexes shows that the energy estimation can be applied to value binding free energies for protein-protein associations.     

The influence of protein-protein interactions on the organization of proteins within thylakoid membranes

The influence of attractive protein-protein interactions on the organization of photosynthetic proteins within the thylakoid membrane was investigated. Protein-protein interactions were simulated using Monte Carlo techniques and the influence of different interaction energies was examined. It was found that weak interactions led to protein clusters whereas strong interactions led to ramified chains. An optimum curve for the relationship between interaction energy and the number of contact sites emerged. With increasing particle densities the effect decreased. In a mixture of interacting and noninteracting particles the distance between the noninteracting particles was increased and there seemed to be much more free space around them. In thylakoids, this could lead to a more homogeneous distribution of the noninteracting but rate-limiting cytochrome bf complexes. Due to the increased free space between cytochrome bf, obstruction of binding sites--occurring unavoidably in a random distribution--may be drastically reduced. Furthermore, protein-protein interactions in thylakoids may lead to a decrease in plastoquinone diffusion.     

Protein-protein interactions as a target for drugs in proteomics

Protein-protein interactions play a central role in numerous processes in the cell and are one of the main fields of functional proteomics. This review highlights the methods of bioinformatics and functional proteomics of protein-protein interaction investigation. The structures and properties of contact surfaces, forces involved in protein-protein interactions, kinetic and thermodynamic parameters of these reactions were considered. The properties of protein contact surfaces depend on their functions. The contact surfaces of permanent complexes resemble domain contacts or the protein core and it is reasonable to consider such complex formation as a continuation of protein folding. Characteristics of contact surfaces of temporary protein complexes share some similarities with active sites of enzymes. The contact surfaces of the temporary protein complexes have unique structure and properties and they are more conservative in comparison with active site of enzymes. So they represent prospective targets for a new generation of drugs. During the last decade, numerous investigations were undertaken to find or design small molecules that block protein dimerization or protein(peptide)-receptor interaction, or, on the contrary, to induce protein dimerization.     

Accounting for conformational entropy in predicting binding free energies of protein-protein interactions

Protein-protein interactions are governed by the change in free energy upon binding, ΔG = ΔH - TΔS. These interactions are often marginally stable, so one must examine the balance between the change in enthalpy, ΔH, and the change in entropy, ΔS, when investigating known complexes, characterizing the effects of mutations, or designing optimized variants. To perform a large-scale study into the contribution of conformational entropy to binding free energy, we developed a technique called GOBLIN (Graphical mOdel for BiomoLecular INteractions) that performs physics-based free energy calculations for protein-protein complexes under both side-chain and backbone flexibility. Goblin uses a probabilistic graphical model that exploits conditional independencies in the Boltzmann distribution and employs variational inference techniques that approximate the free energy of binding in only a few minutes. We examined the role of conformational entropy on a benchmark set of more than 700 mutants in eight large, well-studied complexes. Our findings suggest that conformational entropy is important in protein-protein interactions--the root mean square error (RMSE) between calculated and experimentally measured ΔΔGs decreases by 12% when explicit entropic contributions were incorporated. GOBLIN models all atoms of the protein complex and detects changes to the binding entropy along the interface as well as positions distal to the binding interface. Our results also suggest that a variational approach to entropy calculations may be quantitatively more accurate than the knowledge-based approaches used by the well-known programs FOLDX and Rosetta--GOBLIN's RMSEs are 10 and 36% lower than these programs, respectively.     

Protein-protein binding is often associated with changes in protonation state

pK(a) values of ionizable residues have been calculated using the PROPKA method and structures of 75 protein-protein complexes and their corresponding free forms. These pK(a) values were used to compute changes in protonation state of individual residues, net changes in protonation state of the complex relative to the uncomplexed proteins, and the correction to a binding energy calculated assuming standard protonation states at pH 7. For each complex, two different structures for the uncomplexed form of the proteins were used: the X-ray structures determined for the proteins in the absence of the other protein and the individual protein structures taken from the structure of the complex (referred to as unbound and bound structures, respectively). In 28 and 77% of the cases considered here, protein-protein binding is accompanied by a complete (>95%) or significant (>50%) change in protonation state of at least one residue using unbound structures. Furthermore, in 36 and 61% of the cases, protein-protein binding is accompanied by a complete or significant net change in protonation state of the complex relative to the separated monomers. Using bound structures, the corresponding values are 12, 51, 20, and 48%. Comparison to experimental data suggest that using unbound and bound structures lead to over- and underestimation of binding-induced protonation state changes, respectively. Thus, we conclude that protein-protein binding is often associated with changes in protonation state of amino acid residues and with changes in the net protonation state of the proteins. The pH-dependent correction to the binding energy contributes at least one order of magnitude to the binding constant in 45 and 23%, using unbound and bound structures, respectively.     

Preferential interactions determine protein solubility in three-component solutions: the MgCl2 system

The correlation between protein solubility and the preferential interactions of proteins with solvent components was critically examined with aqueous MgCl2 as the solvent system. Preferential interaction and solubility measurements with three proteins, beta-lactoglobulin, bovine serum albumin, and lysozyme, resulted in similar patterns of interaction. At acid pH (pH 2-3) and lower salt concentrations (less than 2 M), the proteins were preferentially hydrated, while at higher salt concentrations, the interaction was either that of preferential salt binding or low salt exclusion. At pH 4.5-5, all three proteins exhibited either very low preferential hydration or preferential binding of MgCl2. These results were analyzed in terms of the balance between salt binding and salt exclusion attributed to the increase in the surface tension of water by salts, which is invariant with conditions. It was shown that the increase in salt binding at high salt concentration is a reflection of mass action, while its decrease at acid pH is due to the electrostatic repulsion between Mg2+ ions and the high net positive charge on the protein. The preferential interaction pattern was paralleled by the variation of protein solubility with solvent conditions. Calculation of the transfer free energies from water to the salt solutions for proteins in solution and in the precipitate showed dependencies on salt concentration. This indicates that the nature of interactions between proteins and solvent components is the same in solution and in the solid state, which implies no change in protein structure during precipitation. Analysis of the transfer free energies and preferential interaction parameter in terms of the salting-in, salting-out, and weak ion binding contributions has led to the conclusions that, when the weak ion binding contribution is small, the predominant protein-salt interaction must be that of preferential salt exclusion most probably caused by the increase of the surface tension of water by addition of the salt. A necessary consequence of this is salting-out of the protein, if the protein structure is to remain unaltered.     

Effective knowledge‐based potentials

Empirical or knowledge‐based potentials have many applications in structural biology such as the prediction of protein structure, protein–protein, and protein–ligand interactions and in the evaluation of stability for mutant proteins, the assessment of errors in experimentally solved structures, and the design of new proteins. Here, we describe a simple procedure to derive and use pairwise distance‐dependent potentials that rely on the definition of effective atomic interactions, which attempt to capture interactions that are more likely to be physically relevant. Based on a difficult benchmark test composed of proteins with different secondary structure composition and representing many different folds, we show that the use of effective atomic interactions significantly improves the performance of potentials at discriminating between native and near‐native conformations. We also found that, in agreement with previous reports, the potentials derived from the observed effective atomic interactions in native protein structures contain a larger amount of mutual information. A detailed analysis of the effective energy functions shows that atom connectivity effects, which mostly arise when deriving the potential by the incorporation of those indirect atomic interactions occurring beyond the first atomic shell, are clearly filtered out. The shape of the energy functions for direct atomic interactions representing hydrogen bonding and disulfide and salt bridges formation is almost unaffected when effective interactions are taken into account. On the contrary, the shape of the energy functions for indirect atom interactions (i.e., those describing the interaction between two atoms bound to a direct interacting pair) is clearly different when effective interactions are considered. Effective energy functions for indirect interacting atom pairs are not influenced by the shape or the energy minimum observed for the corresponding direct interacting atom pair. Our results suggest that the dependency between the signals in different energy functions is a key aspect that need to be addressed when empirical energy functions are derived and used, and also highlight the importance of additivity assumptions in the use of potential energy functions.     

Selected peptide extension contacts hydrophobic patch on neighboring zinc finger and mediates dimerization on DNA.

Protein-protein interactions often play a crucial role in stabilizing protein-DNA complexes and thus facilitate site-specific DNA recognition. We have worked to incorporate such protein-protein contacts into our design and selection strategies for short peptide extensions that promote cooperative binding of zinc finger proteins to DNA. We have determined the crystal structure of one of these fusion protein-DNA complexes. The selected peptide extension was found to mediate dimerization by reaching across the dyad axis and contacting a hydrophobic patch on the surface of the zinc finger bound to the adjacent DNA site. The peptide-zinc finger protein interactions observed in this structure are similar to those of some homeodomain heterodimers. We also find that the region of the zinc finger surface contacted by the selected peptide extension corresponds to surfaces that also make key interactions in the zinc finger proteins GLI and SWI5.     

Protein stabilization by salt bridges: concepts, experimental approaches and clarification of some misunderstandings

Salt bridges in proteins are bonds between oppositely charged residues that are sufficiently close to each other to experience electrostatic attraction. They contribute to protein structure and to the specificity of interaction of proteins with other biomolecules, but in doing so they need not necessarily increase a protein's free energy of unfolding. The net electrostatic free energy of a salt bridge can be partitioned into three components: charge-charge interactions, interactions of charges with permanent dipoles, and desolvation of charges. Energetically favorable Coulombic charge-charge interaction is opposed by often unfavorable desolvation of interacting charges. As a consequence, salt bridges may destabilize the structure of the folded protein. There are two ways to estimate the free energy contribution of salt bridges by experiment: the pK(a) approach and the mutation approach. In the pK(a) approach, the contribution of charges to the free energy of unfolding of a protein is obtained from the change of pK(a) of ionizable groups caused by altered electrostatic interactions upon folding of the protein. The pK(a) approach provides the relative free energy gained or lost when ionizable groups are being charged. In the mutation approach, the coupling free energy between interacting charges is obtained from a double mutant cycle. The coupling free energy is an indirect and approximate measure of the free energy of charge-charge interaction. Neither the pK(a) approach nor the mutation approach can provide the net free energy of a salt bridge. Currently, this is obtained only by computational methods which, however, are often prone to large uncertainties due to simplifying assumptions and insufficient structural information on which calculations are based. This state of affairs makes the precise thermodynamic quantification of salt bridge energies very difficult. This review is focused on concepts and on the assessment of experimental methods and does not cover the vast literature.     

Molecular principles of the interactions of disordered proteins

Thorough knowledge of the molecular principles of protein-protein recognition is essential to our understanding of protein function at the cellular level. Whereas interactions of ordered proteins have been analyzed in great detail, complexes of intrinsically unstructured/disordered proteins (IUPs) have hardly been addressed so far. Here, we have collected a database of 39 complexes of experimentally verified IUPs, and compared their interfaces with those of 72 complexes of ordered, globular proteins. The characteristic differences found between the two types of complexes suggest that IUPs represent a distinct molecular implementation of the principles of protein-protein recognition. The interfaces do not differ in size, but those of IUPs cover a much larger part of the surface of the protein than for their ordered counterparts. Moreover, IUP interfaces are significantly more hydrophobic relative to their overall amino acid composition, but also in absolute terms. They rely more on hydrophobic-hydrophobic than on polar-polar interactions. Their amino acids in the interface realize more intermolecular contacts, which suggests a better fit with the partner due to induced folding upon binding that results in a better adaptation to the partner. The two modes of interaction also differ in that IUPs usually use only a single continuous segment for partner binding, whereas the binding sites of ordered proteins are more segmented. Probably, all these features contribute to the increased evolutionary conservation of IUP interface residues. These noted molecular differences are also manifested in the interaction energies of IUPs. Our approximation of these by low-resolution force-fields shows that IUPs gain much more stabilization energy from intermolecular contacts, than from folding, i.e. they use their binding energy for folding. Overall, our findings provide a structural rationale to the prior suggestions that many IUPs are specialized for functions realized by protein-protein interactions.     

Protein-protein interactions: mechanisms and modification by drugs

Protein-protein interactions form the proteinaceous network, which plays a central role in numerous processes in the cell. This review highlights the main structures, properties of contact surfaces, and forces involved in protein-protein interactions. The properties of protein contact surfaces depend on their functions. The characteristics of contact surfaces of short-lived protein complexes share some similarities with the active sites of enzymes. The contact surfaces of permanent complexes resemble domain contacts or the protein core. It is reasonable to consider protein-protein complex formation as a continuation of protein folding. The contact surfaces of the protein complexes have unique structure and properties, so they represent prospective targets for a new generation of drugs. During the last decade, numerous investigations have been undertaken to find or design small molecules that block protein dimerization or protein(peptide)-receptor interaction, or on the other hand, induce protein dimerization.     

Hydrophobic forces between protein molecules in aqueous solutions of concentrated electrolyte

Protein-protein interactions have been measured for a mutant (D101F) lysozyme and for native lysozyme in concentrated solutions of ammonium sulfate at pH 7 and sodium chloride at pH 4.5. In the mutant lysozyme, a surface aspartate residue has been replaced with a hydrophobic phenylalanine residue. The protein-protein interactions of D101F lysozyme are more attractive than those of native lysozyme for all conditions studied. The salt-induced attraction is correlated with a solvation potential of mean force given by the work required to desolvate the part of the protein surfaces that is buried by the protein-protein interaction. This work is proportional to the aqueous surface-tension increment of the salt and the fractional non-polar surface coverage of the protein. Experimental measurements of osmotic second virial coefficients validate a proposed potential of mean force that ascribes the salt-induced attraction between protein molecules to an enhancement of the hydrophobic attraction. This model provides a first approximation for predicting the protein-protein potential of mean force in concentrated aqueous electrolyte solutions; this potential is useful for determining solution conditions favorable for protein crystallization.     

Shape complementarity of protein-protein complexes at multiple resolutions

Biological complexes typically exhibit intermolecular interfaces of high shape complementarity. Many computational docking approaches use this surface complementarity as a guide in the search for predicting the structures of protein-protein complexes. Proteins often undergo conformational changes to create a highly complementary interface when associating. These conformational changes are a major cause of failure for automated docking procedures when predicting binding modes between proteins using their unbound conformations. Low resolution surfaces in which high frequency geometric details are omitted have been used to address this problem. These smoothed, or blurred, surfaces are expected to minimize the differences between free and bound structures, especially those that are due to side chain conformations or small backbone deviations. Despite the fact that this approach has been used in many docking protocols, there has yet to be a systematic study of the effects of such surface smoothing on the shape complementarity of the resulting interfaces. Here we investigate this question by computing shape complementarity of a set of 66 protein-protein complexes represented by multiresolution blurred surfaces. Complexed and unbound structures are available for these protein-protein complexes. They are a subset of complexes from a nonredundant docking benchmark selected for rigidity (i.e. the proteins undergo limited conformational changes between their bound and unbound states). In this work, we construct the surfaces by isocontouring a density map obtained by accumulating the densities of Gaussian functions placed at all atom centers of the molecule. The smoothness or resolution is specified by a Gaussian fall-off coefficient, termed "blobbyness." Shape complementarity is quantified using a histogram of the shortest distances between two proteins' surface mesh vertices for both the crystallographic complexes and the complexes built using the protein structures in their unbound conformation. The histograms calculated for the bound complex structures demonstrate that medium resolution smoothing (blobbyness = -0.9) can reproduce about 88% of the shape complementarity of atomic resolution surfaces. Complexes formed from the free component structures show a partial loss of shape complementarity (more overlaps and gaps) with the atomic resolution surfaces. For surfaces smoothed to low resolution (blobbyness = -0.3), we find more consistency of shape complementarity between the complexed and free cases. To further reduce bad contacts without significantly impacting the good contacts we introduce another blurred surface, in which the Gaussian densities of flexible atoms are reduced. From these results we discuss the use of shape complementarity in protein-protein docking.     

Computational protein design with explicit consideration of surface hydrophobic patches

De novo protein design requires the identification of amino-acid sequences that favor the target-folded conformation and are soluble in water. One strategy for promoting solubility is to disallow hydrophobic residues on the protein surface during design. However, naturally occurring proteins often have hydrophobic amino acids on their surface that contribute to protein stability via the partial burial of hydrophobic surface area or play a key role in the formation of protein-protein interactions. A less restrictive approach for surface design that is used by the modeling program Rosetta is to parameterize the energy function so that the number of hydrophobic amino acids designed on the protein surface is similar to what is observed in naturally occurring monomeric proteins. Previous studies with Rosetta have shown that this limits surface hydrophobics to the naturally occurring frequency (～28%), but that it does not prevent the formation of hydrophobic patches that are considerably larger than those observed in naturally occurring proteins. Here, we describe a new score term that explicitly detects and penalizes the formation of hydrophobic patches during computational protein design. With the new term, we are able to design protein surfaces that include hydrophobic amino acids at naturally occurring frequencies, but do not have large hydrophobic patches. By adjusting the strength of the new score term, the emphasis of surface redesigns can be switched between maintaining solubility and maximizing folding free energy.     

Membrane-SPINE: an improved method to identify protein-protein interaction partners of membrane proteins in vivo

Membrane proteins are crucial for many essential cellular processes. As membrane proteins function in complexes, methods to detect and to characterize membrane protein-protein interactions are undoubtedly required. Therefore, we developed the "Membrane-Strep-tagged protein interaction experiment" (Membrane-SPINE) that combines the specific purification of a Strep-tagged membrane protein with the reversible fixation of protein complexes by formaldehyde cross-linking. In combination with MS analysis, we suggest Membrane-SPINE as a powerful tool to identify unknown interaction partners of membrane proteins in vivo.