Predicting the rate enhancement of protein complex formation from the electrostatic energy of interaction

ProMateus--an open research approach to protein-binding sites analysis

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

Immunophysical exploration of C3d-CR2(CCP1-2) interaction using molecular dynamics and electrostatics

The formation of the complex between the d-fragment of the complement component C3 (C3d) and the modular complement receptor-2 (CR2) is important for cross-linking foreign antigens with surface-bound antibodies and C3d on the surface of B cells. The first two modules of CR2, complement control protein modules (CCPs), participate in non-bonded interactions with C3d. We have used computational methods to analyze the dynamic and electrostatic properties of the C3d-CR2(CCP1-2) complex. The interaction between C3d and CR2 is known to depend on pH and ionic strength. Also, the intermodular mobility of the CR2 modules has been questioned before. We performed a 10 ns molecular dynamics simulation to generate a relaxed structure from crystal packing effects for the C3d-CR2(CCP1-2) complex and to study the energetics of the C3d-CR2(CCP1-2) association. The MD simulation suggests a tendency for intermodular twisting in CR2(CCP1-2). We propose a two-step model for recognition and binding of C3d with CR2(CCP1-2), driven by long and short/medium-range electrostatic interactions. We have calculated the matrix of specific short/medium-range pairwise electrostatic free energies of interaction involved in binding and in intermodular communications. Electrostatic interactions may mediate allosteric effects important for C3d-CR2(CCP1-2) association. We present calculations for the pH and ionic strength-dependence of C3d-CR2(CCP1-2) ionization free energies, which are in overall agreement with experimental binding data. We show how comparison of the calculated and experimental data allows for the decomposition of the contributions of electrostatic from other effects in association. We critically compare predicted stabilities for several mutants of the C3d-CR2(CCP1-2) complex with the available experimental data for binding ability. Finally, we propose that CR2(CCP1-2) is capable of assuming a large array of intermodular topologies, ranging from closed V-shaped to open linear states, with similar recognition properties for C3d, but we cannot exclude an additional contact site with C3d.

The immunoglobulin heavy chain constant region affects kinetic and thermodynamic parameters of antibody variable region interactions with antigen

A central dogma in immunology is that antibody specificity is a function of the variable (V) region. However serological analysis of IgG(1), IgG(2a), and IgG(2b) switch variants of murine monoclonal antibody (mAb) 3E5 IgG(3) with identical V domains revealed apparent specificity differences for Cryptococcus neoformans glucuronoxylomannan (GXM). Kinetic and thermodynamic binding properties of mAbs 3E5 to a 12-mer peptide mimetic of GXM revealed differences in the affinity of these mAbs for a monovalent ligand, a result that implied that the constant (C) region affects the secondary structure of the antigen binding site, thus accounting for variations in specificity. Structural models of mAbs 3E5 suggested that isotype-related differences in binding resulted from amino acid sequence polymorphisms in the C region. This study implies that isotype switching is another mechanism for generating diversity in antigen binding and that isotype restriction of certain antibody responses may reflect structural constraints imposed by C region on V region binding. Furthermore, isotype affected the polyreactivity of V region identical antibodies, implying a role for C region in determining self-reactivity.

Dynamic Docking of Cytochrome b5 with Myoglobin and α-Hemoglobin: Heme-Neutralization “Squares” and the Binding of Electron-Transfer-Reactive Configurations

Intracomplex electron transfer (ET) occurs most often in intrinsically transient, low affinity complexes. As a result, the means by which adequate specificity and reactivity are obtained to support effective ET is still poorly understood. We report here on two such ET complexes: cytochrome b5 (cyt b5) in reaction with its physiological partners, myoglobin (Mb) and hemoglobin (Hb). These complexes obey the Dynamic Docking (DD) paradigm: a large ensemble of weakly bound protein-protein configurations contribute to binding in the rapid-exchange limit, but only a few are ET-active. We report the ionic-strength dependence of the second-order rate constant, k2, for photoinitiated ET from within all four combinations of heme-neutralized Zn deuteroporphyrin-substituted Mb/alphaHb undergoing ET with cyt b5, the four "corners" of a "heme-neutralization square". These experiments provide insights into the relative importance of both global and local electrostatic contributions to the binding of reactive configurations, which are too few to be observed directly. To interpret the variations of k2 arising from heme neutralization, we have developed a procedure by which comparisons of the ET rate constants for a heme-neutralization square permit us to decompose the free energy of reactive binding into individual local electrostatic contributions associated with interactions between (i) the propionates of the two hemes and (ii) the heme of each protein with the polypeptide of its partner. Most notably, we find the contribution from the repulsion between propionates of partner hemes to the reactive binding free energy to be surprisingly small, DeltaG(Hb) approximately +1 kcal/mol at ambient temperature, 18 mM ionic strength, and we speculate about possible causes of this observation. To confirm the fundamental assumption of these studies, that the structure of a heme-neutralized protein is unaltered either by substitution of Zn or by heme neutralization, we have obtained the X-ray structure of ZnMb prepared with the porphyrin dimethyl ester and find it to be nearly isostructural with the native protein.

