Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “electrostatics”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3Linked to original sources
PubMed1996

Roles of electrostatic interaction in proteins.

More than 60 years after the analyses by Linderstrom-Lang and Kirkwood of their hypothetical 'protein' structures, we have now a plethora of experimental evidence and computational estimates of the electrostatic forces in proteins, with very many protein 3D structures at atomic resolution. In the mean time, there were in the beginning, many arguments and suggestions about the roles of electrostatics, mainly from empirical findings and tendencies. A few experimental results indicated that the electrostatic contribution is of the order of several kcal/mol, which was theoretically difficult to reproduce correctly, because a large opposing reaction field should be subtracted from a large, direct Coulombic field. Although the importance of the reaction field was recognized even 70 years ago, appropriate applications to protein molecules were made only in this decade, with the development of numerical computation. Now, an electrostatic molecular surface is one of the most popular pictures in journals of structural biology, indicating that the electrostatic force is one of the important components contributing to molecular recognition, which is a major focus of current biology and biochemistry. The development of NMR techniques has made it possible to observe the individual ionizations of ionizable groups in a protein, in addition to the determination of the 3D structure. Since it does not require any additional probe, each charge state can report the very local and heterogeneous electrostatic potentials working in the protein, without disturbing the original field. From the pKa values, the contributions of electrostatic interactions, ion pairs, charge-dipole interactions, and hydrogen bonds to protein stability have been correctly evaluated. Protein engineering also provides much more information than that obtainable from the native proteins, as the residues concerned can now be easily substituted with other amino acid residues having electrostatically different characteristics. Those experimental results have revealed smaller contributions than previously expected, probably because we underestimated the reaction field effects. Especially, a single ion pair stabilizes a protein only slightly, although a cooperative salt-bridge network can contribute significantly to protein stability. Marginal stabilities of proteins arise from small difference between many factors with driving and opposing forces. In spite of the small contribution of each single electrostatic interaction to the protein stability, the sum of their actions works to maintain the specific 3D structure of the protein. The 'negative' roles of electrostatics, which might destabilize protein conformation, should be pointed out. Unpaired buried charges are energetically too expensive to exit in the hydrophobic core. Isolated hydrogen bond donors and acceptors also exert negative effects, but they are not as expensive as the unpaired buried charges, with costs of a few kcal/mol. Therefore, statistical analyses of protein 3D structures reveal only rare instances of isolated hydrogen bond donors and acceptors. This must be the main reason why alpha-helices and beta-sheets are only observed in protein cores as the backbone structures. Such secondary structures do not leave any backbone hydrogen donors or acceptors unpaired, because of their intrinsically regular packing. Otherwise, it might be very difficult to construct a backbone structure, in which all the backbone amide and carbonyl groups had their own hydrogen bond partners in the protein core. There are two theoretical approaches to protein electrostatics, the macroscopic or continuum model, and the microscopic or molecular model. As described in this article, the macroscopic model has inherent problems because the protein-solvent system is very hetergeneous from the physical point of view...

Actins↗
PubMed2002

Electrostatics in protein-protein docking.

A novel geometric-electrostatic docking algorithm is presented, which tests and quantifies the electrostatic complementarity of the molecular surfaces together with the shape complementarity. We represent each molecule to be docked as a grid of complex numbers, storing information regarding the shape of the molecule in the real part and information regarding the electrostatic character of the molecule in the imaginary part. The electrostatic descriptors are derived from the electrostatic potential of the molecule. Thus, the electrostatic character of the molecule is represented as patches of positive, neutral, or negative values. The potential for each molecule is calculated only once and stored as potential spheres adequate for exhaustive rotation/translation scans. The geometric-electrostatic docking algorithm is applied to 17 systems, starting form the structures of the unbound molecules. The results-in terms of the complementarity scores of the nearly correct solutions, their ranking in the lists of sorted solutions, and their statistical uniqueness-are compared with those of geometric docking, showing that the inclusion of electrostatic complementarity in docking is very important, in particular in docking of unbound structures. Based on our results, we formulate several "good electrostatic docking rules": The geometric-electrostatic docking procedure is more successful than geometric docking when the potential patches are large and when the potential extends away from the molecular surface and protrudes into the solvent. In contrast, geometric docking is recommended when the electrostatic potential around the molecules to be docked appears homogeneous, that is, with a similar sign all around the molecule.

