Strength and co-operativity of contributions of surface salt bridges to protein stability

Many of the interactions that stabilize proteins are co-operative and cannot be reduced to a sum of pairwise interactions. Such interactions may be analysed by protein engineering methods using multiple thermodynamic cycles comprising wild-type protein and all combinations of mutants in the interacting residues. There is a triad of charged residues on the surface of barnase, comprising residues Asp8, Asp12 and Arg110, that interact by forming two exposed salt bridges. The three residues have been mutated to alanine to give all the single, double and triple mutants. The free energies of unfolding of wild-type and the seven mutant proteins have been determined and the results analysed to give the contributions of the residues in the two salt bridges to protein stability. It is possible to isolate the energies of forming the salt bridges relative to the solvation of the separated ions by water. In the intact triad, the apparent contribution to the stabilization energy of the protein of the salt bridge between Asp12 and Arg110 is -1.25 kcal mol-1, whereas that of the salt bridge between Asp8 with Arg110 is -0.98 kcal mol-1. The strengths of the two salt bridges are coupled: the energy of each is reduced by 0.77 kcal mol-1 when the other is absent. The salt-linked triad, relative to alanine residues at the same positions, does not contribute to the stability of the protein since the favourable interactions of the salt bridges are more than offset by other electrostatic and non-electrostatic energy terms. Salt-linked triads occur in other proteins, for example, haemoglobin, where the energy of only the salt-bridge term is important and so the coupling of salt bridges could be of general importance to the stability and function of proteins.     

Geometric cooperativity and anticooperativity of three-body interactions in native proteins

Characterizing multibody interactions of hydrophobic, polar, and ionizable residues in protein is important for understanding the stability of protein structures. We introduce a geometric model for quantifying 3-body interactions in native proteins. With this model, empirical propensity values for many types of 3-body interactions can be reliably estimated from a database of native protein structures, despite the overwhelming presence of pairwise contacts. In addition, we define a nonadditive coefficient that characterizes cooperativity and anticooperativity of residue interactions in native proteins by measuring the deviation of 3-body interactions from 3 independent pairwise interactions. It compares the 3-body propensity value from what would be expected if only pairwise interactions were considered, and highlights the distinction of propensity and cooperativity of 3-body interaction. Based on the geometric model, and what can be inferred from statistical analysis of such a model, we find that hydrophobic interactions and hydrogen-bonding interactions make nonadditive contributions to protein stability, but the nonadditive nature depends on whether such interactions are located in the protein interior or on the protein surface. When located in the interior, many hydrophobic interactions such as those involving alkyl residues are anticooperative. Salt-bridge and regular hydrogen-bonding interactions, such as those involving ionizable residues and polar residues, are cooperative. When located on the protein surface, these salt-bridge and regular hydrogen-bonding interactions are anticooperative, and hydrophobic interactions involving alkyl residues become cooperative. We show with examples that incorporating 3-body interactions improves discrimination of protein native structures against decoy conformations. In addition, analysis of cooperative 3-body interaction may reveal spatial motifs that can suggest specific protein functions.     

Photosynthetic water oxidation in Synechocystis sp. PCC6803: mutations D1-E189K, R and Q are without influence on electron transfer at the donor side of photosystem II.

The oxygen-evolving manganese cluster (OEC) of photosynthesis is oxidised by the photochemically generated primary oxidant (P(+*)(680)) of photosystem II via a tyrosine residue (Y(Z), Tyr161 on the D1 subunit of Synechocystis sp. PCC6803). The redox span between these components is rather small and probably tuned by protonic equilibria. The very efficient electron transfer from Y(Z) to P(+*)(680) in nanoseconds requires the intactness of a hydrogen bonded network involving Y(Z), D1-His190, and presumably D1-Glu189. We studied photosystem II core particles from photoautotrophic mutants where the residue D1-E189 was replaced by glutamine, arginine and lysine which were expected to electrostatically differ from the glutamate in the wild-type (WT). Surprisingly, the rates of electron transfer from Y(Z) to P(+*)(680) as well as from the OEC to Y(ox)(Z) were the same as in the WT. With the generally assumed proximity between D1-His190 (and thus D1-Glu189) and Y(Z), the lack of any influence on the electron transfer around Y(Z) straightforwardly implies a strongly hydrophobic environment forcing Glu (acid) and Lys, Arg (basic) at position D1-189 into electro-neutrality. As one alternative, D1-Glu189 could be located at such a large distance from the OEC, Y(Z) and P(+*)(680) that a charge on D1-189X does not influence the electron transfer. This seems less likely in the light of the drastic influence of its direct neighbour, D1-His190, on Y(Z) function. Another alternative is that D1-Glu189 is negatively charged, but is located in a cluster of acid/base groups that compensates for an alteration of charge at position 189, leaving the overall net charge unchanged in the Gln, Lys, and Arg mutants.     

