Salt bridge as a gatekeeper against partial unfolding

Salt bridges are frequently observed in protein structures. Because the energetic contribution of salt bridges is strongly dependent on the environmental context, salt bridges are believed to contribute to the structural specificity rather than the stability. To test the role of salt bridges in enhancing structural specificity, we investigated the contribution of a salt bridge to the energetics of native‐state partial unfolding in a cysteine‐free version of Escherichia coli ribonuclease H (RNase H*). Thermolysin cleaves a protruding loop of RNase H^* through transient partial unfolding under native conditions. Lys86 and Asp108 in RNase H^* form a partially buried salt bridge that tethers the protruding loop. Investigation of the global stability of K86Q/D108N RNase H^* showed that the salt bridge does not significantly contribute to the global stability. However, K86Q/D108N RNase H^* is greatly more susceptible to proteolysis by thermolysin than wild‐type RNase H^* is. The free energy for partial unfolding determined by native‐state proteolysis indicates that the salt bridge significantly increases the energy for partial unfolding by destabilizing the partially unfolded form. Double mutant cycles with single and double mutations of the salt bridge suggest that the partially unfolded form is destabilized due to a significant decrease in the interaction energy between Lys86 and Asp108 upon partial unfolding. This study demonstrates that, even in the case that a salt bridge does not contribute to the global stability, the salt bridge may function as a gatekeeper against partial unfolding that disturbs the optimal geometry of the salt bridge.     

Contribution of surface salt bridges to protein stability: guidelines for protein engineering

The small globular protein, ubiquitin, contains a pair of oppositely charged residues, K11 and E34, that according to the three-dimensional structure are located on the surface of this protein with a spatial orientation characteristic of a salt bridge. We investigated the strength of this salt bridge and its contribution to the global stability of the ubiquitin molecule. Using the "double mutant cycle" analysis, the strength of the pairwise interactions between K11 and E34 was estimated to be favorable by 3.6kJ/mol. Further, the salt bridge of the reverse orientation, i.e. E11/K34, can be formed and is found to have a strength (3.8kJ/mol) similar to that of the K11/E34 pair. However, the global stability of the K11/E34 variant of ubiquitin is 2.2kJ/mol higher than that of the E11/K34 variant. The difference in the contribution of the opposing salt bridge orientations to the overall stability of the ubiquitin molecule is attributed to the difference in the charge-charge interactions between residues forming the salt bridge and the rest of the ionizable groups in this protein. On the basis of these results, we concluded that surface salt bridges are stabilizing, but their contribution to the overall protein stability is strongly context-dependent, with charge-charge interactions being the largest determinant. Analysis of 16 salt bridges from six different proteins, for which detailed experimental data on energetics have been reported, support the conclusions made from the analysis of the salt bridge in ubiquitin. Implications of these findings for engineering proteins with enhanced thermostability are discussed.     

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.     

Surface salt bridges modulate DNA wrapping by the type II DNA-binding protein TF1

The histone-like protein HU is involved in compaction of the bacterial genome. Up to 37 bp of DNA may be wrapped about some HU homologues in a process that has been proposed to depend on a linked disruption of surface salt bridges that liberates cationic side chains for interaction with the DNA. Despite significant sequence conservation between HU homologues, binding sites from 9 to 37 bp have been reported. TF1, an HU homologue that is encoded by Bacillus subtilis bacteriophage SPO1, has nM affinity for 37 bp preferred sites in DNA with 5-hydroxymethyluracil (hmU) in place of thymine. On the basis of electrophoretic mobility shift assays, we show that TF1-DNA complex formation is associated with a net release of only approximately 0.5 cations. The structure of TF1 suggests that Asp13 can form a dehydrated surface salt bridge with Lys23; substitution of Asp13 with Ala increases the net release of cations to approximately 1. These data are consistent with complex formation linked to disruption of surface salt bridges. Substitution of Glu90 with Ala, which would expose Lys87 predicted to contact DNA immediately distal to a proline-mediated DNA kink, causes an increase in affinity and an abrogation of the preference for hmU-containing DNA. We propose that hmU preference is due to finely tuned interactions at the sites of kinking that expose a differential flexibility of hmU- and T-containing DNA. Our data further suggest that the difference in binding site size for HU homologues is based on a differential ability to stabilize the DNA kinks.     