Kinetic and structural studies of specific protein-protein interactions in substrate catalysis by Cdc25B phosphatase

Using a combination of steady-state and single-turnover kinetics, we probe substrate association, dissociation, and chemistry for the reaction of Cdc25B phosphatase with its Cdk2-pTpY/CycA protein substrate. The rate constant for substrate association for the wild-type enzyme is 1.3 x 10(6) M(-1) s(-1). The rate constant for dissociation is slow compared to the rate constant for phosphate transfer to form the phospho-enzyme intermediate (k2 = 1.1 s(-1)), making Cdk2-pTpY/CycA a sticky substrate. Compared to the wild type, all hotspot mutants of residues at the remote docking site that specifically affect catalysis with the protein substrate (Arg488, Arg492, and Tyr497 on Cdc25B and Asp206 on Cdk2) have greatly slowed rate constants of association (70- to 4500-fold), and some mutants have decreased k2 values compared to that of the wild type. Most dramatically, R492L, despite showing no significant changes in a crystal structure at 2.0 A resolution, has an approximately 100-fold decrease in k2 compared to that of wild-type Cdc25B. The active site C473S mutant binds tightly to and dissociates slowly from Cdk2-pTpY/CycA (Kd = 10 nM, k(off) = 0.01 s(-1)). In contrast, the C473D mutant, despite showing only localized perturbations in the active site at 1.6 A resolution, has a much weaker affinity and dissociates rapidly (Kd of 2 microM, k(off) > 2 s(-1)) from the protein substrate. Overall, we demonstrate that the association of Cdc25B with its Cdk2-pTpY/CycA substrate is governed to a significant extent by the interactions of the remote hotspot residues, whereas dissociation is governed by interactions at the active site.

Inhibition of human pancreatic ribonuclease by the human ribonuclease inhibitor protein

The ribonuclease inhibitor protein (RI) binds to members of the bovine pancreatic ribonuclease (RNase A) superfamily with an affinity in the femtomolar range. Here, we report on structural and energetic aspects of the interaction between human RI (hRI) and human pancreatic ribonuclease (RNase 1). The structure of the crystalline hRI x RNase 1 complex was determined at a resolution of 1.95 A, revealing the formation of 19 intermolecular hydrogen bonds involving 13 residues of RNase 1. In contrast, only nine such hydrogen bonds are apparent in the structure of the complex between porcine RI and RNase A. hRI, which is anionic, also appears to use its horseshoe-shaped structure to engender long-range Coulombic interactions with RNase 1, which is cationic. In accordance with the structural data, the hRI.RNase 1 complex was found to be extremely stable (t(1/2)=81 days; K(d)=2.9 x 10(-16) M). Site-directed mutagenesis experiments enabled the identification of two cationic residues in RNase 1, Arg39 and Arg91, that are especially important for both the formation and stability of the complex, and are thus termed "electrostatic targeting residues". Disturbing the electrostatic attraction between hRI and RNase 1 yielded a variant of RNase 1 that maintained ribonucleolytic activity and conformational stability but had a 2.8 x 10(3)-fold lower association rate for complex formation and 5.9 x 10(9)-fold lower affinity for hRI. This variant of RNase 1, which exhibits the largest decrease in RI affinity of any engineered ribonuclease, is also toxic to human erythroleukemia cells. Together, these results provide new insight into an unusual and important protein-protein interaction, and could expedite the development of human ribonucleases as chemotherapeutic agents.

Temperature dependence of binding and catalysis for the Cdc25B phosphatase

Using a combination of steady-state and single-turnover kinetics, we probe the temperature dependence of substrate association and chemistry for the reaction of Cdc25B phosphatase with its Cdk2-pTpY/CycA protein substrate. The transition state for substrate association is dominated by an enthalpic barrier (DeltaH(++) of 13 kcal/mol) and has a favorable entropic contribution of 4 kcal/mol at 298 K. Phosphate transfer from Cdk2-pTpY/CycA to enzyme (DeltaH(++) of 12 kcal/mol) is enthalpically more favorable than for the small molecule substrate p-nitrophenyl phosphate (DeltaH(++) of 18 kcal/mol), yet entropically less favorable (TDeltaS(++) of 2 vs. -6 kcal/mol at 298 K, respectively). By measuring the temperature dependence of binding and catalysis for several hotspot mutants involved in binding of protein substrate, we determine the enthalpy-entropy compensations for changes in rates of association and phosphate transfer compared to the wild type system. We conclude that the transition state for enzyme-substrate association involves tight and specific contacts at the remote docking site and that phospho-transfer from Cdk2-pTpY/CycA to the pre-organized active site of the enzyme is accompanied by unfavorable entropic rearrangements that promote rapid product dissociation.

Electrostatically driven protein aggregation: beta-lactoglobulin at low ionic strength