Understanding antibody–antigen associations by molecular dynamics simulations: Detection of important intra- and inter-molecular salt bridges

1 NSec molecular dynamics (MD) simulation of anti-hen egg white antibody, HyHEL63 (HH63), complexed with HEL reveals important molecular interactions, not revealed in its X-ray crystal structure. These molecular interactions were predicted to be critical for the complex formation, based on structure–function studies of this complex and 3-other anti-HEL antibodies, HH8, HH10 and HH26, HEL complexes. All four antibodies belong to the same structural family, referred to here as HH10 family. Ala scanning results show that they recognize ‘coincident epitopes’. 1 NSec explicit, with periodic boundary condition, MD simulation of HH63-HEL reveals the presence of functionally important salt-bridges. Around 200 ps in vacuo and an additional 20 ps explicit simulation agree with the observations from 1 Nsec simulation. Intra-molecular salt-bridges predicted to play significant roles in the complex formation, were revealed during MD simulation. A very stabilizing salt-bridge network, and another intra-molecular salt-bridge, at the binding site of HEL, revealed during the MD simulation, is proposed to predipose binding site geometry for specific binding. All the revealed salt-bridges are present in one or more of the other three complexes and/or involve “hot-spot” epitope and paratope residues. Most of these charged epitope residues make large contribution to the binding free energy. The “hot spot” epitope residue Lys97_Y, which significantly contributes to the free energy of binding in all the complexes, forms an intermolecular salt-bridge in several MD conformers. Our earlier computations have shown that this inter-molecular salt-bridge plays a significant role in determining specificity and flexibility of binding in the HH8-HEL and HH26-HEL complexes. Using a robust criterion of salt-bridge detection, this inter-molecular salt-bridge was detected in the native structures of the HH8-HEL and HH26-HEL complexes, but was not revealed in the crystal structure of HH63-HEL complex. The electrostatic strength of this revealed salt-bridge was very strong. During 1 Nsec MD simulation this salt-bridge networks with another inter-molecular salt-bridge to form an inter-molecular salt-bridge triad. Participation of Lys97_Y in the formation of inter-molecular triad further validates the functional importance of Lys97_Y in HH63-HEL associations. These results demonstrate that many important structural details of biomolecular interactions can be better understood when studied in a dynamic environment, and that MD simulations can complement and expand information obtained from static X-ray structure. This study also highlights “hot-spot” molecular interactions in HyHEL63-HEL complex. The publisher or recipient acknowledges right of the U.S. Government to retain a non exclusive, royalty-free license in and to any copyright covering the article.     

Dimeric Crystal Structure of Rabbit l-Gulonate 3-Dehydrogenase/λ-Crystallin: Insights into the Catalytic Mechanism