Algorithms↗
PubMed1999

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.

Acetylcholinesterase↗
PubMed2005

PCE: web tools to compute protein continuum electrostatics.

PCE (protein continuum electrostatics) is an online service for protein electrostatic computations presently based on the MEAD (macroscopic electrostatics with atomic detail) package initially developed by D. Bashford [(2004) Front Biosci., 9, 1082-1099]. This computer method uses a macroscopic electrostatic model for the calculation of protein electrostatic properties, such as pK(a) values of titratable groups and electrostatic potentials. The MEAD package generates electrostatic energies via finite difference solution to the Poisson-Boltzmann equation. Users submit a PDB file and PCE returns potentials and pK(a) values as well as color (static or animated) figures displaying electrostatic potentials mapped on the molecular surface. This service is intended to facilitate electrostatics analyses of proteins and thereby broaden the accessibility to continuum electrostatics to the biological community. PCE can be accessed at http://bioserv.rpbs.jussieu.fr/PCE.

Hydrogen-Ion Concentration↗
PubMed1999

Electrostatic factors in DNA intercalation.

The factors that determine the binding of a chromophore between the base pairs in DNA intercalation complexes are dissected. The electrostatic potential in the intercalation plane is calculated using an accurate ab initio based distributed multipole electrostatic model for a range of intercalation sites, involving different sequences of base pairs and relative twist angles. There will be a significant electrostatic contribution to the binding energy for chromophores with a predominantly positive electrostatic potential, but this varies significantly with sequence, and somewhat with twist angle. The usefulness of these potential maps for understanding the binding of intercalators is explored by calculating the electrostatic binding energy for 9-aminoacridine, ethidium, and daunomycin in a variety of model binding sites. The electrostatic forces play a major role in the positioning of an intercalating 9-aminoacridine and a significant stabilizing role in the binding of ethidium in its sterically constrained position, but the intercalation of daunomycin is determined by the side-chain binding. Sequence preferences are likely to be determined by a complex and subtle mixture of effects, with electrostatics being just one component. The electrostatic binding energy is also unlikely to be a major determinant of the twist angle, as its variation with angle is modest for most intercalation sites. Overall, the electrostatic potential maps give guidance on how positively charged chromophores can be chemically adapted by heteroatomic substitution to optimise their binding.

Aminacrine↗
PubMed2001

Microflotation Suppression and Enhancement Caused by Particle/Bubble Electrostatic Interaction.

The processes of attachment and detachment of small or medium-sized particles to relatively large bubbles during microflotation are considered in terms of the heterocoagulation theory. Calculations are made for the conditions that the surface potentials are of similar sign and constant, that one of the surface potentials is small, that hydrophobic attraction is absent, and that there are no surface deformations. Under these conditions bubble-particle aggregates may form as a result of an electrostatic attraction which exceeds the repulsive van der Waals force at intermediate distances. Next to electrostatic and van der Waals forces, hydrodynamic and gravitational forces are considered. These forces may overcome the electrostatic repulsion at large distances and promote particle bubble attachment. Strong electrostatic attraction at small distances, arising at a large difference of the surface potentials of the bubble and the particle and of low electrolyte concentrations, can prevent subsequent detachment by hydrodynamic and gravitational forces. With increasing electrolyte concentration the electrostatic barrier increases and the attractive electrostatic force diminishes. As a result, a critical electrolyte concentration for microflotation exists. Above this concentration attachment may still occur but it is followed by detachment. At lower electrolyte concentrations the electrostatic attractive force prevents the detachment. The dependence of the critical electrolyte concentration on the values of the bubble and particle potentials and the Hamaker constant is calculated. The critical concentration does not depend on particle or bubble size if the absolute values of the total detachment force and the total pressing force coincide, which is the case for Stokes and potential flow. For every electrolyte concentration lower than the critical value there are two critical particle sizes that limit the flotation possibility. For small particle sizes attachment is impossible because the pressing force is smaller than the electrostatic barrier. For large particle sizes detachment cannot be prevented because the detachment force exceeds the maximum electrostatic attraction. A microflotation domain of intermediate particle sizes exists in which irreversible heterocoagulation occurs. Copyright 2001 Academic Press.

