Stabilization of a fibronectin type III domain by the removal of unfavorable electrostatic interactions on the protein surface

It is generally considered that electrostatic interactions on the protein surface, such as ion pairs, contribute little to protein stability, although they may play important roles in conformational specificity. We found that the tenth fibronectin type III domain of human fibronectin (FNfn10) is more stable at acidic pH than neutral pH, with an apparent midpoint of transition near pH 4. Determination of pK(a)'s for all the side chain carboxyl groups of Asp and Glu residues revealed that Asp 23 and Glu 9 have an upshifted pK(a). These residues and Asp 7 form a negatively charged patch on the surface of FNfn10, with Asp 7 centrally located between Asp 23 and Glu 9, suggesting repulsive electrostatic interactions among these residues at neutral pH. Mutant proteins, D7N and D7K, in which Asp 7 was replaced with Asn and Lys, respectively, exhibited a modest but significant increase in stability at neutral pH, compared to the wild type, and they no longer showed pH dependence of stability. The pK(a)'s of Asp 23 and Glu 9 in these mutant proteins shifted closer to their respective unperturbed values, indicating that the unfavorable electrostatic interactions have been reduced in the mutant proteins. Interestingly, the wild-type and mutant proteins were all stabilized to a similar degree by the addition of 1 M sodium chloride at both neutral and acidic pH, suggesting that the repulsive interactions between the carboxyl groups cannot be effectively shielded by 1 M sodium chloride. These results indicate that repulsive interactions between like charges on the protein surface can destabilize a protein, and protein stability can be significantly improved by relieving these interactions.     

Cooperation between a salt bridge and the hydrophobic core triggers fold stabilization in a Trp-cage miniprotein

Miniproteins are adequate models to study various protein-structure modifying effects such as temperature, pH, point mutation(s), H-bonds, salt bridges, molecular packing, etc. Tc5b, a 20-residue Trp-cage protein is one of the smallest of such models with a stable 3D fold (Neidigh J. W. et al. (2002) Nat. Struct. Biol. 9, 425-430). However, Tc5b exhibits considerable heat-sensitivity and is only stable at relatively low temperatures. Here we report a systematic investigation of structural factors influencing the stability of Tc5b by solving its solution structure in different environments, varying temperature, and pH. The key interactions identified are the hydrophobic stacking of the aromatic rings of Tyr3 and Trp6 and the salt bridge formed between Asp9 and Arg18. To verify the importance of these interactions, selected variants (mutated, glycosylated and truncated) of Tc5b were designed, prepared, and investigated by NMR. Indeed, elimination of either of the key interactions highly destabilizes the structure. These observations enabled us to design a new variant, Tc6b, differing only by a methylene group from Tc5b, in which both key interactions are optimized simultaneously. Tc6b exhibits enhanced heat stability and adopts a stable fold at physiological temperature.     

Dynamics and cooperativity of Trp-cage folding

In this study, multiple independent molecular dynamics (MD) simulations on Trp-cage folding were performed at 300, 325 and 375 K using generalized Born (GB) implicit solvent model. The orientational movement of the side-chain of Trp6 to form a hydrophobic core with 3₁₀-helix was observed. The breaking/formation of a salt bridge between Asp9 and Arg16 was proposed to be the prerequisite for Trp-cage folding/refolding. Our results demonstrate that the cooperation between the salt bridge and the Trp6 orientation leads to a stable tertiary structure of Trp-cage. Analyses on backbone concerted motions at different temperatures indicate that interactions between Trp6 and 3₁₀-helix & Pro18 and between Pro12 and Pro17 & Pro18 are weakened at 375 K but strengthened at lower temperatures, suggesting that they could be the potential driving force of hydrophobic collapse.     

A survey of aspartate-phenylalanine and glutamate-phenylalanine interactions in the protein data bank: searching for anion-π pairs