l-Gulonate 3-dehydrogenase (GDH) is a bifunctional dimeric protein that functions not only as an NAD⁺-dependent enzyme in the uronate cycle but also as a taxon-specific λ-crystallin in rabbit lens. Here we report the first crystal structure of GDH in both apo form and NADH-bound holo form. The GDH protomer consists of two structural domains: the N-terminal domain with a Rossmann fold and the C-terminal domain with a novel helical fold. In the N-terminal domain of the NADH-bound structure, we identified 11 coenzyme-binding residues and found 2 distinct side-chain conformers of Ser124, which is a putative coenzyme/substrate-binding residue. A structural comparison between apo form and holo form and a mutagenesis study with E97Q mutant suggest an induced-fit mechanism upon coenzyme binding; coenzyme binding induces a conformational change in the coenzyme-binding residues Glu97 and Ser124 to switch their activation state from resting to active, which is required for the subsequent substrate recruitment. Subunit dimerization is mediated by numerous intersubunit interactions, including 22 hydrogen bonds and 104 residue pairs of van der Waals interactions, of which those between two cognate C-terminal domains are predominant. From a structure/sequence comparison within GDH homologues, a much greater degree of interprotomer interactions (both polar and hydrophobic) in the rabbit GDH would contribute to its higher thermostability, which may be relevant to the other function of this enzyme as λ-crystallin, a constitutive structural protein in rabbit lens. The present crystal structures and amino acid mutagenesis studies assigned the role of active-site residues: catalytic base for His145 and substrate binding for Ser124, Cys125, Asn196, and Arg231. Notably, Arg231 participates in substrate binding from the other subunit of the GDH dimer, indicating the functional significance of the dimeric state. Proper orientation of the substrate-binding residues for catalysis is likely to be maintained by an interprotomer hydrogen-bonding network of residues Asn196, Gln199, and Arg231, suggesting a network-based substrate recognition of GDH.     

SBION: A Program for Analyses of Salt-Bridges from Multiple Structure Files

Salt-bridge and network salt-bridge are specific electrostatic interactions that contribute to the overall stability of proteins. In hierarchical protein folding model, these interactions play crucial role in nucleation process. The advent and growth of protein structure database and its availability in public domain made an urgent need for context dependent rapid analysis of salt-bridges. While these analyses on single protein is cumbersome and time-consuming, batch analyses need efficient software for rapid topological scan of a large number of protein for extracting details on (i) fraction of salt-bridge residues (acidic and basic). (ii) Chain specific intra-molecular salt-bridges, (iii) inter-molecular salt-bridges (protein-protein interactions) in all possible binary combinations (iv) network salt-bridges and (v) secondary structure distribution of salt-bridge residues. To the best of our knowledge, such efficient software is not available in public domain. At this juncture, we have developed a program i.e. SBION which can perform all the above mentioned computations for any number of protein with any number of chain at any given distance of ion-pair. It is highly efficient, fast, error-free and user friendly. Finally we would say that our SBION indeed possesses potential for applications in the field of structural and comparative bioinformatics studies.

Availability

SBION is freely available for non-commercial/academic institutions on formal request to the corresponding author (ni.ca.vinurub.hcetoib@eejrenabka).     

Electrostatic stabilization of a thermophilic cold shock protein.

The cold shock protein Bc-Csp from the thermophile Bacillus caldolyticus differs from its mesophilic homolog Bs-CspB from Bacillus subtilis by 15.8 kJ mol(-1) in the Gibbs free energy of denaturation (DeltaG(D)). The two proteins vary in sequence at 12 positions but only two of them, Arg3 and Leu66 of Bc-Csp, which replace Glu3 and Glu66 of Bs-CspB, are responsible for the additional stability of Bc-Csp. These two positions are near the ends of the protein chain, but close to each other in the three-dimensional structure. The Glu3Arg exchange alone changed the stability by more than 11 kJ mol(-1). Here, we elucidated the molecular origins of the stability difference between the two proteins by a mutational analysis. Electrostatic contributions to stability were characterized by measuring the thermodynamic stabilities of many variants as a function of salt concentration. Double and triple mutant analyses indicate that the stabilization by the Glu3Arg exchange originates from three sources. Improved hydrophobic interactions of the aliphatic moiety of Arg3 contribute about 4 kJ mol(-1). Another 4 kJ mol(-1) is gained from the relief of a pairwise electrostatic repulsion between Glu3 and Glu66, as in the mesophilic protein, and 3 kJ mol(-1) originate from a general electrostatic stabilization by the positive charge of Arg3, which is not caused by a pairwise interaction. Mutations of all potential partners for an ion pair within a radius of 10 A around Arg3 had only marginal effects on stability. The Glu3-->Arg3 charge reversal thus optimizes ionic interactions at the protein surface by both local and global effects. However, it cannot convert the coulombic repulsion with another Glu residue into a corresponding attraction. Avoidance of unfavorable coulombic repulsions is probably a much simpler route to thermostability than the creation of stabilizing surface ion pairs, which can form only at the expense of conformational entropy.     

