In silico study on the effect of surface lysines and arginines on the electrostatic interactions and protein stability

Charged amino acids are mostly exposed on a protein surface, thereby forming a network of interactions with the surrounding amino acids as well as with water. In particular, positively charged arginine and lysine have different side chain geometries and provide a different number of potential electrostatic interactions. This study reports a comparative analysis of the difference in the number of two representative electrostatic interactions, such as salt-bridges and hydrogen bonds, contributed by surface arginine and lysine, as well as their effect on protein stability using molecular modeling and dynamics simulation techniques. Two in silico variants, the R variant with all arginines and the K variant with all lysines on the protein surface, were modeled by mutating all the surface lysines to arginines and the surface arginines to lysines, respectively, for each of the 10 model proteins. A structural comparison of the respective two variants showed that the majority of R variants possessed more salt-bridges and hydrogen bond interactions than the K variants, indicating that arginine provides a higher probability of electrostatic interactions than lysine owing to its side chain geometry. Molecular dynamics simulations of these variants revealed the R variants to be more stable than the K variants at room temperature but this effect was not prominent under protein denaturating conditions, such as 353 and 333 K with 8 M urea. These results suggest that the arginine residues on a protein surface contribute to the protein stability slightly more than lysine by enhancing the electrostatic interactions.     

Deciphering the factors responsible for the stability of a GFP variant resistant to alkaline pH using molecular dynamics simulations

Charged amino acids having ionizable side chains play crucial roles in maintaining the solubility and stability of a protein. These charged amino acids are mostly exposed on protein surface and participate in electrostatic interactions with neighboring charged amino acids as well as with solvent. Therefore, the change in the solvent pH affects the protein stability in most cases. Previously, we reported a GFP variant, GFP14R having 14 surface lysines replaced with arginines, that showed enhanced stability under alkaline pH. Here, we analyzed the factors that contribute to the stability of the GFP14R under alkaline pH quantitatively using molecular dynamics simulations. Protonation state of the charged amino acids of GFP14R and control GFP under neutral pH and alkaline pH were modeled, and molecular dynamics simulations were performed. This comparative analysis revealed that the GFP14R with more arginine frequency on the surface maintained the stability under both pH conditions without much change in their salt-bridge interactions as well as the hydrogen bond interactions with solvent. On the other hand, these interactions were significantly reduced for the control GFP under alkaline pH due to the deprotonated lysine side chains. These results suggest that the advantageous property of arginine over lysine can be considered one of the parameter for the protein stability engineering under alkaline pH conditions.     

Interaction of arginine,lysine, and guanidine with surface residues of lysozyme: implication to protein stability

Additives are widely used to suppress aggregation of therapeutic proteins. However, the molecular mechanisms of effect of additives to stabilize proteins are still unclear. To understand this, we herein perform molecular dynamics simulations of lysozyme in the presence of three commonly used additives: arginine, lysine, and guanidine. These additives have different effects on stability of proteins and have different structures with some similarities; arginine and lysine have aliphatic side chain, while arginine has a guanidinium group. We analyze atomic contact frequencies to study the interactions of the additives with individual residues of lysozyme. Contact coefficient, quantified from contact frequencies, is helpful in analyzing the interactions with the guanidine groups as well as aliphatic side chains of arginine and lysine. Strong preference for contacts to the additives (over water) is seen for the acidic followed by polar and the aromatic residues. Further analysis suggests that the hydration layer around the protein surface is depleted more in the presence of arginine, followed by lysine and guanidine. Molecular dynamics simulations also reveal that the internal dynamics of protein, as indicated by the lifetimes of the hydrogen bonds within the protein, changes depending on the additives. Particularly, we note that the side-chain hydrogen-bonding patterns within the protein differ with the additives, with several side-chain hydrogen bonds missing in the presence of guanidine. These results collectively indicate that the aliphatic chain of arginine and lysine plays a critical role in the stabilization of the protein.     

The effect of denaturants on protein structure

Virtually all studies of the protein-folding reaction add either heat, acid, or a chemical denaturant to an aqueous protein solution in order to perturb the protein structure. When chemical denaturants are used, very high concentrations are usually necessary to observe any change in protein structure. In a solution with such high denaturant concentrations, both the structure of the protein and the structure of the solvent around the protein can be altered. X-ray crystallography is the obvious experimental technique to probe both types of changes. In this paper, we report the crystal structures of dihydrofolate reductase with urea and of ribonuclease A with guanidinium chloride. These two classic denaturants have similar effects on the native structure of the protein. The most important change that occurs is a reduction in the overall thermal factor. These structures offer a molecular explanation for the reduction in mobility. Although the reduction is observed only with the native enzyme in the crystal, a similar decrease in mobility has also been observed in the unfolded state in solution (Makhatadze G, Privalov PL. 1992. Protein interactions with urea and guanidinium chloride: A calorimetric study. J Mol Biol 226:491-505).     

