Lateral interactions among membrane proteins. Valid estimates based on freeze-fracture electron microscopy.

We consider the lateral distribution of intrinsic membrane proteins from the viewpoint of the statistical-mechanical theory of liquids. We connect the information in freeze-fracture electron micrographs--positions of proteins but not lipids or aqueous species--to a well developed theory of liquid mixtures. An algorithm, based on the Born-Green-Yvon integral equation, is presented for deducing forces between proteins from correlations among protein positions that are observed in micrographs. The algorithm is tested on simulated micrographs, obtained by Monte-Carlo methods, where forces between proteins are known analytically. We conclude that valid estimates of such forces, both attractions and repulsions, can be obtained from the positions of a few thousand proteins.     

Fluctuations and membrane heterogeneity

Biological membranes contain many specialized domains, ranging from tens of nanometers to several microns in size and characterized by different concentrations and compositions of protein. Because these domains influence membrane function, considerable attention has focused on understanding their origin. Here it is shown that number fluctuations and nonspecific interprotein interactions can lead to considerable heterogeneity in the distribution of membrane proteins, and to an associated submicron-scale domain structure. Number fluctuations were analyzed by modeling the membrane as a two-dimensional fluid containing interacting protein solutes. The characteristic size and lifetime of a domain in which one would expect to observe a fluctuation of specified magnitude was calculated; snapshots showing fluctuation-induced heterogeneity were generated by Monte Carlo simulation. Domain size was found to depend on the nature of the interprotein force (e.g., attractive or repulsive) and on the average protein concentration. Domain size was largest at low protein concentrations and in the presence of attractive interprotein forces, and was smallest at high protein concentrations and in the presence of repulsive interprotein forces. Domain lifetime was found to depend on domain size and on the diffusion coefficient of the proteins. In a 'typical' membrane containing 5-nm proteins with diffusion coefficient 10(-10) cm(2)/s at a density of 1000 proteins/microm(2), a 30% fluctuation will yield domains characterized by a 2-fold difference in local concentration; these domains persist over a distance of about 100 nm and have a lifetime of about 0.25 s. These results can be used to analyze the domain structure commonly observed in electron micrographs, and have implications for both number fluctuation and Monte Carlo studies of the distribution and dynamics of membrane proteins.     

Suppression of insulin aggregation by heparin

The aggregation of insulin near its isoelectric point and at low ionic strength was suppressed in the presence of heparin. To understand this effect, we used turbidimetry and stopped-flow to study the pH- and ionic strength ( I)-dependence of the aggregation of heparin-free insulin. The results supported the role of interprotein electrostatic interactions, contrary to the commonly held view that such forces are minimized at pH = pI. Electrostatic modeling of insulin (DelPhi) revealed that attractive interactions arise from the marked charge anisotropy of insulin near pI. We show how screening of the interprotein attractions by added salt lead to maximum aggregation near I = 0.01 M, corresponding to a Debye length nearly equal to the diameter of the insulin dimer, consistent with a dipole-like protein charge distribution. This analysis is also consistent with suppression of aggregation by heparin, a strong polyanion that by binding to the positive domain of one protein, inhibits its interaction with the negative domain of another.     

Similar chemistry, but different bond preferences in inter versus intra-protein interactions

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.     

i‐Patch: Interprotein contact prediction using local network information

Biological processes are commonly controlled by precise protein‐protein interactions. These connections rely on specific amino acids at the binding interfaces. Here we predict the binding residues of such interprotein complexes. We have developed a suite of methods, i‐Patch, which predict the interprotein contact sites by considering the two proteins as a network, with residues as nodes and contacts as edges. i‐Patch starts with two proteins, A and B, which are assumed to interact, but for which the structure of the complex is not available. However, we assume that for each protein, we have a reference structure and a multiple sequence alignment of homologues. i‐Patch then uses the propensities of patches of residues to interact, to predict interprotein contact sites. i‐Patch outperforms several other tested algorithms for prediction of interprotein contact sites. It gives 59% precision with 20% recall on a blind test set of 31 protein pairs. Combining the i‐Patch scores with an existing correlated mutation algorithm, McBASC, using a logistic model gave little improvement. Results from a case study, on bacterial chemotaxis protein complexes, demonstrate that our predictions can identify contact residues, as well as suggesting unknown interfaces in multiprotein complexes. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Protein interactions and misfolding analyzed by AFM force spectroscopy