Upward shift of the pH optimum of Acremonium ascorbate oxidase

A gene encoding a thermostable Acremonium ascorbate oxidase (ASOM) was randomly mutated to generate mutant enzymes with altered pH optima. One of the mutants, which exhibited a significantly higher activity in the pH range 4.5-7 compared to ASOM, had a Gln183Arg substitution in the region corresponding to SBR1, one of the substrate binding regions of the zucchini enzyme. The other mutant with almost the same pH profile as Gln183Arg had a Thr527Ala substitution near the type 3 copper center and became more sensitive to azide than ASOM. Site-directed mutagenesis in the substrate binding regions with reference to the amino acid sequences of plant enzymes led to isolation of mutants shifted upward in the pH optimum; Val193Pro and Val193Pro/Pro190Ile increased the pH optimum by 1 and 0.5 units, respectively, while retaining the near-wild-type thermostability and azide sensitivity. The homology model of ASOM constructed from the zucchini enzyme coordinates suggested that replacement of Val193 by Pro could disturb the ion pair networks among Arg309, Glu192, Arg194 and Glu311. This perturbation could affect either the molecular recognition between the substrate and ASOM or the electron transfer from the substrate to the type 1 copper center, leading to the alkaline shift of the catalytic activity of the mutant enzyme. The other mutations, Val193Pro/Pro190Ile, could also induce similar structural perturbations involving the ion pair networks.     

Electrostatic steering and ionic tethering in the formation of thrombin-hirudin complexes: the role of the thrombin anion-binding exosite-I

Electrostatic interactions between the thrombin anion-binding exosite-I (ABE-I) and the hirudin C-terminal tail play an important role in the formation of the thrombin-hirudin inhibitor complex and serves as a model for the interactions of thrombin with its many other ligands. The role of each solvent exposed basic residue in ABE-I (Arg(35), Lys(36), Arg(67), Arg(73), Arg(75), Arg(77a), Lys(81), Lys(109), Lys(110), and Lys(149e)) in electrostatic steering and ionic tethering in the formation of thrombin-hirudin inhibitor complexes was explored by site directed mutagenesis. The contribution to the binding energy (deltaG(degrees)b) by each residue varied from 1.9 kJ mol(-)(1) (Lys(110)) to 15.3 kJ mol(-1) (Arg(73)) and were in general agreement to their observed interactions with hirudin residues in the thrombin-hirudin crystal structure [Rydel, T. J., Tulinsky, A., Bode, W., and Huber, R. (1991) J. Mol. Biol. 221, 583-601]. Coupling energies (delta deltaG(degrees) int) were calculated for the major ion-pair interactions involved in ionic tethering using complementary hirudin mutants (h-D55N, h-E57Q, and h-E58Q). Cooperativity was seen for the h-Asp(55)/Arg(73) ion pair (2.4 kJ mol(-1)); however, low coupling energies for h-Asp(55)/Lys(149e) (deltadeltaG(degrees)int 0.6 kJ mol(-1)) and h-Glu(58)/Arg(77a) (deltadeltaG(degrees)int 0.9 kJ mol(-1)) suggest these are not major interactions, as anticipated by the crystal structure. Interestingly, high coupling energies were seen for the intermolecular ion-pair h-Glu(57)/Arg(75) (deltadeltaG(degrees)int 2.3 kJ mol(-1)) and for the solvent bridge h-Glu(57)/Arg(77a) (deltadeltaG(degrees)int 2.7 kJ mol(-1)) indicating that h-Glu(57) interacts directly with both Arg(75) and Arg(77a) in the thrombin-hirudin inhibitor complex. The remaining ABE-I residues that do not form major contacts in tethering the C-terminal tail of hirudin make small but collectively important contributions to the overall positive electrostatic field generated by ABE-I important in electrostatic steering.     

Energetics of charge-charge interactions between residues adjacent in sequence