Protein structures are stabilized using noncovalent interactions. In addition to the traditional noncovalent interactions, newer types of interactions are thought to be present in proteins. One such interaction, an anion-π pair, in which the positively charged edge of an aromatic ring interacts with an anion, forming a favorable anion-quadrupole interaction, has been previously proposed [Jackson, M. R., et al. (2007) J. Phys. Chem. B111, 8242-8249]. To study the role of anion-π interactions in stabilizing protein structure, we analyzed pairwise interactions between phenylalanine (Phe) and the anionic amino acids, aspartate (Asp) and glutamate (Glu). Particular emphasis was focused on identification of Phe-Asp or -Glu pairs separated by less than 7 ? in the high-resolution, nonredundant Protein Data Bank. Simplifying Phe to benzene and Asp or Glu to formate molecules facilitated in silico analysis of the pairs. Kitaura-Morokuma energy calculations were performed on roughly 19000 benzene-formate pairs and the resulting energies analyzed as a function of distance and angle. Edgewise interactions typically produced strongly stabilizing interaction energies (-2 to -7.3 kcal/mol), while interactions involving the ring face resulted in weakly stabilizing to repulsive interaction energies. The strongest, most stabilizing interactions were identified as preferentially occurring in buried residues. Anion-π pairs are found throughout protein structures, in helices as well as β strands. Numerous pairs also had nearby cation-π interactions as well as potential π-π stacking. While more than 1000 structures did not contain an anion-π pair, the 3134 remaining structures contained approximately 2.6 anion-π pairs per protein, suggesting it is a reasonably common motif that could contribute to the overall structural stability of a protein.     

Importance of hydrogen-bonding interactions involving the side chain of Asp158 in the catalytic mechanism of papain

In a previous study, it was shown that replacing Asp158 in papain by Asn had little effect on activity and that the negatively charged carboxylate of Asp158 does not significantly stabilize the active site thiolate-imidazolium ion pair of papain (Ménard et al., 1990). In this paper, we report the kinetic characterization of three more mutants at this position: Asp158Gly, Asp158Ala, and Asp158Glu. From the pH-activity profiles of these and other mutants of papain, it has been possible to develop a model that enables us to dissect out the contribution of the various mutations toward (i) intrinsic activity, (ii) ion pair stability, and (iii) the electrostatic potential at the active site. Results obtained with mutants that place either Gly or Ala at position 158 indicate that the hydrogen bonds involving the side chain of Asp158 in wild-type papain are indirectly important for enzyme activity. When CBZ-Phe-Arg-MCA is used as a substrate, the (kcat/KM)obs values at pH 6.5 are 3650 and 494 M-1 s-1 for Asp158Gly and Asp158Ala, respectively, as compared to 119,000 M-1 s-1 for papain. Results with the Asp158Glu mutant suggest that the side chain of Glu moves closer to the active site and cannot form hydrogen bonds similar to those involving Asp158 in papain. From the four mutations introduced at position 158 in papain, we can conclude that it is not the charge but the hydrogen-bonding interactions involving the side chain of Asp158 that contribute the most to the stabilization of the thiolate-imidazolium ion pair in papain. However, the charge and the hydrogen bonds of Asp158 both contribute to the intrinsic activity of the enzyme.     

Mutational analysis of the interaction between active site residues and the loop region in mammalian purple acid phosphatases

Mammalian purple acid phosphatases (PAPs) can be divided into two groups, which exhibit distinct spectroscopic and kinetics properties: PAPs that consist of a single 36 kDa polypeptide, and PAPs that have undergone limited proteolysis to give two fragments with masses of 16 and 20 kDa, respectively. Proteolysis results in an increase in enzymatic activity, an increase in the optimal pH for activity, and a change in the g(z)() value of the characteristic EPR spectrum of the mixed-valence binuclear iron center. It has been proposed that these changes are due to the loss of interactions between Asp146 in an exposed loop region and active site residues upon proteolysis. In the present study, site-directed mutagenesis of Asp146 in recombinant rat bone PAP (recRPAP) has confirmed this hypothesis. Conversion of Asp146 into Ala, which eliminates the interaction of the side chain with the active site, resulted in an enzyme with properties typical of PAPs isolated in proteolytically cleaved forms. The Asp146Asn and Asp146Glu mutants were also prepared and examined to assess the effects of altered electrostatic interactions and side-chain length. Limited proteolysis of all three mutant enzymes with cathepsin L resulted in a significant increase in catalytic activity. Thus, although the interaction between Asp146 and (an) active site residue(s) is the major factor responsible for the low catalytic activity of uncleaved PAPs, other interactions are also important. Since both p-nitrophenyl phosphate and osteopontin, a potential in vivo substrate, show the same level of activation, the observed increase in catalytic activity upon proteolysis is likely to be due to electrostatic rather than steric effects. EPR spectra of FeZn-recRPAP before and after cleavage by cathepsin L suggest that cleavage primarily affects the divalent metal site. The observation that pK(es,1) is also sensitive to changes at the divalent site is consistent with the proposal that the nucleophilic hydroxide is that bridging the divalent and trivalent metals.     