The aggregation of beta-lactoglobulin (BLG) at ambient temperature was studied using turbidimetry and dynamic light scattering in the range 3.8<pH<5.2 in 0.0045 M NaCl, and in the ionic strength range 0.0045-0.5 M at fixed pH=5.0. The initial rate of aggregation, taken as the initial slope of turbidity vs time, (dtau/dt)0, indicated maximum aggregation near pH 4.6 (below the isoelectric point of 5.2), but the dependence of the initial rate of aggregation on pH was highly asymmetric. At pH 5.0, (dtau/dt)0 strongly increased with a decrease in ionic strength I from 0.1 to 0.0045 M and was found to be nearly linear with 1/I. DLS measurements revealed an increase in particle size with time, with the appearance of bimodal distributions in which the fast and slow modes corresponded, respectively, to a BLG dimer and to larger aggregates in the 100-800 nm range. At conditions of slower aggregation, DLS revealed the consumption of dimers to form higher order aggregates with no intermediate species. Computer modeling (Delphi) was used to visualize the electrostatic potential around the dimer to elucidate the pH and ionic strength dependence of the initial aggregation rates. The aggregation process appears to comprise an initial fast consumption of the dimer, whose dependence on pH and I arises from the interaction of the positive and negative domains of interacting dimers, followed by the slow formation of much larger aggregates with relatively little sensitivity to pH and I. The open-ended nature of BLG aggregation is thought to arise from the asymmetry of the dimer charge distribution.

Direct Observation of Salt Effects on Molecular Interactions through Explicit-Solvent Molecular Dynamics Simulations: Differential Effects on Electrostatic and Hydrophobic Interactions and Comparisons to Poisson−Boltzmann Theory

Proteins and other biomolecules function in cellular environments that contain significant concentrations of dissolved salts and even simple salts such as NaCl can significantly affect both the kinetics and thermodynamics of macromolecular interactions. As one approach to directly observing the effects of salt on molecular associations, explicit-solvent molecular dynamics (MD) simulations have been used here to model the association of pairs of the amino acid analogues acetate and methylammonium in aqueous NaCl solutions of concentrations 0, 0.1, 0.3, 0.5, 1, and 2 M. By performing simulations of 500 ns duration for each salt concentration properly converged estimates of the free energy of interaction of the two molecules have been obtained for all intermolecular separation distances and geometries. The resulting free energy surfaces are shown to give significant new insights into the way salt modulates interactions between molecules containing both charged and hydrophobic groups and are shown to provide valuable new benchmarks for testing the description of salt effects provided by the simpler but faster Poisson-Boltzmann method. In addition, the complex many-dimensional free energy surfaces are shown to be decomposable into a number of one-dimensional effective energy functions. This decomposition (a) allows an unambiguous view of the qualitative differences between the salt dependence of electrostatic and hydrophobic interactions, (b) gives a clear rationalization for why salt exerts different effects on protein-protein association and dissociation rates, and (c) produces simplified energy functions that can be readily used in much faster Brownian dynamics simulations.

Exploring the charge space of protein-protein association: a proteomic study

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

The electrostatic character of the ribosomal surface enables extraordinarily rapid target location by ribotoxins

Alpha-sarcin ribotoxins comprise a unique family of ribonucleases that cripple the ribosome by catalyzing endoribonucleolytic cleavage of ribosomal RNA at a specific location in the sarcin/ricin loop (SRL). The SRL structure alone is cleaved site-specifically by the ribotoxin, but the ribosomal context enhances the reaction rate by several orders of magnitude. We show that, for the alpha-sarcin-like ribotoxin restrictocin, this catalytic advantage arises from favorable electrostatic interactions with the ribosome. Restrictocin binds at many sites on the ribosomal surface and under certain conditions cleaves the SRL with a second-order rate constant of 1.7 x 10(10) M(-1) s(-1), a value that matches the predicted frequency of random restrictocin-ribosome encounters. The results suggest a mechanism of target location whereby restrictocin encounters ribosomes randomly and diffuses within the ribosomal electrostatic field to the SRL. These studies show a role for electrostatics in protein-ribosome recognition.

Diffusional encounter of barnase and barstar

We present an analysis of trajectories from Brownian dynamics simulations of diffusional protein-protein encounter for the well-studied system of barnase and barstar. This analysis reveals details about the optimal association pathways, the regions of the encounter complex, possible differences of the pathways for dissociation and association, the coupling of translational and rotation motion, and the effect of mutations on the trajectories. We found that a small free-energy barrier divides the energetically most favorable region into a region of the encounter complex above the barnase binding interface and a region around a second energy minimum near the RNA binding loop. When entering the region of the encounter complex from the region near the RNA binding loop, barstar has to change its orientation to increase the electrostatic attraction between the proteins. By concentrating the analysis on the successful binding trajectories, we found that the region of the second minimum is not essential for the binding of barstar to barnase. Nevertheless, this region may be helpful to steer barstar into the region of the encounter complex. When applying the same analysis to several barnase mutants, we found that single mutations may drastically change the free-energy landscape and may significantly alter the population of the two minima. Therefore, certain protein-protein pairs may require careful adaptation of the positions of encounter and transition states when interpreting mutation effects on kinetic rates of association and/or dissociation.

Rational design of intercellular adhesion molecule-1 (ICAM-1) variants for antagonizing integrin lymphocyte function-associated antigen-1-dependent adhesion