Journal Article↗
PubMed1993

Stabilization of protein crystals by electrostatic interactions as revealed by a numerical approach.

We developed a novel algorithm to solve numerically the Poisson-Boltzmann equations under a periodic boundary condition. By employing this algorithm to calculate the electrostatic potentials in two different types of protein crystals, a bovine pancreatic trypsin inhibitor (BPTI) orthorhombic crystal and a pig-insulin cubic crystal, the energy contributions of the electrostatic interactions to the crystals' stability were evaluated. At a high ionic strength, the condensed state of proteins in the crystal was stabilized electrostatically compared with that isolated in dilute solution because the attractive electrostatic interactions between neighboring protein molecules overcame the repulsive forces that originated from the same net charges of the equivalent protein molecules. On the other hand, at a low ionic strength the electrostatic interactions destabilized the crystalline state of both proteins, although a different dependence on the ionic strength was found between them. Here, the insulin crystal was more stable than the BPTI one because of the higher charge density in the BPTI crystal. In all of the solvent ionic strengths investigated, the attractive electrostatic interactions between charge pairs separated by less than 5 A on the respective protein molecules prominently stabilize the protein crystals. Therefore, two protein molecules in the crystals are oriented to compensate each other for their opposite charges on the surfaces. We also found a specific role for bound phosphate ions in the stabilization of the BPTI crystal, based on comparison of the electrostatic energies of the two crystals with and without the ions. By determining the contribution of each atomic charge in the crystals to the electrostatic energy, it was revealed that several electrostatic pairs specifically contributed to the crystal's stability. On the basis of our numerical calculation results, we propose a new method to design protein molecules that adopt stable crystals by replacing destabilizing residues with stabilizing ones and by introducing specific hydrogen bonds or salt bridges between neighboring protein molecules.

Algorithms↗
PubMed2002

Electrostatic free energy of the DNA double helix in counterion condensation theory.

Polyelectrolyte theory based on counterion condensation is extended from the standard line charge model to helical and double helical charge arrays. The number of condensed counterions turns out to be the same as for a line charge with charge density equal to the axial charge density of the helix. Also, the logarithmic salt dependence of the electrostatic free energy is the same in the range of lower salt concentration, so that the limiting laws remain unchanged. However, the internal free energy of the condensed layer of counterions and the overall electrostatic free energy depend on the helical parameters. At higher salt, the free energies of both single and double helix are negative, indicating electrostatic stabilization of the helical charge lattices due to the mixing entropy of the condensed counterions. Except at very low salt, the free energy of a single helix is higher than the free energy of a double helix with twice the charge density. With B-DNA parameters and single strands modeled as single helices, the predicted salt dependence of the free energy of transition from double helix to separated single strands has a maximum at approximately 0.2 M salt, close to the location in the laboratory of this well-known feature of the DNA strand separation transition. We also calculate the electrostatic free energy for the transition of the DNA double helix from the B to the A conformation. The B form is electrostatically stable over most of the salt range, but there is a spontaneous electrostatic transition to A near 1 M salt. The electrostatic free energy values are close to the experimental values of the overall (electrostatic plus non-electrostatic) transition free energies for A-philic base pair sequences. We are led to suggest that the experimentally observed B-to-A transition for A-philic sequences near 1 M salt in water is governed by the polyelectrolyte properties of these two conformations of the DNA double helix. The effect of ethanol, however, cannot be attributed to lowering of the bulk dielectric constant.

DNA↗
PubMed1994

Electrostatic interactions control the parallel and antiparallel orientation of alpha-helical chains in two-stranded alpha-helical coiled-coils.