Protein misfolding is conformational transition dramatically facilitating the assembly of protein molecules into aggregates of various morphologies. Spontaneous formation of specific aggregates, mostly amyloid fibrils, was initially believed to be limited to proteins involved in the development of amyloidoses. However, recent studies show that, depending on conditions, the majority of proteins undergo structural transitions leading to the appearance of amyloidogenic intermediates followed by aggregate formation. Various techniques have been used to characterize the protein misfolding facilitating the aggregation process, but no direct evidence as to how such a conformational transition increases the intermolecular interactions has been obtained as of yet. We have applied atomic force microscopy (AFM) to follow the interaction between protein molecules as a function of pH. These studies were performed for three unrelated and structurally distinctive proteins, alpha-synuclein, amyloid beta-peptide (Abeta) and lysozyme. It was shown that the attractive force between homologous protein molecules is minimal at physiological pH and increases dramatically at acidic pH. Moreover, the dependence of the pulling forces is sharp, suggesting a pH-dependent conformational transition within the protein. Parallel circular dichroism (CD) measurements performed for alpha-synuclein and Abeta revealed that the decrease in pH is accompanied by a sharp conformational transition from a random coil at neutral pH to the more ordered, predominantly beta-sheet, structure at low pH. Importantly, the pH ranges for these conformational transitions coincide with those of pulling forces changes detected by AFM. In addition, protein self-assembly into filamentous aggregates studied by AFM imaging was shown to be facilitated at pH values corresponding to the maximum of pulling forces. Overall, these results indicate that proteins at acidic pH undergo structural transition into conformations responsible for the dramatic increase in interprotein interaction and promoting the formation of protein aggregates.     

SuperStar: a knowledge-based approach for identifying interaction sites in proteins.

An empirical method for identifying interaction sites in proteins is described and validated. The method is based entirely on experimental information about non-bonded interactions occurring in small-molecule crystal structures. These data are used in the form of scatterplots that show the experimentally observed distribution of one functional group (the "contact group" or "probe") around another. A template molecule (e.g. a protein binding site) is broken down into structure fragments and the scatterplots, showing the distribution of a chosen probe around these structure fragments, are superimposed on the corresponding parts of the template. The scatterplots are then translated into a three-dimensional map that shows the propensity of the probe at different positions around the template molecule. The method is illustrated for l -arabinose-binding protein, complexed with l -arabinose and with d -fucose, and for dihydrofolate reductase complexed with methotrexate. The method is validated on 122 X-ray structures of protein-ligand complexes. For all the binding sites of these proteins, propensity maps are generated for four different probes: a charged NH+3nitrogen, a carbonyl oxygen, a hydroxyl oxygen and a methyl carbon atom. Next, the maps are compared with the experimentally observed positions of ligand atoms of these types. For 74% of these ligand atoms (84% of the solvent-inaccessible ones) the calculated propensity of the matching probe at the experimental positions is higher than expected by chance. For 68% of the atoms (82% of the solvent-inaccessible ones) the propensity of the matching probe is higher than that of the other three probes. These results indicate that the approach generally gives good predictions for protein-ligand interactions. The potential applications of the propensity maps range from an aid in manual docking and structure-based drug design to their use in pharmacophore development.     

Identifying optimal lipid raft characteristics required to promote nanoscale protein-protein interactions on the plasma membrane

The dynamic lateral segregation of signaling proteins into microdomains is proposed to facilitate signal transduction, but the constraints on microdomain size, mobility, and diffusion that might realize this function are undefined. Here we interrogate a stochastic spatial model of the plasma membrane to determine how microdomains affect protein dynamics. Taking lipid rafts as representative microdomains, we show that reduced protein mobility in rafts segregates dynamically partitioning proteins, but the equilibrium concentration is largely independent of raft size and mobility. Rafts weakly impede small-scale protein diffusion but more strongly impede long-range protein mobility. The long-range mobility of raft-partitioning and raft-excluded proteins, however, is reduced to a similar extent. Dynamic partitioning into rafts increases specific interprotein collision rates, but to maximize this critical, biologically relevant function, rafts must be small (diameter, 6 to 14 nm) and mobile. Intermolecular collisions can also be favored by the selective capture and exclusion of proteins by rafts, although this mechanism is generally less efficient than simple dynamic partitioning. Generalizing these results, we conclude that microdomains can readily operate as protein concentrators or isolators but there appear to be significant constraints on size and mobility if microdomains are also required to function as reaction chambers that facilitate nanoscale protein-protein interactions. These results may have significant implications for the many signaling cascades that are scaffolded or assembled in plasma membrane microdomains.     

