Abstract
The protein isoelectric point (pI) can be calculated from an amino acid sequence using computational analysis in a good agreement with experimental data. Availability of whole-genome sequences empowers comparative studies of proteome-wide pI distributions. It was found that the whole-proteome distributions of protein pI values are multimodal in different species. It was further hypothesized that the observed multimodality is associated with subcellular localization-specific differences in local pI distributions. Here, we overview the multimodality of proteome-wide pI distributions in different organisms focusing on the relationships between protein pI and subcellular localization. We also discuss the probable factors responsible for variation of the intracellular localization-specific pI profiles.
Introduction
The isoelectric point (pI) of a protein is defined as the pH at which the net charge of a protein molecule is zero. Accordingly, proteins are positively charged at a pH below their pI and negatively charged at a pH above their pI. The protein pI varies greatly from extremely acidic to highly alkaline values ranging from about 4.0 to 12.0. Hence, pI values have long been used to distinguish between proteins in methods for protein isolation, separation, purification, crystallization, etc. Amino acid composition of a protein sequence primarily defines its pI, based on the combination of dissociation constant (pKa) values of the constituent amino acids. Out of twenty common amino acids, two amino acids, aspartic acid, and glutamic acid, are negatively charged and three amino acids, lysine, arginine, and histidine, are positively charged at the neutral pH, as defined by their pKa values. Thus, the integral property of a protein, such as protein pI, was supposed to result from discrete local acidic and basic pKas of amino acid side chains. It was demonstrated that the protein pI can be estimated based on a polypeptide sequence in close agreement with experimentally determined pI values (), and the focusing positions of proteins in immobilized pH gradients and two-dimentional gels can be reliably predicted from their amino acid composition (; ; ). Notably, three-dimensional structure and pH of surrounding environment can influence ionizable groups and affect the net charge on the molecule significantly ().
Various calculative algorithms have been developed for estimating protein pIs in agreement with experiments regardless of structural aspect (; ; ; ; ). Some methods take into account the effect of the amino acids residues adjacent to the charged residues, such as aspartate and glutamate (), effects of posttranslational modifications, such as phosphorylation and N-terminal acetylation (), or effects of the presence of polyelectrolyte chains around proteins (). In addition, the experimentally observed protein pI values were summarized in the experimental databases (; ). Also, a database of protein pIs that were predicted using multiple available methods has been presented ().
Thus, protein pI is an integral property of a protein molecule fundamentally important for its characterization. The great variation of protein pI values brings about the question about the cause of this variation. Availability of whole-genome sequences allows comparative and evolutionary studies of proteome-wide pI distributions in different organisms. These studies have revealed important universal features of the whole proteome pI distributions providing insights into spatial organization of cellular proteomes. The localization- and function-specific differences in subcellular pI distributions have been disclosed. Our present paper overviews proteome-wide pI distributions focusing on the relationships between protein pI and subcellular localization.
Intrinsic Bimodality of Protein pI Distributions
The early studies of proteome-wide pI distributions demonstrated that they are bimodal, with distinct acidic and alkaline peaks, in several bacterial strains (; Urquhart et al., 1997; VanBogelen et al., 1999). The two major protein clusters, centering around pI 5.0 and pI 9.0, were observed in full proteomes of bacteria and archaea (; Figure 1A). It was suggested that the low abundance of sequences with unbiased pIs curtails protein precipitation at a near-neutral physiological pH. Indeed, the pI value affects solubility of a protein molecule at a given pH. Proteins display the least solubility in water-based solutions at the pH that corresponds to their pI, often resulting in protein aggregation (). It was demonstrated experimentally, using cell-free protein synthesis, that protein solubility positively correlates with the content of charged residues in the expressed proteins, and the proteins with pI 7.0–7.5 have the lowest rate of soluble expression (; Tokmakov et al., 2014; Figure 1B). On the other hand, the ratio of high to low cell-free expression levels was found to be stable in the wide range of pI values (Tokmakov et al., 2014), suggesting the absence of correlation between protein pI and expression level. Several studies proposed that the pI multimodality observed in different proteomes could be rooted in discrete pKa values for different amino acids (Weiller et al., 2004; Wu et al., 2006; ). Importantly, it was found that the pI distributions of cytosolic and integral membrane proteins corresponded to the two modes observed in the whole-proteome pI distributions. Cytoplasmic proteins clustered at pI 5.0 to 6.0, and integral membrane proteins exhibited a distinct clustering at pI 8.5 to 9.0 (). Also, investigation of complete predicted proteomes using theoretical 2D gels (MW vs pI) indicated that the membrane proteomes are generally more alkaline than the non-membrane ones (). The alkaline bias of the membrane proteins was attributed to the fact that biomembranes generally bear a negative charge due to the presence of negatively charged phospholipids, thus the positive charge of basic proteins at normal pH would promote favorable electrostatic interactions stabilizing the proteins in the membranes (). These data strongly suggested a link between the whole-proteome pI distributions and subcellular localization.
FIGURE 1
Common Multimodality of pI Distributions
Further investigations revealed that the protein pI profiles are trimodal in many eukaryotic proteomes (Figure 1A), and the presence of the third peak was linked to the appearance of the nuclear compartment in eukaryotes. Nuclear proteins were revealed to have a wide distribution varying from pI 4.5 to pI 10.0 (
Adaptation of pI Patterns to Environmental Constraints and Evolutionary Aspects
The average proteome pI and relative abundance of the acidic and alkaline peaks in bimodal pI distributions were analyzed in connection with organism taxonomy and environment. It was reported that proteome pI adapts to the conditions of bacterial growth; a significant positive correlation was observed between predicted proteome distributions on the theoretical 2D gels (MW vs pI) and the Biolog profile, a measure associated with ecological niche (
Notably, investigation of the relationship of genetic distance between bacterial strains and similarity of their theoretical 2D gels could not reveal a dependency on phylogeny (
Variation of Subcellular Localization-Specific pI Patterns
The proteome-wide relationships between protein pI and subcellular localization were analyzed in several bioinformatics studies of multiple proteomes. Initially, it was found that cytoplasmic proteins form the acidic modality and integral membrane proteins constitute the basic modality of the bimodal bacterial proteomes, whereas nuclear proteins may account for the third modality often observed in eukaryotes (
FIGURE 2

Relationships between protein pI and subcellular localization in the human proteome (see text for detailed explanations). Figure adapted from
Factors Behind the Variation of Subcellular pI Distributions
The variation of the localization-specific pI distributions was linked to the fact that local pH is different in subcellular compartments. It was reported that protein pIs averaged over a subcellular location correspond to experimentally measured intra-organellar pH in different compartments of the yeast cell and further speculated that subcellular protein pI and intra-organelle pH might have co-evolved to optimize protein function (
On the other hand, it was reported that the folded states of proteins are often most stable at pH values near their pI, and these values also correlate with their optimal pH for function (
Furthermore, some bioinformatics studies addressed proteome-wide relationships between protein pI, intracellular localization and functional classification. Using the COG database, which lists gene orthologs present across completed genomes and assigns their functional classification, both the invariant and highly changeable proteins, which occur with a high frequency, have been identified in different regions of proteome-wide pI distributions (
Generalized View of Localization-Specific pI Patterns (Importance of Local pH and Membrane Charge)
Thus, multiple bioinformatics studies converge on the assumption that the whole-proteome pI patterns adapt to environmental constraints and, in particular, the specific pI distribution at a certain subcellular location is defined by local environment. Our recent comprehensive analysis of 32,138 human proteins predicted to reside in 10 subcellular compartments, revealed the existence of strong relationships between protein pI and subcellular localization (
Another important result of the study is the finding that organelle-specific protein pI patterns are physically defined by local pH and membrane charge. Relationships between the local subcellular pH and pI distributions have been explicitly addressed in previous studies; they are discussed in section 6 of the present paper. However, the effect of membrane charge on the pI patterns of local sub-proteomes has not been thoroughly scrutinized. Considering that the membrane composition and content of the negatively charged membrane lipids, such as phosphatidylserine and phosphatidylinositol, vary greatly in intracellular organelles, ranging from 2% in peroxisomes to more than 17% in nuclei and ER (Yang et al., 2003; Van Meer et al., 2008;
Concluding Remarks
Genome sequencing has provided the information about all cellular and organismal proteins in many species. However, comprehension of life processes requires their further investigation at different levels. Uncovering subcellular localization of proteins with various physicochemical, structural and functional traits can reveal intracellular organization of proteomes and provide deeper understanding of their functioning. The recently disclosed relationships between protein pI and subcellular localization, as reviewed in this paper, contribute to spatial characterization of cellular processes. Still, the origin and mechanisms driving diversification of intracellular localization-specific pI patterns remain unknown. Although the possibility of positive evolutionary selection, which can promote beneficial protein pI patterns, seems unlikely (see section 4 for details), it was recently suggested that neutral evolution, i.e., accumulation of random mutations that have minimal impact on fitness and functional selection, might underline potential adjustment of protein pI to subcellular pH. It was revealed that the neutral evolutionary process leading to fixation of titratable residues in the protein core could likely be driven by marginal effects on protein stability (
Statements
Author contributions
AT, AK, and K-IS conceived and designed the article, AT wrote the manuscript, AK and K-IS reviewed and edited the manuscript. All authors have read and agreed to the published version of the manuscript.
Funding
This work was supported in part by the Grant-in-Aid for Scientific Research 15K07083 from the Ministry of Education, Culture, Sports, Science, and Technology of Japan and the Collaboration Research Grant 281027 from the Kobe University, Japan. Publication cost was covered by the institutional funds of the Kindai University and Kyoto Sangyo University.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
AlexovE. (2004). Numerical Calculations of the pH of Maximal Protein Stability. The Effect of the Sequence Composition and Three-Dimensional Structure. Eur. J. Biochem.271 (1), 173–185. 10.1046/j.1432-1033.2003.03917.x
2
ArakawaT.TimasheffS. N. (1985). Theory of Protein Solubility. Methods Enzymol.114, 49–77. 10.1016/0076-6879(85)14005-x
3
AudainE.RamosY.HermjakobH.FlowerD. R.Perez-RiverolY. (2016). Accurate Estimation of Isoelectric point of Protein and Peptide Based on Amino Acid Sequences. Bioinformatics (Oxford, England)32 (6), 821–827. 10.1093/bioinformatics/btv674
4
BaskinE. M.BukshpanS.ZilbersteinG. V. (2006). pH-Induced Intracellular Protein Transport. Phys. Biol.3 (2), 101–106. 10.1088/1478-3975/3/2/002
5
BjellqvistB.BasseB.OlsenE.CelisJ. E. (1994). Reference Points for Comparisons of Two-Dimensional Maps of Proteins from Different Human Cell Types Defined in a pH Scale where Isoelectric Points Correlate with Polypeptide Compositions. Electrophoresis15 (3-4), 529–539. 10.1002/elps.1150150171
6
BjellqvistB.HughesG. J.PasqualiC.PaquetN.RavierF.SanchezJ.-C.et al (1993). The Focusing Positions of Polypeptides in Immobilized pH Gradients Can Be Predicted from Their Amino Acid Sequences. Electrophoresis14 (10), 1023–1031. 10.1002/elps.11501401163
7
BlattnerF. R.PlunkettG.3rdBlochC. A.PernaN. T.BurlandV.RileyM.et al (1997). The Complete Genome Sequence of Escherichia C K-12. Science277 (5331), 1453–1462. 10.1126/science.277.5331.1453
8
BrettC. L.DonowitzM.RaoR. (2006). Does the Proteome Encode Organellar pH?FEBS Lett.580 (3), 717–719. 10.1016/j.febslet.2005.12.103
9
BunkuteE.CumminsC.CroftsF. J.BunceG.NabneyI. T.FlowerD. R. (2015). PIP-DB: the Protein Isoelectric Point Database. Bioinformatics31 (2), 295–296. 10.1093/bioinformatics/btu637
10
CargileB. J.BundyJ. L.FreemanT. W.StephensonJ. L.Jr (2004). Gel Based Isoelectric Focusing of Peptides and the Utility of Isoelectric point in Protein Identification. J. Proteome Res.3 (1), 112–119. 10.1021/pr0340431
11
CargileB. J.SevinskyJ. R.EssaderA. S.EuJ. P.StephensonJ. L.Jr (2008). Calculation of the Isoelectric point of Tryptic Peptides in the pH 3.5-4.5 Range Based on Adjacent Amino Acid Effects. Electrophoresis29 (13), 2768–2778. 10.1002/elps.200700701
12
CarugoO. (2007). Isoelectric Points of Multi-Domain Proteins. Bioinformation2 (3), 101–104. 10.6026/97320630002101
13
ChanP.LovrićJ.WarwickerJ. (2006). Subcellular pH and Predicted pH-Dependent Features of Proteins. Proteomics6 (12), 3494–3501. 10.1002/pmic.200500534
14
ChanP.WarwickerJ. (2009). Evidence for the Adaptation of Protein pH-Dependence to Subcellular pH. BMC Biol.7, 69. 10.1186/1741-7007-7-69
15
ChasapisC. T.KonstantinoudisG. (2020). Protein Isoelectric point Distribution in the Interactomes Across the Domains of Life. Biophysical Chem.256, 106269. 10.1016/j.bpc.2019.106269
16
Garcia-MorenoB. (2009). Adaptations of Proteins to Cellular and Subcellular pH. J. Biol.8 (11), 98. 10.1186/jbiol199
17
GasteigerE.GattikerA.HooglandC.IvanyiI.AppelR. D.BairochA. (2003). ExPASy: The Proteomics Server for In-Depth Protein Knowledge and Analysis. Nucleic Acids Res.31 (13), 3784–3788. 10.1093/nar/gkg563
18
GauciS.van BreukelenB.LemeerS. M.KrijgsveldJ.HeckA. J. R. (2008). A Versatile Peptide pI Calculator for Phosphorylated and N-Terminal Acetylated Peptides Experimentally Tested Using Peptide Isoelectric Focusing. Proteomics8 (23-24), 4898–4906. 10.1002/pmic.200800295
19
HoE.HayenA.WilkinsM. R. (2006). Characterisation of Organellar Proteomes: A Guide to Subcellular Proteomic Fractionation and Analysis. Proteomics6 (21), 5746–5757. 10.1002/pmic.200600241
20
HooglandC.MostaguirK.SanchezJ.-C.HochstrasserD. F.AppelR. D. (2004). SWISS-2DPAGE, Ten Years Later. Proteomics4 (8), 2352–2356. 10.1002/pmic.200300830
21
HuhW.-K.FalvoJ. V.GerkeL. C.CarrollA. S.HowsonR. W.WeissmanJ. S.et al (2003). Global Analysis of Protein Localization in Budding Yeast. Nature425 (6959), 686–691. 10.1038/nature02026
22
KawashimaT.AmanoN.KoikeH.MakinoS.-i.HiguchiS.Kawashima-OhyaY.et al (2000). Archaeal Adaptation to Higher Temperatures Revealed by Genomic Sequence of Thermoplasma Volcanium. Proc. Natl. Acad. Sci.97 (26), 14257–14262. 10.1073/pnas.97.26.14257
23
KiragaJ.MackiewiczP.MackiewiczD.KowalczukM.BiecekP.PolakN.et al (2007). The Relationships between the Isoelectric Point and: Length of Proteins, Taxonomy and Ecology of Organisms. BMC genomics8, 163. 10.1186/1471-2164-8-163
24
KnightC. G.KassenR.HebestreitH.RaineyP. B. (2004). From the Cover: Global Analysis of Predicted Proteomes: Functional Adaptation of Physical Properties. Proc. Natl. Acad. Sci.101 (22), 8390–8395. 10.1073/pnas.0307270101
25
KozlowskiL. P. (2017). Proteome-pI: Proteome Isoelectric Point Database. Nucleic Acids Res.45 (D1), D1112–D1116. 10.1093/nar/gkw978
26
KurotaniA.TakagiT.ToyamaM.ShirouzuM.YokoyamaS.FukamiY.et al (2010). Comprehensive Bioinformatics Analysis of Cell‐Free Protein Synthesis: Identification of Multiple Protein Properties that Correlate with Successful Expression. FASEB j.24 (4), 1095–1104. 10.1096/fj.09-139527
27
KurotaniA.TokmakovA. A.SatoK.-I.StefanovV. E.YamadaY.SakuraiT. (2019). Localization-Specific Distributions of Protein pI in Human Proteome Are Governed by Local pH and Membrane Charge. BMC Mol. Cel Biol20 (1), 36. 10.1186/s12860-019-0221-4
28
LinkA. J.RobisonK.ChurchG. M. (1997). Comparing the Predicted and Observed Properties of Proteins Encoded in the Genome ofEscherichia Coli K-12. Electrophoresis18 (8), 1259–1313. 10.1002/elps.1150180807
29
LoellK.NandaV. (2018). Marginal Protein Stability Drives Subcellular Proteome Isoelectric point. Proc. Natl. Acad. Sci. USA115 (46), 11778–11783. 10.1073/pnas.1809098115
30
MaldonadoA. A.RibeiroJ. M.SilleroA. (2010). Isoelectric point, Electric Charge, and Nomenclature of the Acid-Base Residues of Proteins. Biochem. Mol. Biol. Educ.38 (4), 230–237. 10.1002/bmb.20405
31
MedjahedD.SmythersG. W.PowellD. A.StephensR. M.LemkinP. F.MunroeD. J. (2003). VIRTUAL2D: A Web-Accessible Predictive Database for Proteomics Analysis. Proteomics3 (2), 129–138. 10.1002/pmic.200390021
32
NandiS.MehraN.LynnA. M.BhattacharyaA. (2005). Comparison of Theoretical Proteomes: Identification of COGs with Conserved and Variable pI within the Multimodal pI Distribution. BMC Genomics6, 116. 10.1186/1471-2164-6-116
33
RussellS. T.WarshelA. (1985). Calculations of Electrostatic Energies in proteinsThe Energetics of Ionized Groups in Bovine Pancreatic Trypsin Inhibitor. J. Mol. Biol.185 (2), 389–404. 10.1016/0022-2836(85)90411-5
34
SchwartzR.TingC. S.KingJ. (2001). Whole Proteome pI Values Correlate with Subcellular Localizations of Proteins for Organisms within the Three Domains of Life. Genome Res.11 (5), 703–709. 10.1101/gr.gr-1587r
35
SilleroA.RibeiroJ. M. (1989). Isoelectric Points of Proteins: Theoretical Determination. Anal. Biochem.179 (2), 319–325. 10.1016/0003-2697(89)90136-x
36
SrivastavaD.SantisoE.GubbinsK.Barroso da SilvaF. L. (2017). Computationally Mapping pKa Shifts Due to the Presence of a Polyelectrolyte Chain Around Whey Proteins. Langmuir33 (42), 11417–11428. 10.1021/acs.langmuir.7b02271
37
TalleyK.AlexovE. (2010). On the pH-Optimum of Activity and Stability of Proteins. Proteins78 (12), 2699–2706. 10.1002/prot.22786
38
TokmakovA. A.KurotaniA.ShirouzuM.FukamiY.YokoyamaS. (2014). Bioinformatics Analysis and Optimization of Cell-Free Protein Synthesis. Methods Mol. Biol. (Clifton, N.J.)1118, 17–33. 10.1007/978-1-62703-782-2_2
39
UrquhartB. L.AtsalosT. E.RoachD.BassealD. J.BjellqvistB.BrittonW. L.et al (1997). 'Proteomic Contigs' ofMycobacterium Tuberculosis andMycobacterium Bovis (BCG) Using Novel Immobilised pH Gradients. Electrophoresis18 (8), 1384–1392. 10.1002/elps.1150180813
40
Van MeerG.VoelkerD. R.FeigensonG. W. (2008). Membrane Lipids: Where They Are and How They Behave. Nat. Rev. Mol. Cel Biol9 (2), 112–124. 10.1038/nrm2330
41
VanBogelenR. A.SchillerE. E.ThomasJ. D.NeidhardtF. C. (1999). Diagnosis of Cellular States of Microbial Organisms Using Proteomics. Electrophoresis20 (11), 2149–2159. 10.1002/(sici)1522-2683(19990801)20:11<2149:aid-elps2149>3.0.co;2-n
42
WangT.TangH. (2017). The Physical Characteristics of Human Proteins in Different Biological Functions. PloS one12 (5), e0176234. 10.1371/journal.pone.0176234
43
WeillerG. F.CarauxG.SylvesterN. (2004). The Modal Distribution of Protein Isoelectric Points Reflects Amino Acid Properties Rather Than Sequence Evolution. Proteomics4 (4), 943–949. 10.1002/pmic.200200648
44
WilkinsM. R.WilliamsK. L. (1997). Cross-Species Protein Identification Using Amino Acid Composition, Peptide Mass Fingerprinting, Isoelectric point and Molecular Mass: A Theoretical Evaluation. J. Theor. Biol.186 (1), 7–15. 10.1006/jtbi.1996.0346
45
WuS.WanP.LiJ.LiD.ZhuY.HeF. (2006). Multi-Modality of pI Distribution in Whole Proteome. Proteomics6 (2), 449–455. 10.1002/pmic.200500221
46
YangJ.HanX.GrossR. W. (2003). Identification of Hepatic Peroxisomal Phospholipase A(2) and Characterization of Arachidonic Acid-Containing Choline Glycerophospholipids in Hepatic Peroxisomes. FEBS Lett.546 (2-3), 247–250. 10.1016/s0014-5793(03)00581-7
Summary
Keywords
protein pI, proteome-wide analysis, multimodality, subcellular localization, local environment
Citation
Tokmakov AA, Kurotani A and Sato K-I (2021) Protein pI and Intracellular Localization. Front. Mol. Biosci. 8:775736. doi: 10.3389/fmolb.2021.775736
Received
14 September 2021
Accepted
11 November 2021
Published
29 November 2021
Volume
8 - 2021
Edited by
Jan M. Antosiewicz, University of Warsaw, Poland
Reviewed by
Rajini Rao, Johns Hopkins University, United States
Gleb Zilberstein, Spectrophon Ltd., Rehovot, Israel
Michael Gilson, University of California, San Diego, United States
Updates

Check for updates
Copyright
© 2021 Tokmakov, Kurotani and Sato.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Alexander A. Tokmakov, tokmak@waka.kindai.ac.jp
This article was submitted to Biophysics, a section of the journal Frontiers in Molecular Biosciences
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.