The role of interchain electrostatic interactions in orientating alpha-helical chains to form two-stranded parallel and antiparallel coiled-coils has been investigated. Four disulfide-bridged coiled-coils were designed: parallel coiled-coils with interchain electrostatic attractions (P/A) and repulsions (P/R) and antiparallel coiled-coils with interchain electrostatic attractions (AP/A) and repulsions (AP/R). These coiled-coils were made by air oxidation of two 35-residue peptides with the appropriate heptad repeat (LaEbAcLdEeGfKg or LaAbEcLdKeGfEg) to give the desired interchain electrostatic interactions, and the appropriate position of the cysteine residue (C2 or C33) to give the desired chain orientation. The coiled-coils were characterized by circular dichroism spectroscopy, and their stabilities were assessed by guanidine hydrochloride and urea denaturations. The results indicated that the favored chain orientation, that is, the major disulfide-bridged product formed under benign conditions, was the one that provides interchain electrostatic attractions between oppositely-charged amino acid residues in the e-g' and g-e' positions of the parallel coiled-coil and the g-g' and e-e' positions in the antiparallel coiled-coil. When the electrostatic interactions were similar, the antiparallel coiled-coils were more stable than the parallel coiled-coils. However, the overall stability of the coiled-coils was either increased by interchain electrostatic attractions or decreased by interchain electrostatic repulsions, as determined by urea denaturation. Thus, the order of overall stability of these coiled-coils was AP/A > P/A > AP/R > P/R. This study demonstrates the importance of interchain electrostatic interactions in determining the parallel or antiparallel orientation of alpha-helical chains in two-stranded coiled-coils.

Amino Acid Sequence↗
PubMed2003

Conservation of electrostatic properties within enzyme families and superfamilies.

Electrostatic interactions play a key role in enzyme catalytic function. At long range, electrostatics steer the incoming ligand/substrate to the active site, and at short distances, electrostatics provide the specific local interactions for catalysis. In cases in which electrostatics determine enzyme function, orthologs should share the electrostatic properties to maintain function. Often, electrostatic potential maps are employed to depict how conserved surface electrostatics preserve function. We expand on previous efforts to explain conservation of function, using novel electrostatic sequence and structure analyses of four enzyme families and one enzyme superfamily. We show that the spatial charge distribution is conserved within each family and superfamily. Conversely, phylogenetic analysis of key electrostatic residues provide the evolutionary origins of functionality.

Amino Acid Sequence↗
PubMed1991

Role of electrostatic interaction in the stability of the hexamer of constitutive phycocyanin from Fremyella diplosiphon.

The role of electrostatic interactions between the charges carried by the titratable groups in the different aggregates of constitutive phycocyanin from Fremyella diplosiphon has been studied by using a simple theoretical approach based on the modified Tanford-Kirkwood model. The electrostatic potential has also been calculated by means of a numerical solution of the linearized Poisson-Boltzmann equations using the finite-differences technique. The pH dependence of the electrostatic contribution to free energy suggests an electrostatic stabilization of the alpha- and beta-subunits as well as of the (alpha beta)-monomer over a broad pH interval. The charge distributions in the individual alpha- and beta-subunits produce electrostatic complementarity and promote the assembly of the subunits to the (alpha beta)-monomer, as well as of the monomers to the larger trimeric and hexameric aggregates. Trimer-trimer electrostatic interactions exhibit strong pH dependence, predicting an association/dissociation equilibrium with a midpoint at pH 6. The electrostatic trimer-trimer interactions correspond to the steric fit, suggesting that electrostatic interactions may initially help to orient the trimers during aggregation. The distribution of the electrostatic potential of the monomers and of the higher aggregates suggests that it plays an important role also in phycocyanin-linker protein binding.

Computer Simulation↗
PubMed1996

Catalytic subunit of cAMP-dependent protein kinase: electrostatic features and peptide recognition.