Contributions of a hydrogen bond/salt bridge network to the stability of secondary and tertiary structure in lambda repressor.

In the N-terminal domain of lambda repressor, the Asp 14 side chain forms an intrahelical, hydrogen bond/salt bridge with the Arg 17 side chain and a tertiary hydrogen bond with the Ser 77 side chain. By measuring the stabilities to urea denaturation of the wild-type N-terminal domain and variants containing single, double, and triple alanine substitutions at positions 14, 17, and 77, the side-chain interaction energies, the coupling energy between interactions, and the intrinsic effects of each wild-type side chain on protein stability have been estimated. These studies indicate that the Asp 14-Arg 17 and Asp 14-Ser 77 interactions are stabilizing by roughly 0.8 and 1.5 kcal/mol, respectively, but that Asp 14, by itself, is destabilizing by roughly 0.9 kcal/mol. We also show that a peptide model of alpha-helix 1, which contains Asp 14 and Arg 17, forms a reasonably stable, monomeric helix in solution and responds to alanine mutations at positions 14 and 17 in the fashion expected from the intact protein studies. These studies suggest that it is possible to view the stability effects of mutations in intact proteins in a hierarchical fashion, with the stability of units of secondary structure being distinguishable from the stability of tertiary structure.     

Analysis of the interaction between charged side chains and the alpha-helix dipole using designed thermostable mutants of phage T4 lysozyme

It was shown previously that the introduction of a negatively charged amino acid at the N-terminus of an alpha-helix could increase the thermostability of phage T4 lysozyme via an electrostatic interaction with the "helix dipole" [Nicholson, H., Becktel, W. J., & Matthews, B. W. (1988) Nature 336, 651-656]. The prior report focused on the two stabilizing substitutions Ser 38----Asp (S38D) and Asn 144----Asp (N144D). Two additional examples of stabilizing mutants, T109D and N116D, are presented here. Both show the pH-dependent increase in thermal stability expected for the interaction of an aspartic acid with an alpha-helix dipole. Control mutants were also constructed to further characterize the nature of the interaction with the alpha-helix dipole. High-resolution crystal structure analysis was used to determine the nature of the interaction of the substituted amino acids with the end of the alpha-helix in both the primary and the control mutants. Control mutant S38N has stability essentially the same as that of wild-type lysozyme but hydrogen bonding similar to that of the stabilizing mutant S38D. This confirms that it is the electrostatic interaction between Asp 38 and the helix dipole, rather than a change in hydrogen-bonding geometry, that gives enhanced stability. Structural and thermodynamic analysis of mutant T109N provide a similar control for the stabilizing replacement T109D. In the case of mutant N116D, there was concern that the enhanced stability might be due to a favorable salt-bridge interaction between the introduced aspartate and Arg 119, rather than an interaction with the alpha-helix dipole. The additivity of the stabilities of N116D and R119M seen in the double mutant N116D/R119M indicates that favorable interactions are largely independent of residue 119. As a further control, Asp 92, a presumed helix-stabilizing residue in wild-type lysozyme, was replaced with Asn. This decreased the stability of the protein in the manner expected for the loss of a favorable helix dipole interaction. In total, five mutations have been identified that increase the thermostability of T4 lysozyme and appear to do so by favorable interactions with alpha-helix dipoles. As measured by the pH dependence of stability, the strength of the electrostatic interaction between the charged groups studied here and the helix dipole ranges from 0.6 to 1.3 kcal/mol in 150 mM KCl. In the case of mutants S38D and N144H, NMR titration was used to measure the pKa's of Asp 38 and His 144 in the folded structures.(ABSTRACT TRUNCATED AT 400 WORDS)     

Hemoglobin providence. Functional consequences of two alterations of the 2,3-diphosphoglycerate binding site at position beta 82.