Statistical analysis of the residue separation between a pair of ionizable side chains within 4 ? of each other was performed on a set of 1560 non-homologous PDB structures. We found that the frequency of pairs of like charges (i.e., pairs consisting of acidic residues Asp and Glu or pairs consisting of basic residues Arg and Lys) is two orders of magnitude lower than the pairs of oppositely charged residues (salt-bridges). We also found that for pairs of like charges the distribution is skewed dramatically towards short residue separation (<3). On the basis of these observations, we hypothesize that at short residue separation the repulsion between charges does not contribute much to the protein stability and the effects are largely dominated by the long range charge-charge interactions with other ionizable groups in the protein molecule. To test this hypothesis, we incorporated various pairs of charged residues at position 63 and 64 of ubiquitin and compared the stabilities of these variants. We also performed calculations of the expected changes in the charge-charge interactions. A very good correlation between experimental changes in the stability of ubiquitin variants, and changes in the energy of charge-charge interactions provides support for the hypothesis that a pair of ionizable residues next to each other in sequence modulates protein stability via long range charge-charge interactions with the rest of the protein.     

The interaction between residues 62 and 193 play a key role in activity and structural stability of arginine kinase

The purpose of this study is to clarify that the amino acid residues (Asp62 and Arg193) are responsible for the activity and stability of arginine kinase (AK). The amino acid residues Asp62 (D62) and Arg193 (R193) are strictly conserved in monomeric AKs and form an ion pair in the transition state analogue complex. In this research, we replaced D62 with glutamate (E) or glycine (G) and R193 with lysine (K) or glycine (G). The mutants of D62E and R193K retained almost 90% of the wild-type activity, whereas D62G and R193G had a pronounced loss in activity. A detailed comparison was made between the physic-chemical properties and conformational changes of wild-type AK and the mutants by means of ultraviolet (UV) difference and fluorescence spectra. The results indicated that the conformation of all of the mutants had been changed and the stability in a urea solution was also reduced. We speculated that the hydrogen bond and electrostatic interactions formed between residues 62 and 193 play a key role in stabilizing the structure and mediating the synergism in substrate binding of arginine kinase from greasyback shrimp (Metapenaeus ensis).     

Surface Charges and Regulation of FMN to Heme Electron Transfer in Nitric-oxide Synthase

The nitric-oxide synthases (NOS, EC 1.14.13.39) are modular enzymes containing attached flavoprotein and heme (NOSoxy) domains. To generate nitric oxide (NO), the NOS FMN subdomain must interact with the NOSoxy domain to deliver electrons to the heme for O₂ activation during catalysis. The molecular basis and how the interaction is regulated is unclear. We explored the role of eight positively charged residues that create an electropositive patch on NOSoxy in enabling the electron transfer by incorporating mutations that neutralized or reversed their individual charges. Stopped-flow and steady-state experiments revealed that individual charges at Lys⁴²³, Lys⁶²⁰, and Lys⁶⁶⁰ were the most important in enabling heme reduction in nNOS. Charge reversal was more disruptive than neutralization in all cases, and the effects on heme reduction were not due to a weakening in the thermodynamic driving force for heme reduction. Mutant NO synthesis activities displayed a complex pattern that could be simulated by a global model for NOS catalysis. This analysis revealed that the mutations impact the NO synthesis activity only through their effects on heme reduction rates. We conclude that heme reduction and NO synthesis in nNOS is enabled by electrostatic interactions involving Lys⁴²³, Lys⁶²⁰, and Lys⁶⁶⁰, which form a triad of positive charges on the NOSoxy surface. A simulated docking study reveals how electrostatic interactions of this triad can enable an FMN-NOSoxy interaction that is productive for electron transfer.     

Probing possible downhill folding: native contact topology likely places a significant constraint on the folding cooperativity of proteins with approximately 40 residues

