Calculation of the relative metastabilities of proteins in subcellular compartments of Saccharomyces cerevisiae
© Dick; licensee BioMed Central Ltd. 2009
Received: 31 January 2009
Accepted: 18 July 2009
Published: 18 July 2009
Protein subcellular localization and differences in oxidation state between subcellular compartments are two well-studied features of the the cellular organization of S. cerevisiae (yeast). Theories about the origin of subcellular organization are assisted by computational models that can integrate data from observations of compositional and chemical properties of the system.
Presentation and implications of the hypothesis
I adopt the hypothesis that the state of yeast subcellular organization is in a local energy minimum. This hypothesis implies that equilibrium thermodynamic models can yield predictions about the interdependence between populations of proteins and their subcellular chemical environments.
Testing the hypothesis
Three types of tests are proposed. First, there should be correlations between modeled and observed oxidation states for different compartments. Second, there should be a correspondence between the energy requirements of protein formation and the order the appearance of organelles during cellular development. Third, there should be correlations between the predicted and observed relative abundances of interacting proteins within compartments.
The relative metastability fields of subcellular homologs of glutaredoxin and thioredoxin indicate a trend from less to more oxidizing as mitochondrion – cytoplasm – nucleus. Representing the overall amino acid compositions of proteins in 23 different compartments each with a single reference model protein suggests that the formation reactions for proteins in the vacuole (in relatively oxidizing conditions), ER and early Golgi (in relatively reducing conditions) are relatively highly favored, while that for the microtubule is the most costly. The relative abundances of model proteins for each compartment inferred from experimental data were found in some cases to correlate with the predicted abundances, and both positive and negative correlations were found for some assemblages of proteins in known complexes.
The results of these calculations and tests suggest that a tendency toward a metastable energy minimum could underlie some organizational links between the the chemical thermodynamic properties of proteins and subcellular chemical environments. Future models of this kind will benefit from consideration of additional thermodynamic variables together with more detailed subcellular observations.
A complex interplay of chemical and biological forces is responsible for subcellular structure. There exist in eukaryotic cells gradients between subcellular compartments of chemical properties such as pH, oxidation-reduction (or redox) state and chemical activity of water, among others [1–5]. Different population of proteins are localized within each subcellular compartment [6–8]. Within compartments, the relative abundances or levels of different proteins are not equal , and different proteins predominate in the various subcellular populations depending on growth state of the cell and exposure to environmental stress . Physical separation of key enzymes is thought to be essential in the cytoskeletal network and in regulation of metabolic pathways and other cellular functions [11, 12]. The patterns of subcellular structure persist even though populations of proteins turnover through continual degradation and synthesis in cells .
The biosynthesis and transport of proteins in an energy-demanding process . If cells have evolved to minimize their energy expenditure in the maintenance of biological function, it may be reasonable to expect to find signals of energy minimization in cellular organization. One such example is the finding that the relative abundances of amino acids in proteins correlate inversely with the metabolic cost of amino acid synthesis [15, 16], and that this is a temperature-dependent function . This observation is consistent with the notion that not all proteins are equal in energetic terms. In thermodynamic calculations of chemical affinity , the energy demands of protein formation (including synthesis and transport) are also a function of the local physical chemical environment, which includes variables such as oxidation-reduction potential . It follows that subcellular structures that are characterized by differences in the amino acid composition of proteins and in chemical potentials have distinct energetic consequences.
For the purposes of this study, the hypothesis is made that cellular organization is in a local energy minimum. Energy minimization in biological operations is not a new hypothesis, especially in the context of fitness and adaptation to the environment [20–22]. However, the implications of this hypothesis for subcellular organization have not been investigated from the standpoint of equilibrium chemical thermodynamics. Algorithms for computing the requisite standard molal Gibbs energies of proteins  and the relative abundances of proteins in metastable equilibrium  have recently been reported. The goal of this study is to perform these types of calculations for model systems representative of various levels of subcellular organization and to compare the results with observations and measurements reported in the literature. If successful, this exercise may lead to an enhanced awareness about the chemical forces that shape cellular structure.
The theoretical approach adopted here is based on the description of a chemical system in terms of intensive variables. These variables include temperature, pressure and the chemical potentials of the system. It is convenient to denote the chemical potentials by the chemical activities or fugacities of basis species, for example the activity of H+ (which defines pH) or the fugacity of oxygen. This permits comparison of the parameters of the model with reference systems described in experimental and other theoretical biochemical studies. In the following calculations, temperature and pressure were set to 25°C and 1 bar, respectively, and the logarithm of oxygen fugacity is the primary variable of interest. Below, oxidation-reduction potential and oxygen fugacity are used synonymously, and redox refers specifically to Eh. The oxidation-reduction potential of a system can be expressed in terms of Eh using an equation given in the Methods.
S. cerevisiae was chosen as a model system for the current investigation because there is abundant information about the subcellular distribution of proteins as well as independent measurements of the pH and oxidation state of some subcellular compartments. Also, the cellular development of yeast is extensively documented, which can yield other comparisons for some of the results of the model calculations.
There are two major parts to this paper. In the first part, the reactions corresponding to intercompartmental interactions between subcellular homologs (or isoforms) of particular enzymes and between reference model proteins for different compartments are quantified by calculating the oxygen fugacities for equal chemical activities of the reacting proteins in metastable equilibrium. The relative metastabilities of the reference model proteins are compared with some observations from the literature about reaction progress during the cell cycle. Specific known interactions between compartments are considered in order to derive values of the oxygen fugacity within compartments that best metastabilize the proteins contained within them. In the second part of this paper, the relative abundances of model proteins in metastable equilibrium are calculated and compared with measured abundances. The range of protein abundances in a metastable equilibrium population often approaches that seen in experiments over a narrow window of oxygen fugacity. Positive and negative correlations between the calculated and experimental relative abundances are found in some cases. The paper concludes with a summary of the findings and other implications of the hypothesis.
Presentation of the hypothesis
Part of a cell's expenditure of metabolic fuel is directed toward the formation of proteins, including their synthesis and transport to other compartments. Even when it is normalized to the lengths of the proteins, the energy required for protein formation is not a constant, but depends on the composition and environment of the protein. If these energy differences are quantified, the relative abundances of model proteins in metastable equilibrium can be calculated. The compositions of these metastable assemblages depend on local environmental variables such as oxygen fugacity, which is a scale for oxidation-reduction potential in a system. The major hypothesis adopted for this investigation is that energy minimization is a force contributing to the organization of cells; this implies the possibility of an evolutionary convergence between biomolecular composition and the chemical properties of subcellular compartments. In a first set of tests of this hypothesis, chemical reactions among model proteins in known intercompartmental interactions were used to obtain values of oxygen fugacity for subcellular compartments that can be compared with measured redox values. A second set of calculations presented here shows that the relative abundances of proteins within compartments and of those that form complexes can be correlated in some cases with metastable equilibrium assemblages. These results provide theoretical constraints on the spontaneous generation of order in the distributions of proteins within cells and imply that work done by maintaining oxidation-reduction gradients can selectively alter the degrees of formation of assemblages of proteins.
Testing the hypothesis
Relative metastabilities of subcellular homologs of redoxins
Subcellular isoforms of glutaredoxin, thioredoxin and thioredoxin reductase in yeasta.
To quantify the metastability limits of the proteins in terms of the chemical environment, one can assess the energetics of formation reactions for the proteins. At metastable equilibrium, the predominant protein in a population is the one with the highest chemical activity (in this communication, activity refers to chemical activity rather than enzymatic activity; activity is equivalent to concentration for ideal systems, where activity coefficients are unity). This statement implies that the overall formation reaction from basis species (see Methods) for the predominant protein has a lower Gibbs energy (or higher chemical affinity) than any of the others. The consequences of these relationships can be portrayed on chemical activity diagrams using a previously described procedure that is encoded in the CHNOSZ software package , which was used to perform the calculations reported below (see Methods). Additional File 1 includes the program script and data files that were used to carry out these calculations and generate the tables and figures. So that the results described in this section can be reconstructed at pH = 7, the standard molal Gibbs energies ( ) and net charges of ionized proteins at this pH are listed in Table 1.
The equal-activity lines in these pH diagrams are curved because the ionization states of the proteins depend on pH. The observation apparent in Fig. 1a that increasing log favors formation of the cytoplasmic protein homolog relative to its mitochondrial counterpart is also true for the thioredoxin/thioredoxin reductase system shown in Fig. 1b. In comparing Figs. 1a and 1b note that in the latter figure, predominance fields for a greater number of candidate proteins appear, and that the predominance field boundary between mitochondrial and cytoplasmic proteins occurs at a lower oxidation-reduction potential. The dashed lines shown in each diagram of Fig. 1 are reference lines denoting the reduction stability limit of H2O (log ≈ -83.1 at 25°C and 1 bar ).
Predominance diagrams as a function of Eh and pH for the glutaredoxin and thioredoxin/thioredoxin reductase systems are shown in Figs. 1c and 1d. Like log , Eh and pH together are a measure of the oxidation-reduction potential of the system; the different scales can be converted using Eqn. (5) in the Methods. The trapezoidal areas bounded by dotted lines in Figs. 1c and 1d show the ranges of Eh and pH corresponding to the limits of the log -pH diagrams of Figs. 1a and 1b. It can be deduced from these diagrams that if the upper log limit of Fig. 1a were extended upward, this diagram would include a portion of the predominance field for the nuclear protein GLRX3.
It appears from Figs. 1a–b that increasing increasing log at constant pH, or increasing pH at constant oxidation-reduction potential have similar consequences for the relative metastabilities of the cytoplasmic and mitochondrial homologs. In this analysis, however, pH does not appear to be a very descriptive variable; the magnitude of the effect of changing oxygen fugacity over several log units is greater than the effect of changing pH by several units. In further calculations described below pH was set to 7.
In Figs. 1e and 1f the logarithm of activity of water (log ) appears as a variable. In Fig. 1e it can be seen that the formation of a nuclear homolog of glutaredoxin is favored relative to the cytoplasmic homologs by decreasing activity of water and/or increasing oxygen fugacity, and that increasing relative metastabilities of the mitochondrial proteins are consistent with lower oxidation-reduction potentials. In Fig. 1f it appears that the formation of the thioredoxin reductase relative to thioredoxin is favored by increasing , and that for the thioredoxin the relative metastabilities of the mitochondrial proteins increase with decreasing .
Comparison with subcellular redox measurements
Nominal electrochemical characteristics of subcellular environments in eukaryotes. Values refer to yeast cells unless noted otherwise.
The current understanding of the major trends of redox states in compartments of eukaryotic cells can be summarized as, from most reducing to most oxidizing, mitochondrion – nucleus – cytoplasm – endoplasmic reticulum (ER) – extracellular . Strong redox gradients within the mitochondrion are essential to its function , which is not captured by the single values listed in Table 2. Comparison nevertheless with the computational results shown in Fig. 1 indicates that a relatively reducing environment does favor the mitochondrial homologs over the others shown in the diagram.
Measurements of GSH/GSSG concentrations point to a lower redox state in the nucleus than in the cytoplasm, but the present model has the nuclear proteins favored by relatively oxidizing conditions. Studies using nuclear magnetic resonance (NMR) showing that the hydration state of the nucleus is higher than the cytoplasm [42, 3] also seem to contradict the trend in Fig. 1e that the formation of the nuclear proteins is favored relative to their cytoplasmic counterparts by decreasing activity of water. Finally, mitochondrial pH is somewhat higher than that of the cytoplasm [37, 38], but in Figs. 1a and 1b it appears that the predicted energetic constraints favor the cytoplasmic proteins at higher pHs. These comparisons indicate that the investigated metastable equilibrium constraints are not entirely responsible for the spatial distribution of the isoforms of redoxins in the cell.
Relative metastabilities of reference model proteins
Chemical formulas, standard molal Gibbs energies and other properties for reference model proteins for compartments in yeast cellsa.
Can the relative metastabilities of proteins be linked to the order of their appearance in the cell cycle? It is noteworthy that the reference model proteins representing the two cytoskeletal systems in yeast cells, actin and microtubule, appear near opposite ends of the energy spectrum. This outcome may be consistent with the observation that actin in different forms appears to be present at most stages of the cell cycle , but that the microtubule cytoskeleton grows during anaphase (i.e., the stage of the cell cycle characterized by physical separation of the chromosomes ) and is degraded during other stages of the cell cycle [43, 44]. The outcome of the mitotic cycle in S. cerevisiae is the growth of a new cell in the form of a bud . Not all structures in the bud form simultaneously. Instead, it has been observed that "the endoplasmic reticulum, Golgi, mitochondria, and vacuoles all begin to populate the bud well before anaphase and that their segregation into the bud does not require microtubules". From Fig. 2 it is apparent that the proteins in the vacuole, ER, mitochondria and Golgi are all energetically less costly than many of their counterparts in other subcellular locations. The proteins in actin and lipid particles are also relatively metastable, which could imply that they too have a primary position in the formation of new cells. These and other potential consequences of energetic differences between the biomacromolecules in subcellular compartments have not been fully explored.
Intercompartmental protein interactions
The diagrams in Fig. 2 show the metastability limits for interactions between predominant reference model proteins for different subcellular compartments. However, many subcellular interactions may in fact be meta-metastable with respect to the reaction boundaries shown in Fig. 2. For example, interactions occur between proteins in the cytoplasm and nucleus , but the reference model proteins for these compartments do not share a reaction boundary in Fig. 2. Below, known intercompartmental interactions are combined with the oxygen fugacity requirements for equal activities of the reference model proteins to characterize compartmental oxidation-reduction potentials.
Major intercompartmental protein interactions in yeasta.
nuclear periphery-bud neck
Chemical reactions corresponding to each of the interactions listed in Table 4 are listed in Additional File 4. The values of (coefficient on O2(g)in the reactions) are listed in Table 4 together with the values of log calculated for equal chemical activities of the two reference model proteins in each reaction. Note that there are some reactions where the absolute value of is substantially smaller than the others; these include spindle pole-cytoplasm and mitochondrion-nucleus. Because of the small value of in these reactions, the values of log for equal activities of these proteins tend to be more extreme than for other reactions. The sign of denotes the thermodynamically favored direction of the reaction as log is changed from its equal-activity value; for example, at log = -74.9, the reference model proteins of actin and bud can metastably coexist with equal chemical activities, but at higher values that of actin predominates in a metastable assemblage.
The interactions listed in Table 4 were used to obtain model values of the oxygen fugacity in each compartment that are listed in Table 3. The model value of the oxygen fugacity for each compartment was selected so that in as many cases as possible the reactions listed in Table 4 favor the formation of the reference model protein for this compartment relative to those of interacting compartments. For example, the log listed for the actin compartment is -74.7, which allows this reference model protein to be metastable relative to its interacting partners in the first four reactions listed in Table 4. The limits of the model values of log were set to between ca. -70 and -80, so that some of the more extreme values listed in Table 4 were not considered in the analysis.
It is notable that at log = -75, the reference model protein for microtubule is not metastable with respect to any of its interacting partners except for bud neck. The reference model protein for microtubule only becomes relatively metastable at high oxygen fugacities (w.r.t. bud and cell periphery) or at low oxygen fugacities (w.r.t. cytoplasm and spindle pole). Hence, the value of log -75 taken here for the microtubule compartment is different from the values for other compartments, in that this represents conditions where the formation of its reference model protein is more unfavorable than that of most of its interacting partners.
Calculation of relative abundances of proteins
Above, the interactions between subcellular homologs of enzymes and reference model proteins for subcellular compartments were used to derive oxygen fugacity limits for metastable reactions of proteins in different compartments. In the second part of this study, attention is focused on the relative abundances and intracompartmental interactions of proteins.
Relative abundances of proteins within compartments
Oxygen fugacities, root mean square deviations and correlation coefficients in comparisons of intracompartmental protein interactionsa.
Most abundant proteins
ER to Golgi
The experimental abundances of the proteins reported by  are 21400, 12200, 1840, 1720 and 358, respectively, in relative units. These abundances were scaled to the same total activity of amino acid residues (unity) used in the calculations to generate the experimental relative abundances plotted at the dashed line in Fig. 4 at log = -78. Under these conditions, the metastable equilibrium abundances of the proteins do not occur in exactly the same rank order as the experimental ones, but there is a greater overall correspondence with the experimental relative abundances.
The retrieval of optimal values of log was aided by calculating the root mean square deviation (RMSD) of logarithms of activities using Eqn. (7) and the Spearman rank correlation coefficient (ρ; Eqn. 8) between experimental and calculated logarithms of activities. The dotted lines in Fig. 5 were drawn at one RMSD on either side of the one-to-one correspondence, denoted by the solid lines in this figure. The RMSD values were used to identify outliers that are identified in Fig. 5 by letters that are listed in Additional File 5. To aid in distinguishing the points, they were assigned colors on a red (reduced) – blue (oxidized) scale that reflects the average nominal oxidation state of carbon of the protein (Eqn. 6).
There is a considerable degree of scatter apparent in many of the plots shown in Fig. 5, so a low degree of certainty may be associated with the log values regressed from these comparisons. In specific cases such as peroxisome and nuclear periphery a lower overall deviation is apparent and a positive correlation appears between the calculated and experimental relative abundances. Because they were regressed from intracompartmental protein abundance data, the values of log listed in Table 5 might not be as representative of subcellular oxidation-reduction conditions as those listed in Table 3, which have the additional benefit of being based partly on known subcellular interactions (see above).
The comparisons depicted in Fig. 5 and in Additional File 6 are also significant because they reveal that the range of protein abundances observed in cells is accessible in a metastable equilibrium assemblage at some values of log . For example, the range of experimental abundances of the model proteins in actin covers about 1.6 orders of magnitude, while the calculated abundances vary over about 2.2 orders of magnitude. Extreme values of log tend to weaken this correspondence. The lowest degree of correspondence occurs for the cytoplasmic proteins, where ~5 orders of magnitude separate the predicted relative abundances of the top 50 most abundant proteins, which in the experimental measurements have a dynamic range spanning about 1.2 orders of magnitude. The great degree of scatter apparent in many of the comparisons in Fig. 5 could be partly a consequence of including in the comparisons model proteins that do not actually interact with each other, despite their high relative abundances. To address this concern, a more directed approach was adopted below that takes account of fewer numbers of proteins that are known to interact through the formation of complexes.
Relative abundances of proteins in complexes
The correspondence between the calculated and experimental relative abundances of the five model proteins in ER to Golgi raises the question of what characteristics of the proteins might be responsible for this result. Scanning the functional annotations of these proteins reveals that they are part of the COPII coat complex . The results for this model system suggested that focusing on specific complexes in other compartments could yield interesting results. Because the interactions of proteins to form complexes is essential in cellular structure and regulating the functions of enzymes , factors that affect the relative abundances of the complexing proteins may be fundamental to the control of metabolic processes.
List of selected complexesa.
cyclin-dependent protein kinase complex
actin-associated motor protein complex 2
translation initiation factor eIF3
ESCRT I & II complexes
signal recognition complex
coatomer COPII complex
Golgi transport complex
retrograde protein complex
sterol biosynthesis enzymes
mitochondrial ribosome small subunit
nuclear pore complex
small subunit processome
RNA polymerase I
integral to peroxisomal membrane
proteins localized here and early.Golgi
spindle-pole body complex
VO vacuolar ATPase complex
vacuolar proteases and other canonical proteins
The calculated metastable equilibrium logarithms of activities of the proteins in each complex are shown as a function of log in Additional File 8. The calculated logarithms of activities of the proteins were compared with experimental ones by constructing scatterplots at 0.5 log unit intervals from log = -82 to -70.5, which are shown in Additional File 6. As described above, visual assessment of fit was used in combination with the RMSD and Spearman rank correlation coefficients to obtain values of log that maximize the correspondence with experimental relative abundances. The resulting calculated relative abundances are listed together with the experimental ones in Additional File 9.
The results in Fig. 6 permit an interpretation of the relative energetic requirements for formation of different groups of interacting proteins. Take for example complex 14, which is the DASH complex that associates with the microtubule. A negative correlation between the experimental and calculated relative abundances is apparent for this complex in Fig. 6. The RMSD between calculated and experimental logarithms of activities of proteins is 1.04, which is among the highest listed in Table 5. Note from Eqn. (11) that a ~1 log unit change in the chemical activity of a chemical species corresponds to a Gibbs energy difference equal to 2.303RT. An average difference of ~1 between calculated and experimental logarithms of activity indicates that the formation of the proteins requires 2.303RT = 1364 cal mol-1 beyond what would be needed if the proteins formed in metastable equilibrium relative abundances. On the other hand, the formation in specific oxidation-reduction conditions of proteins making up other assemblages where cellular abundances positively correlate with and span the same range as the metastable equilibrium distribution can proceed close to a local minimum energy required for protein formation.
Because of their relatively high energy demands, proteins in complexes such as the DASH complex and the spindle pole body are likely to be more dynamic in the cell. (Note that although a positive rank correlation coefficient for the latter complex is reported in Table 5, at a lower oxygen fugacity (log = -77.5) an inverse correlation results between experimental abundances and calculated metastable equilibrium relative abundances of the proteins in this complex; see Additional File 6). The finding made elsewhere of some inverse relationships between relative abundance of proteins and corresponding mRNA levels was interpreted as evidence for additional effort on the part of the cell . An inverse relationship that opposes equilibrium may be favored in evolution because of the strategic advantage of incorporating otherwise costly (rare) amino acids that increase enzymatic diversity .
The differences in the numbers of proteins considered in each of the comparisons implies that the values of the correlation coefficients are not directly comparable. The p-values for each of the correlations listed in Table 5 were calculated and are reported in Additional File 10. The p-value is the probability that the value of the observed correlation coefficient can be met or exceeded by a random configuration of the system. The present calculations suggest that the lowest p-values are associated with the collections of greater than ca. 40 proteins listed in Table 5 that have a Spearman rank correlation coefficient greater than ca. 0.4. Using the p-value as a criterion, the most convincing demonstrations of the existence of correlations appear in the most abundant proteins of the vacuolar membrane and cell periphery. By comparison, the smaller systems of proteins making up complexes, which in some cases have higher correlation coefficients, also have relatively high p-values, indicating a greater probability that the same result can be obtained in a random configuration.
Implications of the hypothesis
Subcellular homologs of glutaredoxin and thioredoxin are metastable at different log ranges. The mitochondrial homolog appears to be more reduced, and the nuclear one most oxidized. Reactions within the glutaredoxin system also exhibit sensitivity to hydration state.
Reference model proteins for 23 subcellular locations also have metastability limits in log space. The relationships are consistent with a relatively oxidized nuclear reference model protein, but that for the mitochondrion is intermediate between the nucleus and the cytoplasm. Among the reference model proteins predicted to be most stable (Fig. 2a–b), four are known to be involved in the early stages of formation of the bud. Golgi is predicted to be a reduced compartment while actin and the vacuole are relatively oxidized.
The least stable reference model proteins are those for the microtubule and bud neck (Fig. 2f and Fig. 3). This observation suggests that the proteins in the microtubule are very reactive with other cellular components, and/or that they have a relatively high turnover rate.
Observed trends in the relative abundances of the most abundant proteins in some compartments can be correlated with the relative abundances of proteins predicted using a metastable equilibrium model. Correlations between observed and predicted relative abundances for smaller numbers of proteins that make up complexes can also be documented. In some cases negative correlations may be supported, such as for the DASH complex (microtubule), translation initiation factor, and the early Golgi SNARE complex. The maintenance of these complexes might entail a higher energy demand than for others in the cell.
If the hypothesis adopted for this study was true, it would imply that there are processes that impart an energetic bias on the appearance of proteins in specific compartments. The thermodynamic model described above by itself gives no information about the possible nature of processes involved. Two processes that could be important are the work against diffusional gradients required for active protein transport and the turnover rates of subcellular populations of proteins. Regarding the former, it is the gradient of chemical potential (not concentration) of the biomacromolecule, that appears in statements such as Fick's Law. Differences in the chemical conditions between compartments would be expected to differentially contribute to the activity coefficients of proteins, so that the cost of transport to various compartment is not equal. Although the activity coefficients of proteins were not considered in this study, their values might depend on oxidation or hydration potential, so the the current results could be implicitly influenced by nonideality in the subcellular system. Regarding the latter process, one may expect that the turnover rates of proteins are tied to the local chemical environment. If, for example, the turnover rate of a population of proteins minimizes at a specific oxygen fugacity, then any deviations away from this oxidation potential would increase the turnover rate and cause the cell to expend more energy in maintaining this population.
If chemical energy minimization by the cell results in an increase in fitness, the energetic effects of the physical-chemical processes outlined above may constrain the overall process of natural selection. Therefore, the hypothesis also implies that an expected outcome of evolution is the formation of biomacromolecules with lower energy demands compared with other possible, and otherwise equal, products. How does this connect with the mechanism of protein transport and trafficking, i.e. that from their place of synthesis (ribosomes) proteins are transported to different subcellular locations, often under the influence of specific signal sequences? At one extreme, it is possible that the energetic differences between compartments are not influential in the evolution of the mechanism of protein sorting and trafficking. However, if the signal sequences themselves are chemically reactive to varying degrees depending on their subcellular environment, then selection for mutations in them might be tuned to both function and chemistry. It would not be surprising then to find evidence for the chemical adaptation of signal sequences to specific compartments.
These results and observations support the notion that changing oxidation-reduction potential can selectively alter the potential for reactions leading to formation of proteins and their complexes. Chemical selectivity in the dynamic formation in the cell of high-energy proteins could lead to transient formation of complexes that function only under certain conditions. Because of the different stability limits of the reference model proteins in log space, these results also in principle support the notion that "a fundamental redox attractor underpins ... core cellular processes". In reality, many chemical properties vary spatially in cells, including the hydration state, pH, activities of CO2 and H2S, and temperature and pressure in the extracellular environment. These all factor into the Gibbs energy changes accompanying the chemical transformations between proteins, as do the thermodynamic properties of protein folding reactions and nonideality in protein solutions. Because of its energetic basis, the model used here can be extended in the future to incorporate the effects of these variables. Building these relationships into a multidimensional thermodynamic assessment is a promising avenue for predicting the chemical features of proteomic adaptation in the context of the cellular environment.
The essential steps in the calculations reported here are 1) defining standard states, 2) identifying model proteins for systems of interest, 3) assessing the relative abundances of model proteins in metastable equilibrium, 4) visualizing the results of the calculations on chemical diagrams and 5) comparing the computational results with experimental biochemical and proteomic data.
Standard states and chemical activities
where R and T represent, respectively, the gas constant and the temperature, μ and μ○ stand for the chemical potential and standard chemical potential, respectively, and a denotes activity. No provision for activity coefficients of proteins or other species was used in this study; under this approximation, the activity of an aqueous species is equal to its concentration (molality).
The standard state for aqueous species including proteins specifies unit activity of the aqueous species in hypothetical one molal solution referenced to infinite dilution. The standard molal Gibbs energies of the proteins were calculated with the CHNOSZ software package  using group additivity properties and parameters taken from Ref. .
Reference model proteins for amino acid compositions
The overall amino acid compositions of proteins in 23 subcellular locations in S. cerevisiae were calculated by combining localization  and abundance  data for proteins measured in the YeastGFP project with amino acid compositions of proteins downloaded from the Saccharomyces Genome Database (SGD) . Of 4155 ORF names listed in the YeastGFP dataset, all but 12 are present in SGD (the missing ones are YAR044W, YBR100W, YDR474C, YFL006W, YFR024C, YGL046W, YGR272C, YJL012C-A, YJL017W, YJL018W, YJL021C and YPR090W).
To generate reference model proteins that are most representative of each compartment, proteins that were annotated in the YeastGFP study as being localized to more than one compartment were excluded from this analysis (except for bud; see below), as were those for which no abundance was reported. The names of the open reading frames (ORFs) corresponding to the proteins in the YeastGFP data set were matched against the SGD's protein_properties.tab file downloaded on 2008-08-04. This search yielded a number of model proteins for each compartment, ranging from 5 (ER to Golgi) to 746 (cytoplasm); see Table 3. The names of the compartments used throughout the tables and figures in this paper correspond to the notation used in the YeastGFP data files.
It was found that no proteins with reported abundances and localized to the bud were exclusive to that compartment, hence all of the proteins localized there (which also have localizations in other compartments) were taken as models for the bud reference model protein. The amino acid composition of the reference model protein for each compartment was calculated by taking the sum of the compositions of each model protein for a compartment in proportion to its fractional abundance in the total model protein population of the compartment. The resulting amino acid compositions are listed in Additional File 2. The corresponding chemical formulas of the nonionized reference model proteins and the calculated standard molal Gibbs energies of formation from the elements at 25°C and 1 bar of the ionized reference model proteins are shown in Table 3.
Diagrams showing the predominant proteins and the relative abundances of proteins in metastable equilibrium were generated using the CHNOSZ software package . These calculations take account of formation reactions of the proteins written for their residue equivalents . An example of this approach is described further below for a specific model system.
The basis species appearing in the formation reactions studied here are CO2(aq), H2O, NH3(aq), O2(g), H2S(aq)and H+. The reference activities used for the basis species were 10-3, 100, 10-4, 10-7 and 10-7, respectively, for CO2(aq), H2O, NH3(aq), H2S(aq)and H+. In the case of diagrams showing Eh as a variable, the aqueous electron (e-) was substituted for O2(g)in the basis species. Reference values for or are not listed here because one or the other is used as an independent variable in each of the calculations described above.
Conversion between scales of oxidation-reduction potential
At 25°C and 1 bar, F/2.303RT = 16.903 volt-1 and log K2 = -41.55; for pH = 7 and log = 0, a value of Eh = 0 V corresponds to log = -55. Eqn. (5) permits the conversion between Eh and log as well at other temperatures, pHs, and activities of H2O.
Average nominal oxidation state of carbon
Eqn. (6) is consistent with the electronegativity rules described in  and is compatible with the equation for average oxidation number of carbon used in . For example, Eqn. (6) can be used to calculate the average nominal oxidation states of carbon in CO2 and CH4, which are +4 and -4, respectively. Note that the proportions of oxygen and other covalently bonded heteroatoms contribute to the value of of a protein or other molecule, but that proton ionization does not alter the nominal carbon oxidation state, because of the opposite contributions from Z and nH in Eqn. (6). In the 4143 proteins identified in the YeastGFP subcellular localization study and found in the Saccharomyces Genome Database, the minimum and maximum of are -0.414 and 0.390, respectively. Of the proteins in this dataset, six have < -0.35 (YDR193W, YDR276C, YEL017C-A, YJL097W, YML007C-A, YMR292W) and six have > 0.15 (YCL028W, YHR053C, YHR055C, YKR092C, YMR173W, YPL223C). The points in the scatterplots in this paper (Figs. 5 and 6 and Additional File 6) are colored on a continuous red-blue scale according to the value of of the proteins, where maximum red occurs at = -0.35 and maximum blue occurs at = 0.15.
Comparison with experimental relative abundances
where Xcalc, iand Xexpt, idenote the calculated and experimental logarithms of activities and n stands for the number of proteins. In the calculations described above, experimental abundances of proteins in each model system were scaled so that the total chemical activity of amino acid residues was equal to unity.
Calculating relative abundances of proteins in metastable equilibrium
The following example demonstrates the procedure used to calculate the relative abundances of proteins in metastable equilibrium. The model proteins for ER to Golgi, in order of decreasing abundance in the cell reported by , are YLR208W, YHR098C, YDL195W, YNL049C and YPL085W. (For simplicity, the proteins are identified here by the names of the open reading frames (ORF).) The formula of the uncharged form of the first protein, YLR208W, is C1485H2274N400O449S4, and its amino acid sequence length is 297 residues. The standard molal Gibbs energy of formation from the elements ( ) of this protein at 25°C and 1 bar calculated using group additivity  is -10670 kcal mol-1. At this temperature and pressure and at pH = 7, group additivity can also be used to calculate the charge of the protein (-10.8832) and the standard molal Gibbs energy of formation from the elements of the charged protein (-10880 kcal mol-1). The formula of the protein in this ionization state is C1485H2263.1168N400O449 . Dividing by the length of the protein, we find that the formula and standard molal Gibbs energy of formation from the elements of the residue equivalent of YLR208W are C5.0000H7.6199N1.3468O1.5118 and -36.633 kcal mol-1, respectively.
At 929 residues, YHR098C is over 3 times as long as YLR208W, but in the formation reactions from the basis species of the residue equivalents of the two proteins, the coefficients on the basis species are similar. The difference between the coefficients of the same basis species in the reactions signifies the response of the metastable equilibrium assemblage to changes in the corresponding chemical activity or fugacity. For example, because and increasing , or at constant T, P and chemical activities of the other basis species shifts the metastable equilibrium in favor of YLR208W at the expense of YHR098C. Here, ν i denotes the reaction coefficient of the i th basis species or protein, which is negative for reactants and positive for products as written. Conversely, because and increasing , or (decreasing pH) at constant T, P and chemical activities of the other basis species shifts the metastable equilibrium in favor of YHR098C at the expense of YLR208W. The magnitude of the effect is proportional to the size of the difference between the coefficients of the basis species in the reactions, and it can be quantified for a specific model system using the following calculations.
In Eqn. (12) the values of of O2(g)and H+ are both zero, which are consistent with the standard state conventions for gases and the hydrogen ion convention used in solution chemistry. The values of of the other basis species are taken from the literature [61–63]. The value of log K9 consistent with Eqn. (12) is -392.19.
The values of a i used to write Eqn. (13) are the reference values listed in the Methods for and . The value of used in Eqn. (13) (log = -75.3) is also a reference value that, it will be shown, characterizes a metastable equilibrium distribution of proteins that is rank-identical to the measured relative abundances of the proteins. Finally, the value of a of the residue equivalent of the protein in Eqn. (13) is set to a reference value of unity (log a = 0). If we are only concerned with the relative abundances of the proteins in metastable equilibrium, the actual value used here does not matter so long as it is the same in the analogous calculations for the other proteins.
where a t denotes the total activity of residue equivalents in the system and n stands for the number of proteins in the system. Note regarding the left-hand side of Eqn. (14) that because we are taking activity coefficients of unity, the ratio a i /a t is equal to the ratio of concentrations of residue equivalents in the system. No negative sign appears in front of A/RT in the exponents Eqn. (14) because the chemical affinity is the negative of Gibbs energy change of the reaction. Note in addition that the values of A/2.303RT given above must be multiplied by ln 10 = 2.303 before being substituted in Eqn. (14). By taking a t = 1, we can combine Eqn. (14) with A/RT of each of the formation reactions to calculate chemical activities of the residue equivalents of the proteins equal to 0.0905, 0.2248, 0.2994, 0.1944 and 0.1909, respectively. The lengths of the proteins are 297, 929, 1273, 876 and 2195, so the corresponding logarithms of activities of the proteins are e.g. log (0.0905/297) = -3.52 for YLR208W, and -3.61, -3.63, -3.65 and -4.06 for the remaining proteins, respectively.
The comments of two anonymous reviewers helped to improve this paper. This material is based upon work supported by the National Science Foundation under grant EAR-0309829 and the Department of Energy under grant DE-FG02-03ER151418, both awarded to the author's Ph.D. advisor, Professor Harold C. Helgeson. His enthusiasm for thermodynamics made this project possible.
- Preston RA, Murphy RF, Jones EW: Assay of vacuolar pH in yeast and identification of acidification-defective mutants. Proc Natl Acad Sci USA. 1989, 86 (18): 7027-7031.PubMed CentralView ArticlePubMedGoogle Scholar
- Hwang C, Sinskey AJ, Lodish HF: Oxidized redox state of glutathione in the endoplasmic reticulum. Science. 1992, 257 (5076): 1496-1502.View ArticlePubMedGoogle Scholar
- Morrill GA, Kostellow AB, Osterlow K, Gupta RK: Differences in hydration state of nucleus and cytoplasm of the amphibian oocyte. J Membrane Biol. 1996, 153: 45-51.View ArticleGoogle Scholar
- Al-Habori M: Microcompartmentation, metabolic channelling and carbohydrate metabolism. Int J Biochem Cell Biol. 1995, 27 (2): 123-132.View ArticlePubMedGoogle Scholar
- Aw TY: Intracellular compartmentation of organelles and gradients of low molecular weight species. Microcompartmentation and Phase Separation in Cytoplasm, Int. Rev. Cytol. Edited by: Walter H, Brooks DE, Srere PA. 2000, 192: 223-253. San Diego: Academic PressGoogle Scholar
- Cedano J, Aloy P, Pérez-Pons JA, Querol E: Relation between amino acid composition and cellular location of proteins. J Mol Biol. 1997, 266: 594-600.View ArticlePubMedGoogle Scholar
- Andrade MA, O'Donoghue SI, Rost B: Adaptation of protein surfaces to subcellular location. J Mol Biol. 1998, 276 (2): 517-525.View ArticlePubMedGoogle Scholar
- Huh WK, Falvo JV, Gerke LC, Carroll AS, Howson RW, Weissman JS, O'Shea EK: Global analysis of protein localization in budding yeast. Nature. 2003, 425 (6959): 686-691.View ArticlePubMedGoogle Scholar
- Ghaemmaghami S, Huh W, Bower K, Howson RW, Belle A, Dephoure N, O'Shea EK, Weissman JS: Global analysis of protein expression in yeast. Nature. 2003, 425 (6959): 737-741.View ArticlePubMedGoogle Scholar
- Gasch AP, Spellman PT, Kao CM, Carmel-Harel O, Eisen MB, Storz G, Botstein D, Brown PO: Genomic expression programs in the response of yeast cells to environmental changes. Mol Biol Cell. 2000, 11 (12): 4241-4257.PubMed CentralView ArticlePubMedGoogle Scholar
- Schekman R: Protein localization and membrane traffic in yeast. Annu Rev Cell Biol. 1985, 1: 115-143.View ArticlePubMedGoogle Scholar
- Doxsey S, McCollum D, Theurkauf W: Centrosomes in cellular regulation. Annu Rev Cell Dev Biol. 2005, 21: 411-434.View ArticlePubMedGoogle Scholar
- Halvorson H: Intracellular protein and nucleic acid turnover in resting yeast cells. Biochim Biophys Acta. 1958, 27 (2): 255-266.View ArticlePubMedGoogle Scholar
- Morowitz HJ: Foundations of Bioenergetics. 1978, New York: Academic PressGoogle Scholar
- Seligmann H: Cost-minimization of amino acid usage. J Mol Evol. 2003, 56 (2): 151-161.View ArticlePubMedGoogle Scholar
- Swire J: Selection on synthesis cost affects interprotein amino acid usage in all three domains of life. J Mol Evol. 2007, 64 (5): 558-571.View ArticlePubMedGoogle Scholar
- Berezovsky IN, Zeldovich KB, Shakhnovich EI: Positive and negative design in stability and thermal adaptation of natural proteins. PLoS Comput Biol. 2007, 3 (3): 498-507.View ArticleGoogle Scholar
- Kondepudi DK, Prigogine I: Modern Thermodynamics: From Heat Engines to Dissipative Structures. 1998, New York: John Wiley & SonsGoogle Scholar
- Dick JM: Calculation of the relative metastabilities of proteins using the CHNOSZ software package. Geochem Trans. 2008, 9: 10-PubMed CentralView ArticlePubMedGoogle Scholar
- Wicken JS: A thermodynamic theory of evolution. J Theor Biol. 1980, 87: 9-23.View ArticlePubMedGoogle Scholar
- Aita T, Husimi Y: Fitness spectrum among random mutants on Mt. Fuji-type fitness landscape. J Theor Biol. 1996, 182 (4): 469-485.View ArticlePubMedGoogle Scholar
- Demetrius L, Ziehe M: Darwinian fitness. Theor Popul Biol. 2007, 72 (3): 323-345.View ArticlePubMedGoogle Scholar
- Dick JM, LaRowe DE, Helgeson HC: Temperature, pressure, and electrochemical constraints on protein speciation: Group additivity calculation of the standard molal thermodynamic properties of ionized unfolded proteins. Biogeosciences. 2006, 3 (3): 311-336.View ArticleGoogle Scholar
- Pedrajas JR, Porras P, Martínez-Galisteo E, Padilla CA, Miranda-Vizuete A, Bárcena JA: Two isoforms of Saccharomyces cerevisiae glutaredoxin 2 are expressed in vivo and localize to different subcellular compartments. Biochem J. 2002, 364: 617-623.PubMed CentralView ArticlePubMedGoogle Scholar
- Molina MM, Bellí G, de la Torre MA, Rodríguez-Manzaneque MT, Herrero E: Nuclear monothiol glutaredoxins of Saccharomyces cerevisiae can function as mitochondrial glutaredoxins. J Biol Chem. 2004, 279: 51923-51930.View ArticlePubMedGoogle Scholar
- Herrero E, Ros J, Tamarit J, Belli G: Glutaredoxins in fungi. Photosynth Res. 2006, 89: 127-140.View ArticlePubMedGoogle Scholar
- Pedrajas JR, Kosmidou E, Miranda-Vizuete A, Gustafsson JA, Wright APH, Spyrou G: Identification and functional characterization of a novel mitochondrial thioredoxin system in Saccharomyces cerevisiae. J Biol Chem. 1999, 274 (10): 6366-6373.View ArticlePubMedGoogle Scholar
- Trotter EW, Grant CM: Overlapping roles of the cytoplasmic and mitochondrial redox regulatory systems in the yeast Saccharomyces cerevisiae. Eukaryot Cell. 2005, 4: 392-400.PubMed CentralView ArticlePubMedGoogle Scholar
- Garrels RM: Mineral Equilibria. 1960, New York: Harper & BrothersGoogle Scholar
- Dahm LJ, Jones DP: Rat jejunum controls luminal thiol-disulfide redox. J Nutr. 2000, 130 (11): 2739-2745.PubMedGoogle Scholar
- Trotter EW, Grant CM: Non-reciprocal regulation of the redox state of the glutathione-glutaredoxin and thioredoxin systems. EMBO Rep. 2003, 4: 184-188.PubMed CentralView ArticlePubMedGoogle Scholar
- Drakulic T, Temple MD, Guido R, Jarolim S, Breitenbach M, Attfield PV, Dawes IW: Involvement of oxidative stress response genes in redox homeostasis, the level of reactive oxygen species, and ageing in Saccharomyces cerevisiae. FEMS Yeast Res. 2005, 5 (12): 1215-1228.View ArticlePubMedGoogle Scholar
- Hanson GT, Aggeler R, Oglesbee D, Cannon M, Capaldi RA, Tsien RY, Remington SJ: Investigating mitochondrial redox potential with redox-sensitive green fluorescent protein indicators. J Biol Chem. 2004, 279 (13): 13044-13053.View ArticlePubMedGoogle Scholar
- Macville M, Schröck E, Padilla-Nash H, Keck C, Ghadimi BM, Zimonjic D, Popescu N, Ried T: Comprehensive and definitive molecular cytogenetic characterization of HeLa cells by spectral karyotyping. Cancer Res. 1999, 59: 141-150.PubMedGoogle Scholar
- ATCC: The Global Bioresource Center: Product Description. CRL-1606. 2008, http://www.atcc.org/tabid/452/Default.aspx?ATCCNum=CRL-1606&Template=cellBiologyGoogle Scholar
- Mojaverian P: Evaluation of gastrointestinal pH and gastric residence time via the Heidelberg Radiotelemetry Capsule: Pharmaceutical application. Drug Dev Res. 1996, 38 (2): 73-85.View ArticleGoogle Scholar
- Imai T, Ohno T: Measurement of yeast intracellular pH by image processing and the change it undergoes during growth phase. J Biotech. 1995, 38 (2): 165-172.View ArticleGoogle Scholar
- Llopis J, McCaffery JM, Miyawaki A, Farquhar MG, Tsien RY: Measurement of cytosolic, mitochondrial, and Golgi pH in single living cells with green fluorescent proteins. Proc Natl Acad Sci USA. 1998, 95 (12): 6803-6808.PubMed CentralView ArticlePubMedGoogle Scholar
- Singh A, Kaur N, Kosman DJ: The metalloreductase Fre6p in Fe-Efflux from the yeast vacuole. J Biol Chem. 2007, 282 (39): 28619-28626.View ArticlePubMedGoogle Scholar
- Hansen JM, Go YM, Jones DP: Nuclear and mitochondrial compartmentation of oxidative stress and redox signaling. Annu Rev Pharmacol Toxicol. 2006, 46: 215-234.View ArticlePubMedGoogle Scholar
- Go YM, Jones DP: Redox compartmentalization in eukaryotic cells. Biochim Biophys Acta-Gen Subj. 2008, 1780 (11): 1271-1290.View ArticleGoogle Scholar
- Päuser S, Zschunke A, Khuen A, Keller K: Estimation of water content and water mobility in the nucleus and cytoplasm of Xenopus laevis oocytes by NMR spectroscopy. Magn Reson Imaging. 1995, 13 (2): 269-276.View ArticlePubMedGoogle Scholar
- Botstein D, Amberg D, Mulholland J, Huffaker T, Adams A, Drubin D, Stearns T: The yeast cytoskeleton. The Molecular and Cellular Biology of the Yeast Saccharomyces: Cell Cycle and Cell Biology. Edited by: Pringle JR, Broach JR, Jones EW. 1997, 1-90. New York: Cold Spring Harbor Laboratory PressGoogle Scholar
- Alberts B, Bray D, Lewis J, Raff M, Roberts K, Watson JD: Molecular Biology of the Cell. 1989, New York: Garland Publishing, Inc, 2Google Scholar
- Lew DJ, Weinert T, Pringle JR: Cell cycle control in Saccharomyces cerevisiae. The Molecular and Cellular Biology of the Yeast Saccharomyces: Cell Cycle and Cell Biology. Edited by: Pringle JR, Broach JR, Jones EW. 1997, 607-695. New York: Cold Spring Harbor Laboratory PressGoogle Scholar
- Wente SR, Gasser SM, Caplan AJ: The nucleus and nucleocytoplasmic transport in Saccharomyces cerevisiae. The Molecular and Cellular Biology of the Yeast Saccharomyces: Cell Cycle and Cell Biology. Edited by: Pringle JR, Broach JR, Jones EW. 1997, 471-546. New York: Cold Spring Harbor Laboratory PressGoogle Scholar
- Kaiser CA, Gimeno RE, Shaywitz DA: Protein secretion, membrane biogenesis, and endocytosis. The Molecular and Cellular Biology of the Yeast Saccharomyces: Cell Cycle and Cell Biology. Edited by: Pringle JR, Broach JR, Jones EW. 1997, 91-227. New York: Cold Spring Harbor Laboratory PressGoogle Scholar
- Jones EW, Webb GC, Hiller MA: Biogenesis and function of the yeast vacuole. The Molecular and Cellular Biology of the Yeast Saccharomyces: Cell Cycle and Cell Biology. Edited by: Pringle JR, Broach JR, Jones EW. 1997, 363-470. New York: Cold Spring Harbor Laboratory PressGoogle Scholar
- Lazarow PB, Kunau W: Peroxisomes. The Molecular and Cellular Biology of the Yeast Saccharomyces: Cell Cycle and Cell Biology. Edited by: Pringle JR, Broach JR, Jones EW. 1997, 547-605. New York: Cold Spring Harbor Laboratory PressGoogle Scholar
- Pon L, Schatz G: Biogenesis of yeast mitochondria. The Molecular and Cellular Biology of the Yeast Saccharomyces: Genome Dynamics, Protein Synthesis, and Energetics. Edited by: Broach JR, Pringle JR, Jones EW. 1991, 333-406. New York: Cold Spring Harbor Laboratory PressGoogle Scholar
- SGD Project: Saccharomyces Genome Database. 2007, http://www.yeastgenome.orgGoogle Scholar
- Sarry JE, Chen S, Collum RP, Liang S, Peng M, Lang A, Naumann B, Dzierszinskil F, Yuan CX, Hippler M, Rea PA: Analysis of the vacuolar luminal proteome of Saccharomyces cerevisiae. FEBS J. 2007, 274 (16): 4287-4305.View ArticlePubMedGoogle Scholar
- Mo CQ, Bard M: Erg28p is a key protein in the yeast sterol biosynthetic enzyme complex. J Lipid Res. 2005, 46 (9): 1991-1998.View ArticlePubMedGoogle Scholar
- Tuller T, Kupiec M, Ruppin E: Determinants of protein abundance and translation efficiency in S. cerevisiae. PLoS Comput Biol. 2007, 3 (12): 2510-2519.View ArticleGoogle Scholar
- Wicken JS: Evolution, Thermodynamics, and Information. 1987, Oxford University PressGoogle Scholar
- Murray DB: On the temporal self-organisation of Saccharomyces cerevisae. Curr Genomics. 2004, 5 (8): 665-671.View ArticleGoogle Scholar
- Drever JI: The Geochemistry of Natural Waters. 1997, Upper Saddle River, New Jersey: Prentice Hall, 3Google Scholar
- Hendrickson JB, Cram DJ, Hammond GS: Organic Chemistry. 1970, New York: McGraw-Hill, 3Google Scholar
- Buvet R: General criteria for the fulfillment of redox reactions. Bioelectrochemistry I: Biological Redox Reactions, of Ettore Majorana International Science Series. Edited by: Milazzo G, Blank M. 1983, 11: 15-50. New York: Plenum PressView ArticleGoogle Scholar
- Prigogine I, Defay R: Chemical Thermodynamics. 1954, London: Longmans, Green and CoGoogle Scholar
- Helgeson HC, Kirkham DH: Theoretical prediction of the thermodynamic behavior of aqueous electrolytes at high pressures and temperatures: I. Summary of the thermodynamic/electrostatic properties of the solvent. Am J Sci. 1974, 274 (10): 1089-1198.View ArticleGoogle Scholar
- Wagman DD, Evans WH, Parker VB, Schumm RH, Halow I, Bailey SM, Churney KL, Nuttall RL: The NBS tables of chemical thermodynamic properties. Selected values for inorganic and C1 and C2 organic substances in SI units. J Phys Chem Ref Data. 1982, 11: 1-392.View ArticleGoogle Scholar
- Shock EL, Helgeson HC, Sverjensky DA: Calculation of the thermodynamic and transport properties of aqueous species at high pressures and temperatures: Standard partial molal properties of inorganic neutral species. Geochim Cosmochim Acta. 1989, 53 (9): 2157-2183.View ArticleGoogle Scholar
- Boeckmann B, Bairoch A, Apweiler R, Blatter MC, Estreicher A, Gasteiger E, Martin MJ, Michoud K, O'Donovan C, Phan I, Pilbout S, Schneider M: The SWISS-PROT protein knowledgebase and its supplement TrEMBL in 2003. Nucleic Acids Res. 2003, 31: 365-370.PubMed CentralView ArticlePubMedGoogle Scholar
- Shock EL, Sassani DC, Willis M, Sverjensky DA: Inorganic species in geologic fluids: Correlations among standard molal thermodynamic properties of aqueous ions and hydroxide complexes. Geochim Cosmochim Acta. 1997, 61 (5): 907-950.View ArticlePubMedGoogle Scholar
- Chan P, Lovrić J, Warwicker J: Subcellular pH and predicted pH-dependent features of proteins. Proteomics. 2006, 6 (12): 3494-3501.View ArticlePubMedGoogle Scholar
- Wu MM, Llopis J, Adams S, McCaffery JM, Kulomaa MS, Machen TE, Moore HPH, Tsien RY: Organelle pH studies using targeted avidin and fluorescein-biotin. Chem Biol. 2000, 7 (3): 197-209.View ArticlePubMedGoogle Scholar
- Welch MD, Iwamatsu A, Mitchison TJ: Actin polymerization is induced by Arp2/3 protein complex at the surface of Listeria monocytogenes. Nature. 1997, 385 (6613): 265-269.View ArticlePubMedGoogle Scholar
- Mullins RD, Pollard TD: Structure and function of the Arp2/3 complex. Curr Opin Struct Biol. 1999, 9 (2): 244-249.View ArticlePubMedGoogle Scholar
- Schmid M, Jaedicke A, Du TG, Jansen RP: Coordination of endoplasmic reticulum and mRNA localization to the yeast bud. Curr Biol. 2006, 16 (15): 1538-1543.View ArticlePubMedGoogle Scholar
- Frazier JA, Wong ML, Longtine MS, Pringle JR, Mann M, Mitchison TJ, Field C: Polymerization of purified yeast septins: Evidence that organized filament arrays may not be required for septin function. J Cell Biol. 1998, 143 (3): 737-749.PubMed CentralView ArticlePubMedGoogle Scholar
- Burri L, Lithgow T: A complete set of SNAREs in yeast. Traffic. 2004, 5: 45-52.View ArticlePubMedGoogle Scholar
- Kostelansky MS, Schluter C, Tam YYC, Lee S, Ghirlando R, Beach B, Conibear E, Hurley JH: Molecular architecture and functional model of the complete yeast ESCRT-I heterotetramer. Cell. 2007, 129 (3): 485-498.PubMed CentralView ArticlePubMedGoogle Scholar
- Hierro A, Sun J, Rusnak AS, Kim J, Prag G, Emr SD, Hurley JH: Structure of the ESCRT-II endosomal trafficking complex. Nature. 2004, 431 (7005): 221-225.View ArticlePubMedGoogle Scholar
- Conibear E, Cleck JN, Stevens TH: Vps51p mediates the association of the GARP (Vps52/53/54) complex with the late Golgi t-SNARE Tlg1p. Mol Biol Cell. 2003, 14 (4): 1610-1623.PubMed CentralView ArticlePubMedGoogle Scholar
- Miranda JJL, De Wulf P, Sorger PK, Harrison SC: The yeast DASH complex forms closed rings on microtubules. Nat Struct Mol Biol. 2005, 12 (2): 138-143.View ArticlePubMedGoogle Scholar
- Rout MP, Aitchison JD, Suprapto A, Hjertaas K, Zhao YM, Chait BT: The yeast nuclear pore complex: Composition, architecture, and transport mechanism. J Cell Biol. 2000, 148 (4): 635-651.PubMed CentralView ArticlePubMedGoogle Scholar
- Bernstein KA, Gallagher JEG, Mitchell BM, Granneman S, Baserga SJ: The small-subunit processome is a ribosome assembly intermediate. Eukaryot Cell. 2004, 3 (6): 1619-1626.PubMed CentralView ArticlePubMedGoogle Scholar
- Vinh DBN, Kern JW, Hancock WO, Howard J, Davis TN: Reconstitution and characterization of budding yeast γ-tubulin complex. Mol Biol Cell. 2002, 13 (4): 1144-1157.PubMed CentralView ArticlePubMedGoogle Scholar
- Gavin AC, Aloy P, Grandi P, Krause R, Boesche M, Marzioch M, Rau C, Jensen LJ, Bastuck S, Dumpelfeld B, Edelmann A, Heurtier MA, Hoffman V, Hoefert C, Klein K, Hudak M, Michon AM, Schelder M, Schirle M, Remor M, Rudi T, Hooper S, Bauer A, Bouwmeester T, Casari G, Drewes G, Neubauer G, Rick JM, Kuster B, Bork P, Russell RB, Superti-Furga G: Proteome survey reveals modularity of the yeast cell machinery. Nature. 2006, 440 (7084): 631-636.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.