The interaction between integrin lymphocyte function-associated antigen-1 (LFA-1) and its ligand intercellular adhesion molecule-1 (ICAM-1) is critical in immunological and inflammatory reactions but, like other adhesive interactions, is of low affinity. Here, multiple rational design methods were used to engineer ICAM-1 mutants with enhanced affinity for LFA-1. Five amino acid substitutions 1) enhance the hydrophobicity and packing of residues surrounding Glu-34 of ICAM-1, which coordinates to a Mg2+ in the LFA-1 I domain, and 2) alter associations at the edges of the binding interface. The affinity of the most improved ICAM-1 mutant for intermediate- and high-affinity LFA-1 I domains was increased by 19-fold and 22-fold, respectively, relative to wild type. Moreover, potency was similarly enhanced for inhibition of LFA-1-dependent ligand binding and cell adhesion. Thus, rational design can be used to engineer novel adhesion molecules with high monomeric affinity; furthermore, the ICAM-1 mutant holds promise for targeting LFA-1-ICAM-1 interaction for biological studies and therapeutic purposes.

IL-1beta epitope mapping using site-directed mutagenesis and hydrogen-deuterium exchange mass spectrometry analysis

Hu007, a humanized IgG1 monoclonal antibody, binds and neutralizes human, cynomolgus, and rabbit IL-1beta but only weakly binds to mouse and rat IL-1beta. Biacore experiments demonstrated that Hu007 and the type-I IL-1 receptor competed for binding to IL-1beta. Increasing salt concentrations decrease the association rate with only moderate effects on the dissociation rate, suggesting that long-range electrostatics are critical for formation of the initial complex. To understand the ligand-binding specificity of Hu007, we have mapped the critical residues involved in the recognition of IL-1beta. Selected residues in cynomolgus IL-1beta were mutated to the corresponding residues in mouse IL-1beta, and the effects of the changes on binding were evaluated by surface plasmon resonance measurements using Biacore. Specifically, substitution of F150S decreased binding affinity by 100-fold, suggesting the importance of hydrophobic interactions in stabilizing the antibody/antigen complex. Substitution of three amino acids near the N- and C-terminal regions of cIL-1beta with those found in mouse IL-1beta (V3I/S5Q/F150S) decreased the binding affinity of Hu007 to IL-1beta by about 1000-fold. Conversely, mutating the corresponding residues in mouse IL-1beta to the human sequence resulted in an increase in binding affinity of about 1000-fold. Hydrogen-deuterium exchange/mass spectrometry analysis confirmed that these regions of IL-1beta were protected from exchange because of antibody binding. The results from this study demonstrate that Hu007 binds to a region located in the open end of the beta-barrel structure of IL-1beta and blocks binding of IL-1beta to its receptor.

Prediction of water and metal binding sites and their affinities by using the Fold-X force field

The empirical force field Fold-X was developed previously to allow rapid free energy calculations in proteins. Here, we present an enhanced version of the force field allowing prediction of the position of structural water molecules and metal ions, together called single atom ligands. Fold-X picks up 76% of water molecules found to interact with two or more polar atoms of proteins in high-resolution crystal structures and predicts their position to within 0.8 A on average. The prediction of metal ion-binding sites have success rates between 90% and 97% depending on the metal, with an overall standard deviation on the position of binding of 0.3-0.6 A. The following metals were included in the force field: Mg2+, Ca2+, Zn2+, Mn2+, and Cu2+. As a result, the current version of Fold-X can accurately decorate a protein structure with biologically important ions and water molecules. Additionally, the free energy of binding of Ca2+ and Zn2+ (i.e., the natural logarithm of the dissociation constant) and its dependence on ionic strength correlate reasonably well with the experimental data available in the literature, allowing one to discriminate between high- and low-affinity binding sites. Importantly, the accuracy of the energy prediction presented here is sufficient to efficiently discriminate between Mg2+, Ca2+, and Zn2+ binding.

Two different T cell receptors use different thermodynamic strategies to recognize the same peptide/MHC ligand

A6 and B7 are two alphabeta T cell receptors (TCRs) that recognize the Tax peptide presented by the class I major histocompatibility molecule HLA-A2 (Tax/HLA-A2). Despite the fact that the two TCRs have different CDR loops and use different amino acid residues to contact their ligand, both receptors bind ligand with similar diagonal orientations. Here we show that they also bind with very similar binding affinities and kinetics (the DeltaDeltaG degrees for binding is approximately 0.3kcal/mol at 25 degrees C). The two receptors respond similarly to alterations in the MHC molecule, yet differ dramatically in their responses to ionic strength and temperature. The different responses to temperature indicate markedly different binding thermodynamics, which are not predictable from the surface area buried in the interfaces. A6 and B7 thus represent two TCRs that are both compatible with Tax/HLA-A2, although compatibility has been achieved through the use of different thermodynamic strategies. Finally, neither A6 nor B7 are predicted to undergo large conformational adaptations upon binding, distinguishing them from a number of other TCRs whose structure, thermodynamics, and kinetics have been characterized.

A detailed thermodynamic analysis of ras/effector complex interfaces

