Role of charges on cytochrome f from the cyanobacterium Phormidium laminosum in its interaction with plastocyanin

The role of charge on the surface of cytochrome f from the cyanobacterium Phormidium laminosum in the reaction with plastocyanin was investigated in vitro using site-directed mutagenesis. Charge was neutralized at five acidic residues individually and introduced at a residue close to the interface between the two proteins. The effects on the kinetics of the reaction were measured using stopped-flow spectrophotometry, and the midpoint potentials of the mutant proteins were determined. The dependence of the bimolecular rate constant of reaction, k(2), on ionic strength was determined for the reactions of the cytochrome f mutants with wild-type and mutant forms of plastocyanin. Double mutant cycle analysis was carried out to probe for the presence of specific electrostatic interactions. The effects of mutations on Cyt f were smaller than those seen previously for mutants of plastocyanin [Schlarb-Ridley, B. G. et al. (2002) Biochemistry 41, 3279-3285]. One specific short-range interaction between charged residues of wild-type plastocyanin (Arg93) and wild-type cytochrome f (Asp63) was identified. The kinetic evidence from this study and that of Schlarb-Ridley et al., 2002, appears to conflict with the NMR structure of the P. laminosum complex, which suggests the absence of electrostatic interactions in the final complex [Crowley, P. et al. (2001) J. Am. Chem. Soc. 123, 10444-10453]. The most likely explanation of the apparent paradox is that the overall rate is diffusion controlled and that electrostatics specifically influence the encounter complex and not the reaction complex.     

Electron-transfer mechanism of the triplet state quenching of aluminium tetrasulfonated phthalocyanine by cytochrome c

The mechanism of electron-transfer from aluminium tetrasulfonated phthalocyanine triplet state to cytochrome c was investigated in this work. This reaction successfully quenches the dye triplet state due to the formation of complexes between the solute and the protein at the active site. The electron-transfer rate constant is around 3x10(7) s(-1), and is in accordance with previous results for the singlet excited state quenching [C.A.T. Laia, S.M.B. Costa, D. Phillips, A. Beeby. Electron-transfer kinetics in sulfonated aluminum phthalocyanines/cytochrome c complexes, J. Phys. Chem. B 108 (2004) 7506-7514.] in the framework of the Marcus theory, with a reorganization energy equal to 0.94 eV. The complex formation is diffusion controlled, but heterogeneities of the protein surface charge distribution lead to quenching rate constants smaller than predicted on a hard-spheres model with electrostatic interactions. Also the binding equilibrium constant is strongly affected by this phenomenon. Ionic strength plays an important role on the complex formation, but its effect on the unimolecular electron-transfer rate constant is negligible within experimental error.     

Direct simulation of plastocyanin and cytochrome f interactions in solution

Most biological functions, including photosynthetic activity, are mediated by protein interactions. The proteins plastocyanin and cytochrome f are reaction partners in a photosynthetic electron transport chain. We designed a 3D computer simulation model of diffusion and interaction of spinach plastocyanin and turnip cytochrome f in solution. It is the first step in simulating the electron transfer from cytochrome f to photosystem 1 in the lumen of thylakoid. The model is multiparticle and it can describe the interaction of several hundreds of proteins. In our model the interacting proteins are represented as rigid bodies with spatial fixed charges. Translational and rotational motion of proteins is the result of the effect of stochastic Brownian force and electrostatic force. The Poisson-Boltzmann formalism is used to determine the electrostatic potential field generated around the proteins. Using this model we studied the kinetic characteristics of plastocyanin-cytochrome f complex formation for plastocyanin mutants at pH 7 and a variety of ionic strength values.     

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.     

Concerted simulations reveal how peroxidase compound III formation results in cellular oscillations

