Salt bridge interactions: Stability of the ionic and neutral complexes in the gas phase,in solution,and in proteins

A theoretical study on the stability of the salt bridges in the gas phase, in solution, and in the interior of proteins is presented. The study is mainly focused on the interaction between acetate and methylguanidinium ions, which were used as model compounds for the salt bridge between Asp (Glu) and Arg. Two different solvents (water and chloroform) were used to analyze the effect of varying the dielectric constant of the surrounding media on the salt bridge interaction. Calculations in protein environments were performed by using a set of selected protein crystal structures. In all cases attention was paid to the difference in stability between the ion pair and neutral hydrogen-bonded forms. Comparison of the results determined in the gas phase and in solution allows us to stress the large influence of the environment on the binding process, as well as on the relative stability between the ionic and neutral complexes. The high anisotropy of proteins and the local microenvironment in the interior of proteins make a decisive contribution in modulating the energetics of the salt bridge. In general, the formation of salt bridges in proteins is not particularly favored, with the ion pair structure being preferred over the interaction between neutral species. Proteins 32:67–79, 1998. © 1998 Wiley-Liss, Inc.     

Electrostatic interactions in the association of proteins: an analysis of the thrombin-hirudin complex.

The role of electrostatic interactions in stabilization of the thrombin-hirudin complex has been investigated by means of two macroscopic approaches: the modified Tanford-Kirkwood model and the finite-difference method for numerical solution of the Poisson-Boltzmann equations. The electrostatic potentials around the thrombin and hirudin molecules were asymmetric and complementary, and it is suggested that these fields influence the initial orientation in the process of the complex formation. The change of the electrostatic binding energy due to mutation of acidic residues in hirudin has been calculated and compared with experimentally determined changes in binding energy. In general, the change in electrostatic binding energy for a particular mutation calculated by the modified Tanford-Kirkwood approach agreed well with the experimentally observed change. The finite-difference approach tended to overestimate changes in binding energy when the mutated residues were involved in short-range electrostatic interactions. Decreases in binding energy caused by mutations of amino acids that do not make any direct ionic interactions (e.g., Glu 61 and Glu 62 of hirudin) can be explained in terms of the interaction of these charges with the positive electrostatic potential of thrombin. Differences between the calculated and observed changes in binding energy are discussed in terms of the crystal structure of the thrombin-hirudin complex.     

Electrostatic stabilization in four-helix bundle proteins.

Charge substitutions generated by site-directed mutagenesis at the termini of adjacent anti-parallel alpha-helices in a four-helix bundle protein were used to determine a precise value for the contribution of indirect charge-charge interactions to overall protein stability, and to simulate the electrostatic effects of alpha-helix macrodipoles. Thermodynamic double mutant cycles were constructed to measure the interaction energy between such charges on adjacent anti-parallel helices in the four-helix bundle cytochrome b562 from Escherichia coli. Previously, theoretical calculations of helix macrodipole interactions using modeled four-helix bundle proteins have predicted values ranging over an order of magnitude from 0.2 to 2.5 kcal/mol. Our system represents the first experimental evidence for electrostatic interactions such as those between partial charges due to helix macrodipole charges. At the positions mutated, we have measured a favorable interaction energy of 0.6 kcal/mol between opposite charges simulating an anti-parallel helix pair. Pairs of negative or positive charges simulating a parallel orientation of helices produce an unfavorable interaction of similar magnitude. The interaction energies show a strong dependence upon ionic strength, consistent with an electrostatic effect. Indirect electrostatic contacts do appear to confer a limited stabilization upon the association of anti-parallel packing of helices, favoring this orientation by as much as 1 kcal/mol at 20 mM K phosphate.     

Optimization of the electrostatic interactions between ionized groups and peptide dipoles in proteins.

