New insights into the mechanism of protein-protein association

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

How FMN binds to anabaena apoflavodoxin: a hydrophobic encounter at an open binding site

Molecular recognition begins when two molecules approach and establish interactions of certain strength. The mechanisms of molecular recognition reactions between biological molecules are not well known, and few systems have been analyzed in detail. We investigate here the reaction between an apoprotein and its physiological cofactor (apoflavodoxin and flavin mononucleotide) that binds reversibly to form a non-covalent complex (flavodoxin) involved in electron transfer reactions. We have analyzed the fast binding reactions between the FMN cofactor (and shorter analogs) and wild type (and nine mutant apoflavodoxins where residues interacting with FMN in the final complex have been replaced). The x-ray structures of two such mutants are reported that show the mutations are well tolerated by the protein. From the calculated microscopic binding rate constants we have performed a Phi analysis of the transition state of complex formation that indicates that the binding starts by interaction of the isoalloxazine-fused rings in FMN with residues of its hydrophobic binding site. In contrast, the phosphate in FMN, known to contribute most to the affinity of the final holoflavodoxin complex, is not bound in the transition state complex. Both the effects of ionic strength and of phosphate concentration on the wild type complex rate constant agree with this scenario. As suggested previously by nmr data, it seems that the isoalloxazine-binding site may be substantially open in solution. Interestingly, although FMN is a charged molecule, electrostatic interactions seem not to play a role in directing the binding, unlike what has been reported for other biological complexes. The binding can thus be best described as a hydrophobic encounter at an open binding site.     

On the dynamic nature of the transition state for protein-protein association as determined by double-mutant cycle analysis and simulation

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

Fruitful and Futile Encounters along the Association Reaction between Proteins

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

Introducing intrinsic disorder reduces electrostatic steering in protein-protein interactions

Protein-protein interactions underlie many critical biology functions, such as cellular signaling and gene expression, in which electrostatic interactions can play a critical role in mediating the specificity and stability of protein complexes. A substantial portion of proteins are intrinsically disordered, and the influences of structural disorder on binding kinetics and thermodynamics have been widely investigated. However, whether the effect of electrostatic steering depends on structural disorder remains unexplored. In this work, we addressed the consequence of introducing intrinsic disorder in the electrostatic steering of the E3/Im3 complex using molecular dynamics simulation. Our results recapitulated the experimental observations that the responses of stability and kinetics to salt concentration for the ordered E3/Im3 complex were larger than those for the disordered E3/Im3 complex. Mechanistic analysis revealed that the native contact interactions involved in the encounter state and the transition state were essentially identical for both ordered and disordered E3. Therefore, the observed difference in electrostatic steering between ordered E3 and disordered E3 may result from their difference in conformation rather than their difference in binding mechanism. Because charged residues are frequently involved in protein-protein interactions, our results suggest that increasing structural disorder is expected to generally modulate the effect of electrostatic steering.     

Mutations in the S4 region isolate the final voltage-dependent cooperative step in potassium channel activation

Charged residues in the S4 transmembrane segment play a key role in determining the sensitivity of voltage-gated ion channels to changes in voltage across the cell membrane. However, cooperative interactions between subunits also affect the voltage dependence of channel opening, and these interactions can be altered by making substitutions at uncharged residues in the S4 region. We have studied the activation of two mutant Shaker channels that have different S4 amino acid sequences, ILT (V369I, I372L, and S376T) and Shaw S4 (the S4 of Drosophila Shaw substituted into Shaker), and yet have very similar ionic current properties. Both mutations affect cooperativity, making a cooperative transition in the activation pathway rate limiting and shifting it to very positive voltages, but analysis of gating and ionic current recordings reveals that the ILT and Shaw S4 mutant channels have different activation pathways. Analysis of gating currents suggests that the dominant effect of the ILT mutation is to make the final cooperative transition to the open state of the channel rate limiting in an activation pathway that otherwise resembles that of Shaker. The charge movement associated with the final gating transition in ILT activation can be measured as an isolated component of charge movement in the voltage range of channel opening and accounts for 13% ( approximately 1.8 e0) of the total charge moved in the ILT activation pathway. The remainder of the ILT gating charge (87%) moves at negative voltages, where channels do not open, and confirms the presence of Shaker-like conformational changes between closed states in the activation pathway. In contrast to ILT, the activation pathway of Shaw S4 seems to involve a single cooperative charge-moving step between a closed and an open state. We cannot detect any voltage-dependent transitions between closed states for Shaw S4. Restoring basic residues that are missing in Shaw S4 (R1, R2, and K7) rescues charge movement between closed states in the activation pathway, but does not alter the voltage dependence of the rate-limiting transition in activation.     