Surface salt bridges stabilize the GCN4 leucine zipper.

We present a study of the role of salt bridges in stabilizing a simplified tertiary structural motif, the coiled-coil. Changes in GCN4 sequence have been engineered that introduce trial patterns of single and multiple salt bridges at solvent exposed sites. At the same sites, a set of alanine mutants was generated to provide a reference for thermodynamic analysis of the salt bridges. Introduction of three alanines stabilizes the dimer by 1.1 kcal/mol relative to the wild-type. An arrangement corresponding to a complex type of salt bridge involving three groups stabilizes the dimer by 1.7 kcal/ mol, an apparent elevation of the melting temperature relative to wild type of about 22 degrees C. While identifying local from nonlocal contributions to protein stability is difficult, stabilizing interactions can be identified by use of cycles. Introduction of alanines for side chains of lower helix propensity and complex salt bridges both stabilize the coiled-coil, so that combining the two should yield melting temperatures substantially higher than the starting species, approaching those of thermophilic sequences.     

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.     

Do salt bridges stabilize proteins? A continuum electrostatic analysis

The electrostatic contribution to the free energy of folding was calculated for 21 salt bridges in 9 protein X-ray crystal structures using a continuum electrostatic approach with the DELPHI computer-program package. The majority (17) were found to be electrostatically destabilizing; the average free energy change, which is analogous to mutation of salt bridging side chains to hydrophobic isosteres, was calculated to be 3.5 kcal/mol. This is fundamentally different from stability measurements using pKa shifts, which effectively measure the strength of a salt bridge relative to 1 or more charged hydrogen bonds. The calculated effect was due to a large, unfavorable desolvation contribution that was not fully compensated by favorable interactions within the salt bridge and between salt-bridge partners and other polar and charged groups in the folded protein. Some of the salt bridges were studied in further detail to determine the effect of the choice of values for atomic radii, internal protein dielectric constant, and ionic strength used in the calculations. Increased ionic strength resulted in little or no change in calculated stability for 3 of 4 salt bridges over a range of 0.1-0.9 M. The results suggest that mutation of salt bridges, particularly those that are buried, to "hydrophobic bridges" (that pack at least as well as wild type) can result in proteins with increased stability. Due to the large penalty for burying uncompensated ionizable groups, salt bridges could help to limit the number of low free energy conformations of a molecule or complex and thus play a role in determining specificity (i.e., the uniqueness of a protein fold or protein-ligand binding geometry).     

Thermodynamic stability of a cold-adapted protein, type III antifreeze protein, and energetic contribution of salt bridges

A thermodynamic analysis of a cold-adapted protein, type III anti-freeze protein (AFP), was carried out. The results indicate that the folding equilibrium of type III AFP is a reversible, unimolecular, two-state process with no populated intermediates. Compared to most mesophilic proteins whose folding is two-state, the psychrophilic type III AFP has a much lower thermodynamic stability at 25 degrees C, approximately 3 kcal/mol, and presents a remarkably downshifted stability-temperature curve, reaching a maximum of 5 kcal/mol around 0 degrees C. Type III AFPs contain few and non-optimally distributed surface charges relative to their mesophilic homologs, the C-terminal domains of sialic acid synthases. We used thermodynamic double mutant cycles to evaluate the energetic role of every surface salt bridge in type III AFP. Two isolated salt bridges provided no contribution to stability, while the Asp36-Arg39 salt bridge, involved in a salt bridge network with the C-terminal carboxylate, had a substantial contribution (approximately 1 kcal/mol). However, this contribution was more than counteracted by the destabilizing effect of the Asp36 carboxylate itself, whose removal led to a net 30% increase in stability at 25 degrees C. This study suggests that type III AFPs may have evolved for a minimally acceptable stability at the restricted, low temperature range (around 0 degrees C) at which AFPs must function. In addition, it indicates that salt bridge networks are used in nature also for the stability of psychrophilic proteins, and has led to a type III AFP variant of increased stability that could be used for biotechnological purposes.     