Many cellular functions are based on the interaction and crosstalk of various signaling proteins. Among these, members of the Ras family of small GTP-binding proteins are important for communicating signals into different pathways. In order to answer the question of how binding affinity and specificity is achieved, we analyzed binding energetics on the molecular level, with reference to the available structural data. The interaction of two members of the Ras subfamily with two different effector proteins, namely Raf and RalGDS, were investigated using isothermal titration calorimetry and a fluorescence-based method. Experiments with alanine mutants, located in the complex interfaces, yielded an energy map for the contact areas of the Ras/effector complexes, which could be differentiated into enthalpy and entropy contributions. In addition, by using double mutant cycle analysis, we probed the energetic contribution of selected pairs of amino acid residues. The resulting energy landscapes of the Ras/effector interface areas show a highly different topology when comparing the two effectors, Raf and RalGDS, demonstrating the specificity of the respective interactions. Particularly, we observe a high degree of compensating effects between enthalpy and entropy; differences between these components are much greater than the overall free energy differences. This is observed also when using the software FOLD-X to predict the effect of point mutations on the crystal structures of the different complexes. Prediction of the free energy changes shows a very good correlation with the experimentally observed energies. Furthermore, in line with experimental data, energy decomposition indicates that many different components of large magnitude counteract each other to produce a smaller change in overall free energy, illustrating the importance of long-range electrostatic forces in complex formation.

A perspective on mechanisms of protein tetramer formation

Homotetrameric proteins can assemble by several different pathways, but have only been observed to use one, in which two monomers associate to form a homodimer, and then two homodimers associate to form a homotetramer. To determine why this pathway should be so uniformly dominant, we have modeled the kinetics of tetramerization for the possible pathways as a function of the rate constants for each step. We have found that competition with the other pathways, in which homotetramers can be formed either by the association of two different types of homodimers or by the successive addition of monomers to homodimers and homotrimers, can cause substantial amounts of protein to be trapped as intermediates of the assembly pathway. We propose that this could lead to undesirable consequences for an organism, and that selective pressure may have caused homotetrameric proteins to evolve to assemble by a single pathway.

A Brownian dynamics study: the effect of a membrane environment on an electron transfer system

During the past few years, three-dimensional crystal structures of many of the important integral membrane proteins responsible for the bioenergetic processes of photosynthesis and respiration have been determined. Moreover, a few crystal structures of protein-protein complexes have become available that characterize the interaction between those membrane proteins and the electron carrier protein cytochrome c. Here, we address the association kinetics for binding of cytochrome c to cytochrome c oxidase (COX) from Paracoccus denitrificans by Brownian dynamics simulations. The effects of ionic strength and protein mutations were studied for two different cytochrome c species: the positively charged, dipolar horse heart cytochrome c and the negatively charged physiological electron transfer partner cytochrome c(552). We studied association toward "naked" COX and toward membrane-embedded COX where the membrane is represented as an uncharged DPPC bilayer modeled in atomistic detail. For the nonnatural association toward "naked" COX, the association rates are >100 times larger for horse heart cytochrome c than for cytochrome c(552). Interestingly, the presence of the lipid bilayer leads to a dramatic decrease of the association rate of horse heart cytochrome c, but slightly enhances association of cytochrome c(552), leading to very similar association rates of both proteins to membrane-embedded COX. This finding from computational modeling studies may reflect the optimization of surface patches and of the total net charge on electron transfer pairs in nature.

Realistic protein-protein association rates from a simple diffusional model neglecting long-range interactions, free energy barriers, and landscape ruggedness

We develop a simple but rigorous model of protein-protein association kinetics based on diffusional association on free energy landscapes obtained by sampling configurations within and surrounding the native complex binding funnels. Guided by results obtained on exactly solvable model problems, we transform the problem of diffusion in a potential into free diffusion in the presence of an absorbing zone spanning the entrance to the binding funnel. The free diffusion problem is solved using a recently derived analytic expression for the rate of association of asymmetrically oriented molecules. Despite the required high steric specificity and the absence of long-range attractive interactions, the computed rates are typically on the order of 10(4)-10(6) M(-1) sec(-1), several orders of magnitude higher than rates obtained using a purely probabilistic model in which the association rate for free diffusion of uniformly reactive molecules is multiplied by the probability of a correct alignment of the two partners in a random collision. As the association rates of many protein-protein complexes are also in the 10(5)-10(6) M(-1) sec(-1) range, our results suggest that free energy barriers arising from desolvation and/or side-chain freezing during complex formation or increased ruggedness within the binding funnel, which are completely neglected in our simple diffusional model, do not contribute significantly to the dynamics of protein-protein association. The transparent physical interpretation of our approach that computes association rates directly from the size and geometry of protein-protein binding funnels makes it a useful complement to Brownian dynamics simulations.

Association and dissociation kinetics for CheY interacting with the P2 domain of CheA