Cation-π Interactions and their Functional Roles in Membrane Proteins

Cation-π interactions arise as a result of strong attractive forces between positively charged entities and the π-electron cloud of aromatic groups. The physicochemical characteristics of cation-π interactions are particularly well-suited to the dual hydrophobic/hydrophilic environment of membrane proteins. As high-resolution structural data of membrane proteins bring molecular features into increasingly sharper view, cation-π interactions are gaining traction as essential contributors to membrane protein chemistry, function, and pharmacology. Here we review the physicochemical properties of cation-π interactions and present several prominent examples which demonstrate significant roles for this specialized biological chemistry.     

Protein surface-distribution and protein-protein interactions in the binding of peripheral proteins to charged lipid membranes.

The binding of native cytochrome c to negatively charged lipid dispersions of dioleoyl phosphatidylglycerol has been studied over a wide range of ionic strengths. Not only is the strength of protein binding found to decrease rapidly with increasing ionic strength, but also the binding curves reach an apparent saturation level that decreases rapidly with increasing ionic strength. Analysis of the binding isotherms with a general statistical thermodynamic model that takes into account not only the free energy of the electrostatic double layer, but also the free energy of the surface distribution of the protein, demonstrates that the apparent saturation effects could arise from a competition between the out-of-plane binding reaction and the lateral in-plane interactions between proteins at the surface. It is found that association with nonlocalized sites results in binding isotherms that display the apparent saturation effect to a much more pronounced extent than does the Langmuir adsorption isotherm for binding to localized sites. With the model for nonlocalized sites, the binding isotherms of native cytochrome c can be described adequately by taking into account only the entropy of the surface distribution of the protein, without appreciable enthalpic interactions between the bound proteins. The binding of cytochrome c to dioleoyl phosphatidylglycerol dispersions at a temperature at which the bound protein is denatured on the lipid surface, but is nondenatured when free in solution, has also been studied. The binding curves for the surface-denatured protein differ from those for the native protein in that the apparent saturation at high ionic strength is less pronounced. This indicates the tendency of the denatured protein to aggregate on the lipid surface, and can be described by the binding isotherms for nonlocalized sites only if attractive interactions between the surface-bound proteins are included in addition to the distributional entropic terms. Additionally, it is found that the binding capacity for the native protein is increased at low ionic strength to a value that is greater than that for complete surface coverage, and that corresponds more closely to neutralization of the effective charge (determined from the ionic strength dependence), rather than of the total net charge, on the protein. Electron spin resonance experiments with spin-labeled lipids indicate that this different mode of binding arises from a penetration or disturbance of the bilayer surface by the protein that may alleviate the effects of in-plane interactions under conditions of strong binding.     

Membrane Diffusion Occurs by Continuous-Time Random Walk Sustained by Vesicular Trafficking

Diffusion in cellular membranes is regulated by processes that occur over a range of spatial and temporal scales. These processes include membrane fluidity, interprotein and interlipid interactions, interactions with membrane microdomains, interactions with the underlying cytoskeleton, and cellular processes that result in net membrane movement. The complex, non-Brownian diffusion that results from these processes has been difficult to characterize, and moreover, the impact of factors such as membrane recycling on membrane diffusion remains largely unexplored. We have used a careful statistical analysis of single-particle tracking data of the single-pass plasma membrane protein CD93 to show that the diffusion of this protein is well described by a continuous-time random walk in parallel with an aging process mediated by membrane corrals. The overall result is an evolution in the diffusion of CD93: proteins initially diffuse freely on the cell surface but over time become increasingly trapped within diffusion-limiting membrane corrals. Stable populations of freely diffusing and corralled CD93 are maintained by an endocytic/exocytic process in which corralled CD93 is selectively endocytosed, whereas freely diffusing CD93 is replenished by exocytosis of newly synthesized and recycled CD93. This trafficking not only maintained CD93 diffusivity but also maintained the heterogeneous distribution of CD93 in the plasma membrane. These results provide insight into the nature of the biological and biophysical processes that can lead to significantly non-Brownian diffusion of membrane proteins and demonstrate that ongoing membrane recycling is critical to maintaining steady-state diffusion and distribution of proteins in the plasma membrane.     

Distinguishing structural and functional restraints in evolution in order to identify interaction sites