A major problem in mathematical modeling of the dynamics of complex biological systems is the frequent lack of knowledge of kinetic parameters. Here, we apply Brownian dynamics simulations, based on protein three-dimensional structures, to estimate a previously undetermined kinetic parameter, which is then used in biochemical network simulations. The peroxidase-oxidase reaction involves many elementary steps and displays oscillatory dynamics important for immune response. Brownian dynamics simulations were performed for three different peroxidases to estimate the rate constant for one of the elementary steps crucial for oscillations in the peroxidase-oxidase reaction, the association of superoxide with peroxidase. Computed second-order rate constants agree well with available experimental data and permit prediction of rate constants at physiological conditions. The simulations show that electrostatic interactions depress the rate of superoxide association with myeloperoxidase, bringing it into the range necessary for oscillatory behavior in activated neutrophils. Such negative electrostatic steering of enzyme-substrate association presents a novel control mechanism and lies in sharp contrast to the electrostatically-steered fast association of superoxide and Cu/Zn superoxide dismutase, which is also simulated here. The results demonstrate the potential of an integrated and concerted application of structure-based simulations and biochemical network simulations in cellular systems biology.     

Direct electrochemistry of protein-protein complexes involving cytochrome c, cytochrome b5, and plastocyanin

The direct electrochemistry of the cytochrome c/cytochrome b5 and cytochrome c/plastocyanin complexes has been investigated at edge-plane graphite and modified gold electrode surfaces, which are selective for one of the two components of the complex. Electrochemical response of one protein at an otherwise electrostatically unfavorable electrode surface was achieved in the presence of the other protein, and the calculated heterogeneous electron-transfer rate constant and diffusion coefficient were found to be in good agreement with the values determined previously from the electrochemistry of the individual proteins [Armstrong, F. A., Hill, H. A. O., & Walton, N. J. (1988) Acc. Chem. Res. 21, 407 and references therein]. A dynamic model of the protein-protein-electrode ternary complex is proposed to explain the promotion effect, and this model is supported by a study comparing the electrochemical responses of covalent and electrostatic cytochrome c/plastocyanin complexes. It is also suggested that the behavior of protein-protein complexes at electrode surfaces could be related to that of the complexes associated with biological membranes.     

Computer simulation of complex formation between plastocyanine and cytochrome f in thylakoid lumen

The diffusion of the protein plastocyanine and complex formation between plastocyanine and cytochrome f (a subunit of a cytochrome b6/f complex) in the chloroplast thylakoid lumen has been studied. A 3D computer simulation model of diffusion and binding of plastocyanine and cytochrome f has been constructed, which considers their electrostatic interaction. Based on the experimental data, the parameters of the model for complex formation between plastocyanine and cytochrome f in solution have been estimated. The dependence of the rate of plastocyanine-cytochrome f reaction on the size of the luminal space has been studied. It was shown that the contraction of the luminal space leads to a decrease in the reaction rate, which is in agreement with the experimental data on the inhibition of the reaction under hyperosmotic stress.     

Electron transfer between flavodoxin semiquinone and c-type cytochromes: correlations between electrostatically corrected rate constants, redox potentials, and surface topologies

We have measured the ionic strength dependence of the rate constants for electron transfer from the semiquinone of Clostridium pasteurianum flavodoxin to 12 c-type cytochromes and several inorganic oxidants using stopped-flow methodology. The experimental data were fit quite well by an electrostatic model that represents the interaction domains as parallel disks with a point charge equal to the charge within this region of the protein. The analysis provides an evaluation of the electrostatic interaction energy and the rate constant at infinite ionic strength (k affinity). The electrostatic charge on the oxidant within the interaction site can be obtained from the electrostatic energy, and for most of those reactants for which structures are available, the results are in good agreement with expectation. The k affinity values were found to correlate with redox potential differences, as expected from the theory of adiabatic (or nonadiabatic) outer-sphere electron-transfer reactions. Deviations from the theoretical curves are interpreted in terms of the influence of surface topology on reaction rate constants. In general, we find that electrostatic effects, steric influences, and redox potential all exert a much larger effect on reaction rate constants for the flavodoxin-cytochrome system than has been previously observed for free flavin-cytochrome interactions. The implications of this for determining biological specificity are discussed.     

Laser flash-induced kinetic analysis of cytochrome f oxidation by wild-type and mutant plastocyanin from the cyanobacterium Nostoc sp. PCC 7119