On the protein–protein diffusional encounter complex

When two proteins diffuse together to form a bound complex, an intermediate is formed at the end‐point of diffusional association which is called the encounter complex. Its characteristics are important in determining association rates, yet its structure cannot be directly observed experimentally. Here, we address the problem of how to construct the ensemble of three‐dimensional structures which constitute the protein–protein diffusional encounter complex using available experimental data describing the dependence of protein association rates on mutation and on solvent ionic strength and viscosity. The magnitude of the association rates is fitted well using a variety of definitions of encounter complexes in which the two proteins are located at up to about 17 Å root‐mean‐squared distance from their relative arrangement in the bound complex. Analysis of the ionic strength dependence of bimolecular association rates shows that this is determined to a greater extent by the (protein charge) – (salt ion) separation distance than by the protein–protein charge separation distance. Consequently, ionic strength dependence of association rates provides little information about the geometry of the encounter complex. On the other hand, experimental data on electrostatic rate enhancement, mutation and viscosity dependence suggest a model of the encounter complex in which the two proteins form a subset of the contacts present in the bound complex and are significantly desolvated. Copyright © 1999 John Wiley & Sons, Ltd.     

Incomplete ion dissociation underlies the weakened attraction between DNA helices at high spermidine concentrations

We have investigated the salt sensitivity of the hexagonal-to-cholesteric phase transition of spermidine-condensed DNA. This transition precedes the resolubilization of precipitated DNA that occurs at high spermidine concentration. The sensitivity of the critical spermidine concentration at the transition point to the anion species and the NaCl concentration indicates that ion pairing of this trivalent ion underlies this unusual transition. Osmotic pressure measurements of spermidine salt solutions are consistent with this interpretation. Spermidine salts are not fully dissociated at higher concentrations. The competition for DNA binding among the fully charged trivalent ion and the lesser charged complex species at higher concentrations significantly weakens attraction between DNA helices in the condensed state. This is contrary to the suggestion that the binding of spermidine at higher concentrations causes DNA overcharging and consequent electrostatic repulsion.     

Role of ionic interactions in ligand binding and catalysis of R67 dihydrofolate reductase

R67 dihydrofolate reductase (DHFR), which catalyzes the NADPH dependent reduction of dihydrofolate to tetrahydrofolate, belongs to a type II family of R-plasmid encoded DHFRs that confer resistance to the antibacterial drug trimethoprim. Crystal structure data reveals this enzyme is a homotetramer that possesses a single active site pore. Only two charged residues in each monomer are located near the pore, K32 and K33. Site-directed mutants were constructed to probe the role of these residues in ligand binding and/or catalysis. As a result of the 222 symmetry of this enzyme, mutagenesis of one residue results in modification at four related sites. All mutants at K32 affected the quaternary structure, producing an inactive dimer. The K33M mutant shows only a 2-4-fold effect on K(m) values. Salt effects on ligand binding and catalysis for K33M and wildtype R67 DHFRs were investigated to determine if these lysines are involved in forming ionic interactions with the negatively charged substrates, dihydrofolate (overall charge of -2) and NADPH (overall charge of -3). Binding studies indicate that two ionic interactions occur between NADPH and R67 DHFR. In contrast, the binding of folate, a poor substrate, to R67 DHFR.NADPH appears weak as a titration in enthalpy is lost at low ionic strength. Steady-state kinetic studies for both wild type (wt) and K33M R67 DHFRs also support a strong electrostatic interaction between NADPH and the enzyme. Interestingly, quantitation of the observed salt effects by measuring the slopes of the log of ionic strength versus the log of k(cat)/K(m) plots indicates that only one ionic interaction is involved in forming the transition state. These data support a model where two ionic interactions are formed between NADPH and symmetry related K32 residues in the ground state. To reach the transition state, an ionic interaction between K32 and the pyrophosphate bridge is broken. This unusual scenario likely arises from the constraints imposed by the 222 symmetry of the enzyme.     

