Estimating the contribution of engineered surface electrostatic interactions to protein stability by using double-mutant cycles

Coulombic interactions between charges on the surface of proteins contribute to stability. It is difficult, however, to estimate their importance by protein engineering methods because mutation of one residue in an ion pair alters the energetics of many interactions in addition to the coulombic energy between the two components. We have estimated the interaction energy between two charged residues, Asp-12 and Arg-16, in an alpha-helix on the surface of a barnase mutant by invoking a double-mutant cycle involving wild-type enzyme (Asp-12, Thr-16), the single mutants Thr----Arg-16 and Asp----Ala-12, and the double mutant Asp----Ala-12, Thr----Arg-16. The changes in free energy of unfolding of the single mutants are not additive because of the coulombic interaction energy. Additivity is restored at high concentrations of salt that shield electrostatic interactions. The geometry of the ion pair in the mutant was assumed to be the same as that in the highly homologous ribonuclease from Bacillus intermedius, binase, which has Asp-12 and Arg-16 in the native enzyme. The ion pair does not form a hydrogen-bonded salt bridge, but the charges are separated by 5-6 A. The mutant barnase containing the ion pair Asp-12/Arg-16 is more stable than wild type by 0.5 kcal/mol, but only a part of the increased stability is attributable to the electrostatic interaction. We present a formal analysis of how double-mutant cycles can be used to measure the energetics of pairwise interactions.     

Effect of surface electrostatic interactions on the stability and folding of formate dehydrogenase from Candida methylica

NAD⁺-dependent formate dehydrogenase (FDH-EC 1.2.1.2) is an important enzyme to regenerate valuable NADH required by NAD⁺-dependent oxidoreductases in enzyme catalysis. The limitation in the thermostability of FDH enzyme is a crucial problem for development of biotechnological and industrial processes, despite of its advantages. In this study, to investigate the contribution of surface electrostatic interaction to the thermostability of FDH from Candida methylica (cmFDH) N187E, H13E, Q105R, N300E, N147R N300E/N147R, N187E/Q105R, N187E/N147R,Y160R, Y302R, Y160E and Y302E mutants were designed using a homology model of cmFDH based on Candida boidinii (cb) by considering electrostatic interactions on the protein surface. The effects of site-specific engineering on the stability of this molecule was analyzed according to minimal model of folding and assembly reaction and deduced equilibrium properties of the native system with respect to its thermal and denaturant sensitivities. It was observed that mutations did not change the unfolding pattern of native cmFDH and increased numbers of electrostatic interactions can cause either stabilizing or destabilizing effect on the thermostability of this protein. The thermodynamic and kinetic results suggested that except relatively improved mutants, three out of the nine single mutations increased the melting temperature of cmFDH enzyme.     

The effect of protein stability on protein-monoglyceride interactions

We have previously shown that proteins such as beta-lactoglobulin and lysozyme insert into monoglyceride monolayers and are able to induce an L(beta) to coagel phase transition in monoglyceride bilayers. These studies gave a first indication that protein stability could be an important factor for these interactions. This study therefore aims at further investigating the potential role of protein stability on protein-monoglyceride interactions. To this end we studied the interaction of stable and destabilized alpha-lactalbumin with monostearoylglycerol. Our results show that protein stability is important for the insertion of proteins into a monostearoylglycerol monolayer, such that the lower the stability of the protein the better the protein inserts. In marked contrast to beta-lactoglobulin and lysozyme we found that destabilized alpha-lactalbumin does not induce the L(beta) to coagel phase transition in monoglyceride bilayers. We propose that this is due to an increased surface coverage by the protein which could result from the unfolding of the protein upon binding to the interface.     

Role of electrostatic interactions for the stability and folding behavior of cold shock protein

The impacts of three charged‐residue‐involved mutations, E46A, R3E, and R3E/L66E, on the thermostability and folding behavior of the cold shock protein from the themophile Bacillus caldolyticus (Bc‐Csp) were investigated by using a modified Gō‐like model, in which the nonspecific electrostatic interactions of charged residues were taken into account. Our simulation results show that the wild‐type Bc‐Csp and its three mutants are all two‐sate folders, which is consistent with the experimental observations. It is found that these three mutations all lead to a decrease of protein thermodynamical stability, and the effect of R3E mutation is the strongest. The lower stability of these three mutants is due to the increase of the enthalpy of the folded state and the entropy of the unfolded state. Using this model, we also studied the folding kinetics and the folding/unfolding pathway of the wild‐type Bc‐Csp as well as its three mutants and then discussed the effects of electrostatic interactions on the folding kinetics. The results indicate that the substitutions at positions 3 and 46 largely decrease the folding kinetics, whereas the mutation of residue 66 only slightly decreases the folding rate. This result agrees well with the experimental observations. It is also found that these mutations have little effects on the folding transition state and the folding pathway, in which the N‐terminal β sheet folds earlier than the C‐terminal region. We also investigated the detailed unfolding pathway and found that it is really the reverse of the folding pathway, providing the validity of our simulation results. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

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

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