Oxidation of the soluble, truncated form of cytochrome f by wild-type and mutant species of plastocyanin has been analyzed by laser flash absorption spectroscopy in the cyanobacterium Nostoc (formerly, Anabaena) sp. PCC 7119. At low ionic strengths, the apparent electron transfer rate constant of cytochrome f oxidation by wild-type plastocyanin is 1.34 x 10(4) s(-)(1), a value much larger than those determined for the same proteins from other organisms. Upon site-directed mutagenesis of specific residues at the plastocyanin interaction area, the rate constant decreases in all cases yet to varying extents. The only exception is the D54K variant, which exhibits a higher reactivity toward cytochrome f. In most cases, the reaction rate constant decreases monotonically with an increase in ionic strength. The observed changes in the reaction mechanism and rate constants are in agreement with the location of the mutated residues at the interface area, as well as with the peculiar orientation of the two partners within the Nostoc plastocyanin-cytochrome f transient complex, whose NMR structure has been determined recently. Furthermore, the experimental data herein reported match well the kinetic behavior exhibited by the same set of plastocyanin mutants when acting as donors of electrons to photosystem I [Molina-Heredia, F. P., et al. (2001) J. Biol. Chem. 276, 601-605], thus indicating that the copper protein uses the same surface areas-one hydrophobic and the other electrostatic-to interact with both cytochrome f and photosystem I.     

Brownian dynamics study of the interaction between plastocyanin and cytochrome f.

The electrostatic interaction between plastocyanin (PC) and cytochrome f (cyt f), electron transfer partners in photosynthesis was studied using Brownian dynamics (BD) simulations. By using the software package MacroDox, which implements the BD algorithm of Northrup et al. (Northrup, S. H., J. O. Boles, and J. C. L. Reynolds. 1987. J. Phys. Chem. 91:5991-5998), we have modeled the interaction of the two proteins based on crystal structures of poplar PC and turnip cyt f at pH 7 and a variety of ionic strengths. We find that the electrostatic attraction between positively charged residues (K58, K65, K187, and R209, among others) on cyt f and negatively charged residues (E43, D44, E59, and E60, among others) on PC steers PC into a single dominant orientation with respect to cyt f, and furthermore, that the single dominant orientation that we observe is one that we had predicted in our previous work (Pearson, D. C., E. L. Gross, and E. S. David. 1996. Biophys. J. 71:64-76). This dominant orientation permits the formation of hydrophobic interactions, which are not implemented in the MacroDox algorithm. This proposed complex between PC and cyt f implicates H87, a copper ligand on PC, as the residue that accepts electrons from the heme on cyt f (and possibly through Y1 as we proposed previously). We argue for the existence of this single dominant complex on the basis of observations that the most favorable orientations of the interaction between PC and cyt f, as determined by grouping successful BD trajectories on the basis of closest contacts of charged residues, tend to overlap one another and have very close distances between the metal centers on the two proteins (copper on PC, iron on cyt f). We use this knowledge to develop a model for PC/cyt f interaction that places a reaction between the two proteins occurring when the copper-to-iron distance is between 16 and 17 A. This reaction distance gives a good estimate of the experimentally observed rate constant for PC-cyt f interaction. Analysis of BD results as a function of ionic strength predicts an interaction that happens less frequently and becomes less specific as ionic strength increases.     

Mapping electrostatic interactions in macromolecular associations.