Differences in electrostatic properties at antibody-antigen binding sites: implications for specificity and cross-reactivity

Antibodies HyHEL8, HyHEL10, and HyHEL26 (HH8, HH10, and HH26, respectively) recognize highly overlapping epitopes on hen egg-white lysozyme (HEL) with similar affinities, but with different specificities. HH8 binding to HEL is least sensitive toward mutations in the epitope and thus is most cross-reactive, HH26 is most sensitive, whereas the sensitivity of HH10 lies in between HH8 and HH26. Here we have investigated intra- and intermolecular interactions in three antibody-protein complexes: theoretical models of HH8-HEL and HH26-HEL complexes, and the x-ray crystal structure of HH10-HEL complex. Our results show that HH8-HEL has the lowest number and HH26-HEL has the highest number of intra- and intermolecular hydrogen bonds. The number of salt bridges is lowest in HH8-HEL and highest in HH26-HEL. The binding site salt bridges in HH8-HEL are not networked, and are weak, whereas, in HH26-HEL, an intramolecular salt-bridge triad at the binding site is networked to an intermolecular triad to form a pentad. The pentad and each salt bridge of this pentad are exceptionally stabilizing. The number of binding-site salt bridges and their strengths are intermediate in HH10-HEL, with an intramolecular triad. Our further calculations show that the electrostatic component contributes the most to binding energy of HH26-HEL, whereas the hydrophobic component contributes the most in the case of HH8-HEL. A "hot-spot" epitope residue Lys-97 forms an intermolecular salt bridge in HH8-HEL, and participates in the intermolecular pentad in the HH26-HEL complex. Mutant modeling and surface plasmon resonance (SPR) studies show that this hot-spot epitope residue contributes significantly more to the binding than an adjacent epitope residue, Lys-96, which does not form a salt bridge in any of the three HH-HEL complexes. Furthermore, the effect of mutating Lys-97 is most severe in HH26-HEL. Lys-96, being a charged residue, also contributes the most in HH26-HEL among the three complexes. The SPR results on these mutants also highlight that the apparent "electrostatic steering" on net on rates actually act at post-collision level stabilization of the complex. The significance of this work is the observed variations in electrostatic interactions among the three complexes. Our work demonstrates that higher electrostatics, both as a number of short-range electrostatic interactions and their contributions, leads to higher binding specificity. Strong salt bridges, their networking, and electrostatically driven binding, limit flexibilities through geometric constrains. In contrast, hydrophobic driven binding and low levels of electrostatic interactions are associated with conformational flexibility and cross-reactivity.     

Salt bridges in prion proteins are necessary for high-affinity binding to the monoclonal antibody T2

We studied the role of the 2 salt bridges (Asp143-Arg147 and Asp146-Arg150) in helix 1 of mouse prion protein (PrP) on the formation of the complex between PrP and the monoclonal antibody T2. We introduced 6 charge-changing mutations to the amino acid residues associated with the salt bridges. Analysis of the circular dichroism spectra of the mutant PrPs showed that the salt bridge mutations did not change the secondary structures. We analyzed the kinetics of the association and dissociation of the PrPs with the T2 antibody. The results showed that the association kinetics were not significantly different among the variants except Arg150Lys, while the dissociation rate of the neutralized-charge variants was 2 orders of magnitude higher than that of the wild type. These results indicate that salt bridges make the interaction of PrP with T2 tighter by slowing down dissociation.     

Specificity and cooperativity at β-lactamase position 104 in TEM-1/BLIP and SHV-1/BLIP interactions