Involvement of basic amino acids in the activity of a nucleic acid helix-destabilizing protein

Under conditions of low ionic strength, ribonuclease A, which binds more tightly to single- than to double-stranded DNA, lowers the melting temperature of DNA helices (Jensen and von Hippel (1976) J. Biol. Chem. 251, 7198-7214). The effects of chemical modification of lysine and arginine residues on the helix-destabilizing properties of this protein have been examined. Removal of the positive charge on the lysine epsilon-amino group, either by maleylation or acetylation, destroys the ability of RNAase A to lower the Tm of poly[d(A-T)]. However, reductive alkylation of these residues, which has not effect on charge, yields derivatives which lower the Tm by only about one-half that seen with unmodified controls. Phenylglyoxalation of arginines can largely remove the Tm-depressing activity of RNAase A. RNAase S, which is produced by cleavage of RNAase A between amino acids 20 and 21, possesses DNA helix-destabilizing activity comparable to that of the parent protein, whereas S-protein (residues 21-124) increases poly[d(A-T)] Tm and S-peptide (1-20) has no effect on Tm. These results suggest that specific location of several basic amino acids situated on the surface of RNAase A is largely responsible for this protein's DNA melting activity.     

Contributions of basic residues to ribosomal protein L11 recognition of RNA

The C-terminal domain of ribosomal protein L11, L11-C76, binds in the distorted minor groove of a helix within a 58 nucleotide domain of 23 S rRNA. To study the electrostatic component of RNA recognition in this protein, arginine and lysine residues have been individually mutated to alanine or methionine residues at the nine sequence positions that are conserved as basic residues among bacterial L11 homologs. In measurements of the salt dependence of RNA-binding, five of these mutants have a reduced value of - partial differentiallog(K(obs))/ partial differentiallog[KCl] as compared to the parent L11-C76 sequence, indicating that these residues interact with the RNA electrostatic field. These five residues are located at the perimeter of the RNA-binding surface of the protein; all five of them form salt bridges with phosphates in the crystal structure of the complex. A sixth residue, Lys47, was found to make an electrostatic contribution to binding when measurements were made at pH 6.0, but not at pH 7.0; its pK in the free protein must be <6.5. The unusual behavior of Lys47 is explained by its burial in the hydrophobic core of the free protein, and unburial in the RNA-bound protein, where it forms a salt bridge with a phosphate. The contributions of these six residues to the electrostatic component of binding are not additive; thus the magnitude of the salt dependence cannot be used to count the number of ionic interactions in this complex. By interacting with irregular features of the RNA backbone, including an S-turn, these basic residues contribute to the specificity of L11 for its target site.     

Critical functional requirement for the guanidinium group of the arginine 41 side chain of human epidermal growth factor as revealed by mutagenic inactivation and chemical reactivation.

In a preliminary study we demonstrated that the formation of the epidermal growth factor (EGF) receptor-ligand complex requires the participation of the highly conserved arginine 41 side chain of the growth factor peptide (Engler, D.A., Montelione, G.T., and Niyogi, S.K. (1990) FEBS Lett. 271, 47-50). In an attempt to gain further insight into the nature of this interaction(s), we used both site-directed mutagenesis and chemical modification reagents to produce human EGF (hEGF) analogues with altered chemical properties of the residue 41 side chain. Eight mutant analogues of hEGF were generated, substituting arginine 41 with lysine, glutamine, isoleucine, tyrosine, glycine, alanine, aspartate, or glutamate. Although each of the mutant analogues was able to displace wild-type hEGF fully in receptor competition binding assays, affinity of the receptor for the mutants was substantially reduced, varying from 0.4 to less than 0.01% of that observed for wild-type growth factor. At sufficiently high concentrations these mutants were able to stimulate DNA synthesis in mouse keratinocytes. Substitution of lysine for arginine 41 reduced the receptor affinity 250-fold from that observed for wild type, despite retention of the positive electrostatic charge. The lysine substitution leaves a reactive amine at position 41 and made it possible, using amine-specific chemical modification reagents, to produce selected arginine homologues that were tested for their effects on receptor binding, receptor tyrosine kinase activation, and stimulation of DNA synthesis in mouse keratinocytes. The reaction of lysine 41 with methyl acetimidate resulted in a lysineacetamidine product which only partially restored activity of the lysine hEGF mutant. However, reaction with O-methylisourea resulted in generation of an arginine 41 homologue (homoarginine) which restored full activity. The results indicate that the chemical properties inherent in the guanidinium group of the arginine 41 side chain of hEGF are responsible for optimal receptor-ligand association.     