Structural genomics projects are producing many three-dimensional structures of proteins that have been identified only from their gene sequences. It is therefore important to develop computational methods that will predict sites involved in productive intermolecular interactions that might give clues about functions. Techniques based on evolutionary conservation of amino acids have the advantage over physiochemical methods in that they are more general. However, the majority of techniques neither use all available structural and sequence information, nor are able to distinguish between evolutionary restraints that arise from the need to maintain structure and those that arise from function. Three methods to identify evolutionary restraints on protein sequence and structure are described here. The first identifies those residues that have a higher degree of conservation than expected: this is achieved by comparing for each amino acid position the sequence conservation observed in the homologous family of proteins with the degree of conservation predicted on the basis of amino acid type and local environment. The second uses information theory to identify those positions where environment-specific substitution tables make poor predictions of the overall amino acid substitution pattern. The third method identifies those residues that have highly conserved positions when three-dimensional structures of proteins in a homologous family are superposed. The scores derived from these methods are mapped onto the protein three-dimensional structures and contoured, allowing identification clusters of residues with strong evolutionary restraints that are sites of interaction in proteins involved in a variety of functions. Our method differs from other published techniques by making use of structural information to identify restraints that arise from the structure of the protein and differentiating these restraints from others that derive from intermolecular interactions that mediate functions in the whole organism.     

Extreme C terminus of bacterial cytoskeletal protein FtsZ plays fundamental role in assembly independent of modulatory proteins

Bacterial cell division typically requires assembly of the cytoskeletal protein FtsZ into a ring (Z-ring) at the nascent division site that serves as a foundation for assembly of the division apparatus. High resolution imaging suggests that the Z-ring consists of short, single-stranded polymers held together by lateral interactions. Several proteins implicated in stabilizing the Z-ring enhance lateral interactions between FtsZ polymers in vitro. Here we report that residues at the C terminus of Bacillus subtilis FtsZ (C-terminal variable region (CTV)) are both necessary and sufficient for stimulating lateral interactions in vitro in the absence of modulatory proteins. Swapping the 6-residue CTV from B. subtilis FtsZ with the 4-residue CTV from Escherichia coli FtsZ completely abolished lateral interactions between chimeric B. subtilis FtsZ polymers. The E. coli FtsZ chimera readily formed higher order structures normally seen only in the presence of molecular crowding agents. CTV-mediated lateral interactions are important for the integrity of the Z-ring because B. subtilis cells expressing the B. subtilis FtsZ chimera had a low frequency of FtsZ ring formation and a high degree of filamentation relative to wild-type cells. Site-directed mutagenesis of the B. subtilis CTV suggests that electrostatic forces are an important determinant of lateral interaction potential.     

Use of protein-protein interactions in affinity chromatography.

Biospecific recognition between proteins is a phenomenon that can be exploited for designing affinity-chromatographic purification systems for proteins. In principle, the approach is straightforward, and there are usually many alternative ways, since a protein can be always found which binds specifically enough to the desired protein. Routine immunoaffinity chromatography utilizes the recognition of antigenic epitopes by antibodies. However, forces involved in protein-protein interactions as well the forces keeping the three-dimensional structures of proteins intact are complicated, and proteins are easily unfolded by various factors with unpredictable results. Because of this and because of the generally high association strength between proteins, the correct adjustment of binding forces between an immobilized protein and the protein to be purified as well as the release of bound proteins in biologically active form from affinity complexes are the main problem. Affinity systems involving interactions like enzyme-enzyme, subunit-oligomer, protein-antibody, protein-chaperone and the specific features involved in each case are presented as examples. This article also aims to sketch prospects for further development of the use of protein-protein interactions for the purification of proteins.     

Depletion-mediated red blood cell aggregation in polymer solutions

Polymer-induced red blood cell (RBC) aggregation is of current basic science and clinical interest, and a depletion-mediated model for this phenomenon has been suggested; to date, however, analytical approaches to this model are lacking. An approach is thus described for calculating the interaction energy between RBC in polymer solutions. The model combines electrostatic repulsion due to RBC surface charge with osmotic attractive forces due to polymer depletion near the RBC surface. The effects of polymer concentration and polymer physicochemical properties on depletion layer thickness and on polymer penetration into the RBC glycocalyx are considered for 40 to 500 kDa dextran and for 18 to 35 kDa poly (ethylene glycol). The calculated results are in excellent agreement with literature data for cell-cell affinities and with RBC aggregation-polymer concentration relations. These findings thus lend strong support to depletion interactions as the basis for polymer-induced RBC aggregation and suggest the usefulness of this approach for exploring interactions between macromolecules and the RBC glycocalyx.     