Position beta 82 in human hemoglobin (Hb) is normally occupied by lysine, a positively charged residue that is involved in the binding of anionic cofactors. This residue is substituted by a neutral residue in Hb Providence Asn and by a negatively charged residue in Hb Providence Asp. Hb Providence Asp shows more differences from Hb A than does Hb Providence Asn in studies of the kinetics and equilibria of ligand binding. For both forms, homotropic (cooperative) interactions are normal with n values of 2.5 to 2.7, while heterotropic (pH and anion) interactions are reduced greatly. The reduction in anion sensitivity is attributed to the absence of a positive residue at position beta 82. Reduction in pH sensitivity may be due to a ligand-linked change in the pK of a neighboring residue, beta 143 histidine, which normally is not a Bohr group. This change in pK would act in opposition to the normal Bohr effect. Reduction in the net positive charge of the central cavity has a further consequence. Relative to Hb A, both Hb Providence Asn and Hb Providence Asp show decreased oxygen affinities at neutral pH in the absence of cofactors. This suggests that in Hb A the binding of anionic cofactors directly influences the oxygen affinity by neutralizing the charged groups of the diphosphoglycerate binding site and thus stabilizing the low affinity (T) conformation. From pH 6 to 9 in the presence of 1 M NaCl, where all the charged groups may be masked, the oxygen-binding properties of Hb A and the Hb Providence mutants are identical. Moreover, subunit dissociation of the liganded Hb Providence mutants appears to be increased, as is known to occur for Hb A in the presence of high salt. The results obtained with Hb Providence Asn and Hb Providence Asp illustrate how single amino acid substitutions can modify hemoglobins' pH and anion interactions without altering cooperative interactions between subunits. The alteration in cofactor effects observed with these mutants also illustrates differences between the allosteric effects induced by organic and inorganic anions.     

Interactions of trifluroethanol with the Trp-cage peptide

It has been suggested that aggregation of fluorinated alcohols in water solutions is involved with the abilities of these alcohols to provoke conformational changes in peptides and proteins. The extent of fluoroalcohol aggregation depends on the degree of fluorination: hexafluoroisopropanol (HFIP) is more extensively aggregated than is TFE. We previously described a study of the interactions of HFIP with the peptide Trp-cage and provided evidence for the formation of long-lived complexes between this fluoroalcohol and the peptide. In the present work, we have examined the interactions of the less-fluorinated TFE with Trp-cage, in order to probe the role of fluoroalcohol aggregation in the phenomena observed. Intermolecular (1)H{(19)F} nuclear Overhauser effects arising from interactions of TFE with the hydrogens of the peptide in a solution containing 42% TFE were determined at sample temperatures from 5 to 45 degrees C. It is shown that the folded state of the peptide under these conditions is essentially the same as that observed in water and in 30% HFIP-water. The observed peptide-solvent NOEs indicate formation of complexes of Trp-cage with TFE that persist for times of the order of 1 ns. The interactions leading to complexes with TFE are somewhat weaker than those involved in complex formation with HFIP. There are no indications that the aggregation of fluoroalcohol is a necessary concomitant of the interactions of TFE or HFIP with Trp-cage. Rather, the stronger and more long-lived interactions of HFIP with Trp-cage appear to be primarily the result of the greater hydrogen-bonding ability and hydrophobicity of this fluoroalcohol.     

Effect of a single aspartate on helix stability at different positions in a neutral alanine-based peptide.

A single aspartate residue has been placed at various positions in individual peptides for which the alanine-based reference peptide is electrically neutral, and the helix contents of the peptides have been measured by circular dichroism. The dependence of peptide helix content on aspartate position has been used to determine the helix propensity (s-value). Both the charged (Asp-) and uncharged (Asp0) forms of the aspartate residue are strong helix breakers and have identical s-values of 0.29 at 0 degree C. The interaction of Asp- with the helix dipole affects helix stability at positions throughout the helix, not only near the N-terminus, where the interaction is helix stabilizing, and the C-terminus, where it is destabilizing. Comparison of the helix contents at acidic pH (Asp0) and at neutral pH (Asp-) shows that the charge-helix dipole interaction is screened slowly with increasing NaCl concentration, and screening is not complete even at 4.8 M NaCl. Lastly, a helix-stabilizing hydrogen-bond interaction between glutamine and aspartate (spacing i, i + 4) has been found. This side-chain interaction is specific for both the orientation and spacing of the glutamine and aspartate residues and is resistant to screening by NaCl.     

Characterisation of the effect of electrostatic interaction on the structure of Trp-cage using molecular dynamics simulation