Charge-charge interactions influence the denatured state ensemble and contribute to protein stability

Several recent studies have shown that it is possible to increase protein stability by improving electrostatic interactions among charged groups on the surface of the folded protein. However, the stability increases are considerably smaller than predicted by a simple Coulomb's law calculation, and in some cases, a charge reversal on the surface leads to a decrease in stability when an increase was predicted. These results suggest that favorable charge-charge interactions are important in determining the denatured state ensemble, and that the free energy of the denatured state may be decreased more than that of the native state by reversing the charge of a side chain. We suggest that when the hydrophobic and hydrogen bonding interactions that stabilize the folded state are disrupted, the unfolded polypeptide chain rearranges to compact conformations with favorable long-range electrostatic interactions. These charge-charge interactions in the denatured state will reduce the net contribution of electrostatic interactions to protein stability and will help determine the denatured state ensemble. To support this idea, we show that the denatured state ensemble of ribonuclease Sa is considerably more compact at pH 7 where favorable charge-charge interactions are possible than at pH 3, where unfavorable electrostatic repulsion among the positive charges causes an expansion of the denatured state ensemble. Further support is provided by studies of the ionic strength dependence of the stability of charge-reversal mutants of ribonuclease Sa. These results may have important implications for the mechanism of protein folding.     

In-vitro selection of highly stabilized protein variants with optimized surface

Thermostable proteins are of prime importance in protein science, but it has remained difficult to develop general strategies for stabilizing a protein. Site-directed mutagenesis based on comparisons with thermophilic homologs is rarely successful because the sequence differences are too numerous and dominated by neutral mutations. Here we used a method of directed evolution to increase the stability of a mesophilic protein, the cold shock protein Bs-CspB from Bacillus subtilis. It differs from its thermophilic counterpart Bc-Csp from Bacillus caldolyticus at 12 surface-exposed positions. To elucidate the stabilizing potential of exposed amino acid residues, six of these variant positions were randomized by saturation mutagenesis, the corresponding library of sequences was inserted into the gene-3-protein of the filamentous phage fd, and stabilized variants were selected by the Proside technique. Proside links the increased protease resistance of stabilized protein variants with the infectivity of the phage. Many strongly stabilized variants of Bs-CspB were identified in two selections, one in the presence of a denaturant and the other at elevated temperature. Several of them are significantly more stable than the naturally thermostable homolog Bc-Csp, and the best variant reaches Tm-Csp (the homolog from the hyperthermophile Thermotoga maritima) in stability. Remarkably, this variant differs from Tm-Csp at five and from Bc-Csp at all six randomized positions. This indicates that proteins can be strongly stabilized by many different sets of surface mutations, and Proside selects them efficiently from large libraries. The course of the selection could be directed by the conditions. In an ionic denaturant non-polar surface interactions were optimized, whereas at elevated temperature variants with improved electrostatics were selected, pointing to two different strategies for stabilization at protein surfaces.     

Citrate and phosphate influence on green fluorescent protein thermal stability

Protein structure and function can be regulated by no specific interactions, such as ionic interactions in the presence of salts. Green fluorescent protein (GFP) shows remarkable structural stability and high fluorescence; its stability can be directly related to its fluorescence output, among other characteristics. GFP is stable under increasing temperatures, and its thermal denaturation is highly reproducible. The aim of this study was to evaluate the thermal stability of GFP in the presence of different salts at several concentrations and exposed to constant temperatures, in a range of 70–95°C. Thermal stability was expressed in decimal reduction time. It was observed that the D‐values obtained were higher in the presence of citrate and phosphate, when compared with that obtained in their absence, indicating that these salts stabilized the protein against thermal denaturation. © 2010 American Institute of Chemical Engineers Biotechnol. Prog., 2011     

The stability of cell surface protein to surfactants and denaturants

The effects of several denaturants and detergents on the structure and stability of cell surface protein have been evaluated by circular dichroism and fluorescence measurements. Cell surface protein undergoes a single broad transition in both urea and guanidinium chloride. Although guanidinium chloride is twice as effective as urea on a molar basis, both appear to eliminate all of the organized structure present in the native molecule. Nonionic surfactants and lysolecithin have little effect on cell surface protein. However, sodium dodecyl sulfate increases the alpha helical content and cetyltrimethylammonium bromide increases the beta structure of cell surface protein. The reorganization of the polypeptide backbone requires the loss of certain restraints imposed by tertiary interactions as evidenced by a decrease in ellipticity in the far ultraviolet and in the polarization of tryptophanyl fluorescence. These results along with the data of a previous paper (Alexander, S. S., Jr., Colonna, G., Yamada, K. M., Pastan, I., and Edelhoch, H. (1978) J. Biol. Chem. 253, 5820--5824) suggest the presence of structural domains distributed along the flexible polypeptide chain of cell surface protein.     