The three-dimensional optimization of the electrostatic interactions between the charged amino acid residues and the peptide partial charges was studied by Monte-Carlo simulations on a set of 127 nonhomologous protein structures with known atomic coordinates. It was shown that this type of interaction is very well optimized for all proteins in the data set, which suggests that they are a valuable driving force, at least for the native side-chain conformations. Similar to the optimization of the charge-charge interactions (Spassov VZ, Karshikoff AD, Ladenstein R, 1995, Protein Sci 4:1516-1527), the optimization effect was found more pronounced for enzymes than for proteins without enzymatic function. The asymmetry in the interactions of acidic and basic groups with the peptide dipoles was analyzed and a hypothesis was proposed that the properties of peptide dipoles are a factor contributing to the natural selection of the basic amino acids in the chemical composition of proteins.     

The optimization of protein-solvent interactions: thermostability and the role of hydrophobic and electrostatic interactions.

Protein-solvent interactions were analyzed using an optimization parameter based on the ratio of the solvent-accessible area in the native and the unfolded protein structure. The calculations were performed for a set of 183 nonhomologous proteins with known three-dimensional structure available in the Protein Data Bank. The dependence of the total solvent-accessible surface area on the protein molecular mass was analyzed. It was shown that there is no difference between the monomeric and oligomeric proteins with respect to the solvent-accessible area. The results also suggested that for proteins with molecular mass above some critical mass, which is about 28 kDa, a formation of domain structure or subunit aggregation into oligomers is preferred rather than a further enlargement of a single domain structure. An analysis of the optimization of both protein-solvent and charge-charge interactions was performed for 14 proteins from thermophilic organisms. The comparison of the optimization parameters calculated for proteins from thermophiles and mesophiles showed that the former are generally characterized by a high degree of optimization of the hydrophobic interactions or, in cases where the optimization of the hydrophobic interactions is not sufficiently high, by highly optimized charge-charge interactions.     

Simple theoretical models for biochemical systems,with applications to DNA

A set oflogically connected models, to study chemical systems ofbiological interest, is presented. The sequence in the set is dictated by a progressive reduction of details with a corresponding enlargement of the field of application. The exposition starts with models suitable for interactions among a finite number of molecules, passes then to models considering also solvent effects and ends with models specialized for DNA containing systems.     

Hydrophobic vs coulombic interactions in the binding of steroidal polyamines to DNA

Seven new steroidal polyamines derived from bile acids, either lithocholic or deoxycholic acid, have been studied as DNA-binding agents using four complimentary methods: an ethidium displacement assay, observed changes in the thermal denaturation of poly[d(AT)], effects on hyperchromicity of DNA, and circular dichroism. In addition, modelling studies were conducted to examine the electrostatic surface potential of the polycations. The results point to a key role for a large hydrophobic surface area on the steroid in addition to the Coulombic attraction by ammonium and guanidinium groups on the steroid interacting with the polyphosphate backbone.     

Brownian dynamics of interactions between aldolase mutants and F-actin

Previous Brownian dynamics (BD) simulations (Ouporov IG, Knull HR and Thomasson KA 1999. Biophys. J. 76: 17-27) of complex formation between rabbit aldolase and F-actin have identified three lysine residues (K288, K293 and K341) on aldolase and acidic residues (DEDE) at the N-terminus of actin as important to binding. BD simulations of computer models of aldolase mutants with any of these lysine residues replaced by alanine show reduced binding energy; the greatest effect of a single substitution is for K341A, and replacement of all three lysines greatly reduces binding. BD simulations of wild-type rabbit aldolase vs altered F-actin show that binding is decreased if any one of the four N-terminal acidic residues is replaced by alanine and binding is greatly reduced if three or more of the N-terminal acidic residues are replaced; none of the four actin residues appear more critical for binding than the others.     

A probabilistic method for calculating the energy of hydrophobic interactions

A theoretical approach is developed for estimating the hydrophobic interaction energy on the molecular level. The underlying idea of the model is the fundamental relationship between the probability of the appearance of a solvent-excluded volume in fluid water as a result of density fluctuations, on one hand, and the free energy of primitive hydrophobic solutes, on the other. This probability is estimated basing on an information theory model in combination with experimental data on the distribution density and radial distribution function. The free energy of hydrophobic interactions is determined for a complex that consists of several hard spheres. The critical number of particles in the complex is estimated.     

Intrinsically unstructured proteins by design—electrostatic interactions can control binding,folding, and function of a helix‐loop‐helix heterodimer