The electrostatic field was calculated for the mammalian cAMP-dependent protein kinase (PKA) catalytic subunit (C-subunit) complexed with a 20-residue peptide from a heat stable protein kinase inhibitor (PKI: 5-24). The electrostatic field was also calculated for the C-subunit complexed with a modeled heptapeptide substrate that has been used extensively in structure/function studies for the C-subunit. Perturbations in the electrostatic free energy were calculated when single ionizable active site residues were mutated to alanine. These perturbations in electrostatic free energy were correlated to changes in the binding energy measured in a charge-to-alanine scan of the homologous yeast C-subunit by M. J. Zoller and C. S. Gibbs [(1991) Journal of Biological Chemistry, Vol. 266, pp. 8923-8931; C. S. Gibbs and M. J. Zoller (1991) Biochemistry, Vol. 30, p. 22]. This analysis indicated that the substrate binding parameters primarily depend on electrostatic interactions between a substrate or inhibitor and the C-subunit. Amino acid replacements that led to large perturbations in the electrostatic field are listed in the text. pKa shifts were also calculated for the substrate's phosphate accepting atom, the serine hydroxyl oxygen, when the active site ionizable residues were changed to structurally similar uncharged amino acids. The theoretical mutation of three active site residues caused large shifts in this parameter: E91Q, D166N, and D184N. The calculated pKa shifts for these mutants indicate that the rate of phosphotransfer should be markedly reduced in these cases. This prediction has been experimentally confirmed for the D166N mutant. The correlation between calculated electrostatic free energy changes and measured binding energy, and pKa shifts with phosphotransfer for C-subunit mutants were within experimental error of the measurements. The calculations of electrostatic energy and delta pKa have identified previously unconsidered active site residues in the mammalian C-subunit that contribute to binding energy and phosphotransfer.

Animals↗
PubMed2004

Calculation of electrostatic interaction energies in molecular dimers from atomic multipole moments obtained by different methods of electron density partitioning.

Accurate and fast evaluation of electrostatic interactions in molecular systems is still one of the most challenging tasks in the rapidly advancing field of macromolecular chemistry, including molecular recognition, protein modeling and drug design. One of the most convenient and accurate approaches is based on a Buckingham-type approximation that uses the multipole moment expansion of molecular/atomic charge distributions. In the mid-1980s it was shown that the pseudoatom model commonly used in experimental X-ray charge density studies can be easily combined with the Buckingham-type approach for calculation of electrostatic interactions, plus atom-atom potentials for evaluation of the total interaction energies in molecular systems. While many such studies have been reported, little attention has been paid to the accuracy of evaluation of the purely electrostatic interactions as errors may be absorbed in the semiempirical atom-atom potentials that have to be used to account for exchange repulsion and dispersion forces. This study is aimed at the evaluation of the accuracy of the calculation of electrostatic interaction energies with the Buckingham approach. To eliminate experimental uncertainties, the atomic moments are based on theoretical single-molecule electron densities calculated at various levels of theory. The electrostatic interaction energies for a total of 11 dimers of alpha-glycine, N-acetylglycine and L-(+)-lactic acid structures calculated according to Buckingham with pseudoatom, stockholder and atoms-in-molecules moments are compared with those evaluated with the Morokuma-Ziegler energy decomposition scheme. For alpha-glycine a comparison with direct "pixel-by-pixel" integration method, recently developed Gavezzotti, is also made. It is found that the theoretical pseudoatom moments combined with the Buckingham model do predict the correct relative electrostatic interactions energies, although the absolute interaction energies are underestimated in some cases. The good agreement between electrostatic interaction energies computed with Morokuma-Ziegler partitioning, Gavezzotti's method, and the Buckingham approach with atoms-in-molecules moments demonstrates that reliable and accurate evaluation of electrostatic interactions in molecular systems of considerable complexity is now feasible.

Algorithms↗
PubMed2003

Electrostatic properties in the catalytic site of papain: A possible regulatory mechanism for the reactivity of the ion pair.

