Skip to main content

MINI REVIEW article

Front. Mol. Biosci., 29 November 2021
Sec. Biophysics
This article is part of the Research Topic Intracellular Molecular Processes Affected by pH View all 6 articles

Protein pI and Intracellular Localization

  • 1Department of Genetic Engineering, Faculty of Biology-Oriented Science and Technology, Kindai University, Wakayama, Japan
  • 2Center for Sustainable Resource Science, RIKEN Yokohama Institute, Yokohama, Japan
  • 3Laboratory of Cell Signaling and Development, Faculty of Life Sciences, Kyoto Sangyo University, Kyoto, Japan

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 (Sillero and Ribeiro, 1989), and the focusing positions of proteins in immobilized pH gradients and two-dimentional gels can be reliably predicted from their amino acid composition (Bjellqvist et al., 1993; Bjellqvist et al., 1994; Link et al., 1997). Notably, three-dimensional structure and pH of surrounding environment can influence ionizable groups and affect the net charge on the molecule significantly (Russell and Warshel, 1985).

Various calculative algorithms have been developed for estimating protein pIs in agreement with experiments regardless of structural aspect (Gasteiger et al., 2003; Cargile et al., 2004; Gauci et al., 2008; Maldonado et al., 2010; Audain et al., 2016). Some methods take into account the effect of the amino acids residues adjacent to the charged residues, such as aspartate and glutamate (Cargile et al., 2008), effects of posttranslational modifications, such as phosphorylation and N-terminal acetylation (Gauci et al., 2008), or effects of the presence of polyelectrolyte chains around proteins (Srivastava et al., 2017). In addition, the experimentally observed protein pI values were summarized in the experimental databases (Hoogland et al., 2004; Bunkute et al., 2015). Also, a database of protein pIs that were predicted using multiple available methods has been presented (Kozlowski, 2017).

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 (Blattner et al., 1997; 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 (Schwartz et al., 2001; 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 (Arakawa and Timasheff, 1985). 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 (Kurotani et al., 2010; 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; Garcia-Moreno, 2009). 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 (Schwartz et al., 2001). 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 (Knight et al., 2004). 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 (Schwartz et al., 2001). These data strongly suggested a link between the whole-proteome pI distributions and subcellular localization.

FIGURE 1
www.frontiersin.org

FIGURE 1. (A) Smoothened distributions of protein pI in the proteomes of Escherichia coli (solid line), Drosophila melanogaster (dotted line) and Caenorhabditis elegans (dashed line), as based on the histogram data presented by Schwartz et al. (Schwartz et al., 2001) (B) Distributions of protein pI and solubility in the human proteome. The thick blue line and thin black line show smoothened and un-smoothened whole-proteome distributions of protein pI values, respectively. The red line presents solubility of cell free-expressed human proteins. Arrowheads indicate the positions of shoulders and sub-peaks observed in the whole-proteome pI distribution.

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 (Schwartz et al., 2001). Several additional modes, such as a minor peak at the pI above 11.0, were distinguished in the whole-proteome pI distributions of eukaryotic proteins (Wu et al., 2006; Carugo, 2007), further suggesting the presence of divergent subcellular protein pI profiles. Markedly, the trimodality of proteome-wide pI distributions is not conserved across eukaryotic species. Although trimodal distributions of protein pI have been observed in some proteomes of eukaryotic species, such as Saccharomyces cerevisiae, Caenorhabditis elegans, and Drosophila melanogaster (Schwartz et al., 2001; Figure 1A), bimodal distributions of pI were witnessed in the proteomes of human, mouse, and malaria plasmodium (Medjahed et al., 2003). In addition, contrary to the earlier study, it was reported that the global pI distribution of C. elegans and S. cerevisiae proteins are bimodal (Medjahed et al., 2003; Ho et al., 2006), as explained by difference in the algorithms employed for calculation of protein pI. Also, our recent study demonstrated that the profile of protein pI values determined for the complete human proteome is essentially bimodal with the major acidic and alkaline peaks at around pI 6.0 and pI 8.25 (Kurotani et al., 2019, Figure 1B). Notably, the two major peaks of the pI distribution are not Gaussian and not well-resolved, leaving open the possibility that the broad modality corresponding to nuclear proteins may be obscured by the two major overlapping peaks. Moreover, the distribution of human proteins displayed some additional statistical features, such as minor peaks and peak shoulders (Kurotani et al., 2019, Figure 1B). Protein localization patterns were further analyzed throughout the whole-proteome pI distribution, and it was found that the observed major and minor peaks of the distribution were associated with specific subcellular localizations (Kurotani et al., 2019).

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 (Knight et al., 2004). It was noted that smaller proteomes of intracellular parasites are more alkaline because of their adaptation to elevated host pH (Knight et al., 2004). It was also reported that, proteome pI adjusts to high-temperature environmental conditions of Thermoplasma volcanium growth (Kawashima et al., 2000). A later bioinformatics study confirmed significant relationships between pI and habitat, such as salinity and host environments, in prokaryotic proteomes, but it could not reveal significant correlations with oxygen and temperature requirements (Kiraga et al., 2007).

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 (Knight et al., 2004). The most closely related organisms displayed very different proteome distributions as those typically observed between the organisms from different domains of life. Other study reported, based on analysis of pI distribution of 115 fully sequenced genomes, that the modal distributions do not reflect phylogeny or sequence evolution, but rather the chemical properties of amino acids (Weiller et al., 2004). Similarly, more recent investigation could not reveal any relation between pI bias and taxonomy both in prokaryotic and eukaryotic proteomes, however a phylogenetic signal was observed in mitochondrial proteomes (Kiraga et al., 2007). These findings are consistent with other observations that the pI values of protein orthologs are poorly conserved from species to species (Wilkins and Williams 1997; Nandi et al., 2005), further challenging the possibility of phylogenic pI adaptation to evolutionary constraints.

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 (Schwartz et al., 2001). Furthermore, it was demonstrated, using the experimental data of protein localization based on GFP tagging and microscopic detection of about 4,000 yeast proteins in 22 subcellular compartments, that the distributions of protein pI differ significantly in subcellular compartments (Huh et al., 2003; Ho et al., 2006). Although both the global and local intracellular pI values showed a bimodal distribution, the ratio between proteins of acidic and basic pI varied significantly among individual compartments. It was found that the proteomes of the cytoplasm, Golgi apparatus and vacuole are highly biased towards acidic pI, whereas the mitochondrial sub-proteome has a bias towards proteins of basic pI (Ho et al., 2006). Similarly, it was reported that yeast proteins localized in the organelles with alkaline pH, such as peroxisomes, endoplasmic reticulum and mitochondria, had relatively high pI values, whereas the proteins contained in the acidic organelles, such as vacuoles, Golgi and endosomes, tended to have rather low pIs (Brett et al., 2006). A detailed study of multiple proteomes from different biological species also confirmed that the proteomes of the cytoplasm, lysosomes, vacuoles and cytoskeleton are acidic, whereas those of mitochondria and the plasma membrane tend to be basic (Kiraga et al., 2007). Our recent study using one of the latest updates of human genome data disclosed a plethora of strong statistically significant correlations between protein pI and subcellular localization. Protein pI was found to correlate positively with mitochondrial and nuclear locations and negatively with lysosomal, cytoskeletal, peroxisomal and cytoplasmic ones (Kurotani et al., 2019, Figure 2). The most recent analysis of protein pI distributions in the interactomes across life domains has largely confirmed the above relationships between protein pI and subcellular localization (Chasapis and Konstantinoudis, 2020). The study also revealed that acidic proteins have the highest average number of interactions, whereas basic proteins have the lowest number of interactions in both prokaryotic and eukaryotic proteomes. A rationale behind these relationships remains unknown. Of note, the difference in the intracellular spatial distributions of proteins was proposed to be driven by a non-uniform distribution of intracellular pH (Baskin et al., 2006). This phenomenon based on the mechanism of pH-induced protein trapping was witnessed both in artificial systems and in living cells.

FIGURE 2
www.frontiersin.org

FIGURE 2. Relationships between protein pI and subcellular localization in the human proteome (see text for detailed explanations). Figure adapted from Kurotani et al., 2019 (CC BY 4 license, authors retain the copyright).

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 (Brett et al., 2006). However, this finding is difficult to reconcile with the notion that proteins are least soluble at the pH that corresponds to their pI. Indeed, a tendency has been observed for the averaged values of local pI distributions to differ from local pH (Chan et al., 2006; Chan and Warwicker, 2009). Furthermore, some analyses of multiple bacterial and eukaryotic proteomes failed to detect any statistically significant relationship between local pI distributions and subcellular intra-organelle pH (Wu et al., 2006; Kiraga et al., 2007).

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 (Alexov, 2004; Talley and Alexov, 2010; Loell and Nanda, 2018). The evidence has been presented for adaptation of the protein pH dependence, rather than protein pI, to local subcellular pH. The average pH of maximal stability, but not the average pI of proteins in a subcellular compartment, was demonstrated to correlate with subcellular pH (Chan et al., 2006; Chan and Warwicker, 2009; Garcia-Moreno, 2009). In this connection, it was shown that the pH optimum for protein stability and activity can differ significantly from the pI value (Alexov, 2004; Talley and Alexov, 2010). The recent bioinformatics analysis of the human proteome confirmed that the specific pI distributions at different subcellular locations are governed by local physicochemical environment and further suggested that the local pH and organelle membrane charge are the main factors responsible for variation of the intracellular localization-specific pI profiles (Kurotani et al., 2019; see next section for details). Notably, the study failed to detect a statistically significant correlation between the mean values of local pI distributions and intra-organelle pH alone, however, it was observed that the proteins in alkaline compartments tended to have higher mean pI values than those in acidic organelles.

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 (Nandi et al., 2005). In addition, a significant pI distribution bias, acidic or alkaline, was reported for certain protein functional classes localized in specific subcellular compartments (Wang and Tang, 2017).

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 (Kurotani et al., 2019). Particularly, a robust positive correlation was witnessed between protein pI and propensity for mitochondrial and nuclear localization, and a negative correlation was observed for cytoskeletal, cytoplasmic, peroxisomal, lysosomal and endoplasmic reticulum proteins. These findings are broadly consistent with the data obtained by previous analyses of multiple prokaryotic and eukaryotic proteomes (Schwartz et al., 2001; Brett et al., 2006; Ho et al., 2006; Kiraga et al., 2007). The proteome-wide relationships between protein pI and subcellular localization are summarized in Figure 2.

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; Kurotani et al., 2019), the membrane charge could be regarded as a likely factor related to the variation of intracellular localization-specific patterns. Although the correlation between organelle membrane charge and mean local pI was not statistically significant, a composite function of the two variables, compartment pH and membrane charge, could approximate localization-specific mean pI with a statistically significant coefficient of determination (Kurotani et al., 2019). The result indicates that local pH and membrane charge jointly define intracellular localization-specific pI patterns. In a practical sense, the finding that membrane charge affects organelle-specific protein pI patterns can be useful when considering intracellular targeting of both endogenous and ectopically expressed exogenous proteins.

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 (Loell and Nanda, 2018). Further proteomics and evolutionary studies are necessary to elucidate the factors that define subcellular localization of proteins with different physicochemical and functional traits.

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

Alexov, E. (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. doi:10.1046/j.1432-1033.2003.03917.x

CrossRef Full Text | Google Scholar

Arakawa, T., and Timasheff, S. N. (1985). Theory of Protein Solubility. Methods Enzymol. 114, 49–77. doi:10.1016/0076-6879(85)14005-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Audain, E., Ramos, Y., Hermjakob, H., Flower, D. R., and Perez-Riverol, Y. (2016). Accurate Estimation of Isoelectric point of Protein and Peptide Based on Amino Acid Sequences. Bioinformatics (Oxford, England) 32 (6), 821–827. doi:10.1093/bioinformatics/btv674

PubMed Abstract | CrossRef Full Text | Google Scholar

Baskin, E. M., Bukshpan, S., and Zilberstein, G. V. (2006). pH-Induced Intracellular Protein Transport. Phys. Biol. 3 (2), 101–106. doi:10.1088/1478-3975/3/2/002

PubMed Abstract | CrossRef Full Text | Google Scholar

Bjellqvist, B., Basse, B., Olsen, E., and Celis, J. 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. Electrophoresis 15 (3-4), 529–539. doi:10.1002/elps.1150150171

PubMed Abstract | CrossRef Full Text | Google Scholar

Bjellqvist, B., Hughes, G. J., Pasquali, C., Paquet, N., Ravier, F., Sanchez, J.-C., et al. (1993). The Focusing Positions of Polypeptides in Immobilized pH Gradients Can Be Predicted from Their Amino Acid Sequences. Electrophoresis 14 (10), 1023–1031. doi:10.1002/elps.11501401163

PubMed Abstract | CrossRef Full Text | Google Scholar

Blattner, F. R., Plunkett, G., Bloch, C. A., Perna, N. T., Burland, V., Riley, M., et al. (1997). The Complete Genome Sequence of Escherichia C K-12. Science 277 (5331), 1453–1462. doi:10.1126/science.277.5331.1453

PubMed Abstract | CrossRef Full Text | Google Scholar

Brett, C. L., Donowitz, M., and Rao, R. (2006). Does the Proteome Encode Organellar pH? FEBS Lett. 580 (3), 717–719. doi:10.1016/j.febslet.2005.12.103

PubMed Abstract | CrossRef Full Text | Google Scholar

Bunkute, E., Cummins, C., Crofts, F. J., Bunce, G., Nabney, I. T., and Flower, D. R. (2015). PIP-DB: the Protein Isoelectric Point Database. Bioinformatics 31 (2), 295–296. doi:10.1093/bioinformatics/btu637

PubMed Abstract | CrossRef Full Text | Google Scholar

Cargile, B. J., Bundy, J. L., Freeman, T. W., and Stephenson, J. L. (2004). Gel Based Isoelectric Focusing of Peptides and the Utility of Isoelectric point in Protein Identification. J. Proteome Res. 3 (1), 112–119. doi:10.1021/pr0340431

CrossRef Full Text | Google Scholar

Cargile, B. J., Sevinsky, J. R., Essader, A. S., Eu, J. P., and Stephenson, J. L. (2008). Calculation of the Isoelectric point of Tryptic Peptides in the pH 3.5-4.5 Range Based on Adjacent Amino Acid Effects. Electrophoresis 29 (13), 2768–2778. doi:10.1002/elps.200700701

PubMed Abstract | CrossRef Full Text | Google Scholar

Carugo, O. (2007). Isoelectric Points of Multi-Domain Proteins. Bioinformation 2 (3), 101–104. doi:10.6026/97320630002101

PubMed Abstract | CrossRef Full Text | Google Scholar

Chan, P., Lovrić, J., and Warwicker, J. (2006). Subcellular pH and Predicted pH-Dependent Features of Proteins. Proteomics 6 (12), 3494–3501. doi:10.1002/pmic.200500534

PubMed Abstract | CrossRef Full Text | Google Scholar

Chan, P., and Warwicker, J. (2009). Evidence for the Adaptation of Protein pH-Dependence to Subcellular pH. BMC Biol. 7, 69. doi:10.1186/1741-7007-7-69

PubMed Abstract | CrossRef Full Text | Google Scholar

Chasapis, C. T., and Konstantinoudis, G. (2020). Protein Isoelectric point Distribution in the Interactomes Across the Domains of Life. Biophysical Chem. 256, 106269. doi:10.1016/j.bpc.2019.106269

PubMed Abstract | CrossRef Full Text | Google Scholar

Garcia-Moreno, B. (2009). Adaptations of Proteins to Cellular and Subcellular pH. J. Biol. 8 (11), 98. doi:10.1186/jbiol199

CrossRef Full Text | Google Scholar

Gasteiger, E., Gattiker, A., Hoogland, C., Ivanyi, I., Appel, R. D., and Bairoch, A. (2003). ExPASy: The Proteomics Server for In-Depth Protein Knowledge and Analysis. Nucleic Acids Res. 31 (13), 3784–3788. doi:10.1093/nar/gkg563

PubMed Abstract | CrossRef Full Text | Google Scholar

Gauci, S., van Breukelen, B., Lemeer, S. M., Krijgsveld, J., and Heck, A. J. R. (2008). A Versatile Peptide pI Calculator for Phosphorylated and N-Terminal Acetylated Peptides Experimentally Tested Using Peptide Isoelectric Focusing. Proteomics 8 (23-24), 4898–4906. doi:10.1002/pmic.200800295

PubMed Abstract | CrossRef Full Text | Google Scholar

Ho, E., Hayen, A., and Wilkins, M. R. (2006). Characterisation of Organellar Proteomes: A Guide to Subcellular Proteomic Fractionation and Analysis. Proteomics 6 (21), 5746–5757. doi:10.1002/pmic.200600241

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoogland, C., Mostaguir, K., Sanchez, J.-C., Hochstrasser, D. F., and Appel, R. D. (2004). SWISS-2DPAGE, Ten Years Later. Proteomics 4 (8), 2352–2356. doi:10.1002/pmic.200300830

PubMed Abstract | CrossRef Full Text | Google Scholar

Huh, W.-K., Falvo, J. V., Gerke, L. C., Carroll, A. S., Howson, R. W., Weissman, J. S., et al. (2003). Global Analysis of Protein Localization in Budding Yeast. Nature 425 (6959), 686–691. doi:10.1038/nature02026

PubMed Abstract | CrossRef Full Text | Google Scholar

Kawashima, T., Amano, N., Koike, H., Makino, S.-i., Higuchi, S., Kawashima-Ohya, Y., et al. (2000). Archaeal Adaptation to Higher Temperatures Revealed by Genomic Sequence of Thermoplasma Volcanium. Proc. Natl. Acad. Sci. 97 (26), 14257–14262. doi:10.1073/pnas.97.26.14257

PubMed Abstract | CrossRef Full Text | Google Scholar

Kiraga, J., Mackiewicz, P., Mackiewicz, D., Kowalczuk, M., Biecek, P., Polak, N., et al. (2007). The Relationships between the Isoelectric Point and: Length of Proteins, Taxonomy and Ecology of Organisms. BMC genomics 8, 163. doi:10.1186/1471-2164-8-163

PubMed Abstract | CrossRef Full Text | Google Scholar

Knight, C. G., Kassen, R., Hebestreit, H., and Rainey, P. B. (2004). From the Cover: Global Analysis of Predicted Proteomes: Functional Adaptation of Physical Properties. Proc. Natl. Acad. Sci. 101 (22), 8390–8395. doi:10.1073/pnas.0307270101

PubMed Abstract | CrossRef Full Text | Google Scholar

Kozlowski, L. P. (2017). Proteome-pI: Proteome Isoelectric Point Database. Nucleic Acids Res. 45 (D1), D1112–D1116. doi:10.1093/nar/gkw978

PubMed Abstract | CrossRef Full Text | Google Scholar

Kurotani, A., Takagi, T., Toyama, M., Shirouzu, M., Yokoyama, S., Fukami, Y., 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. doi:10.1096/fj.09-139527

PubMed Abstract | CrossRef Full Text | Google Scholar

Kurotani, A., Tokmakov, A. A., Sato, K.-I., Stefanov, V. E., Yamada, Y., and Sakurai, T. (2019). Localization-Specific Distributions of Protein pI in Human Proteome Are Governed by Local pH and Membrane Charge. BMC Mol. Cel Biol 20 (1), 36. doi:10.1186/s12860-019-0221-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Link, A. J., Robison, K., and Church, G. M. (1997). Comparing the Predicted and Observed Properties of Proteins Encoded in the Genome ofEscherichia Coli K-12. Electrophoresis 18 (8), 1259–1313. doi:10.1002/elps.1150180807

PubMed Abstract | CrossRef Full Text | Google Scholar

Loell, K., and Nanda, V. (2018). Marginal Protein Stability Drives Subcellular Proteome Isoelectric point. Proc. Natl. Acad. Sci. USA 115 (46), 11778–11783. doi:10.1073/pnas.1809098115

PubMed Abstract | CrossRef Full Text | Google Scholar

Maldonado, A. A., Ribeiro, J. M., and Sillero, A. (2010). Isoelectric point, Electric Charge, and Nomenclature of the Acid-Base Residues of Proteins. Biochem. Mol. Biol. Educ. 38 (4), 230–237. doi:10.1002/bmb.20405

PubMed Abstract | CrossRef Full Text | Google Scholar

Medjahed, D., Smythers, G. W., Powell, D. A., Stephens, R. M., Lemkin, P. F., and Munroe, D. J. (2003). VIRTUAL2D: A Web-Accessible Predictive Database for Proteomics Analysis. Proteomics 3 (2), 129–138. doi:10.1002/pmic.200390021

PubMed Abstract | CrossRef Full Text | Google Scholar

Nandi, S., Mehra, N., Lynn, A. M., and Bhattacharya, A. (2005). Comparison of Theoretical Proteomes: Identification of COGs with Conserved and Variable pI within the Multimodal pI Distribution. BMC Genomics 6, 116. doi:10.1186/1471-2164-6-116

PubMed Abstract | CrossRef Full Text | Google Scholar

Russell, S. T., and Warshel, A. (1985). Calculations of Electrostatic Energies in proteinsThe Energetics of Ionized Groups in Bovine Pancreatic Trypsin Inhibitor. J. Mol. Biol. 185 (2), 389–404. doi:10.1016/0022-2836(85)90411-5

CrossRef Full Text | Google Scholar

Schwartz, R., Ting, C. S., and King, J. (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. doi:10.1101/gr.gr-1587r

PubMed Abstract | CrossRef Full Text | Google Scholar

Sillero, A., and Ribeiro, J. M. (1989). Isoelectric Points of Proteins: Theoretical Determination. Anal. Biochem. 179 (2), 319–325. doi:10.1016/0003-2697(89)90136-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Srivastava, D., Santiso, E., Gubbins, K., and Barroso da Silva, F. L. (2017). Computationally Mapping pKa Shifts Due to the Presence of a Polyelectrolyte Chain Around Whey Proteins. Langmuir 33 (42), 11417–11428. doi:10.1021/acs.langmuir.7b02271

PubMed Abstract | CrossRef Full Text | Google Scholar

Talley, K., and Alexov, E. (2010). On the pH-Optimum of Activity and Stability of Proteins. Proteins 78 (12), 2699–2706. doi:10.1002/prot.22786

PubMed Abstract | CrossRef Full Text | Google Scholar

Tokmakov, A. A., Kurotani, A., Shirouzu, M., Fukami, Y., and Yokoyama, S. (2014). Bioinformatics Analysis and Optimization of Cell-Free Protein Synthesis. Methods Mol. Biol. (Clifton, N.J.) 1118, 17–33. doi:10.1007/978-1-62703-782-2_2

PubMed Abstract | CrossRef Full Text | Google Scholar

Urquhart, B. L., Atsalos, T. E., Roach, D., Basseal, D. J., Bjellqvist, B., Britton, W. L., et al. (1997). 'Proteomic Contigs' ofMycobacterium Tuberculosis andMycobacterium Bovis (BCG) Using Novel Immobilised pH Gradients. Electrophoresis 18 (8), 1384–1392. doi:10.1002/elps.1150180813

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Meer, G., Voelker, D. R., and Feigenson, G. W. (2008). Membrane Lipids: Where They Are and How They Behave. Nat. Rev. Mol. Cel Biol 9 (2), 112–124. doi:10.1038/nrm2330

CrossRef Full Text | Google Scholar

VanBogelen, R. A., Schiller, E. E., Thomas, J. D., and Neidhardt, F. C. (1999). Diagnosis of Cellular States of Microbial Organisms Using Proteomics. Electrophoresis 20 (11), 2149–2159. doi:10.1002/(sici)1522-2683(19990801)20:11<2149:aid-elps2149>3.0.co;2-n

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, T., and Tang, H. (2017). The Physical Characteristics of Human Proteins in Different Biological Functions. PloS one 12 (5), e0176234. doi:10.1371/journal.pone.0176234

PubMed Abstract | CrossRef Full Text | Google Scholar

Weiller, G. F., Caraux, G., and Sylvester, N. (2004). The Modal Distribution of Protein Isoelectric Points Reflects Amino Acid Properties Rather Than Sequence Evolution. Proteomics 4 (4), 943–949. doi:10.1002/pmic.200200648

PubMed Abstract | CrossRef Full Text | Google Scholar

Wilkins, M. R., and Williams, K. 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. doi:10.1006/jtbi.1996.0346

CrossRef Full Text | Google Scholar

Wu, S., Wan, P., Li, J., Li, D., Zhu, Y., and He, F. (2006). Multi-Modality of pI Distribution in Whole Proteome. Proteomics 6 (2), 449–455. doi:10.1002/pmic.200500221

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, J., Han, X., and Gross, R. 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. doi:10.1016/s0014-5793(03)00581-7

PubMed Abstract | CrossRef Full Text | Google Scholar

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.

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

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, dG9rbWFrQHdha2Eua2luZGFpLmFjLmpw

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.