Experiments point to appreciable variations in folding cooperativity among natural proteins with approximately 40 residues, indicating that the behaviors of these proteins are valuable for delineating the contributing factors to cooperative folding. To explore the role of native topology in a protein's propensity to fold cooperatively and how native topology might constrain the degree of cooperativity achievable by a given set of physical interactions, we compared folding/unfolding kinetics simulated using three classes of native-centric C^α chain models with different interaction schemes. The approach was applied to two homologous 45-residue fragments from the peripheral subunit-binding domain family and a 39-residue fragment of the N-terminal domain of ribosomal protein L9. Free-energy profiles as functions of native contact number were computed to assess the heights of thermodynamic barriers to folding. In addition, chevron plots of folding/unfolding rates were constructed as functions of native stability to facilitate comparison with available experimental data. Although common Gō-like models with pairwise Lennard-Jones-type interactions generally fold less cooperatively than real proteins, the rank ordering of cooperativity predicted by these models is consistent with experiment for the proteins investigated, showing increasing folding cooperativity with increasing nonlocality of a protein's native contacts. Models that account for water-expulsion (desolvation) barriers and models with many-body (nonadditive) interactions generally entail higher degrees of folding cooperativity indicated by more linear model chevron plots, but the rank ordering of cooperativity remains unchanged. A robust, experimentally valid rank ordering of model folding cooperativity independent of the multiple native-centric interaction schemes tested here argues that native topology places significant constraints on how cooperatively a protein can fold.     

Aromatic stacking as a determinant of the thermal stability of CYP119 from Sulfolobus solfataricus

Two notable features of the thermophilic CYP119, an Arg154-Glu212 salt bridge between the F-G loop and the I helix and an extended aromatic cluster, were studied to determine their contributions to the thermal stability of the enzyme. Site-specific mutants of the salt bridge (Arg154, Glu212) and aromatic cluster (Tyr2, Trp4, Trp231, Tyr250, Trp281) were expressed and purified. The substrate-binding and kinetic constants for lauric acid hydroxylation are little affected in most mutants, but the E212D mutant is inactive and the R154Q mutant has higher K(s),K(m), and k(cat) values. The salt bridge mutants, like wild-type CYP119, melt at 91+/-1 degrees C, whereas mutation of individual residues in the extended aromatic cluster lowers the T(m) by 10-15 degrees C even though no change is observed on mutation of an unrelated aromatic residue. The extended aromatic cluster, but not the Arg154-Glu212 salt bridge, contributes to the thermal stability of CYP119.     

Double-mutant cycle scanning of the interaction of a peptide ligand and its G protein-coupled receptor

The interaction between the yeast G protein-coupled receptor (GPCR), Ste2p, and its alpha-factor tridecapeptide ligand was subjected to double-mutant cycle scanning analysis by which the pairwise interaction energy of each ligand residue with two receptor residues, N205 and Y266, was determined. The mutations N205A and Y266A were previously shown to result in deficient signaling but cause only a 2.5-fold and 6-fold decrease, respectively, in the affinity for alpha-factor. The analysis shows that residues at the amine terminus of alpha-factor interact strongly with N205 and Y266 whereas residues in the center and at the carboxyl terminus of the peptide interact only weakly if at all with these receptor residues. Multiple-mutant thermodynamic cycle analysis was used to assess whether the energies of selected pairwise interactions between residues of the alpha-factor peptide changed upon binding to Ste2p. Strong positive cooperativity between residues 1 through 4 of alpha-factor was observed during receptor binding. In contrast, no thermodynamic evidence was found for an interaction between a residue near the carboxyl terminus of alpha-factor (position 11) and one at the N-terminus (position 3). The study shows that multiple-mutant cycle analyses of the binding of an alanine-scanned peptide to wild-type and mutant GPCRs can provide detailed information on contributions of inter- and intramolecular interactions to the binding energy and potentially prove useful in developing 3D models of ligand docked to its receptor.     

Substrate specificity at the P1' site of Escherichia coli OmpT under denaturing conditions

Though OmpT has been reported to mainly cleave the peptide bond between consecutive basic amino acids, we identified more precise substrate specificity by using a series of modified substrates, termed PRX fusion proteins, consisting of 184 residues. The cleavage site of the substrate PRR was Arg140-Arg141 and the modified substrates PRX substituted all 19 natural amino acids at the P1' site instead of Arg141. OmpT under denaturing conditions (in the presence of 4 M urea) cleaved not only between two consecutive basic amino acids but also at the carboxyl side of Arg140 except for the Arg140-Asp141, -Glu141, and -Pro141 pairs. In addition to Arg140 at the P1 site, similar results were obtained when Lys140 was substituted into the P1 site. In the absence of urea, an aspartic acid residue at the P1' site was unfavorable for OmpT cleavage of synthetic decapeptides, the enzyme showed a preference for a dibasic site.     