A mutated variant of 20 amino acid miniprotein Trp-cage (TC5b), called TC5c (Asp9 replaced by Asn9), was designed to demonstrate the effect of a salt bridge. As a result of strong electrostatic interaction, the distance distribution between Asp9 and Arg16 exhibited a larger probability in the range of the salt bridge for TC5b compared to TC5c. The probability of α-helix formation for residues 3–8, as well as for residues 11–14, was high for TC5b. The salt bridge formation between Asp9 and Arg16 in TC5b was indicated by (a) a strong correlation of their distance of separation with the subtended angle with the centre and (b) a step decrease in the distance between Gly11O and Arg16H at 12 ns. Replica exchange molecular dynamics simulation at different temperatures in the range of 270–590 K indicated that the average distance between Asp9 and Arg16, end-to-end distance, root mean square deviation with respect to a reference NMR structure of TC5b did not change significantly with temperature below 370 K for TC5b and increased at higher temperatures. These values were higher for TC5c for the whole temperature range, with their rate of increase with temperature being higher below 370 K.     

Aspartate-histidine interaction in the retinal schiff base counterion of the light-driven proton pump of Exiguobacterium sibiricum

One of the distinctive features of eubacterial retinal-based proton pumps, proteorhodopsins, xanthorhodopsin, and others, is hydrogen bonding of the key aspartate residue, the counterion to the retinal Schiff base, to a histidine. We describe properties of the recently found eubacterium proton pump from Exiguobacterium sibiricum (named ESR) expressed in Escherichia coli, especially features that depend on Asp-His interaction, the protonation state of the key aspartate, Asp85, and its ability to accept a proton from the Schiff base during the photocycle. Proton pumping by liposomes and E. coli cells containing ESR occurs in a broad pH range above pH 4.5. Large light-induced pH changes indicate that ESR is a potent proton pump. Replacement of His57 with methionine or asparagine strongly affects the pH-dependent properties of ESR. In the H57M mutant, a dramatic decrease in the quantum yield of chromophore fluorescence emission and a 45 nm blue shift of the absorption maximum with an increase in the pH from 5 to 8 indicate deprotonation of the counterion with a pK(a) of 6.3, which is also the pK(a) at which the M intermediate is observed in the photocycle of the protein solubilized in detergent [dodecyl maltoside (DDM)]. This is in contrast with the case for the wild-type protein, for which the same experiments show that the major fraction of Asp85 is deprotonated at pH >3 and that it protonates only at low pH, with a pK(a) of 2.3. The M intermediate in the wild-type photocycle accumulates only at high pH, with an apparent pK(a) of 9, via deprotonation of a residue interacting with Asp85, presumably His57. In liposomes reconstituted with ESR, the pK(a) values for M formation and spectral shifts are 2-3 pH units lower than in DDM. The distinctively different pH dependencies of the protonation of Asp85 and the accumulation of the M intermediate in the wild-type protein versus the H57M mutant indicate that there is strong Asp-His interaction, which substantially lowers the pK(a) of Asp85 by stabilizing its deprotonated state.     

Optimal salt bridge for Trp-cage stabilization

Gai and co-workers [Bunagan, M. R., et al. (2006) J. Phys. Chem. B 110, 3759-3763] reported computational design studies suggesting that a D9E mutation would stabilize the Trp-cage. Experimental studies for this mutation were reported in 2008 [Hudaky, P., et al. (2008) Biochemistry 47, 1007-1016]; the authors suggested that [D9E]-TC5b presented a more compact and melting resistant structure because of the "optimal distance between the two sides of the molecule". Nonetheless, the authors reported essentially the same circular dichroism (CD) melting temperature, 38 ± 0.3 °C, for TC5b and its [D9E] mutant. In this study, a more stable Trp-cage, DAYAQ WLKDG GPSSG RPPPS, was examined by nuclear magnetic resonance and CD with the following mutations: [D9E], [D9R,R16E], [R16O], [D9E,R16O], [R16K], and [D9E,R16K]. Of these, the [D9E] mutant displayed the smallest acidification-induced change in the apparent T(m). In analogy to the prior study, the CD melts of TC10b and its [D9E] mutant were, however, very similar; all of the other mutations were significantly fold destabilizing by all measures. A detailed analysis indicates that the original D9-R16 salt bridge is optimal with regard to fold cooperativity and fold stabilization. Evidence of salt bridge formation is also provided for a swapped pair, the [D9R,R16E] mutant. Model systems reveal that an ionized aspartate at the C-terminus of a helix significantly decreases intrinsic helicity, a requirement for Trp-cage fold stability. The CD evidence that was cited as supporting increased fold stability for [D9E]-TC5b at higher temperatures appears to be a reflection of increased helix stability in both the folded and unfolded states rather than a more favorable salt bridge. Our study also provides evidence of other Trp-cage stabilizing roles of the R16 side chain.     