We present an analysis of the electrostatic properties in the catalytic site of papain (EC 3.4.22.2), an archetype enzyme of the C1 cysteine proteinase family, and we investigate their possible role in the formation, stabilization and regulation of the Cys25((-))...His159((+)) catalytic ion pair. The electrostatic properties were computed using a reassociation method based in multicentered multipolar expansions obtained from ab initio quantum calculations of overlapping protein fragments. Solvent effects were introduced by coupling the use of multicentered multipolar expansions to two continuum boundary element methods to solve the Poisson and the linearized Poisson-Boltzmann equations. The electrostatic profile found in the proton transfer region of papain showed that this enzyme has a well-defined electrostatic environment to favor the formation and stabilization of the catalytic ion pair. The papain catalytic site electrostatic profile can be considered as an electrostatic fingerprint of the papain family with the following characteristics: (i) the presence of a net electric field highly aligned in the (Cys25)-SG-->(His159)-ND1 direction; (ii) the electrostatic profile has a saddle-point character; (iii) it is basically a local environmental effect. Furthermore, our analysis describes a possible regulatory mechanism (the E(SG-->ND1) attenuation effect) controlling the ion pair reactivity and permits to infer the Asp57 acidic residue as the most probable candidate to act as the electrostatic modulator.

Amino Acids↗
PubMed1995

The matching of electrostatic extrema: a useful method in drug design? A study of phosphodiesterase III inhibitors.

Ligands which bind to a specific protein binding site are often expected to have a similar electrostatic environment which complements that of the binding site. One method of assessing molecular electrostatic similarity is to examine the possible overlay of the maxima and minima in the electrostatic potential outside the molecules and thereby match the regions where strong electrostatic interactions, including hydrogen bonds, with the residues of the binding site may be possible. This approach is validated with accurate calculations of the electrostatic potential, derived from a distributed multiple analysis of an ab initio charge density of the molecule, so that the effects of lone pair and pi-electron density are correctly included. We have applied this method to the phosphodiesterase (PDE) III substrate adenosine-3',5'-cyclic monophosphate (cAMP) and a range of nonspecific and specific PDE III inhibitors. Despite the structural variation between cAMP and the inhibitors, it is possible to match three or four extrema to produce relative orientations in which the inhibitors are sufficiently sterically and electrostatically similar to the natural substrate to account for their affinity for PDE III. This matching of extrema is more apparent using the accurate electrostatic models than it was when this approach was first applied, using semiempirical point charge models. These results reinforce the hypothesis of electrostatic similarity and give weight to the technique of extrema matching as a useful tool in drug design.

3',5'-Cyclic-AMP Phosphodiesterases↗
PubMed2002

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.

Antigen-Antibody Reactions↗
PubMed2003

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.

Binding Sites↗
PubMed2001

Optimization of binding electrostatics: charge complementarity in the barnase-barstar protein complex.

Theoretical and experimental studies have shown that the large desolvation penalty required for polar and charged groups frequently precludes their involvement in electrostatic interactions that contribute strongly to net stability in the folding or binding of proteins in aqueous solution near room temperature. We have previously developed a theoretical framework for computing optimized electrostatic interactions and illustrated use of the algorithm with simplified geometries. Given a receptor and model assumptions, the method computes the ligand-charge distribution that provides the most favorable balance of desolvation and interaction effects on binding. In this paper the method has been extended to treat complexes using actual molecular shapes. The barnase-barstar protein complex was investigated with barnase treated as a target receptor. The atomic point charges of barstar were varied to optimize the electrostatic binding free energy. Barnase and natural barstar form a tight complex (K(d) approximately 10(-14) M) with many charged and polar groups near the interface that make this a particularly relevant system for investigating the role of electrostatic effects on binding. The results show that sets of barstar charges (resulting from optimization with different constraints) can be found that give rise to relatively large predicted improvements in electrostatic binding free energy. Principles for enhancing the effect of electrostatic interactions in molecular binding in aqueous environments are discussed in light of the optima. Our findings suggest that, in general, the enhancements in electrostatic binding free energy resulting from modification of polar and charged groups can be substantial. Moreover, a recently proposed definition of electrostatic complementarity is shown to be a useful tool for examining binding interfaces. Finally, calculational results suggest that wild-type barstar is closer to being affinity optimized than is barnase for their mutual binding, consistent with the known roles of these proteins.

Algorithms↗