Peptide Binding to a PDZ Domain by Electrostatic Steering via Nonnative Salt Bridges

We have captured the binding of a peptide to a PDZ domain by unbiased molecular dynamics simulations. Analysis of the trajectories reveals on-pathway encounter complex formation, which is driven by electrostatic interactions between negatively charged carboxylate groups in the peptide and positively charged side chains surrounding the binding site. In contrast, the final stereospecific complex, which matches the crystal structure, features completely different interactions, namely the burial of the hydrophobic side chain of the peptide C-terminal residue and backbone hydrogen bonds. The simulations show that nonnative salt bridges stabilize kinetically the encounter complex during binding. Unbinding follows the inverse sequence of events with the same nonnative salt bridges in the encounter complex. Thus, in contrast to protein folding, which is driven by native interactions, the binding of charged peptides can be steered by nonnative interactions, which might be a general mechanism, e.g., in the recognition of histone tails by bromodomains.     

GroEL stability and function. Contribution of the ionic interactions at the inter-ring contact sites

The chaperonin GroEL consists of a double ring structure made of identical subunits that display different modes of allosteric communication. The protein folding cycle requires the simultaneous positive intra-ring and negative inter-ring cooperativities of ATP binding. This ensures GroES binding to one ring and release of the ligands from the opposite one. To better characterize inter-ring allosterism, the thermal stability as well as the temperature dependence of the functional and conformational properties of wild type GroEL, a single ring mutant (SR1) and two single point mutants suppressing one interring salt bridge (E434K and E461K) were studied. The results indicate that ionic interactions at the two interring contact sites are essential to maintain the negative cooperativity for protein substrate binding and to set the protein thermostat at 39 degrees C. These electrostatic interactions contribute distinctly to the stability of the inter-ring interface and the overall protein stability, e.g. the E434K thermal inactivation curve is shifted to lower temperatures, and its unfolding temperature and activation energy are also lowered. An analysis of the ionic interactions at the inter-ring contact sites reveals that at the so called "left site" a network of electrostatic interactions involving three charged residues might be established, in contrast to what is found at the "right site" where only two oppositely charged residues interact. Our data suggest that electrostatic interactions stabilize protein-protein interfaces depending on both the number of ionic interactions and the number of residues engaged in each of these interactions. In the case of GroEL, this combination sets the thermostat of the protein so that the chaperonin distinguishes physiological from stress temperatures.     

Ion Specificity and Nonmonotonic Protein Solubility from Salt Entropy

The addition of salt to protein solutions can either increase or decrease the protein solubility, and the magnitude of this effect depends on the salt used. We show that these effects can be captured using a theory that includes attractive and repulsive electrostatic interactions, nonelectrostatic protein-ion interactions, and ion-solvent interactions via an effective solvated ion radius. We find that the ion radius has significant effects on the translational entropy of the salt, which leads to salt specificity in the protein solubility. At low salt, the dominant effect comes from the entropic cost of confining ions within the aggregate, whereas at high concentrations, the salt drives a depletion attraction that favors aggregation. Our theory explains the reversal in the Hofmeister series observed in lysozyme cloud point measurements and semi-quantitatively describes the solubility of lysozyme and chymosin crystals. We present a comparison of the contributions to the free energy and give guidelines for when salting in or salting out should be expected.     

Molecular Interactions and Residues Involved in Force Generation in the T4 Viral DNA Packaging Motor