Acetylation of C-terminal lysines modulates protein turnover and stability of Connexin-32

Background

The gap junction protein, Connexin32 (Cx32), is expressed in various tissues including liver, exocrine pancreas, gastrointestinal epithelium, and the glia of the central and peripheral nervous system. Gap junction-mediated cell-cell communication and channel-independent processes of Cx32 contribute to the regulation of physiological and cellular activities such as glial differentiation, survival, and proliferation; maintenance of the hepatic epithelium; and axonal myelination. Mutations in Cx32 cause X-linked Charcot–Marie–Tooth disease (CMT1X), an inherited peripheral neuropathy. Several CMT1X causing mutations are found in the cytoplasmic domains of Cx32, a region implicated in the regulation of gap junction assembly, turnover and function. Here we investigate the roles of acetylation and ubiquitination in the C-terminus on Cx32 protein function. Cx32 protein turnover, ubiquitination, and response to deacetylase inhibitors were determined for wild-type and C-terminus lysine mutants using transiently transfected Neuro2A (N2a) cells.

Results

We report here that Cx32 is acetylated in transfected N2a cells and that inhibition of the histone deacetylase, HDAC6, results in an accumulation of Cx32. We identified five lysine acetylation targets in the C-terminus. Mutational analysis demonstrates that these lysines are involved in the regulation of Cx32 ubiquitination and turnover. While these lysines are not required for functional Cx32 mediated cell-cell communication, BrdU incorporation studies demonstrate that their relative acetylation state plays a channel-independent role in Cx32-mediated control of cell proliferation.

Conclusion

Taken together these results highlight the role of post translational modifications and lysines in the C-terminal tail of Cx32 in the fine-tuning of Cx32 protein stability and channel-independent functions.

     

A study on the effect of surface lysine to arginine mutagenesis on protein stability and structure using green fluorescent protein

Two positively charged basic amino acids, arginine and lysine, are mostly exposed to protein surface, and play important roles in protein stability by forming electrostatic interactions. In particular, the guanidinium group of arginine allows interactions in three possible directions, which enables arginine to form a larger number of electrostatic interactions compared to lysine. The higher pKa of the basic residue in arginine may also generate more stable ionic interactions than lysine. This paper reports an investigation whether the advantageous properties of arginine over lysine can be utilized to enhance protein stability. A variant of green fluorescent protein (GFP) was created by mutating the maximum possible number of lysine residues on the surface to arginines while retaining the activity. When the stability of the variant was examined under a range of denaturing conditions, the variant was relatively more stable compared to control GFP in the presence of chemical denaturants such as urea, alkaline pH and ionic detergents, but the thermal stability of the protein was not changed. The modeled structure of the variant indicated putative new salt bridges and hydrogen bond interactions that help improve the rigidity of the protein against different chemical denaturants. Structural analyses of the electrostatic interactions also confirmed that the geometric properties of the guanidinium group in arginine had such effects. On the other hand, the altered electrostatic interactions induced by the mutagenesis of surface lysines to arginines adversely affected protein folding, which decreased the productivity of the functional form of the variant. These results suggest that the surface lysine mutagenesis to arginines can be considered one of the parameters in protein stability engineering.     

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.     

State-dependent electrostatic interactions of S4 arginines with E1 in S2 during Kv7.1 activation

The voltage-sensing domain of voltage-gated channels is comprised of four transmembrane helices (S1–S4), with conserved positively charged residues in S4 moving across the membrane in response to changes in transmembrane voltage. Although it has been shown that positive charges in S4 interact with negative countercharges in S2 and S3 to facilitate protein maturation, how these electrostatic interactions participate in channel gating remains unclear. We studied a mutation in Kv7.1 (also known as KCNQ1 or KvLQT1) channels associated with long QT syndrome (E1K in S2) and found that reversal of the charge at E1 eliminates macroscopic current without inhibiting protein trafficking to the membrane. Pairing E1R with individual charge reversal mutations of arginines in S4 (R1–R4) can restore current, demonstrating that R1–R4 interact with E1. After mutating E1 to cysteine, we probed E1C with charged methanethiosulfonate (MTS) reagents. MTS reagents could not modify E1C in the absence of KCNE1. With KCNE1, (2-sulfonatoethyl) MTS (MTSES)⁻ could modify E1C, but [2-(trimethylammonium)ethyl] MTS (MTSET)⁺ could not, confirming the presence of a positively charged environment around E1C that allows approach by MTSES⁻ but repels MTSET⁺. We could change the local electrostatic environment of E1C by making charge reversal and/or neutralization mutations of R1 and R4, such that MTSET⁺ modified these constructs depending on activation states of the voltage sensor. Our results confirm the interaction between E1 and the fourth arginine in S4 (R4) predicted from open-state crystal structures of Kv channels and reveal an E1–R1 interaction in the resting state. Thus, E1 engages in electrostatic interactions with arginines in S4 sequentially during the gating movement of S4. These electrostatic interactions contribute energetically to voltage-dependent gating and are important in setting the limits for S4 movement.     