Peptide tags for enhanced cellular and protein adhesion to single-crystalline sapphire

In order to facilitate a novel means for coupling proteins to metal oxides, peptides were identified from a dodecamer peptide yeast surface display library that bound a model metal oxide material, the C, A, and R crystalline faces of synthetic sapphire (alpha-Al(2)O(3)). Seven rounds of screening yielded peptides enriched in basic amino acids compared to the naive library. While the C-face had a high background of endogenous yeast cell binding, the A- and R faces displayed clear peptide-mediated cell adhesion. Cell detachment assays showed that cell adhesion strength correlated positively with increasing basicity of expressed peptides. Cell adhesion was also shown to be sensitive to buffer ionic strength as well as incubation with soluble peptide (with half maximal inhibition of cell binding at approximately 5 microM peptide). Next, dodecamer peptides cloned into yeast showed that lysine led to stronger interactions than arginine, and that charge distribution affected adhesion strength. We postulate binding to arise from peptide geometries that permit conformation alignment of the basic amino acids towards the surface so that the charged groups can undergo local electrostatic interactions with the surface oxide. Lastly, peptide K1 (-(GK)(6)) was cloned onto the c-terminus of maltose binding protein (MBP) and the resultant mutant protein showed a half-maximal binding at approximately 10(-7)-10(-6) M, which marked a approximately 500- to 1,000-fold binding improvement to sapphire's A-face as compared with wild-type MBP. Targeting proteins to metal oxide surfaces with peptide tags may provide a facile one-step alternative coupling chemistry for the formation of protein bioassays and biosensors.     

Effect of protein denaturants on the structure of water: Electron spin resonance spin probe study

The effect of the urea-class of protein denaturants on the structure of liquid water was studied. The method chosen was to monitor the microviscosity of the medium by estimating the reorientational correlation time, τ₈ and proton hyperfine linewidth, WH of an inert stable organic free radical (2, 2, 6, 4-tetramethyl-piperid 4-one-l-oxyl) by measuring its electron spin resonance spectra in aqueous denaturant solutions. Urea, thiourea and guanidinium chloride were found to disrupt water structure efficiently and continuously, with the effect significant at low molarities. Dimethyl urea was found to be somewhat less efficient. The relation between the structure-breaking tendency of the denaturants and their ability to weaken hydrophobic interactions is rationalised.     

Probing electrostatic channeling in protozoal bifunctional thymidylate synthase-dihydrofolate reductase using site-directed mutagenesis

In this study we used site-directed mutagenesis to test the hypothesis that substrate channeling in the bifunctional thymidylate synthase-dihydrofolate reductase enzyme from Leishmania major occurs via electrostatic interactions between the negatively charged dihydrofolate produced at thymidylate synthase and a series of lysine and arginine residues on the surface of the protein. Accordingly, 12 charge reversal or charge neutralization mutants were made, with up to 6 putative channel residues changed at once. The mutants were assessed for impaired channeling using two criteria: a lag in product formation at dihydrofolate reductase and an increase in dihydrofolate accumulation. Surprisingly, none of the mutations produced changes consistent with impaired channeling, so our findings do not support the electrostatic channeling hypothesis. Burst experiments confirmed that the mutants also did not interfere with intermediate formation at thymidylate synthase. One mutant, K282E/R283E, was found to be thymidylate synthase-dead because of an impaired ability to form the covalent enzyme-methylene tetrahydrofolate-deoxyuridate complex prerequisite for chemical catalysis.     

Analysis of a nuclear localization signal of simian virus 40 major capsid protein Vp1.

The nuclear localization signal of the major structural protein, Vp1, of simian virus 40 was further defined by mutagenesis. The targeting activity was examined in cells microinjected with SV-Vp1 variant viral DNAs bearing either an initiation codon mutation of the agnoprotein or mutations in the Vp1 coding sequence or microinjected with pSG5-Vp1 and pSG5-Vp1 mutant DNAs in which Vp1 or mutant Vp1 is expressed from simian virus 40 early promoter. The Vp1 nuclear localization signal functioned autonomously without agno-protein once the Vp1 protein was synthesized in the cytoplasm. The targeting activity was localized to the amino-terminal 19 residues. While replacement of cysteine 10 with glycine, alanine, or serine did not affect the activity, replacement of arginine 6 with glycine caused the cytoplasmic phenotype. When multiple mutations were introduced among residue 5, 6, 7, 16, 17, or 19, the targeting activity was found to reside in two clusters of basic residues, a cluster of lysine 5, arginine 6, and lysine 7 and a cluster of lysine 16, lysine 17, and lysine 19. The clusters are independently important for nuclear localization activity.     