Establishing a quantitative understanding of the determinants of affinity in protein–protein interactions remains challenging. For example, TEM‐1/β‐lactamase inhibitor protein (BLIP) and SHV‐1/BLIP are homologous β‐lactamase/β‐lactamase inhibitor protein complexes with disparate K_d values (3 nM and 2 μM, respectively), and a single substitution, D104E in SHV‐1, results in a 1000‐fold enhancement in binding affinity. In TEM‐1, E104 participates in a salt bridge with BLIP K74, whereas the corresponding SHV‐1 D104 does not in the wild type SHV‐1/BLIP co‐structure. Here, we present a 1.6 Å crystal structure of the SHV‐1 D104E/BLIP complex that demonstrates that this point mutation restores this salt bridge. Additionally, mutation of a neighboring residue, BLIP E73M, results in salt bridge formation between SHV‐1 D104 and BLIP K74 and a 400‐fold increase in binding affinity. To understand how this salt bridge contributes to complex affinity, the cooperativity between the E/K or D/K salt bridge pair and a neighboring hot spot residue (BLIP F142) was investigated using double mutant cycle analyses in the background of the E73M mutation. We find that BLIP F142 cooperatively stabilizes both interactions, illustrating how a single mutation at a hot spot position can drive large perturbations in interface stability and specificity through a cooperative interaction network. Proteins 2011. © 2011 Wiley‐Liss, Inc.     

pH stability of penicillin acylase from <Emphasis Type="Italic">Escherichia coli</Emphasis>

The inactivation kinetics of penicillin acylase from Escherichia coli have been investigated over a wide pH range at 25 and 50 degrees C. The enzyme was very stable in neutral solutions and quickly lost its catalytic activity in acidic and alkaline solutions. In all cases, the inactivation proceeded according to first order reaction kinetics. Analysis of the pH dependence of enzyme stability provides evidence that stable penicillin acylase conformation is maintained by salt bridges. Destruction of the salt bridges due to protonation/deprotonation of the amino acid residues forming these ion pairs causes inactivation by formation of the unstable "acidic" EH(4)(3+), EH(3)(2+), EH(2)(+) and "alkaline" E(-) enzyme forms. At temperatures above 35 degrees C penicillin acylase apparently undergoes a conformational change that is accompanied by destruction of one of these salt bridges and change in the catalytic properties.     

Molecular Dynamics Simulations of Peptides from the Central Domain of Smooth Muscle Caldesmon

Abstract

The central domain of smooth muscle caldesmon contains a highly charged region consisting of ten 13-residue repeats. Experimental evidence obtained from the intact protein and fragments thereof suggests that this entire region forms a single stretch of stable α-helix. We have carried out molecular dynamics simulations on peptides consisting of one, two and three repeats to examine the mechanism of α-helical stability of the central domain at the atomic level. All three peptides show high helical stability on the timescale of the MD simulations. Deviations from α-helical structure in all the simulations arise mainly from the formation of long stretches of π-helix. Interconversion between α-helical and π-helical conformations occurs through insertion of water molecules into α-helical hydrogen bonds and subsequent formation of reverse turns. The α-helical structure is stabilized by electrostatic interactions (salt bridges) between oppositely charged sidechains with i,i+4 spacings, while the π-helix is stabilized by i,i+5 salt bridge interactions. Possible i,i+3 salt bridges are of minor importance. There is a strong preference for salt bridges with a Glu residue N-terminal to a basic sidechain as compared to the opposite orientation. In the double and triple repeat peptides, strong i,i+4 salt bridges exist between the last Glu residue of one repeat and the first Lys residue of the next. This demonstrates a relationship between the repetitive nature of the central domain sequence and its ability to form very long stretches of α-helical structure.     

Effects of engineered salt bridges on the stability of subtilisin BPN'

Variants designed using PROTEUS have been produced in an attempt to engineer stabilizing salt bridges into subtilisin BPN'. All the mutants constructed by site-directed mutagenesis were secreted by Bacillus subtilis, except L75K. Q19E, expressed as a single variant and also in a double variant, Q19E/Q271E, appears to form a stabilizing salt bridge based on X-ray crystal structure determination and differential scanning calorimeter measurements. Although the double mutant was found to be less thermodynamically stable than the wild-type, it did exhibit an autolytic stability about two-fold greater under hydrophobic conditions. Four variants, A98K, S89E, V26R and L235R, were found to be nearly identical to wild-type in thermal stability, indicative of stable structures without evidence of salt bridge formation. Variants Q271E, V51K and T164R led to structures that resulted in varying degrees of thermodynamic and autolytic instability. A computer-modeling analysis of the PROTEUS predictions reveals that the low percentage of salt bridge formation is probably due to an overly simplistic electrostatic model, which does not account for the geometry of the pairwise interactions.     