Molecular interactions in protein crystals: solvent accessible surface and stability

The accessible surface areas of 58 monomeric and dimeric proteins, when measured in the crystalline environment, are found to be simply related to molecular weight. The loss of accessible surface when the proteins go from a free to their crystalline environment is well defined, implying that the hydrophobic interaction, which has been found to contribute to protein folding and stability in living systems, also contributes to protein crystal stability.     

Mechanism of thermostabilization in a designed cold shock protein with optimized surface electrostatic interactions

Using computational and sequence analysis of bacterial cold shock proteins, we designed a protein (CspB-TB) that has the core residues of mesophilic protein from Bacillus subtilis(CspB-Bs) and altered distribution of surface charged residues. This designed protein was characterized by circular dichroism spectroscopy, and found to have secondary and tertiary structure similar to that of CspB-Bs. The activity of the CspB-TB protein as measured by the affinity to a single-stranded DNA (ssDNA) template at 25 degrees C is somewhat higher than that of CspB-Bs. Furthermore, the decrease in the apparent binding constant to ssDNA upon increase in temperature is much more pronounced for CspB-Bs than for CspB-TB. Temperature-induced unfolding (as monitored by differential scanning calorimetry and circular dichroism spectroscopy) and urea-induced unfolding experiments were used to compare the stabilities of CspB-Bs and CspB-TB. It was found that CspB-TB is approximately 20 degrees C more thermostable than CspB-Bs. The thermostabilization of CspB-TB relative to CspB-Bs is achieved by decrease in the enthalpy and entropy of unfolding without affecting their temperature dependencies, i.e. these proteins have similar heat capacity changes upon unfolding. These changes in the thermodynamic parameters result in the global stability function, i.e. Gibbs energy, deltaG(T), that is shifted to higher temperatures with only small changes in the maximum stability. Such a mechanism of thermostabilization, although predicted from the basic thermodynamic considerations, has never been identified experimentally.     

Surface electrostatic interactions contribute little of stability of barnase

Electrostatic interactions are believed to play an important role in stabilizing the native structure of proteins. We have quantified the contribution to stability of an interaction between two oppositely charged side-chains on the surface of barnase. Using site-directed mutagenesis, glutamate 28 and lysine 32 were introduced onto the solvent-accessible side of the second alpha-helix in barnase. These two residues are separated by one turn of the helix, and so are ideally situated for their opposite charges to interact. Double mutant cycle analysis reveals that the interaction between Glu28 and Lys32 contributes only approximately 0.2 kcal/mol to stability of the protein. All other interactions between exposed charged side-chains in barnase examined so far also contribute little to stability. We explain this low value by their location on the surface, rather than in the interior, of the protein.     

In silico discovery and experimental validation of new protein-protein interactions

We introduce a framework for predicting novel protein-protein interactions (PPIs), based on Fisher's method for combining probabilities of predictions that are based on different data sources, such as the biomedical literature, protein domain and mRNA expression information. Our method compares favorably to our previous method based on text-mining alone and other methods such as STRING. We evaluated our algorithms through the prediction of experimentally found protein interactions underlying Muscular Dystrophy, Huntington's Disease and Polycystic Kidney Disease, which had not yet been recorded in protein-protein interaction databases. We found a 1.74-fold increase in the mean average prediction precision for dysferlin and a 3.09-fold for huntingtin when compared to STRING. The top 10 of predicted interaction partners of huntingtin were analysed in depth. Five were identified previously, and the other five were new potential interaction partners. The full matrix of human protein pairs and their prediction scores are available for download. Our framework can be extended to predict other types of relationships such as proteins in a complex, pathway or related disease mechanisms.     

In silico identification of functional protein interfaces