The chemotaxis system of Escherichia coli makes use of an extended two-component sensory response pathway in which CheA, an autophosphorylating protein histidine kinase (PHK) rapidly passes its phosphoryl group to CheY, a phospho-accepting response regulator protein (RR). The CheA-->CheY phospho-transfer reaction is 100-1000 times faster than the His-->Asp phospho-relays that operate in other (non-chemotaxis) two-component regulatory systems, suggesting that CheA and CheY have unique features that enhance His-->Asp phospho-transfer kinetics. One such feature could be the P2 domain of CheA. P2 encompasses a binding site for CheY, but an analogous RR-binding domain is not found in other PHKs. In previous work, we removed P2 from CheA, and this decreased the catalytic efficiency of CheA-->CheY phospho-transfer by a factor of 50-100. Here we examined the kinetics of the binding interactions between CheY and P2. The rapid association reaction (k(assn) approximately 10(8)M(-1)s(-1) at 25 degrees C and micro=0.03 M) exhibited a simple first-order dependence on P2 concentration and appeared to be largely diffusion-limited. Ionic strength (micro) had a moderate effect on k(assn) in a manner predictable based on the calculated electrostatic interaction energy of the protein binding surfaces and the expected Debye-Hückel shielding. The speed of binding reflects, in part, electrostatic interactions, but there is also an important contribution from the inherent plasticity of the complex and the resulting flexibility that this allows during the process of complex formation. Our results support the idea that the P2 domain of CheA contributes to the overall speed of phospho-transfer by promoting rapid association between CheY and CheA. However, this alone does not account for the ability of the chemotaxis system to operate much more rapidly than other two-component systems: k(cat) differences indicate that CheA and CheY also achieve the chemical events of phospho-transfer more rapidly than do PHK-RR pairs of slower systems.

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.

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.

Redesigning an antibody fragment for faster association with its antigen

Traditional approaches for increasing the affinity of protein-protein complexes focus on constructing highly complementary binding surfaces. Recent theoretical simulations and experimental results suggest that electrostatic steering forces can also be manipulated to increase association rates while leaving dissociation rates unchanged, thus increasing affinity. Here we demonstrate that electrostatic attraction can be enhanced between an antibody fragment and its cognate antigen through application of a few simple rules to identify potential on-rate amplification sites that lie at the periphery of the antigen-antibody interface.

Kinetic analysis of sequence-specific recognition of ssDNA by an autoantibody

11F8 is a pathogenic monoclonal anti-ssDNA autoantibody isolated from a lupus prone mouse. Previous studies established that 11F8 is sequence-specific and identified the thermodynamic and kinetic basis for the specific recognition of ssDNA, and binding site mutations of a single-chain construct reveal that (Y32)LCDR1, (R31)HCDR1, (W33)HCDR1, (R98)HCDR3, (L97)HCDR3, and (Y100)HCDR3 are responsible for approximately 80% of the binding free energy. Here we evaluate the role of these residues along with a group of basic residues (K62, K64, R24, K52) within the context of the binding mechanism. Binding of 11F8 takes place in two steps. In the first step, the overall positive charge of the antigen binding site attracts the negatively charged DNA to form an encounter complex that is stabilized by two salt bridges and a hydrogen bond. The second step is a slow process in which minor conformational changes occur. During this step, aromatic side chains become desolvated, presumably through stacking interactions involving two thymine bases within the DNA recognition epitope. Although the stability of the complex arises primarily from interactions formed in the second step, sequence specificity results from interactions with residues involved in both steps. These studies also show that the way in which 11F8 achieves high affinity sequence-specific binding is more closely related to RNA binding proteins than those that bind DNA and point to strategies for disrupting DNA binding that could prove to be therapeutically useful.

Multiple metal ions drive DNA association by PvuII endonuclease

Restriction enzymes serve as important model systems for understanding the role of metal ions in phosphodiester hydrolysis. To this end, a number of laboratories have reported dramatic differences between the metal ion-dependent and metal ion-independent DNA binding behaviors of these systems. In an effort to illuminate the underlying mechanistic details which give rise to these differences, we have quantitatively dissected these equilibrium behaviors into component association and dissociation rates for the representative PvuII endonuclease and use these data to assess the stoichiometry of metal ion involvement in the binding process. The dependence of PvuII cognate DNA on Ca(II) concentration binding appears to be cooperative, exhibiting half-saturation at 0.6 mM metal ion and yielding an n(H) of 3.5 +/- 0.2 per enzyme homodimer. Using both nitrocellulose filter binding and fluorescence assays, we observe that the cognate DNA dissociation rate (k(-)(1) or k(off)) is very slow (10(-)(3) s(-)(1)) and exhibits a shallow dependence on metal ion concentration. DNA trap cleavage experiments with Mg(II) confirm the general irreversibility of DNA binding relative to cleavage, even at low metal ion concentrations. More dramatically, the association rate (k(1) or k(on)) also appears to be cooperative, increasing more than 100-fold between 0.2 and 10 mM Ca(II), with an optimum value of 2.7 x 10(7) M(-)(1) s (-)(1). Hill analysis of the metal ion dependence of k(on) indicates an n(H) of 3.6 +/- 0.2 per enzyme dimer. This value is consistent with the involvement in DNA association of two metal ions per subunit active site, a result which lends new strength to arguments for two-metal ion mechanisms in restriction enzymes.

Glycosylation of erythropoietin affects receptor binding kinetics: role of electrostatic interactions