Interactions of Arg2 in the Mnt N-terminal arm with the central and flanking regions of the Mnt operator

Arg2, in the N-terminal arm of the Mnt repressor, plays an important role in determining operator-binding specificity. In the complex of the Mnt tetramer with the 21 base-pair mnt operator, there are four potential sites for Arg2 interactions, two in the central region of the operator and two on the outer flanks of the operator. Single-chain variants of the dimeric N-terminal domain of Mnt containing one Arg2 residue and one Lys2 or Met2 residue were constructed and interactions with operator DNA were probed using Fe. EDTA affinity cleavage. The results of these orientation studies show that the majority of the energetically significant interactions mediated by Arg2 occur in the central region of the mnt operator. The RK2, RA2, and RM2 mutations reduce the free energy of operator binding by 1.7 kcal/mol, 3.3 kcal/mol, and 4.9 kcal/mol, respectively. Double-mutant thermodynamic cycle analyses using the RA2, RM2, and operator variants also reveal interaction free energies between Arg2 and operator base-pairs 9, 10, 11, 12 and 13, which in aggregate account for most of the Arg2 contribution to operator binding.     

Why is the Arg5-Glu13 salt bridge conserved in mammalian alpha-defensins?

Mammalian alpha-defensins, expressed primarily in leukocytes and epithelia, kill a broad range of microbes, constituting one of the first lines of innate immune defense against infection. Nine amino acid residues, including six cysteines, one glycine, and a pair of oppositely charged residues Arg/Glu, are conserved in the otherwise diverse sequences of all known mammalian alpha-defensins. Structural analysis indicates that the two charged residues form a salt bridge, likely stabilizing a protruding loop in the molecule. To investigate the structural and functional roles of the conserved Arg5-Glu13 salt bridge in alpha-defensins, we chemically prepared human neutrophil alpha-defensin 2 (HNP2) and five HNP2 analogs, R5E/E13R, E13Q, E13R, R5T/E13Y, and R14A. In contrast to HNP2 and R14A-HNP2, none of the four salt bridge analogs was capable of folding into a native conformation in the context of isolated defensin domains. However, when covalently attached to the 45-residue pro-HNP2 propeptide, the salt bridge analogs of HNP2 in their pro-forms all folded productively, suggesting that the Arg5-Glu13 salt bridge is not required for correct pro-alpha-defensin folding. When assayed against both Escherichia coli and Staphylococcus aureus, the six alpha-defensins showed bactericidal activity that correlated with the number of net positive charges carried by individual molecules in the panel, irrespective of whether or not the Arg5-Glu13 salt bridge was decimated, suggesting that Arg5 and Glu13 are not functionally conserved. Proteolytic resistance analysis with human neutrophil elastase, one major protease contained in azurophils with HNPs, revealed that destabilization of the salt bridge dramatically accelerated defensin degradation by the enzyme. Thus, we propose that the Arg5-Glu13 salt bridge found in most mammalian alpha-defensins is conserved for defensin in vivo stability.     

Cooperative helix stabilization by complex Arg-Glu salt bridges

Among the interactions that stabilize the native state of proteins, the role of electrostatic interactions has been difficult to quantify precisely. Surface salt bridges or ion pairs between acidic and basic side chains have only a modest stabilizing effect on the stability of helical peptides or proteins: estimates are roughly 0.5 kcal/mol or less. On the other hand, theoretical arguments and the occurrence of salt bridge networks in thermophilic proteins suggest that multiple salt bridges may exert a stronger stabilizing effect. We show here that triads of charged side chains, Arg(+)-Glu(-)-Arg(+) spaced at i,i+4 or i,i+3 intervals in a helical peptide stabilize alpha helix by more than the additive contribution of two single salt bridges. The free energy of the triad is more than 1 kcal/mol in excess of the sum of the individual pairs, measured in low salt concentration (10 mM). The effect of spacing the three groups is severe; placing the charges at i,i+4 or i,i+3 sites has a strong effect on stability relative to single bridges; other combinations are weaker. A conservative calculation suggests that interactions of this kind between salt bridges can account for much of the stabilization of certain thermophilic proteins.