Proteins perform many of their biological roles through protein-protein, protein-DNA or protein-ligand interfaces. The identification of the amino acids comprising these interfaces often enhances our understanding of the biological function of the proteins. Many methods for the detection of functional interfaces have been developed, and large-scale analyses have provided assessments of their accuracy. Among them are those that consider the size of the protein interface, its amino acid composition and its physicochemical and geometrical properties. Other methods to this effect use statistical potential functions of pairwise interactions, and evolutionary information. The rationale of the evolutionary approach is that functional and structural constraints impose selective pressure; hence, biologically important interfaces often evolve at a slower pace than do other external regions of the protein. Recently, an algorithm, Rate4Site, and a web-server, ConSurf (http://consurf.tau.ac.il/), for the identification of functional interfaces based on the evolutionary relations among homologous proteins as reflected in phylogenetic trees, were developed in our laboratory. The explicit use of the tree topology and branch lengths makes the method remarkably accurate and sensitive. Here we demonstrate its potency in the identification of the functional interfaces of a hypothetical protein, the structure of which was determined as part of the international structural genomics effort. Finally, we propose to combine complementary procedures, in order to enhance the overall performance of methods for the identification of functional interfaces in proteins.     

Perturbation of the stability of amyloid fibrils through alteration of electrostatic interactions

The self-assembly of proteins and peptides into polymeric amyloid fibrils is a process that has important implications ranging from the understanding of protein misfolding disorders to the discovery of novel nanobiomaterials. In this study, we probe the stability of fibrils prepared at pH 2.0 and composed of the protein insulin by manipulating electrostatic interactions within the fibril architecture. We demonstrate that strong electrostatic repulsion is sufficient to disrupt the hydrogen-bonded, cross-β network that links insulin molecules and ultimately results in fibril dissociation. The extent of this dissociation correlates well with predictions for colloidal models considering the net global charge of the polypeptide chain, although the kinetics of the process is regulated by the charge state of a single amino acid. We found the fibrils to be maximally stable under their formation conditions. Partial disruption of the cross-β network under conditions where the fibrils remain intact leads to a reduction in their stability. Together, these results support the contention that a major determinant of amyloid stability stems from the interactions in the structured core, and show how the control of electrostatic interactions can be used to characterize the factors that modulate fibril stability.     

In silico studies of urease inhibitors to explore ligand-enzyme interactions

In continuation of our previous study on the urease inhibition by a number of chalcones, 2,3-dihydro-1,5-benzothiazepines and 2,3,4,5-tetrahydro-1,5-benzothiazepines, FlexX docking has been exploited to get a deeper insight into the mechanism of their inhibitory action. A comparison of the IC(50) values of the active compounds reveals that, of the three classes of compounds studied, 2,3-dihydro-1,5-benzothiazepines were the most potent urease inhibitors. An in silico examination of these compounds showed that the activity is related to the interaction of ligand with the nickel metallocentre, its interaction with two amino acid residues, Asp224 and Cys322, in addition to the orientation of rings A and B in the catalytic core of the enzyme. The most active compound 2,3-dihydro-1,5-benzothiazepine (4) anchor tightly through a network of interactions with Ni701 and Ni702. This includes a number of hydrogen bonds and hydrophobic contacts with the amino acid residues in its vicinity. For their reduced analogs, the difference in the activity of different diastereomers has been observed to be configuration-dependent. This may be ascribed mainly to the difference in the orientation of ring B of the two stereoisomers and the extent of their interaction with Asp224 and Cys322 present in the catalytic core of the enzyme.     

In silico studies of urease inhibitors to explore ligand-enzyme interactions

In continuation of our previous study on the urease inhibition by a number of chalcones, 2,3-dihydro-1,5-benzothiazepines and 2,3,4,5-tetrahydro-1,5-benzothiazepines, FlexX docking has been exploited to get a deeper insight into the mechanism of their inhibitory action. A comparison of the IC₅₀ values of the active compounds reveals that, of the three classes of compounds studied, 2,3-dihydro-1,5-benzothiazepines were the most potent urease inhibitors. An in silico examination of these compounds showed that the activity is related to the interaction of ligand with the nickel metallocentre, its interaction with two amino acid residues, Asp224 and Cys322, in addition to the orientation of rings A and B in the catalytic core of the enzyme. The most active compound 2,3-dihydro-1,5-benzothiazepine (4) anchor tightly through a network of interactions with Ni701 and Ni702. This includes a number of hydrogen bonds and hydrophobic contacts with the amino acid residues in its vicinity. For their reduced analogs, the difference in the activity of different diastereomers has been observed to be configuration-dependent. This may be ascribed mainly to the difference in the orientation of ring B of the two stereoisomers and the extent of their interaction with Asp224 and Cys322 present in the catalytic core of the enzyme.