Identifying stabilizing key residues in proteins using interresidue interaction energy matrix

We are proposing an interresidue interaction energy map (IEM)--a new tool for protein structure analysis and protein bioinformatics. This approach employs the sum of pair-wise interaction energies of a particular residue as a measure of its structural importance. We will show that the IEM can serve as a means for identifying key residues responsible for the stability of a protein. Our method can be compared with the interresidue contact map but has the advantage of weighting the contacts by the stabilization energy content which they bring to the protein structure. For the theoretical adjustment of the proposed method, we chose the Trp-cage mini protein as a model system to compare a spectrum of computational methods ranging from the ab initio MP2 level through the DFT method to empirical force-field methods. The IEM method correctly identifies Tryptophane 6 as the key residue in the Trp-cage. The other residues with the highest stabilizing contributions correspond to the structurally important positions in the protein. We have further tested our method on the Trp2Cage miniprotein--a P12W mutant of the Trp-cage and on two proteins from the rubredoxin family that differ in their thermostability. Our method correctly identified the thermodynamically more stable variants in both cases and therefore can also be used as a tool for the relative measurement of protein stability. Finally, we will point out the important role played by dispersion energy, which contributes significantly to the total stabilization energy and whose role in aromatic pairs is clearly dominant. Surprisingly, the dispersion energy plays an even more important role in the interaction of prolines with aromatic systems.     

Effect of the N1 residue on the stability of the alpha-helix for all 20 amino acids

N1 is the first residue in an alpha-helix. We have measured the contribution of all 20 amino acids to the stability of a small helical peptide CH(3)CO-XAAAAQAAAAQAAGY-NH(2) at the N1 position. By substituting every residue into the N1 position, we were able to investigate the stabilizing role of each amino acid in an isolated context. The helix content of each of the 20 peptides was measured by circular dichroism (CD) spectroscopy. The data were analyzed by our modified Lifson-Roig helix-coil theory, which includes the n1 parameter, to find free energies for placing a residue into the N1 position. The rank order for free energies is Asp(-), Ala > Glu(-) > Glu(0) > Trp, Leu, Ser > Asp(0), Thr, Gln, Met, Ile > Val, Pro > Lys(+), Arg, His(0) > Cys, Gly > Phe > Asn, Tyr, His(+). N1 preferences are clearly distinct from preferences for the preceding N-cap and alpha-helix interior. pK(a) values were measured for Asp, Glu, and His, and protonation-free energies were calculated for Asp and Glu. The dissociation of the Asp proton is less favorable than that of Glu, and this reflects its involvement in a stronger stabilizing interaction at the N terminus. Proline is not energetically favored at the alpha-helix N terminus despite having a high propensity for this position in crystal structures. The data presented are of value both in rationalizing mutations at N1 alpha-helix sites in proteins and in predicting the helix contents of peptides.     

Diverse gating in K channels: Differential role of the pore-helix glutamate in stabilizing the channel pore

The selectivity filter and adjacent regions in the bacterial KcsA and inwardly rectifying K⁺ (Kir) channels reveal significant conformational changes that cause the channel pore to transition from an activated to inactive state (C-type inactivation) once the channel is open. The meshwork of residues stabilizing the pore of KcsA involves Glu71–Asp80 carboxyl–carboxylate interaction ‘behind’ the selectivity filter. Interestingly, the Kir channels do not have this exact interaction, but instead have a Glu–Arg salt bridge where the Glu is in the same position but the Arg is one position N-terminal compared to the Asp in KcsA. Also, the Kir channels lack the Trp that hydrogen bonds to Asp80 in KcsA. Here, the sequence and structural information are combined to understand the dissimilarity in the role of the pore-helix Glu in stabilizing the pore structure in KcsA and Kir channels. This review illustrates that although Glu is quite conserved among both types of channels, the network of interactions is not translatable from one channel to the other; thereby suggesting a unique phenomenon of diverse gating patterns in K⁺ channels.     

Polarity Alteration of a Calcium Site Induces a Hydrophobic Interaction Network and Enhances Cel9A Endoglucanase Thermostability