Intrinsically disordered proteins that exist as unordered monomeric structures in aqueous solution at pH 7 but fold into four‐helix bundles upon binding to recognized polypeptide targets have been designed. NMR and CD spectra of the monomeric polypeptides show the hallmarks of unordered structures, whereas in the bound state they are highly helical. Analytical ultracentrifugation data shows that the polypeptides bind to their targets to form exclusively heterodimers at neutral pH. To demonstrate the relationship between binding, folding, and function, a catalytic site for ester hydrolysis was introduced into an unordered and largely inactive monomer, but that was structured and catalytically active in the presence of a specific polypeptide target. Electrostatic interactions between surface‐exposed residues inhibited the binding and folding of the monomers at pH 7. Charge–charge repulsion between ionizable amino acids was thus found to be sufficient to disrupt binding between polypeptide chains despite their inherent propensities for structure formation and may be involved in the folding and function of inherently disordered proteins in biology. Copyright © 2013 European Peptide Society and John Wiley & Sons, Ltd.     

Intra-Helical Salt Bridge Contribution to Membrane Protein Insertion

Salt bridges between negatively (D, E) and positively charged (K, R, H) amino acids play an important role in protein stabilization. This has a more prevalent effect in membrane proteins where polar amino acids are exposed to a hydrophobic environment. In transmembrane (TM) helices the presence of charged residues can hinder the insertion of the helices into the membrane. It is possible that the formation of salt bridges could decrease the cost of membrane integration. However, the presence of intra-helical salt bridges in TM domains and their effect on insertion has not been properly studied yet. In this work, we show that potentially salt-bridge forming pairs are statistically over-represented in TM-helices. We then selected some candidates to experimentally determine the contribution of these electrostatic interactions to the translocon-assisted membrane insertion process. Using both in vitro and whole cell systems, we confirm the presence of intra-helical salt bridges in TM segments during biogenesis and determined that they contribute ～0.5 kcal/mol to the apparent free energy of membrane insertion (ΔG_app). Our observations suggest that salt bridge interactions can be stabilized during translocon-mediated insertion and thus could be relevant to consider for the future development of membrane protein prediction software.     

Electrostatic contributions to protein-protein interactions: fast energetic filters for docking and their physical basis

The methods of continuum electrostatics are used to calculate the binding free energies of a set of protein-protein complexes including experimentally determined structures as well as other orientations generated by a fast docking algorithm. In the native structures, charged groups that are deeply buried were often found to favor complex formation (relative to isosteric nonpolar groups), whereas in nonnative complexes generated by a geometric docking algorithm, they were equally likely to be stabilizing as destabilizing. These observations were used to design a new filter for screening docked conformations that was applied, in conjunction with a number of geometric filters that assess shape complementarity, to 15 antibody-antigen complexes and 14 enzyme-inhibitor complexes. For the bound docking problem, which is the major focus of this paper, native and near-native solutions were ranked first or second in all but two enzyme-inhibitor complexes. Less success was encountered for antibody-antigen complexes, but in all cases studied, the more complete free energy evaluation was able to identify native and near-native structures. A filter based on the enrichment of tyrosines and tryptophans in antibody binding sites was applied to the antibody-antigen complexes and resulted in a native and near-native solution being ranked first and second in all cases. A clear improvement over previously reported results was obtained for the unbound antibody-antigen examples as well. The algorithm and various filters used in this work are quite efficient and are able to reduce the number of plausible docking orientations to a size small enough so that a final more complete free energy evaluation on the reduced set becomes computationally feasible.     

The role of electrostatic interaction in the molecular recognition of selective agonists to metabotropic glutamate receptors

The influence of electrostatic interactions in determining selectivity for individual subtypes of metabotropic glutamate receptors (mGluRs) is evaluated for a small set of agonists by using the program Delphi and the information thus obtained is compared with docking experiments carried out with AutoDock. The evaluation of the electrostatic component of the free energy of binding for L-Glu, L-AP4, or S-PPG to mGluR1, mGluR2, and mGluR4 subtypes allowed for the detection of subtle differences in the electronic properties of the three subtypes, differences that can account for the observed agonist selectivity.     