Many viruses utilize molecular motors to package their genomes into preformed capsids. A striking feature of these motors is their ability to generate large forces to drive DNA translocation against entropic, electrostatic, and bending forces resisting DNA confinement. A model based on recently resolved structures of the bacteriophage T4 motor protein gp17 suggests that this motor generates large forces by undergoing a conformational change from an extended to a compact state. This transition is proposed to be driven by electrostatic interactions between complementarily charged residues across the interface between the N- and C-terminal domains of gp17. Here we use atomistic molecular dynamics simulations to investigate in detail the molecular interactions and residues involved in such a compaction transition of gp17. We find that although electrostatic interactions between charged residues contribute significantly to the overall free energy change of compaction, interactions mediated by the uncharged residues are equally if not more important. We identify five charged residues and six uncharged residues at the interface that play a dominant role in the compaction transition and also reveal salt bridging, van der Waals, and solvent hydrogen-bonding interactions mediated by these residues in stabilizing the compact form of gp17. The formation of a salt bridge between Glu309 and Arg494 is found to be particularly crucial, consistent with experiments showing complete abrogation in packaging upon Glu309Lys mutation. The computed contributions of several other residues are also found to correlate well with single-molecule measurements of impairments in DNA translocation activity caused by site-directed mutations.     

Continuum electrostatic model for the binding of cytochrome c2 to the photosynthetic reaction center from Rhodobacter sphaeroides

Electrostatic interactions are important for protein-protein association. In this study, we examined the electrostatic interactions between two proteins, cytochrome c(2) (cyt c(2)) and the reaction center (RC) from the photosynthetic bacterium Rhodobacter sphaeroides, that function in intermolecular electron transfer in photosynthesis. Electrostatic contributions to the binding energy for the cyt c(2)-RC complex were calculated using continuum electrostatic methods based on the recent cocrystal structure [Axelrod, H. L., et al. (2002) J. Mol. Biol. 319, 501-515]. Calculated changes in binding energy due to mutations of charged interface residues agreed with experimental results for a protein dielectric constant epsilon(in) of 10. However, the electrostatic contribution to the binding energy for the complex was close to zero due to unfavorable desolvation energies that compensate for the favorable Coulomb attraction. The electrostatic energy calculated as a function of displacement of the cyt c(2) from the bound position showed a shallow minimum at a position near but displaced from the cocrystal configuration. These results show that although electrostatic steering is present, other short-range interactions must be present to contribute to the binding energy and to determine the structure of the complex. Calculations made to model the experimental data on association rates indicate a solvent-separated transition state for binding in which the cyt c(2) is displaced approximately 8 A above its position in the bound complex. These results are consistent with a two-step model for protein association: electrostatic docking of the cyt c(2) followed by desolvation to form short-range van der Waals contacts for rapid electron transfer.     

Theoretical study of structural transition in a nucleosome at low ionic strength

The theoretical analysis of nucleosome stability at low ionic strength has been performed on the basis of consideration of different contributions to the free energy of compact state of the nucleosome DNA terminal regions. The proposed model explains: the fact of low-salt structural change; the transition point (approximately 1.7 mM NaCl) and width (approximately 1 mM); the shift of the transition to the higher salt concentrations in the case of histones tails removal by trypsin. According to the model the increase of electrostatic repulsion between neighbouring turns of DNA superhelix is the main cause of the unwinding of nucleosomal DNA terminal regions in the course of low-salt structural change. The interactions between histone (H2A-H2B) dimer and (H3-H4)2 tetramer provide the compact state of the nucleosomal DNA terminal regions. The existence of electrostatic interactions of nucleosomal DNA terminal regions with tetramer was suggested. These interactions can provide the compact state of nucleosomal DNA at physiological ionic strength even in the absence of (H2A-H2B) dimer.     

Apparent Debye-Huckel electrostatic effects in the folding of a simple, single domain protein