Structural calcium sites control protein thermostability and activity by stabilizing native folds and changing local conformations. Alicyclobacillus acidocaldarius survives in thermal-acidic conditions and produces an endoglucanase Cel9A (AaCel9A) which contains a calcium-binding site (Ser465 to Val470) near the catalytic cleft. By superimposing the Ca²⁺-free and Ca²⁺-bounded conformations of the calcium site, we found that Ca²⁺ induces hydrophobic interactions between the calcium site and its nearby region by driving a conformational change. The hydrophobic interactions at the high-B-factor region could be enhanced further by replacing the surrounding polar residues with hydrophobic residues to affect enzyme thermostability and activity. Therefore, the calcium-binding residue Asp468 (whose side chain directly ligates Ca²⁺), Asp469, and Asp471 of AaCel9A were separately replaced by alanine and valine. Mutants D468A and D468V showed increased activity compared with those of the wild type with 0 mM or 10 mM Ca²⁺ added, whereas the Asp469 or Asp471 substitution resulted in decreased activity. The D468A crystal structure revealed that mutation D468A triggered a conformational change similar to that induced by Ca²⁺ in the wild type and developed a hydrophobic interaction network between the calcium site and the neighboring hydrophobic region (Ala113 to Ala117). Mutations D468V and D468A increased 4.5°C and 5.9°C, respectively, in melting temperature, and enzyme half-life at 75°C increased approximately 13 times. Structural comparisons between AaCel9A and other endoglucanases of the GH9 family suggested that the stability of the regions corresponding to the AaCel9A calcium site plays an important role in GH9 endoglucanase catalysis at high temperature.     

Effect of lysine side chain length on intra-helical glutamate--lysine ion pairing interactions

Ion-pairing interactions are important for protein stabilization. Despite the apparent electrostatic nature of these interactions, natural positively charged amino acids Lys and Arg have multiple methylenes linking the charged functionality to the backbone. Interestingly, the amino acids Lys and Orn have positively charged side chains that differ by only one methylene. However, only Lys is encoded and incorporated into proteins. To investigate the effect of side chain length of Lys on ion-pairing interactions, a series of 12 monomeric alpha-helical peptides containing potential Glu-Xaa (i, i+3), (i, i+4) and (i, i+5) (Xaa = Lys, Orn, Dab, Dap) interactions were studied by circular dichroism (CD) spectroscopy at pH 7 and 2. At pH 7, no Glu-Xaa (i, i+5) interaction was observed, regardless of the Xaa side chain length. Furthermore, only Lys was capable of supporting Glu-Xaa (i, i+3) interactions, whereas any Xaa side chain length supported Glu-Xaa (i, i+4) interactions. Side chain conformational analysis by molecular mechanics calculations showed that the side chain length of Lys enables the Glu-Xaa (i, i+3) interaction with lower energy conformations compared to residues with side chain lengths shorter than that of Lys. Furthermore, these calculated low energy conformers were consistent with conformations of intra-helical Glu-Lys salt bridges in a non-redundant protein structure database. Importantly, the CD spectra for peptides with Glu-Lys interactions did not alter significantly upon changing the pH because of a greater contribution to these interactions by forces other than electrostatics. Incorporating side chains just one methylene shorter (Orn) resulted in significant pH dependence or lack of interaction, suggesting that nature has chosen Lys to form durable interactions with negatively charged functional groups.     

Exploring the environmental preference of weak interactions in (alpha/beta)8 barrel proteins

The environmental preference for the occurrence of noncanonical hydrogen bonding and cation-pi interactions, in a data set containing 71 nonredundant (alpha/beta)(8) barrel proteins, with respect to amino acid type, secondary structure, solvent accessibility, and stabilizing residues has been performed. Our analysis reveals some important findings, which include (a) higher contribution of weak interactions mediated by main-chain atoms irrespective of the amino acids involved; (b) domination of the aromatic amino acids among interactions involving side-chain atoms; (c) involvement of strands as the principal secondary structural unit, accommodating cross strand ion pair interaction and clustering of aromatic amino acid residues; (d) significant contribution to weak interactions occur in the solvent exposed areas of the protein; (e) majority of the interactions involve long-range contacts; (f) the preference of Arg is higher than Lys to form cation-pi interaction; and (g) probability of theoretically predicted stabilizing amino acid residues involved in weak interaction is higher for polar amino acids such as Trp, Glu, and Gln. On the whole, the present study reveals that the weak interactions contribute to the global stability of (alpha/beta)(8) TIM-barrel proteins in an environment-specific manner, which can possibly be exploited for protein engineering applications.