Electrostatic interactions with histone tails may bend linker DNA in chromatin

Is linker DNA bent in the 30‐nm chromatin fiber at physiological conditions? We show here that electrostatic interactions between linker DNA and histone tails including salt condensation and release may bend linker DNA, thus affecting the higher order organization of chromatin. © 2005 Wiley Periodicals, Inc. Biopolymers 81: 20–28, 2006 This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Calculation of ligand-nucleic acid binding free energies with the generalized-born model in DOCK

The calculation of ligand-nucleic acid binding free energies is investigated by including solvation effects computed with the generalized-Born model. Modifications of the solvation module in DOCK, including introduction of all-atom parameters and revision of coefficients in front of different terms, are shown to improve calculations involving nucleic acids. This computing scheme is capable of calculating binding energies, with reasonable accuracy, for a wide variety of DNA-ligand complexes, RNA-ligand complexes, and even for the formation of double-stranded DNA. This implementation of GB/SA is also shown to be capable of discriminating strong ligands from poor ligands for a series of RNA aptamers without sacrificing the high efficiency of the previous implementation. These results validate this approach to screening large databases against nucleic acid targets.     

Use of protein charge ladders to study electrostatic interactions during protein ultrafiltration

Recent studies have demonstrated the importance of electrostatic interactions in membrane systems, but there is still controversy about the underlying phenomena. Protein charge ladders, consisting of a set of chemical derivatives of a given protein that differ by single charge groups, were used to quantify the electrostatic interactions during protein ultrafiltration. Myoglobin charge ladders were generated by acylation, with the different derivatives analyzed simultaneously by capillary electrophoresis. Filtration experiments were performed using polyethersulfone and composite regenerated cellulose membranes, with the membrane charge determined from the streaming potential. As expected, the rejection increased as the protein became more heavily charged due to the increase in electrostatic repulsion. However, the transmission of the weakly charged myoglobin species increased dramatically at very low ionic strength. This increase in transmission was attributed to a shift in pH within the pore caused by hydrogen ion partitioning into the charged membrane. The sieving data were in good agreement with theoretical calculations accounting for the effects of this pH shift on the electrostatic interactions.     

The influence of molecular variation on protein interactions

We investigated the influence of solvation forces on protein-protein interactions for two forms of lysozyme: hen egg white (HEWL) and turkey egg white (TEWL). Turkey egg white has more surface exposed hydrophobic residues than HEWL and the protein-protein interactions of TEWL are shown to be more attractive than those of HEWL, for the conditions studied. The importance of including a solvation term in the potential of mean force model, to account for molecular variation in protein surface characteristics, is highlighted. We also show that the magnitude of this solvation term can be estimated using readily available data.     

Energetics of cell-cell and cell-biopolymer interactions

The energy vs distance balance of cell suspensions (in the presence and in the absence of extracellular biopolymer solutions) is studied, not only in the light of the classical Derjaguin-Landau-Verwey-Over-beek (DLVO) theory (which considered just the electrostatic (EL) and Lifshitz-van der Waals (LW) interactions), but also by taking electron-acceptor/electron-donor, or Lewis acid-base (AB) and osmotic (OS) interactions into account. Since cell surfaces, as well as many biopolymers tend to have strong monopolar electron-donor properties, they are able to engage in a strong mutual AB repulsion when immersed in a polar liquid such as water. The effects of that repulsion have been observed earlier in the guise of hydration pressure. The AB repulsion is, at close range, typically one or two orders of magnitude stronger than the EL repulsion, but its rate of decay is much steeper. In most cases, AB interactions are quantitatively the dominant factor in cell stability (when repulsive) and in “hydrophobic interactions” (when attractive). OS interactions exerted by extracellularly dissolved biopolymers are weak, but their rate of decay is very gradual, so OS repulsions engendered by biopolymer solutions may be of importance in certain long-range interactions. OS interactions exerted by biopolymers attached to cells or particles (e.g., by glycocalix glycoproteins), are very short-ranged and usually are negligibly small in comparison with the other interaction forces, in aqueous media.