A structural basis of Equisetum arvense ferredoxin isoform II producing an alternative electron transfer with ferredoxin-NADP+ reductase

We have determined the crystal structure, at 1.2-A resolution, of Equisetum arvense ferredoxin isoform II (FdII), which lacks residues equivalent to Arg(39) and Glu(28) highly conserved among other ferredoxins (Fds). In other Fds these residues form an intramolecular salt bridge crucial for stabilization of the [2Fe-2S] cluster, which is disrupted upon complex formation with Fd-NADP(+) oxidoreductase (FNR) to form two intermolecular salt bridges. The overall structure of FdII resembles the known backbone structures of E. arvense isoform I (FdI) and other plant-type Fds. Dramatically, in the FdII structure a unique, alternative salt bridge is formed between Arg(22) and Glu(58). This results in a different relative orientation of the alpha-helix formed by Leu(23)-Glu(29) and eliminates the possibility of forming three of the five intermolecular salt bridges identified on formation of a complex between maize FdI and maize FNR. Mutation of FdII, informed by structural differences with FdI, showed that the alternative salt bridge and the absence of an otherwise conserved Tyr residue are important for the alternative stabilization of the FdII [2Fe-2S] cluster. We also investigated FdI and FdII electron transfer to FNR on chloroplast thylakoid membranes. The K(m) and V(max) values of FdII are similar to those of FdI, contrary to previous measurements of the reverse reaction, from FNR to Fd. The affinity between reduced FdI and oxidized FNR is much greater than that between oxidized FdI and reduced FNR, whereas this is not the case with FdII. The pH dependence of electron transfer by FdI, FdII, and an FdII mutant with FdI features was measured and further indicated that the binding mode to FNR differs between FdI and FdII. Based on this evidence, we hypothesize that binding modes with other Fd-dependent reductases may also vary between FdI and FdII. The structural differences between FdI and FdII therefore result in functional differences that may influence partitioning of electrons into different redox metabolic pathways.     

Molecular dynamics simulations of peptides from the central domain of smooth muscle caldesmon

The central domain of smooth muscle caldesmon contains a highly charged region consisting of ten 13-residue repeats. Experimental evidence obtained from the intact protein and fragments thereof suggests that this entire region forms a single stretch of stable alpha-helix. We have carried out molecular dynamics simulations on peptides consisting of one, two and three repeats to examine the mechanism of alpha-helical stability of the central domain at the atomic level. All three peptides show high helical stability on the timescale of the MD simulations. Deviations from alpha-helical structure in all the simulations arise mainly from the formation of long stretches of pi-helix. Interconversion between alpha-helical and pi-helical conformations occurs through insertion of water molecules into alpha-helical hydrogen bonds and subsequent formation of reverse turns. The alpha-helical structure is stabilized by electrostatic interactions (salt bridges) between oppositely charged sidechains with i,i+4 spacings, while the pi-helix is stabilized by i,i+5 salt bridge interactions. Possible i,i+3 salt bridges are of minor importance. There is a strong preference for salt bridges with a Glu residue N-terminal to a basic sidechain as compared to the opposite orientation. In the double and triple repeat peptides, strong i,i+4 salt bridges exist between the last Glu residue of one repeat and the first Lys residue of the next. This demonstrates a relationship between the repetitive nature of the central domain sequence and its ability to form very long stretches of alpha-helical structure.     

Statistical characterization of salt bridges in proteins

The structure and folding mechanism of a given protein are determined by many factors, including the electrostatic interactions between charged residues of protein molecules known in general as salt bridges. In this study, analyses were conducted on 10,370 salt bridges in 2017 proteins and the results compared to previous statistical surveys of 36 protein structures. Although many of the general trends remained consistent with other studies, more detailed information was illuminated by the larger dataset. In particular, it was shown that there is a strong correlation between secondary structure and salt bridge formation, and that salt bridges display preferential formation in an environment of about 30% solvent accessible surface area.     