In the association of electron transfer proteins, electrostatics has been proposed to play a role in maintaining the stability and specificity of the biomolecular complexes formed. An excellent model system is the interaction between mammalian cytochrome b5 and cytochrome c, in which the X-ray structures of the individual components reveal a complementary asymmetry of charges surrounding their respective redox centers. Determining the exact extent of the electrostatic interactions and identifying the specific residues involved in the formation of the electron transfer complex has proved more elusive. We report herein the utilization of high-pressure techniques, together with site-directed mutagenesis, to provide a map of the interaction domains in biomolecular complex formation. The application of high pressure disrupts macromolecular associations since dissociation of the complex results in a decreased volume of the system due to the solvation of charges that had been previously sequestered in the interface region and force solvation of hydrophobic surfaces. Site-directed mutagenesis of a totally synthetic gene for rat liver cytochrome b5, which expresses this mammalian protein in Escherichia coli as a hemecontaining soluble component, was used to selectively alter negatively charged residues of cytochrome b5 to neutral amide side-chains. We have demonstrated that the interaction domain of cytochrome b5 with cytochrome c can be mapped from a comparison of dissociation volumes of these modified cytochrome b5-cytochrome c complexes with the native complex. Using these techniques we can specifically investigate the role of particular residues in the equilibrium association of these two electron transfer proteins. Single-point mutations in the interaction domain give nearly identical effects on the measured dissociation volumes, yet removal of acidic residues outside the recognition surface yield volumes similar to wild-type protein. Multiple mutations in the proposed protein-protein interaction site are found to allow greater solvent-accessibility of the interface as reflected in a diminution in the volume changes on subsequent charge removal. This is indicative that the interprotein salt-bridges in this complex provide a mechanism for a greater exclusion of solvent from the interfacial domain of the complex, resulting in a more stable association.     

Study of the cytochrome c2-reaction center interaction by site-directed mutagenesis.

Photooxidation of Rhodobacter capsulatus cytochrome c2 and four site-directed mutants by detergent solubilized Rhodobacter sphaeroides reaction centers was studied as a function of ionic strength at pH 8.0. Mutants of cytochrome c2 included K12D (lysine 12 substituted by aspartate), K14E (lysine 14 substituted by glutamate), K32E (lysine 32 substituted by glutamate), and K14E/K32E (lysines 14 and 32 substituted by glutamates). With respect to the wild-type, the mutants exhibited decreased second-order rate constants, indicating perturbation of their electrostatic interaction with the reaction center. In the transient complex, the interaction domain charges of the reaction center and wild-type cytochrome c2 were estimated to be -4.8 and +4.8, respectively. In contrast, the interaction domain charges of mutants K12D, K14E, K32E, and K14E/K32E were estimated to be +2.8, +3.7, +3.6 and +1.3, respectively. At infinite ionic strength, the second-order rate constant of the wild-type cytochrome c2 photooxidation (k infinity) was estimated to be 8.7 x 10(6) M-1 s-1. In the case of K32E, k infinity was not changed significantly (8.2 x 10(6) m-1 s-1), suggesting that the electrostatic perturbation of this mutant was largely overcome at high ionic strength. In contrast, the k infinity for K12D, K14E, and K14E/K32E were estimated to be decreased 2-7-fold. Consequently, mutations to R. capsulatus lysines 12 and 14 appear to perturb the distance and/or orientation of the cytochrome c2 relative to the reaction center in the reactive complex, as well as alter electrostatic interactions. Based upon the kinetic results presented here, the cytochrome c2-reaction center transient complex has been modeled.     

Transient kinetics of electron transfer reactions of flavodoxin: ionic strength dependence of semiquinone oxidation by cytochrome c, ferricyanide, and ferric ethylenediaminetetraacetic acid and computer modeling of reaction complexes