Effects of arginine on heat-induced aggregation of concentrated protein solutions

Arginine is one of the commonly used additives to enhance refolding yield of proteins, to suppress aggregation of proteins, and to increase solubility of proteins, and yet the molecular interactions that contribute to the role of arginine are unclear. Here, we present experiments, using bovine serum albumin (BSA), lysozyme (LYZ), and β-lactoglobulin (BLG) as model proteins, to show that arginine can enhance heat-induced aggregation of concentrated protein solutions, contrary to the conventional belief that arginine is a universal suppressor of aggregation. Results show that the enhancement in aggregation is caused only for BSA and BLG, but not for LYZ, indicating that arginine's preferential interactions with certain residues over others could determine the effect of the additive on aggregation. We use this previously unrecognized behavior of arginine, in combination with density functional theory calculations, to identify the molecular-level interactions of arginine with various residues that determine arginine's role as an enhancer or suppressor of aggregation of proteins. The experimental and computational results suggest that the guanidinium group of arginine promotes aggregation through the hydrogen-bond-based bridging interactions with the acidic residues of a protein, whereas the binding of the guanidinium group to aromatic residues (aggregation-prone) contributes to the stability and solubilization of the proteins. The approach, we describe here, can be used to select suitable additives to stabilize a protein solution at high concentrations based on an analysis of the amino acid content of the protein.     

Structural characterization of protein-denaturant interactions: crystal structures of hen egg-white lysozyme in complex with DMSO and guanidinium chloride

A variety of physico-chemical methods employ chemical denaturants to unfold proteins, and study different biophysical processes involved therein. Chemical denaturants are believed to induce unfolding by stabilizing the unfolded state of proteins over the folded state, either macroscopically or through specific interactions. In order to characterize the nature of specific interactions between proteins and denaturants, we have solved crystal structures of hen egg-white lysozyme complexed with denaturants, and report here dimethyl sulfoxide and guanidinium chloride complexes. The dimethyl sulfoxide molecules and guanidinium ions were seen to bind the protein at specific sites and were involved in characteristic interactions. They share a major binding site between them, the C site in the sugar binding cleft of the enzyme. Although the overall conformations of the complexes were very similar to the native structure, spectacular conformational changes were seen to occur locally. Temperature factors were also seen to drop dramatically in the local regions close to the denaturant binding sites. An interesting observation of the present study was the generation of a sodium ion binding site in hen egg-white lysozyme in the presence of denaturants, which was hitherto unknown in any of the other lysozyme structures solved so far. Loss of some of the crucial side chain-main chain interactions may form the initial events in lysozyme unfolding.     

Folding Study of Venus Reveals a Strong Ion Dependence of Its Yellow Fluorescence under Mildly Acidic Conditions

Venus is a yellow fluorescent protein that has been developed for its fast chromophore maturation rate and bright yellow fluorescence that is relatively insensitive to changes in pH and ion concentrations. Here, we present a detailed study of the stability and folding of Venus in the pH range from 6.0 to 8.0 using chemical denaturants and a variety of spectroscopic probes. By following hydrogen-deuterium exchange of ¹⁵N-labeled Venus using NMR spectroscopy over 13 months, residue-specific free energies of unfolding of some highly protected amide groups have been determined. Exchange rates of less than one per year are observed for some amide groups. A super-stable core is identified for Venus and compared with that previously reported for green fluorescent protein. These results are discussed in terms of the stability and folding of fluorescent proteins. Under mildly acidic conditions, we show that Venus undergoes a drastic decrease in yellow fluorescence at relatively low concentrations of guanidinium chloride. A detailed study of this effect establishes that it is due to pH-dependent, nonspecific interactions of ions with the protein. In contrast to previous studies on enhanced green fluorescence protein variant S65T/T203Y, which showed a specific halide ion-binding site, NMR chemical shift mapping shows no evidence for specific ion binding. Instead, chemical shift perturbations are observed for many residues primarily located in both lids of the β-barrel structure, which suggests that small scale structural rearrangements occur on increasing ionic strength under mildly acidic conditions and that these are propagated to the chromophore resulting in fluorescence quenching.