Investigation on the mechanism of crystallization of soluble protein in the presence of nonionic surfactant

The mechanism of crystallization of soluble, globular protein (lysozyme) in the presence of nonionic surfactant C8E4 (tetraoxyethylene glycol monooctyl ether) was examined using both static and dynamic light scattering. The interprotein interaction was found to be attractive in solution conditions that yielded crystals and repulsive in the noncrystallizing solution conditions. The validity of the second virial coefficient as a criterion for predicting protein crystallization could be established even in the presence of nonionic surfactants. Our experiments indicate that the origin of the change in interactions can be attributed to the adsorption of nonionic surfactant monomers on soluble proteins, which is generally assumed to be the case with only membrane proteins. This adsorption screens the hydrophobic attractive force and enhances the hydration and electrostatic repulsive forces between protein molecules. Thus at low surfactant concentration, the effective protein-protein interaction remains repulsive. Large surfactant concentrations promote protein crystallization, possibly due to the attractive depletion force caused by the intervening free surfactant micelles.     

A theoretical study of lipid-protein interactions in bilayers.

We present a theoretical study of the effect of different types of lipid-protein interactions on the thermodynamic properties of protein-containing lipid bilayers. The basis of this work is a theoretical model for pure lipid bilayer phase transitions developed earlier by Scott. Simple assumptions on the nature of the lipid conformations near a protein strongly affect the predicted properties of the model. Here we consider (a) random protein-lipid contacts, (b) enhanced contact between protein and lipid with a number of gauche bonds, and (c) enhanced contact between protein and all-trans but tilted lipid chains. Comparison of predicted results with experimental data seems to favor point c above but, by itself point c does not work well at larger protein concentrations. The results are discussed in the light of spectroscopic data, lipid-protein (plus annular lipid) miscibility, and interprotein forces.     

Bioactive proteinaceous hydrogels from designed bifunctional building blocks

Stimulus-responsive, or "smart" protein-based hydrogels are of interest for many bioengineering applications, but have yet to include biological activity independent of structural functionality. We have genetically engineered bifunctional building blocks incorporating fluorescent proteins that self-assemble into robust and active hydrogels. Gelation occurs when protein building blocks are cross-linked through native protein-protein interactions and the aggregation of alpha-helical hydrogel-forming appendages. Building blocks constructed from different fluorescent proteins can be mixed to enable tuning of fluorescence loading and hydrogel strength with a high degree of independence. FRET experiments suggest a macro-homogeneous structure and that intragel and interprotein reactions can be engineered. This design approach will enable the facile construction of complex hydrogels with broad applicability.     

Starling forces that oppose filtration after tissue oncotic pressure is increased

We tested the hypothesis that the effective oncotic force that opposes fluid filtration across the microvessel wall is the local oncotic pressure difference across the endothelial surface glycocalyx and not the global difference between the plasma and tissue. In single frog mesenteric microvessels perfused and superfused with solutions containing 50 mg/ml albumin, the effective oncotic pressure exerted across the microvessel wall was not significantly different from that measured when the perfusate alone contained albumin at 50 mg/ml. Measurements were made during transient and steady-state filtration at capillary pressures between 10 and 35 cmH(2)O. A cellular-level model of coupled water and solute flows in the interendothelial cleft showed water flux through small breaks in the junctional strand limited back diffusion of albumin into the protected space on the tissue side of the glycocalyx. Thus oncotic forces opposing filtration are larger than those estimated from blood-to-tissue protein concentration differences, and transcapillary fluid flux is smaller than estimated from global differences in oncotic and hydrostatic pressures.     

Formation of the gap junction nexus: binding partners for connexins.

Gap junctions are the morphological correlates of direct cell-cell communication and are formed of hexameric assemblies of gap junction proteins (connexins) into hemichannels (or connexons) provided by each coupled cell. Gap junction channels formed by each of the connexin subtypes (of which there are as many as 20) display different properties, which have been attributed to differences in amino acid sequences of gating domains of the connexins. Recent studies additionally indicate that connexin proteins interact with other cellular components to form a protein complex termed the Nexus. This review summarizes current knowledge regarding the protein-protein interactions involving of connexin proteins and proposes hypothesized functions for these interactions.