Erythropoietin (EPO) is a cytokine produced by the kidney whose function is to stimulate red blood cell production in the bone marrow. Previously, it was shown that the affinity of EPO for its receptor, EPOR, is inversely related to the sialylation of EPO carbohydrate. To better understand the properties of EPO that modulate its receptor affinity, various glycoforms were analyzed using surface plasmon resonance. The system used has been well characterized and is based on previous reports employing an EPOR-Fc chimera captured on a Protein A surface. Using three variants of EPO containing different levels of sialylation, we determined that sialic acid decreased the association rate constant (k(on)) about 3-fold. Furthermore, glycosylated EPO had a 20-fold slower k(on) than nonglycosylated EPO, indicating that the core carbohydrate also negatively impacted k(on). The effect of electrostatic forces on EPO binding was studied by measuring binding kinetics in varying NaCl concentrations. Increasing NaCl concentration resulted in a slower k(on) while having little impact on k(off), suggesting that long-range electrostatic interactions are primarily important in determining the rate of association between EPO and EPOR. Furthermore, the glycosylation content (i.e., nonglycosylated vs glycosylated, sialylated vs desialylated) affected the overall sensitivities of k(on) to [NaCl], indicating that sialic acid and the glycan itself each impact the overall effect of these electrostatic forces.

MIAX: a new paradigm for modeling biomacromolecular interactions and complex formation in condensed phases

A new paradigm is proposed for modeling biomacromolecular interactions and complex formation in solution (protein-protein interactions so far in this report) that constitutes the scaffold of the automatic system MIAX (acronym for Macromolecular Interaction Assessment X). It combines in a rational way a series of computational methodologies, the goal being the prediction of the most native-like protein complex that may be formed when two isolated (unbound) protein monomers interact in a liquid environment. The overall strategy consists of first inferring putative precomplex structures by identification of binding sites or epitopes on the proteins surfaces and a simultaneous rigid-body docking process using geometric instances alone. Precomplex configurations are defined here as all those decoys the interfaces of which comply substantially with the inferred binding sites and whose free energy values are lower. Retaining all those precomplex configurations with low energies leads to a reasonable number of decoys for which a flexible treatment is amenable. A novel algorithm is introduced here for automatically inferring binding sites in proteins given their 3-D structure. The procedure combines an unsupervised learning algorithm based on the self-organizing map or Kohonen network with a 2-D Fourier spectral analysis. To model interaction, the potential function proposed here plays a central role in the system and is constituted by empirical terms expressing well-characterized factors influencing biomacromolecular interaction processes, essentially electrostatic, van der Waals, and hydrophobic. Each of these procedures is validated by comparing results with observed instances. Finally, the more demanding process of flexible docking is performed in MIAX embedding the potential function in a simulated annealing optimization procedure. Whereas search of the entire configuration hyperspace is a major factor precluding hitherto systems from efficiently modeling macromolecular interaction modes and complex structures, the paradigm presented here may constitute a step forward in the field because it is shown that a rational treatment of the information available from the 3-D structure of the interacting monomers combined with conveniently selected computational techniques can assist to elude search of regions of low probability in configuration space and indeed lead to a highly efficient system oriented to solve this intriguing and fundamental biologic problem.

On the role of electrostatic interactions in the design of protein-protein interfaces

Here, the methods of continuum electrostatics are used to investigate the contribution of electrostatic interactions to the binding of four protein-protein complexes; barnase-barstar, human growth hormone and its receptor, subtype N9 influenza virus neuraminidase and the NC41 antibody, the Ras binding domain (RBD) of kinase cRaf and a Ras homologue Rap1A. In two of the four complexes electrostatics are found to strongly oppose binding (hormone-receptor and neuraminidase-antibody complexes), in one case the net effect is close to zero (barnase-barstar) and in one case electrostatics provides a significant driving force favoring binding (RBD-Rap1A). In order to help understand the wide range of electrostatic contributions that were calculated, the electrostatic free energy was partitioned into contributions of individual charged and polar residues, salt bridges and networks involving salt bridges and hydrogen bonds. Although there is no one structural feature that accounts for the differences between the four interfaces, the extent to which the desolvation of buried charges is compensated by the formation of hydrogen bonds and ion pairs appears to be an important factor. Structural features that are correlated with contribution of an individual residue to stability are also discussed. These include partial burial of a charged group in the free monomer, the formation of networks involving charged and polar amino acids, and the formation of partially exposed ion-pairs. The total electrostatic contribution to binding is found to be inversely correlated with buried total and non-polar surface area. This suggests that different interfaces can be designed to exploit electrostatic and hydrophobic forces in very different ways.

Biomolecular diffusional association

The primary method for the computational study of biomolecular diffusional association is Brownian dynamics. Recent work has seen advances in the efficiency of computing association rates and in the accuracy of simulation models. New areas to which Brownian dynamics has been applied include protein polymerisation and protein adsorption to a surface. There has recently been particularly intense study of protein-protein association, and Brownian dynamics, together with other theoretical and experimental approaches, has led to new insights into the determinants of protein-protein binding kinetics.

Kinetic studies of protein-protein interactions

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

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.

Mutational screen identifies critical amino acid residues of beta-actin mediating interaction between its folding intermediates and eukaryotic cytosolic chaperonin CCT