Molecular mechanisms of functional rescue mediated by P53 tumor suppressor mutations

We have utilized both molecular dynamics simulations and solution biophysical measurements to investigate the rescue mechanism of mutation N235K, which plays a key role in the recently identified global suppressor motif of K235/Y239/R240 in the human p53 DNA-binding domain (DBD). Previous genetic analysis indicates that N235K alone rescues five out of six destabilized cancer mutants. However, the solution biophysical measurement shows that N235K generates only a slight increase to the stability of DBD, implying a rescue mechanism that is not a simple additive contribution to thermodynamic stability. Our molecular simulations show that the N235K substitution generates two non-native salt bridges with residues D186 and E198. We find that the nonnative salt bridges, D186-K235 and E198-K235, and a native salt bridge, E171-R249, are mutually exclusive, thus resulting in only a marginal increase in stability as compared to the wild type protein. When a destabilized V157F is paired with N235K, the native salt bridge E171-R249 is retained. In this context, the non-native salt bridges, D186-K235 and E198-K235, produce a net increase in stability as compared to V157F alone. A similar rescue mechanism may explain how N235K stabilize other highly unstable β-sandwich cancer mutants.     

Salt bridge stability in monomeric proteins.

Here, we present the results of continuum electrostatic calculations on a dataset of 222 non-equivalent salt bridges derived from 36 non-homologous high-resolution monomeric protein crystal structures. Most of the salt bridges in our dataset are stabilizing, regardless of whether they are buried or exposed, isolated or networked, hydrogen bonded or non-hydrogen bonded. One-third of the salt bridges in our dataset are buried in the protein core, with the remainder exposed to the solvent. The difference in the dielectric properties of water versus the hydrophobic protein interior cost buried salt bridges large desolvation penalties. However, the electrostatic interactions both between the salt-bridging side-chains, and between the salt bridges and charges in their protein surroundings, are also stronger in the interior, due to the absence of solvent screening. Even large desolvation penalties for burying salt bridges are frequently more than compensated for, primarily by the electrostatic interactions between the salt-bridging side-chains. In networked salt bridges both types of electrostatic interactions, those between the salt-bridging side-chains, and those between the salt bridge and its protein environment, are of similar magnitudes. In particular, a major finding of this work is that salt bridge geometry is a critical factor in determining salt bridge stability. Salt bridges with favorable geometrical positioning of the interacting side-chain charged groups are likely to be stabilizing anywhere in the protein structure. We further find that most of the salt bridges are formed between residues that are relatively near each other in the sequence.     

Electrostatic interactions in a peptide--RNA complex

Parallel experimental measurements and theoretical calculations have been used to investigate the energetics of electrostatic interactions in the complex formed between a 22 residue, alpha-helical peptide from the N protein of phage lambda and its cognate 19 nucleotide box B RNA hairpin. Salt-dependent free energies were measured for both peptide folding from coil to helix and peptide binding to RNA, and from these the salt-dependence of binding pre-folded, helical peptide to RNA was determined ( partial differential (DeltaG degrees (dock))/ partial differential log[KCl]=5.98(+/-0.21)kcal/mol). (A folding transition taking place in the RNA hairpin loop was shown to have a negligible dependence on salt concentration.) The non-linear Poisson-Boltzmann equation was used to calculate the same salt dependence of the binding free energy as 5.87(+/-0.22)kcal/mol, in excellent agreement with the measured value. Close agreement between experimental measurements and calculations was also obtained for two variant peptides in which either a basic or acidic residue was replaced with an uncharged residue, and for an RNA variant with a deletion of a single loop nucleotide. The calculations suggest that the strength of electrostatic interactions between a peptide residue and RNA varies considerably with environment, but that all 12 positive and negative N peptide charges contribute significantly to the electrostatic free energy of RNA binding, even at distances up to 11A from backbone phosphate groups. Calculations also show that the net release of ions that accompanies complex formation originates from rearrangements of both peptide and RNA ion atmospheres, and includes accumulation of ions in some regions of the complex as well as displacement of cations and anions from the ion atmospheres of the RNA and peptide, respectively.