We have monitored the effects of salts and denaturants on the folding of the simple, two-state protein FynSH3. As predicted by Debye-Huckel limiting law, both the stability and (log) folding rate of FynSH3 increase nearly perfectly linearly (r(2)> 0.99) with the square root of ionic strength upon increasing concentrations of the relatively nonchaotropic salt sodium chloride. The stability of FynSH3 is also linear in square root ionic strength when the relatively nonchaotropic salts sodium bromide, potassium bromide, and potassium chloride are employed. Comparison of the kinetic and equilibrium effects of sodium chloride suggests that the electrostatic interactions formed in the folding transition state are approximately 50% as destabilizing as those formed in the native state, presumably reflecting the more compact nature of the latter. In contrast, the relationship between concentration and folding kinetics is more complex when the highly chaotropic salt guanidine hydrochloride (GuHCl) is employed. At moderate to high GuHCl concentrations the net effect of the linear, presumably chaotrope-induced deceleration and the presumed, square root-dependent ionic strength-induced acceleration is well approximated as linear, thus accounting for the observation of "chevron behavior" (log folding rate linear in denaturant concentration) typically reported for the folding of single domain proteins. At very low GuHCl concentrations, however, significant kinetic rollover is observed. This rollover is reasonably well fitted as a sum of a linear, presumably chaotropic effect and a square root-dependent, presumably electrostatic effect. These results thus not only provide insight into the nature of the folding transition state but also suggest that caution is in order when extrapolating GuHCl-based chevrons to estimate folding rates in the absence of denaturant and in interpreting deviations from chevron linearity as evidence for non-two-state kinetics.     

Influence of the pyruvate content of xanthan on macromolecular association in solution

The unusual increase in viscosity and pseudoplasticity often observed when salts are added to moderately concentrated aqueous solutions of xanthan gum is shown to arise from an increase in the extent of macromolecular association. The fractional change in viscosity on addition of KCl to salt-free 1% (w/v) solutions of purified polysaccharide in the K⁺ salt form is found to be positive only when the degree of pyruvate substitution (fraction of side chains which carry pyruvate ketal substituents) exceeds ≈0.31. Above this value, the fractional change in viscosity increases with further increase in the degree of pyruvate substitution. These differences cannot arise from different degrees of conformational ordering, since the magnitude of the thermally induced order disorder transition (monitored by optical rotation at low ionic strength) is independent of pyruvate content. The temperature at the transition midpoint, however, falls with increasing degree of pyruvate substitution. This is attributed to destabilization of the ordered structure by intramolecular electrostatic repulsion between pyruvate groups, and stabilization through apolar interactions of acetate methyl groups. Viscosity-concentration relationships show changes of slope which mark the onset of macromolecular association. Association commences at lower concentrations when ionic strength and degree of pyruvate substitution are high. It is suggested that once electrostatic repulsions have been diminished at high ionic strength, association is promoted by intermolecular apolar interactions of pyruvate methyl groups, which are suitably near the periphery of the helical conformation.     

Ionic-strength-dependent effects in protein folding: analysis of rate equilibrium free-energy relationships and their interpretation

The traditional approach to studying protein folding involves applying a perturbation, usually denaturant or mutation, and determining the effect upon the free energy of folding, DeltaG0, and the activation free energy, DeltaG(not equal). Data collected as a function of the perturbation can be used to construct rate equilibrium free-energy relationships, which report on the development of interactions in the transition state for folding. We examine the use of the ionic-strength-dependent rate equilibrium free-energy relationship in protein folding using the N-terminal domain of L9, a small alpha-beta protein, as a model system. Folding is two-state for the range of ionic strength examined, 0.045-1.52 M. The plot of DeltaG(not equal) versus DeltaG0 is linear (r2= 0.918), with a slope equal to 0.45. The relatively low value of the slope indicates that the ionic-strength-dependent interactions are modestly developed in the transition state. The slope is, however, greater than that of a plot of DeltaG(not equal) versus DeltaG0 constructed by varying pH, thus demonstrating directly that ionic-strength-dependent studies probe more than simple electrostatic interactions. Potential transition movement was probed by analysis of the denaturant, ionic strength cross-interaction parameters. The values are small but nonzero and positive, suggesting a small shift of the transition state toward the native state as the protein is destabilized, i.e., Hammond behavior. The complications that arise in the interpretation of ionic-strength-dependent rate equilibrium free-energy relationships are discussed, and it is concluded that the ionic-strength-dependent studies do not provide a reliable indicator of the role of electrostatic interactions. Complications include incomplete screening of electrostatic interactions, specific ion binding, Hofmeister effects, and the potential presence of electrostatic interactions in the denatured state ensemble.