The three-dimensional reconstruction of apo-CCT-alpha-actin by cryoelectron microscopy shows that actin binds either the CCTbeta-CCTdelta or the CCTepsilon-CCTdelta subunit pairs of the chaperonin in an open and apparently quasi-native conformation. The CCT-binding sites are seen located at the tips of the two arms of actin and these same regions of actin have been implicated in CCT binding through beta-actin peptide-array screening. Three main CCT binding regions exist: actin Sites I, II, and III, which are composed of loops that are surface-exposed in native actin. Sixty-eight amino acid residues on beta-actin have been screened by mutagenesis for effects on CCT interaction in quantitative in vitro translation assays in rabbit reticulocyte lysate. Actin seems to be folding cooperatively on chaperonin, since certain mutants discriminate CCT binding from processing. Actin Site II, located at the tip of actin subdomain 4, is the major determinant for CCT binding. Site II is composed of two anti-parallel extended beta-strands, with F200-T203 and D244 contributing substantially to the binding site. The substrate recognition chemistry of CCT thus seems different from that of Group I chaperonins and probably reflects the fact that it needs to be highly specific to enable capture and folding of the actins and tubulins.

Experimental assignment of the structure of the transition state for the association of barnase and barstar

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

Protein-protein association: investigation of factors influencing association rates by brownian dynamics simulations

The rate of protein-protein association limits the response time due to protein-protein interactions. The bimolecular association rate may be diffusion-controlled or influenced, and in such cases, Brownian dynamics simulations of protein-protein diffusional association may be used to compute association rates. Here, we report Brownian dynamics simulations of the diffusional association of five different protein-protein pairs: barnase and barstar, acetylcholinesterase and fasciculin-2, cytochrome c peroxidase and cytochrome c, the HyHEL-5 antibody and hen egg lysozyme (HEL), and the HyHEL-10 antibody and HEL. The same protocol was used to compute the diffusional association rates for all the protein pairs in order to assess, by comparison to experimentally measured rates, whether the association of these proteins can be explained solely on the basis of diffusional encounter. The simulation protocol is similar to those previously derived for simulation of the association of barnase and barstar, and of acetylcholinesterase and fasciculin-2; these produced results in excellent agreement with experimental data for these protein pairs, with changes in association rate due to mutations reproduced within the limits of expected computational and modeling errors. Here, we find that for all protein pairs, the effects of mutations can be well reproduced by the simulations, even though the degree of the electrostatic translational and orientational steering varies widely between the cases. However, the absolute values of association rates for the acetylcholinesterase: fasciculin-2 and HyHEL-10 antibody: HEL pairs are overestimated. Comparison of bound and unbound protein structures shows that this may be due to gating resulting from protein flexibility in some of the proteins. This may lower the association rates compared to their bimolecular diffusional encounter rates.

Substrate specificity of the integral membrane protease OmpT determined by spatially addressed peptide libraries

Escherichia coli outer membrane protease T (OmpT) is an endopeptidase that specifically cleaves between two consecutive basic residues. In this study we have investigated the substrate specificity of OmpT using spatially addressed SPOT peptide libraries. The peptide acetyl-Dap(dnp)-Ala-Arg/Arg-Ala-Lys(Abz)-Gly was synthesized directly onto cellulose membrane. The peptide contained the aminobenzoyl (Abz) fluorophore, which was internally quenched by the dinitrophenyl (dnp) moiety. Treatment of the SPOT membrane with the small, water-soluble protease trypsin resulted in highly fluorescent peptide SPOTs. However, no peptide cleavage was observed after incubation with detergent-solubilized OmpT, a macromolecular complex with an estimated molecular mass of 180 kDa. This problem could be solved by the introduction of a long, polar polyoxyethylene glycol linker between the membrane support and the peptide. Peptide libraries for the P(2), P(1), P(1)', and P(2)' positions in the substrate were screened with OmpT, and peptides of positive SPOTs were resynthesized and subjected to kinetic measurements in solution. The best substrate Abz-Ala-Lys-Lys-Ala-Dap(dnp)-Gly had a turnover number k(cat) of 40 s(-)(1), which is 12-fold higher than the starting substrate. Peptides containing an acidic residue at P(2) or P(2)' were not substrates for OmpT, suggesting that long-range electrostatic interactions are important for the formation of the enzyme-substrate complex. OmpT was highly selective toward L-amino acids at P(1) but was less so at P(1)' where a peptide with D-Arg at P(1)' was a competitive inhibitor (K(i) of 19 microM). An affinity chromatography resin based on these findings was developed, which allowed for the one-step purification of OmpT from a bacterial lysate. The implications of the determined consensus substrate sequence (Arg/Lys)/(Arg/Lys)-Ala for the proposed biological function of OmpT in defense against antimicrobial peptides are discussed.

New structural and functional aspects of the type I interferon-receptor interaction revealed by comprehensive mutational analysis of the binding interface

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

Evaluation of direct and cooperative contributions towards the strength of buried hydrogen bonds and salt bridges

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

Mutational and structural analysis of the binding interface between type I interferons and their receptor Ifnar2

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

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 A 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.

Biophysical analysis of the interaction of human ifnar2 expressed in E. coli with IFNalpha2

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