Electron transfer reactions between Clostridum pasteurianum flavodoxin semiquinone and various oxidants [horse heart cytochrome c, ferricyanide, and ferric ethylenediaminetetraacetic [horse heart cytochrome c, ferricyanide, and ferric ethylenediaminetetraacetic acid (EDTA)] have been studied as a function of ionic strength by using stopped-flow spectrophotometry. The cytochrome c reaction is complicated by the existence of two cytochrome species which react at different rates and whose relative concentrations are ionic strength dependent. Only the faster of these two reactions is considered here. At low ionic strength, complex formation between cytochrome c and flavodoxin is indicated by a leveling off of the pseudo-first-order rate constant at high cytochrome c concentration. This is not observed for either ferricyanide or ferric EDTA. For cytochrome c, the rate and association constants for complex formation were found to increase with decreasing ionic strength, consistent with negative charges on flavodoxin interacting with the positively charged cytochrome electron transfer site. Both ferricyanide and ferric EDTA are negatively charged oxidants, and the rate data respond to ionic strength changes as would be predicted for reactants of the same charge sign. These results demonstrate that electrostatic interactions involving negatively charged groups are important in orienting flavodoxin with respect to oxidants during electron transfer. We have also carried out computer modeling studies of putative complexes of flavodoxin with cytochrome c and ferricyanide, which relate their structural properties to both the observed kinetic behavior and some more general features of physiological electron transfer processes. The results of this study are consistent with the ionic strength behavior described above.     

Effect of Macromolecular Crowding on Reaction Rates: A Computational and Theoretical Study

The effect of macromolecular crowding on the rates of association reactions are investigated using theory and computer simulations. Reactants and crowding agents are both hard spheres, and when two reactants collide they form product with a reaction probability, p_rxn. A value of p_rxn < 1 crudely mimics the fact that proteins must be oriented properly for an association reaction to occur. The simulations show that the dependence of the reaction rate on the volume fraction of crowding agents varies with the reaction probability. For reaction probabilities close to unity where most of encounters between reactants lead to a reaction, the reaction rate always decreases as the volume fraction of crowding agents is increased due to the reduced diffusion coefficient of reactants. On the other hand, for very small reaction probabilities where, in most of encounters, the reaction does not occur, the reaction rate increases with the volume fraction of crowding agents—in this case, due to the increase probability of a recollision. The Smoluchowski theory refined with the radiation boundary condition and the radial distribution function at contact is in quantitative agreement with simulations for the reaction rate constant and allows the quantitative analysis of both effects separately.     

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

The rate of association of proteins is dictated by diffusion, but can be enhanced by favorable electrostatic forces. Here the relationship between the electrostatic energy of interaction, and the kinetics of protein-complex formation was analyzed for the protein pairs of: hirudin-thrombin, acetylcholinesterase-fasciculin and barnase-barstar, and for a panel of point mutants of these proteins. Electrostatic energies of interaction were calculated as the difference between the electrostatic energy of the complex and the sum of the energies of the two individual proteins, using the computer simulation package DelPhi. Calculated electrostatic energies of interaction were compared to experimentally determined rates of association. One kcal/mol of Coulombic interaction energy increased the rate of association by a factor of 2.8, independent of the protein-complex or mutant analyzed. Electrostatic energies of interaction were also determined from the salt dependence of the association rate constant, using the same basic equation as for the theoretical calculation. A Br?nsted analysis of the electrostatic energies of interactions plotted versus experimentally determined ln(rate)s of association shows a linear relation between the two, with a beta value close to 1. This is interpreted as the energy of the transition state varies according to the electrostatic interaction energy, fitting a two state model for the association reaction. Calculating electrostatic rate enhancement from the electrostatic interaction energy can be used as a powerful tool to design protein complexes with altered rates of association and affinities.     

Brownian dynamics simulations of molecular recognition in an antibody-antigen system.

The crystal structure for an antibody-antigen system, that of the anti-hen egg lysozyme monoclonal antibody HyHEL-5 complexed to lysozyme, is used as the starting point for computer simulations of diffusional encounters between the two proteins. The investigation consists of two parts: first, the linearized Poisson-Boltzmann equation is solved to determine the long-range electrostatic forces between antibody and antigen, and then, the relative motion as influenced by these forces is modeled within Brownian motion theory. The effects of various point mutations on the calculated reaction rate are considered. It is found that charged residues close to the binding site exert the greatest influence in steering the proteins into a configuration favorable for their binding, while more distant mutations are qualitatively described by the Smoluchowski model for the mutual diffusion of two uniformly charged spheres. The antibody residues involved in forming salt links with the lysozyme, Glu-H35 and Glu-H50, appear to be particularly important in electrostatic steering, as neutralization of both of them yields reaction rates that are two to three orders of magnitude below those of wild-type rates. The relative rates obtained from the simulations can be tested through kinetic measurements on mutant protein complexes. Kinetically efficient partners can also be designed and constructed through directed mutagenesis.     

The reaction of cytochrome aa3 with (porphyrin) cytochrome c as studied by pulse radiolysis

(1) Using the pulse-radiolysis and stopped-flow techniques, the reactions of iron-free (porphyrin) cytochrome c and native cytochrome c with cytochrome aa3 were investigated. The porphyrin cytochrome c anion radical (generated by reduction of porphyrin cytochrome c by the hydrated electron) can transfer its electron to cytochrome aa3. The bimolecular rate constant for this reaction is 2 x 10(7) M-1 . s-1 (5 mM potassium phosphate, 0.5% Tween 20, pH 7.0, 20 degrees C). (2) The ionic strength dependence of the cytochrome c-cytochrome aa3 interaction was measured in the ionic strength range between 40 and 120 mM. At ionic strengths below 30 mM, a cytochrome c-cytochrome aa3 complex is formed in which cytochrome c is no longer reducible by the hydrated electron. A method is described by which the contributions of electrostatic forces to the reaction rate can be determined. (3) Using the stopped-flow technique, the effect of the dielectric constant (epsilon) of the reaction medium on the reaction of cytochrome C with cytochrome aa3 was investigated. With increasing epsilon the second-order rate constant decreased.     

Electrostatic analysis and Brownian dynamics simulation of the association of plastocyanin and cytochrome f.

The oxidation of cytochrome f by the soluble cupredoxin plastocyanin is a central reaction in the photosynthetic electron transfer chain of all oxygenic organisms. Here, two different computational approaches are used to gain new insights into the role of molecular recognition and protein-protein association processes in this redox reaction. First, a comparative analysis of the computed molecular electrostatic potentials of seven single and multiple point mutants of spinach plastocyanin (D42N, E43K, E43N, E43Q/D44N, E59K/E60Q, E59K/E60Q/E43N, Q88E) and the wt protein was carried out. The experimentally determined relative rates (k(2)) for the set of plastocyanin mutants are found to correlate well (r(2) = 0.90 - 0.97) with the computed measure of the similarity of the plastocyanin electrostatic potentials. Second, the effects on the plastocyanin/cytochrome f association rate of these mutations in the plastocyanin "eastern site" were evaluated by simulating the association of the wild type and mutant plastocyanins with cytochrome f by Brownian dynamics. Good agreement between the computed and experimental relative rates (k(2)) (r(2) = 0.89 - 0.92) was achieved for the plastocyanin mutants. The results obtained by applying both computational techniques provide support for the fundamental role of the acidic residues at the plastocyanin eastern site in the association with cytochrome f and in the overall electron-transfer process.     

Kinetics of intracomplex electron transfer and of reduction of the components of covalent and noncovalent complexes of cytochrome c and cytochrome c peroxidase by free flavin semiquinones

The kinetics of reduction of free flavin semiquinones of the individual components of 1:1 covalent and electrostatic complexes of yeast ferric and ferryl cytochrome c peroxidase and ferric horse cytochrome c have been studied. Covalent cross-linking between the peroxidase and cytochrome c at low ionic strength results in a complex that has kinetic properties both similar to and different from those of the electrostatic complex. Whereas the cytochrome c heme exposure to exogenous reductants is similar in both complexes, the apparent electrostatic environment near the cytochrome c heme edge is markedly different. In the electrostatic complex, a net positive charge is present, whereas in the covalent complex, an essentially neutral electrostatic charge is found. Intracomplex electron transfer within the two complexes is also different. For the covalent complex, electron transfer from ferrous cytochrome c to the ferryl peroxidase has a rate constant of 1560 s-1, which is invariant with respect to changes in the ionic strength. The rate constant for intracomplex electron transfer within the electrostatic complex is highly ionic strength dependent. At mu = 8 mM a value of 750 s-1 has been obtained [Hazzard, J. T., Poulos, T. L., & Tollin, G. (1987) Biochemistry 26, 2836-2848], whereas at mu = 30 mM the value is 3300 s-1. This ionic strength dependency for the electrostatic complex has been interpreted in terms of the rearrangement of the two proteins comprising the complex to a more favorable orientation for electron transfer. In the case of the covalent complex, such reorientation is apparently impeded.(ABSTRACT TRUNCATED AT 250 WORDS)     

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.