Literature DB >> 36147951

Regulation of yeast Snf1 (AMPK) by a polyhistidine containing pH sensing module.

Kobi J Simpson-Lavy1, Martin Kupiec1.   

Abstract

Cellular regulation of pH is crucial for internal biological processes and for the import and export of ions and nutrients. In the yeast Saccharomyces cerevisiae, the major proton pump (Pma1) is regulated by glucose. Glucose is also an inhibitor of the energy sensor Snf1/AMPK, which is conserved in all eukaryotes. Here, we demonstrate that a poly-histidine (polyHIS) tract in the pre-kinase region (PKR) of Snf1 functions as a pH-sensing module (PSM) and regulates Snf1 activity. This regulation is independent from, and unaffected by, phosphorylation at T210, the major regulatory control of Snf1, but is controlled by the Pma1 plasma-membrane proton pump. By examining the PKR from additional yeast species, and by varying the number of histidines in the PKR, we determined that the polyHIS functions progressively. This regulation mechanism links the activity of a key enzyme with the metabolic status of the cell at any given moment.
© 2022 The Authors.

Entities:  

Keywords:  Biological sciences; biochemistry; molecular biology

Year:  2022        PMID: 36147951      PMCID: PMC9486060          DOI: 10.1016/j.isci.2022.105083

Source DB:  PubMed          Journal:  iScience        ISSN: 2589-0042


Introduction

pH homeostasis is vital for cellular processes and life. Whereas vertebrates regulate both intracellular and plasma pH with a bicarbonate/carbon dioxide buffer system (Burggren and Bautista, 2019), plants (Cosse and Seidel, 2021) and fungi utilize proton pumps to export protons from the cytoplasm into the vacuole and extracellular medium (Kane, 2016). Cytoplasmic pH regulates many processes, such as glycolysis (van Leemputte et al., 2020), membrane biogenesis (Young et al., 2010), storage of metabolic enzymes as filaments (Petrovska et al., 2014) and even solidification of the cytoplasm with a concomitant reduction of macromolecular diffusion (Munder et al., 2016) and in metazoans, apoptosis (Lagadic-Gossmann et al., 2004) . The proton gradient established is used to drive transport of nutrients and minerals such as monovalent cations (Ariño et al., 2019), amino acids (Saliba et al., 2018), glycerol (Ferreira et al., 2005), hexoses (Lagunas, 1993) and nucleosides (Loewen et al., 2003). Yeasts maintain an internal pH of 7.2 during growth on glucose by the pumping of protons from the cytoplasm across the plasma membrane by the Pma1 pump (which is the most abundant plasma membrane protein and the highest consumer of ATP in S. cerevsiae (Serrano et al., 1986; Zhao et al., 2021)), and into the vacuole by Vma1. On glucose starvation, these pumps are inactivated, leading to a rapid acidification of the cytoplasm to pH 5.7 (Martínez-Muñoz and Kane, 2008; Orij et al., 2009; Dechant et al., 2010). Glucose causes a slow increase in expression of Pma1 (Rao et al., 1993; García-Arranz et al., 1994) and activates its phosphorylation at S899 (increasing affinity of Pma1 for ATP) and at S911/T912 (increasing Vmax of Pma1) (Goossens et al., 2000; Eraso et al., 2006). Altogether, addition of glucose causes an increase in Pma1 activity by up to ten-fold. Trafficking of Pma1 to the plasma membrane involves Exp1 and Psg1. In the absence of Exp1, Pma1 is retained in the endoplasmic reticulum, whereas deletion of PSG1 increases Pma1 degradation in the vacuole (Geva et al., 2017). On glucose deprivation the phosphorylation of Pma1 at S899 is reversed by the Glc7 phosphatase (Mazón et al., 2015). Glucose deprivation also induces Hsp30 expression (Hahn and Thiele, 2004) to inhibit Pma1 activity (Panaretou and Piper, 1992; Piper et al., 1997). The switch from fermentation of hexoses (such as glucose or fructose) to produce ethanol, to the respiration of poor carbon sources (glycerol, ethanol, etc.) is promoted by the Snf1 kinase (AMP-activated Protein Kinase or AMPK in metazoans) (Ghillebert et al., 2011). Snf1 induces expression of genes involved in the respiration of alternative carbon sources and gluconeogenesis by activating transcription factors such as Adr1 (e.g., to express ADH2, which is used as a reporter for Snf1 activity in this study) (Ratnakumar et al., 2009) and Cat8 (e.g., to express FBP1, PCK1) (Charbon et al., 2004), and by inhibiting repressors such as Mig1 and Nrg1 [to express SUC2, GAL1 (Treitel et al., 1998; Zhou and Winston, 2001)]. In addition to regulating transcription, Snf1 also phosphorylates metabolic enzymes such as Acc1 (as a downregulation mechanism) (Woods et al., 1994), and the arrestins Rod1 and Rog3 to regulate the stability of plasma membrane carbon transporters (Llopis-Torregrosa et al., 2016). When glucose is not present, Snf1 is activated by phosphorylation at T210 (Hong et al., 2003); this site is dephosphorylated on addition of glucose (Ruiz et al., 2013). Snf1 activity is also downregulated by SUMOylation at lysine 549 in response to glucose (Simpson-Lavy and Johnston, 2013). In the absence of glucose, Snf1 associates with a gamma-protein activator (Snf4), and one of three beta-localizing proteins (Sip1, Sip2, and Gal83) (Elbing et al., 2006); this association increases Snf1 localization at the vacuolar membrane, plasma membrane and nucleus respectively. Snf1 localization to the nucleus by Gal83 is inhibited by glucose (Vincent et al., 2001). The Snf1 protein can be divided into an N-terminal domain (NTD) and a C-terminal domain (CTD, also called regulatory domain-RD). The NTD is composed of a pre-kinase domain (PKR – amino acids 1-53) and a kinase domain (KD, aa 54-391) (Figure 1). In the presence of glucose, the regulatory domain interacts with, and inhibits, the kinase domain. On phosphorylation of Snf1 at T210, Snf4 prevents the association of the kinase and regulatory domain by interacting with aa 460-498 of Snf1, preventing the auto-inhibition of Snf1 activity (Jiang and Carlson, 1996) (Amodeo et al., 2007) ; thus the region between amino acids 392-513 is referred to as RD-γ. The beta proteins (Sip1, Sip2, and Gal83) interact with aa 515-633 of Snf1 (Jiang and Carlson, 1997) (this region is thus referred to as RD-β). Snf1 must associate with Snf4 and one of the beta proteins for a stable, active complex to form (Elbing et al., 2006). Whereas the kinase domain is conserved across species, the PKR (aa1-53) shows sequence diversity. This region of the yeast protein contains an unusual polyHIS motif, comprising 13 histidines followed by a glycine and one more histidine (Celenza and Carlson, 1989) (Figure 1). In this work, we demonstrate that the polyHIS motif in the PKR acts as a progressive pH-sensing module (PSM) that controls Snf1 activity by histidine protonation/deprotonation as a consequence of, and in response to, glucose deprivation. This regulatory mechanism controls Snf1 activity independently from phosphorylation at T210, by modulating both internal interactions within Snf1 and between Snf1 and the β-subunits Gal83, Sip1, and Sip2.
Figure 1

Schematic Diagram of S. cerevisiae Snf1

Schematic diagram showing domains and regions of Snf1, with important residues and motifs. Aa1-53 is referred to as the pre-kinase region (PKR) containing a polyHIS tract at Aa18-32. Aa54-391 comprises the kinase domain (KD) – mutation of K84R prevents ATP-binding (kinase dead) and the activation loop phosphorylation site is at T210. Together the PKR and KD comprise the N-terminal domain (NTD). The regulatory domain (RD) comprises aa392-633 and comprises the entire C-terminal domain (CTD); The γ-subunit Snf4 binds to the auto-inhibitory sequence (AIS) (RD-γ) to prevent it from inhibiting the kinase domain – mutation of L470 inactivates the auto-inhibitory sequence. The β-subunits Sip1, Sip2 or Gal83 bind the far C-terminus (RD-β). K549 is SUMOylated in the presence of glucose and the SUMO interacts with the kinase domain.

Schematic Diagram of S. cerevisiae Snf1 Schematic diagram showing domains and regions of Snf1, with important residues and motifs. Aa1-53 is referred to as the pre-kinase region (PKR) containing a polyHIS tract at Aa18-32. Aa54-391 comprises the kinase domain (KD) – mutation of K84R prevents ATP-binding (kinase dead) and the activation loop phosphorylation site is at T210. Together the PKR and KD comprise the N-terminal domain (NTD). The regulatory domain (RD) comprises aa392-633 and comprises the entire C-terminal domain (CTD); The γ-subunit Snf4 binds to the auto-inhibitory sequence (AIS) (RD-γ) to prevent it from inhibiting the kinase domain – mutation of L470 inactivates the auto-inhibitory sequence. The β-subunits Sip1, Sip2 or Gal83 bind the far C-terminus (RD-β). K549 is SUMOylated in the presence of glucose and the SUMO interacts with the kinase domain.

Results

The polyHIS tract of Snf1 regulates SNF1 activity

Whereas the kinase domain of Snf1 (aa54-391) shows remarkable conservation across fungi and indeed eukaryotes, the PKR (aa1-53) is heterogeneous (Figure S1). This region of the Snf1 of many Saccharomyces, Candida, and Kluyveromyces species, but not fission yeasts such as Schizosaccharomyces pombe or vertebrates, contains a polyhistidine tract of up to 13 consecutive residues which is highly conserved in the Saccharomyces species. This is the maximal number of consecutive histidines in any S. cerevisiae protein, and indeed only 12 other proteins of the S. cerevisiae proteome contain more than 6 His (and thus all would be expected to be potential contaminants of a nickel pulldown). An early report into Snf1 structure and function declared that deletion of this tract had no negative effect on invertase activity (Celenza and Carlson, 1989), although no data were shown. Amino acids 2-52 (the PKR) is required for Snf1 function (Figure 2A). We deleted the 13 histidines, the next glycine and a fourteenth histidine (snf1) and found that although there is no effect on Snf1 function, there is an in ADH2 expression (Figure 2A). Lest this phenotype be because of the absence of amino acids (i.e., joining amino acid 17 to 33) we mutated all of the fourteen histidines to either alanine (A), lysine (K), arginine (R), phenylalanine (F), tryptophan (W) or tyrosine (Y). The rationale behind these changes is as follows: The imidazole ring of histidine is deprotonated at physiological pH (pKa = 6.5) but is protonated at lower pH values (acquiring a positive charge). Deprotonated histidine, as well as phenylalanine, tryptophan and tyrosine also engage in π-bond stacking (Liao et al., 2013). In contrast, the aromatic ring of phenylalanine cannot be protonated and thus at low pH lacks a positive charge. Tyrosine introduces a negative polarity (by the phenol group), whereas the indole group of tryptophan is aromatic, contains nitrogen but is neutral and bulkier than imidazole. Arginine and lysine were used to create a strongly positively charged motif, similar to protonated polyHIS. Substitution of the poly-histidine to poly-alanine results in the same phenotype as a deletion (elevated activity), whereas phenylalanine and tryptophan substitutions eliminate, and tyrosine substitution lowers, Snf1 function (Figure 2A). Surprisingly, poly-lysine and poly-arginine mutations also eliminated Snf1 activity. Because Snf1 regulates the transcription of many genes, we examined the expression of ACS1 (Figure S2A) and JEN1 (Figure S2B) and obtained similar results to ADH2 (Figure 2A). We then asked whether the increased Snf1 function of the snf1mutation renders Snf1 activity independent from phosphorylation at T210. A double snf1ΔH, T210A mutant, however, does not permit expression of ADH2 (Figure 2A); therefore, even in the absence of the polyHIS tract, phosphorylation at T210 is required for Snf1 function. Because the tyrosine mutant still retains about 30% Snf1 activity, we wondered whether hyperphosphorylation of T210 by deletion of Reg1 (Sanz et al., 2000) could suppress this phenotype. Indeed, deletion of reg1 enhances Snf1 activity to the same extent whether wild-type, alanine (A) or tyrosine (Y) mutant (Figure 2B).
Figure 2

The polyHIS tract of Snf1 regulates ADH2 expression

(A and B) ADH2 expression in snf1Δ cells (and snf1Δreg1Δ cells (B)) expressing prSNF1::Snf1-GFP from plasmids with the indicated mutations [deletion of polyHIS (ΔH), substitution to polyA (A), polyF (F), polyW (W), polyY (Y), T210A or empty vector]. N = 3. T-test was performed against WT cells (far left, blue). NS = not significant, ∗ = p ≤ 0.05, ∗∗ = p ≤ 0.01, ∗∗∗ = p ≤ 0.001.

(C) Localization of Snf1-GFP was determined as described in the Star Methods section.. Quantification of these data is provided in Figure S1.

(D) Growth of snf1Δ cells with different mutations of the polyHIS tract in the PKR on different carbon sources. The snf1Δ cells do not grow on galactose or glycerol medium, and snf1Δmig1Δ cells grow on galactose but not glycerol medium.

(E) Rod1-6HA and Snf1 T210 phosphorylation following transfer from 4% glucose medium to 2% ethanol medium for 15 min. Cells were grown and processed as described in materials and methods. Molecular weight markers are indicated to the right of the image in kDa.

(F) Rod1-6HA and Snf1 T210 phosphorylation following addition of glucose to 2% in ethanol media grown cells. Cells were grown and processed as described in materials and methods. Molecular weight markers are indicated to the right of the image in kDa.

(G) Phosphorylation of Rod1-6HA and Snf1 T210 following transfer from 4% glucose to 2% ethanol media. Cells were grown and processed as described in materials and methods. The red dotted line indicates maximal phosphorylation of Rod1. The two images (for WT/ΔH and F/W cells) were aligned using the protein ladder (visible on the left of each image). Molecular weight markers are indicated to the right of the image in kDa.

The polyHIS tract of Snf1 regulates ADH2 expression (A and B) ADH2 expression in snf1Δ cells (and snf1Δreg1Δ cells (B)) expressing prSNF1::Snf1-GFP from plasmids with the indicated mutations [deletion of polyHIS (ΔH), substitution to polyA (A), polyF (F), polyW (W), polyY (Y), T210A or empty vector]. N = 3. T-test was performed against WT cells (far left, blue). NS = not significant, ∗ = p ≤ 0.05, ∗∗ = p ≤ 0.01, ∗∗∗ = p ≤ 0.001. (C) Localization of Snf1-GFP was determined as described in the Star Methods section.. Quantification of these data is provided in Figure S1. (D) Growth of snf1Δ cells with different mutations of the polyHIS tract in the PKR on different carbon sources. The snf1Δ cells do not grow on galactose or glycerol medium, and snf1Δmig1Δ cells grow on galactose but not glycerol medium. (E) Rod1-6HA and Snf1 T210 phosphorylation following transfer from 4% glucose medium to 2% ethanol medium for 15 min. Cells were grown and processed as described in materials and methods. Molecular weight markers are indicated to the right of the image in kDa. (F) Rod1-6HA and Snf1 T210 phosphorylation following addition of glucose to 2% in ethanol media grown cells. Cells were grown and processed as described in materials and methods. Molecular weight markers are indicated to the right of the image in kDa. (G) Phosphorylation of Rod1-6HA and Snf1 T210 following transfer from 4% glucose to 2% ethanol media. Cells were grown and processed as described in materials and methods. The red dotted line indicates maximal phosphorylation of Rod1. The two images (for WT/ΔH and F/W cells) were aligned using the protein ladder (visible on the left of each image). Molecular weight markers are indicated to the right of the image in kDa. In the presence of glucose, Snf1 is cytoplasmic. Snf1 localization to the nucleus, vacuole, and cytoplasm in the absence of glucose is regulated by the beta-subunits Gal83, Sip1 and Sip2 respectively (Vincent et al., 2001). We found the distribution of Snf1ΔH and Snf1A to be similar to that of the wild-type both in the presence and absence of glucose (Figures 2C and S2C). However, Snf1R Snf1F, Snf1W and Snf1T210A are excluded from the nucleus in the absence of glucose, with Snf1Y exhibiting an intermediate phenotype. This suggests that, as phosphorylation of T210, the N-terminal PKR affects the interaction of the far-C-terminal of Snf1 with Gal83. Snf1Y does not exhibit a growth defect on glycerol or galactose media, whereas Snf1W does not grow on glycerol yet only has a slight growth defect on galactose medium (Figure 2D). The mutations of the polyHIS tract did not prevent phosphorylation at T210 following glucose withdrawal (Figure 2E); even when Snf1 protein levels were decreased compared to wild-type, T210 phosphorylation was normal. We noted that Snf1R appears to undergo N-terminal truncation (the GFP is at the C-terminus) but is phosphorylated in ethanol media at T210 (Figure 2E) – this results in a phenotype similar to Snf1Δ2-52. The mutations of the polyHIS tract did not prevent T210 phosphorylation following glucose withdrawal (Figure 2E). Rod1 is an arrestin that targets substrates such as Jen1 (Fujita et al., 2018) for ubiquitylation by Rsp5. Rod1 is phosphorylated and inhibited by Snf1 (Becuwe et al., 2012). Whereas Rod1 is maximally phosphorylated in Snf1WT, Snf1ΔH and Snf1A cells, Rod1 is incompletely phosphorylated in cells with the aromatic mutations or the Snf1Rmutation (Figure 2E), further showing that Snf1 activity is compromised by these mutations. We next performed a time-course for Snf1WT and Snf1ΔH inactivation. Close time points were used because we suspected that any differences might only be apparent during the first 15 min following glucose addition. On addition of glucose to ethanol grown cells, Snf1ΔH is inactivated faster than Snf1WT (as seen both by phosphorylation at T210 and Rod1 phosphorylation) (Figure 2F) which may be because of this mutant having less control mechanisms regulating its activity. On removal of glucose from media (and subsequent growth in 2% ethanol), Rod1 in both WT and Snf1ΔH cells was already maximally phosphorylated just after the washes, whereas Rod1 was not maximally phosphorylated in Snf1F or Snf1W cells even after 30 min. However, phosphorylation of Snf1 at T210 was not affected by these mutations (Figure 2G).

The polyHIS tract regulates internal and external interactions of Snf1

We examined the interactions between Snf1 and the beta and gamma subunits of SNF1 by yeast two-hybrid. For comparison, and as verification of the constructs, we included the published interaction between the kinase domain (1-391) or the regulatory domain (392-633) of Snf1 (Jiang and Carlson, 1996) and the beta or gamma proteins. Snf4 interacts with the RDγ (aa460-498) region of Snf1 (Jiang and Carlson, 1996; Amodeo et al., 2007), far away from the histidine tract (Figure 1). The interaction between Snf1 and Snf4 is unaffected by changes in the polyHIS tract (Figure 3A). In contrast, deletion (ΔH) or substitution to alanine (A) of the polyHIS tract enhances the interaction between Snf1 and Gal83 (Figure 3B) which interacts with RDβ at aa515-633 (Jiang and Carlson, 1997). Substitution to tyrosine (Y) lowers the interaction, and conversion of the histidines to phenylalanine (F), tryptophan (W) or arginine (R), abolishes the Snf1/Gal83 interaction altogether, similarly to the T210A mutation. Interactions of Snf1 with the other beta subunits, Sip1 (Figure S3A) or Sip2 (Figure S3B) were not enhanced by deletion of the polyHIS tract (ΔH) or its mutation to alanine (A) but interaction of Snf1 with Sip1 or Sip2 was lowered to levels seen in the T210A mutant when the polyHIS was substituted with arginine, phenylalanine, tryptophan or tyrosine. Because histidine is more deprotonated in glucose medium (because of high Pma1 pumping of protons into the extracellular medium), it seems that the neutral aromatic amino acids phenylalanine and tryptophan (and to lesser extent tyrosine) mimic the deprotonated, neutral form of histidine.
Figure 3

The polyHIS tract regulates Snf1 interactions

(A) Yeast two-hybrid experiment showing interaction of Snf4 with indicated Snf1 polyHIS mutants in glucose and ethanol media. N = 3.

(B) As A. but with pGAD-Gal83 in lieu of pGAD-Snf4. N = 3.

(C) As A. but with pGAD-Snf1392−633 (CTD) and pGBD-Snf11−391 (NTD) with indicated mutations [deletion of polyHIS (ΔH), substitution to polyA (A), polyF (F), polyW (W), polyY (Y), T210A]. N = 3. T-tests were conducted against WT ethanol (far left, orange) for A and B, and against WT glucose (far left blue) for C. NS = not significant, ∗ = p ≤ 0.05, ∗∗ = p ≤ 0.01, ∗∗∗ = p ≤ 0.001.

(D and E) As (C) but with indicated pGBD and pGAD Snf1 fragments. PKR = amino acids 1-53, KD = 54-391, NTD = 1-391, RD/CTD = 392-633, RD-α = 392-513, RD-β = 514-633. N = 3. T-tests were conducted against NTD/CTD glucose (far right blue brick) for D and E. NS = not significant, ∗ = p ≤ 0.05, ∗∗ = p ≤ 0.01, ∗∗∗ = p ≤ 0.001.

The polyHIS tract regulates Snf1 interactions (A) Yeast two-hybrid experiment showing interaction of Snf4 with indicated Snf1 polyHIS mutants in glucose and ethanol media. N = 3. (B) As A. but with pGAD-Gal83 in lieu of pGAD-Snf4. N = 3. (C) As A. but with pGAD-Snf1392−633 (CTD) and pGBD-Snf11−391 (NTD) with indicated mutations [deletion of polyHIS (ΔH), substitution to polyA (A), polyF (F), polyW (W), polyY (Y), T210A]. N = 3. T-tests were conducted against WT ethanol (far left, orange) for A and B, and against WT glucose (far left blue) for C. NS = not significant, ∗ = p ≤ 0.05, ∗∗ = p ≤ 0.01, ∗∗∗ = p ≤ 0.001. (D and E) As (C) but with indicated pGBD and pGAD Snf1 fragments. PKR = amino acids 1-53, KD = 54-391, NTD = 1-391, RD/CTD = 392-633, RD-α = 392-513, RD-β = 514-633. N = 3. T-tests were conducted against NTD/CTD glucose (far right blue brick) for D and E. NS = not significant, ∗ = p ≤ 0.05, ∗∗ = p ≤ 0.01, ∗∗∗ = p ≤ 0.001. How can mutations at the far N-terminus of Snf1 affect interactions in the C-terminal domain? One possibility is that that the polyHIS tract is involved in the intramolecular interaction of the NTD and CTD (Jiang and Carlson, 1996). Indeed, deletion (ΔH) or substitution to alanine (A) decreases this inhibitory interaction even in glucose and severely in ethanol medium, whereas interaction between the NTD and CTD remains high in ethanol when the polyHIS tract is substituted with phenylalanine, tyrosine or tryptophan. Substitution to arginine (R) results in the lowest interaction between the NTD and CTD, possibly because of the strong positive charge in the PKR of the NTD (Figure 3C). Because SUMOylation at K549 within RD-β has been found to be involved in the inhibition of Snf1 by interaction with SUMO-Interacting-Motifs at I129 in the kinase domain (Simpson-Lavy and Johnston, 2013), we examined the effects of a K549R mutation in the regulatory domain on the interaction with the kinase domain. SUMOylation at K549 is a prerequisite for interaction between the kinase and regulatory domains in glucose media, even when the polyHIS is substituted to aromatic amino acid. (Figure S3C). Therefore, regulation of Snf1 by SUMOylation is independent from regulation of Snf1 through the polyHIS motif in the PKR. Together, these data suggest that the polyHIS tract is a regulator of the kinase/regulatory domain interaction at the β-subunit binding site, thereby controlling accessibility of β-subunits and thus Snf1 localization and activity. We next examined interactions of the PKR (1-53) with the kinase domain (KD) (54-391) and the regulatory domain (CTD) (392-633) of Snf1. In glucose medium, the PKR interacts with the regulatory domain. This interaction is mediated via the C-terminal half of the regulatory domain (aa514-633) that contains the binding sites of the β-subunits (RD-β); and is absent in ethanol-based medium. Whereas interaction of the NTD with the CTD is SUMO dependent, interaction of the PKR with either the whole CTD or RDβ is SUMO-independent (i.e., not affected by K549R) (Figure 3D). The PKR contacts the kinase domain in the absence of glucose (Figure 3D). Because this suggests that the PKR interacts with different areas of Snf1 depending on the glucose status of the environment, we determined the requirements for this interaction. This interaction requires histidine or arginine: deletion or substitution to alanine or the aromatic amino acids all prevent this interaction from occurring, whereas substitution to arginine results in some interaction between the PKR and KD (Figure S3D). This may be because of the diffuse positive charge of the protonated imidazole ring. In addition to the polyHIS tract, this interaction also requires Snf1 kinase activity because the K84R mutation abolishes the interaction, and phosphorylation at T210. Indeed, this interaction is elevated (even in glucose-grown cells) when T210 is replaced with a phospho-mimicking glutamate (T210E, Figure 3E). These results suggest that one of the effects of phosphorylation at T210 may be to attract the PKR via the polyHIS tract, thereby unlocking the beta subunit binding site.

The polyHIS tract of Snf1 functions as a pH sensing module (PSM)

Because the interactions between the kinase and the regulatory domains mediated by the polyHIS tract are mimicked by substitution with neutral aromatic amino acids, and are high in glucose and low in ethanol, we considered how this could be regulated. One possibility is that respiration-derived reactive oxygen species (ROS) could oxidize histidine to produce 2-oxo-histidine as occurs in the PerR peroxide sensor in B. subtilis (Traoré et al., 2009) triggering its proteolysis (Ahn and Baker, 2016). Incubation of cells with the ROS scavenger N-acetyl cysteine does not affect Snf1 activity. Although treatment of cells with diamide does increase Snf1H activity (Figure S4A), diamide treatment also further elevates Snf1ΔH activity, and therefore the oxidative stress induced increase in Snf1 activity is polyHIS-independent. Cytoplasmic alkalization as a consequence of Pma1 pumping protons into the external media in response to glucose functions as a second messenger signal (Dechant et al., 2010). We therefore investigated whether mutations or drugs that affect Pma1 activity regulate Snf1 by modulating the charge of the histidines of the polyHIS tract. We increased Pma1 activity (glucose mimicking, resulting in cytoplasmic alkalization) by two methods simultaneously. Truncation of the final 18 amino acids of Pma1 (PMA1-Δ901) increases Pma1 activity to that observed in the presence of glucose (Mason et al., 1998). The Hsp30 chaperone is expressed on glucose exhaustion (via Snf1 activation of the Hsf1 transcription factor (Hahn and Thiele, 2004) and inhibits Pma1 activity (Panaretou and Piper, 1992; Piper et al., 1997). Deletion of HSP30 thus maintains Pma1 activity even in the absence of glucose. Either PMA1-Δ901 truncation or deletion of HSP30 results in a 20–30% decrease in ADH2 expression (Figure 4A). In combination, a dramatic decrease in Snf1 function is observed, as both pathways of Pma1 downregulation are now disabled. These decreases in ADH2 expression are suppressed by Snf1ΔH, indicating that Pma1 upregulation is signaling via the polyHIS tract (Figure 4A). Similar results were obtained for ACS1 (Figure S4B) and JEN1 (Figure S4C) expression. Upregulating Pma1 activity did not diminish phosphorylation of Snf1 at T210 (Figure 4B). However, on shift to ethanol medium Rod1 is degraded in hsp30Δ PMA1-Δ901 cells (rather than being phosphorylated) – this phenotype is suppressed and normal maximal phosphorylation restored by introduction of Snf1ΔH (Figure 4B), suggesting that Snf1ΔH is not affected by high Pma1 activity.
Figure 4

The polyHIS tract of Snf1 is a pH sensing module

(A) ADH2 expression of snf1Δ cells and snf1 hsp30 cells with prSNF1::Snf1-GFP or prSNF1::Snf1ΔH-GFP, and with PMA1 or PMA1-Δ901. N = 3. T-tests were conducted against WT (far left, blue). NS = not significant, ∗ = p ≤ 0.05, ∗∗ = p ≤ 0.01, ∗∗∗ = p ≤ 0.001.

(B) Rod1-6HA and Snf1 T210 phosphorylation following transfer of cells from 4% glucose to 2% ethanol media for 15 min. Cells were grown and processed as described in materials and methods. Molecular weight markers are indicated to the right of the image in kDa.

(C) Yeast two-hybrid experiment showing interaction of GBD-Snf1 and GAD-Snf1 in glucose and ethanol media N = 3. T-test is against WT cells in glucose (far left, blue). NS = not significant, ∗ = p ≤ 0.05, ∗∗ = p ≤ 0.01, ∗∗∗ = p ≤ 0.001.

(D) Rod1-6HA phosphorylation following addition of glucose to 2% to ethanol grown cells with indicated concentrations of ebselen. Cells were grown and processed as described in materials and methods.

(E) Rod1-6HA dephosphorylation following addition of glucose to 2% to ethanol grown cells. Cells were grown and processed as described in materials and methods. Molecular weight markers are indicated to the right of the image in kDa.

(F) Yeast two-hybrid experiment to determine interaction between GBD-Snf1 and GAD-Snf1 in the presence of the indicated concentration of ebselen. N = 3. T-tests were performed against WT glucose without ebselen (far left, blue). NS = not significant, ∗ = p ≤ 0.05, ∗∗ = p ≤ 0.01, ∗∗∗ = p ≤ 0.001.

The polyHIS tract of Snf1 is a pH sensing module (A) ADH2 expression of snf1Δ cells and snf1 hsp30 cells with prSNF1::Snf1-GFP or prSNF1::Snf1ΔH-GFP, and with PMA1 or PMA1-Δ901. N = 3. T-tests were conducted against WT (far left, blue). NS = not significant, ∗ = p ≤ 0.05, ∗∗ = p ≤ 0.01, ∗∗∗ = p ≤ 0.001. (B) Rod1-6HA and Snf1 T210 phosphorylation following transfer of cells from 4% glucose to 2% ethanol media for 15 min. Cells were grown and processed as described in materials and methods. Molecular weight markers are indicated to the right of the image in kDa. (C) Yeast two-hybrid experiment showing interaction of GBD-Snf1 and GAD-Snf1 in glucose and ethanol media N = 3. T-test is against WT cells in glucose (far left, blue). NS = not significant, ∗ = p ≤ 0.05, ∗∗ = p ≤ 0.01, ∗∗∗ = p ≤ 0.001. (D) Rod1-6HA phosphorylation following addition of glucose to 2% to ethanol grown cells with indicated concentrations of ebselen. Cells were grown and processed as described in materials and methods. (E) Rod1-6HA dephosphorylation following addition of glucose to 2% to ethanol grown cells. Cells were grown and processed as described in materials and methods. Molecular weight markers are indicated to the right of the image in kDa. (F) Yeast two-hybrid experiment to determine interaction between GBD-Snf1 and GAD-Snf1 in the presence of the indicated concentration of ebselen. N = 3. T-tests were performed against WT glucose without ebselen (far left, blue). NS = not significant, ∗ = p ≤ 0.05, ∗∗ = p ≤ 0.01, ∗∗∗ = p ≤ 0.001. Upregulating Pma1 activity by deletion of HSP30 and/or by PMA1-Δ901 increases the interaction between the NTD (1-391) and CTD (392-633) domains of Snf1. (Figure 4C). However, Pma1 upregulation does not change the low interaction between the NTD and CTD when the polyHIS tract in the NTD is deleted (ΔH), nor the high interaction observed by substitution to tryptophan (W) (Figure 4C), demonstrating that Pma1, through pH regulation, governs the intra-molecular interactions of Snf1. We also inhibited Pma1 using the organo-selenium compound ebselen (Chan et al., 2007; Young et al., 2010). A low dose treatment (10 μM) does not affect Rod1 dephosphorylation on addition of 2% glucose to ethanol grown cells. However, a high dose (100 μM) of ebselen prevents dephosphorylation of Rod1 in the presence of glucose (Figure 4D). This is suppressed by Snf1ΔH (Figure 4D) showing that Snf1 is not regulated by Pma1 activity. Genetically lowering Pma1 activity by deletion of EXP1 (resulting in cytoplasmic acidification) results in a marked deceleration of Rod1 dephosphorylation on addition of glucose to ethanol grown cells (Figure 4E). Ebselen was added to yeast-two hybrid cells on transition to glucose and the interaction between the kinase and regulatory domains of Snf1 monitored. Pma1 inhibition prevents association of Snf1NTD and Snf1CTD. However, Pma1 inhibition does not prevent association of Snf1CTD with Snf1NTD−W, in which the polyHIS tract is substituted with neutral deprotonated tryptophan (Figure 4F). Together, these results demonstrate that the polyhistdine tract in the PKR functions as a PSM. Glucose thus inhibits Snf1 via upregulation of Pma1 activity, causing deprotonation of the polyHIS tract, resulting in its association with the β-subunit binding site.

The polyHIS tract elicits a progressive response

Although the kinase domain of AMPK/Snf1 is highly conserved from yeasts to man, the PKR shows immense variation both in length (53 amino acids in S. cerevisiae, 25 amino acids in Homo sapiens, 74 amino acids in Candida parapsilosis) and composition, with metazoans typically having short unenriched sequences, whereas yeasts and fungi often possess a PKR with enrichment or tracts of histidine (H), asparagine (N) and glutamine (Q) (Figure S1). To investigate the effects of sequence heterogeneity on Snf1 activity, we replaced the first 52 acids of S. cerevisiae Snf1 with the equivalent PKR from Kluyveromyces, Saccharomyces and Candida yeasts, all of which contain either histidine enrichment or polyHIS tracts of various lengths. We also created a S. cerevisiae PKR with only 4 or 8 polyhistidines. Removal of the PKR or replacement with the H. sapiens PKR eliminates Snf1 activity (Figure 5A). Most heterologous PKR chimeras did not drastically affect Snf1 activity, but the PKR from Kluyveromyces lactis Fog2 increased Snf1 activity to levels near that of Snf1ΔH. The K. lactis PKR contains only two consecutive histidines (with three more scattered); this could explain why it is ineffective in regulating S. cerevisiae Snf1. The S. cerevisiae 4-histidine construct has a Snf1 activity between that of the wild-type Snf1 (14 histidines) and Snf1ΔH, and the 8-histidine construct has a Snf1 activity between WT (14 histidines) and the 4-histidine mutant. We therefore plotted Snf1 activity as a function of the number of histidines in the PKR and found a strong correlation (Figure 5B). This demonstrates that Snf1 activity progressively decreases on gain of even two consecutive histidines (such as in the PKR of K. lactis Fog2).
Figure 5

The polyHIS tract elicits a progressive response

(A) ADH2 expression of snf1Δ cells transformed with Snf1-GFP plasmids bearing the indicated PKR (see the legend of Figure S1 and Table 2 for details). N = 3. Results which are of note are in turquoise. T-tests were performed against WT glucose (far left, blue). NS = not significant, ∗ = p ≤ 0.05, ∗∗ = p ≤ 0.01, ∗∗∗ = p ≤ 0.001.

(B) Number of consecutive histidines in the PKR was plotted against ADH2 expression using data from (A). R2 value was calculated by Excel.

The polyHIS tract elicits a progressive response (A) ADH2 expression of snf1Δ cells transformed with Snf1-GFP plasmids bearing the indicated PKR (see the legend of Figure S1 and Table 2 for details). N = 3. Results which are of note are in turquoise. T-tests were performed against WT glucose (far left, blue). NS = not significant, ∗ = p ≤ 0.05, ∗∗ = p ≤ 0.01, ∗∗∗ = p ≤ 0.001.
Table 2

Plasmids

NumberGenotypeBackboneSource/Notes
152Snf1-GFPPrs313M. Carlson (Vincent et al., 2001)
728Snf1-GFPPrs315Marker switch of 152
847Snf1-ΔH-GFPPrs315polyHIS of 728 deleted
1244Snf1-4H-GFPPrs3154 polyhistdines in PKR (by PCR from 847)
1248Snf1-8H-GFPPr3158 polyhistdines in PKR (by PCR from 847)
742Snf1-A-GFPPrs315Alanine substituted for histidine in PKR of 728
852Snf1-F-GFPPrs315Phenylalanine substituted for histidine in PKR of 728
744Snf1-K-GFPPrs315Lysine substituted for histidine in PKR of 728
1257Snf1-R-GFPPrs315Arginine substituted for histidine in PKR of 728
745Snf1-W-GFPPrs315Trptophan substituted for histidine in PKR of 728
731Snf1-Y-GFPPrs315Tyrosine substituted for histidine in PKR of 728
982Snf1-Δ2-52-GFPPrs315Pre-kinase region deleted
846Snf1-T210A-GFPPrs315T210A mutation in 728
1 or 919prADH2::LacZPrs316E. Young (Ratnakumar et al., 2009)
63prADH2::LacZPrs313Marker switch of 1
164prADH2::LacZPrs315Marker switch of 1
240prACS1::LacZ2μ URA3H.J. Schüller (Kratzer and Schüller, 1995)
706prJEN1::LacZYep358B. Guiard (Lodi et al., 2002)
526Nup49-CherryPrs314M. Lisby (Germann et al., 2014)
589Nup49-CherryPrs316Marker switch of 526
590Nup49-CherryPrs315Marker switch of 526
816Snf1-3HAPrs313M. Carlson (Liu et al., 2011)
817Snf1- K84R-3HAPrs313M. Carlson (Liu et al., 2011)
818Snf1- T210A 3HAPrs313M. Carlson (Liu et al., 2011)
819Snf1-ΔH-3HAPrs313polyHIS of 816 deleted
820Snf1-A-3HAPrs313Alanine substituted for histidine in 816
835Snf1-ΔH-T210A-3HAPrs313polyHIS of 818 deleted
1234Pma1Prs314Pma1 inserted into prs314
1238Pma1-Δ901Prs314Final 18aa of Pma1 deleted

Yeast 2-hybrid plasmids

153prADH1::GADLEU2, 2μpACT2 Yeast 2-hybrid empty vector
154prADH1::GBDTRP1, 2μpGBT9 Yeast 2-hybrid empty vector
155prADH1::GBDURA3, 2μpGBU9 Yeast 2-hybrid empty vector
753GAD-Snf4LEU2, 2μSnf4 in 153
755GAD-Sip1LEU2, 2μSip1 in 153
756GAD-Sip2LEU2, 2μSip2 in 153
757GAD-Gal83LEU2, 2μGal83 in 153
762GBD-Snf1URA3, 2μSnf1 in 155
763GBD-Snf1 K84RURA3, 2μSnf1 K84R in 155
764GBD-Snf1 T210AURA3, 2μSnf1 T210A in 155
789GBD-Snf1 ΔHURA3, 2μpolyHIS deletion of 762
779GBD-Snf1 AURA3, 2μAlanine substituted for polyHIS in 762
801GBD-Snf1 FURA3, 2μPhenylalanine substituted for polyHIS in 762
787GBD-Snf1 RURA3, 2μArginine substituted for polyHIS in 762
765GBD-Snf1 WURA3, 2μTryptophan substituted for polyHIS in 762
768GBD-Snf1 YURA3, 2μTyrosine substituted for polyHIS in 762
781GBD-Snf1 1-391URA3, 2μSnf1 1-391 in 155
786GBD-Snf1 392-633URA3, 2μSnf1 392-633 in 155
810GBD-Snf1 392-633 K549RURA3, 2μSnf1 392-633 in 155
791GAD-Snf1 392-633LEU2, 2μSnf1 392-633 in 153
822GAD-Snf1 392-633 K549RLEU2, 2μSnf1 392-633 K549R in 153
795GBD-Snf1 1-391 ΔHURA3, 2μSnf1 1-391 ΔH in 155
807GBD-Snf1 1-391 AURA3, 2μSnf1 1-391 A in 155
805GBD-Snf1 1-391 FURA3, 2μSnf1 1-391 F in 155
798GBD-Snf1 1-391 RURA3, 2μSnf1 1-391 R in 155
785GBD-Snf1 1-391 WURA3, 2μSnf1 1-391 W in 155
783GBD-Snf1 1-391 YURA3, 2μSnf1 1-391 Y in 155
828GBD-Snf1 1-53URA3, 2μSnf1 1-53 in 155
836GBD-Snf1 1-53 ΔHURA3, 2μSnf1 1-53 ΔH in 155
829GBD-Snf1 1-53 AURA3, 2μSnf1 1-53 A in 155
830GBD-Snf1 1-53 FURA3, 2μSnf1 1-53 F in 155
831GBD-Snf1 1-53 RURA3, 2μSnf1 1-53 R in 155
837GBD-Snf1 1-53 WURA3, 2μSnf1 1-53 W in 155
836GBD-Snf1 1-53 YURA3, 2μSnf1 1-53 Y in 155
844GAD-Snf1 54-391LEU2, 2μSnf1 54-391 in 153
853GAD-Snf1 54-391 K84RLEU2, 2μSnf1 54-391 K84R in 153
860GAD-Snf1 54-391 T210ALEU2, 2μSnf1 54-391 T210A in 153
954GAD-Snf1 54-391 T210ELEU2, 2μMutagenesis of 844
827GAD-Snf1 392-513LEU2, 2μSnf1 392-513 in 153
826GAD-Snf1 514-633LEU2, 2μSnf1 514-633 in 153
825GAD-Snf1 514-633 K549RLEU2, 2μSnf1 514-633 K549R in 153

Pre-kinase regions from other species. Amino acids 1–52 removed from plasmid 728 prScerSNF1::ScerSNF1-GFP and replaced by:

893CtroPrs3151-52 from C. tropicalis
894KlacPrs3151-32 from K. lactis Fog2
895KwalPrs3151-28 from K. waltii
896LeloPrs3151-71 from L. elongisporans
898CparPrs3151-73 from C. parapsilosis
904CglaPrs3151-34 from C. glabrata
945HsapPrs3151-24 from H. sapiens
981DhanPrs3151-52 from D. hanseni
985SmikPrs3151-40 from S. mikatae
986SkluPrs3151-29 from S. kluyveri
(B) Number of consecutive histidines in the PKR was plotted against ADH2 expression using data from (A). R2 value was calculated by Excel.

Discussion

The importance of pH as a second messenger for carbon metabolism

Protons are an abundant diffusible signaling molecule that act as a second messenger for highly fermentable carbon sources, both via the plasma membrane ATPase (Pma1) and the vacuolar ATPase (Vma1), controlling many processes, including membrane biogenesis (Young et al., 2010) and metabolic enzyme aggregation (Petrovska et al., 2014). It is unclear how glucose activates these proton pumps. Pma1 activity decreases extracellular pH, thus activating Gpa2 (which functions as a pH sensor as well as a glucose sensor). Only the fermentable, non-Snf1 activating sugars glucose, fructose and mannose cause this effect, and other hexoses or disaccharides do not affect pH or PKA activity (Isom et al., 2018). Increased glycolytic flux beyond glucose-6-phosphate increases Vma1 activity leading to activation of PKA (Dechant et al., 2010). Recently, valproate has been shown to decrease glycolytic flux resulting in both increased phosphorylation at T210, and decreased Pma1 activity (Salsaa et al., 2021). Here we have shown that decreasing Pma1 activity is a requirement for Snf1 activation. Because PKA and Snf1 reciprocally regulate each other through multiple mechanisms (PKA regulates Snf1 via inhibitory SUMOylation (Simpson-Lavy et al., 2015), inhibiting Sak1 activity (Barrett et al., 2012), and by controlling Sip1 localization (Hedbacker et al., 2004), whereas Snf1 inhibits adenyl cyclase (Cyr1) (Nicastro et al., 2015) we can now state that PKA and Snf1 are also oppositely regulated by pH, ensuring co-ordination of metabolic signaling across the cell. There are two mechanisms by which Pma1 mediation proton export (and thus cytoplasmic alkalization) is attenuated by glucose deprivation – dephosphorylation of the C-terminal tail (Mason et al., 1998) and Pma1 protein destabilization by Hsp30 (Panaretou and Piper, 1992; Piper et al., 1997). We found that both Pma1 inactivation pathways need to be inhibited (Pma1-Δ901 and hsp30Δ) to obtain the maximal inhibition of Snf1 in ethanol media (Figure 4A). This may suggest that the reduction of Pma1 activity on glucose deprivation may be achieved by either mechanism.

A regulation module that senses pH

The PKR (aa 2-53) is essential for Snf1 activity (Figures 1 and 2A) despite being heterogenous between species, whereas the kinase domain (aa 54-391) is considerably conserved. We found that substitution with the human PKR (aa 2-26 of AMPKα) renders this fusion inactive. This is probably because of the distal part of the PKR and AMPKβ (Snf4 in S. cerevisiae) forming the (allosteric drug and metabolite) ADaM pocket (Calabrese et al., 2014; Jørgensen et al., 2021). This pocket centered around G53 may also explain why the hyperactive ScSnf1G53Rmutant suppresses snf4 (Estruch et al., 1992). It is probable that the human PKR is incompatible with S. cerevisiae Snf4. Indeed, HsAMPKα does not complement snf1Δ in S. cerevisiae – although replacement of all of the Snf1 complex in S. cerevisiae with human AMPKα,β,γ can permit growth on non-fermentable carbon sources (Ye et al., 2014). In this study we have identified a neglected motif in Snf1, the polyHIS tract of the PKR, as a Pma1 activity/pH sensing module (PSM) that regulates Snf1 activity, thus directly linking glucose activation of Pma1 (and thus intracellular alkalization) and inhibition of Snf1. Increasing Pma1 activity by deletion of hsp30 in combination with the PMA1-Δ901 allele lowers Snf1 activity (Figures 4A, 4B, S4B, and S4C) and decreasing Pma1 activity in glucose results in a deceleration of Snf1 inactivation (Figures 4D and 4E), with deletion of the PSM suppressing these phenotypes. The polyHIS tract acts by interacting with, and blocking, the β-subunit binding site in the far C-terminal of Snf1 (Figures 3B, 3C, and 3D); preventing Snf1 activation and entry into the nucleus (Figures 2A and 2C). By substituting other poly amino-acids for the polyHIS we found that Snf1A mimics deletion of the polyHIS tract, whereas aromatic, neutral amino acids (Snf1F and Snf1W) cause the kinase and regulatory domains of Snf1 to maintain their interaction even in the absence of glucose, mimicking deprotonated polyHIS as would be found in glucose-grown cells. Because the aromatic amino acids mimic deprotonated histidine, but mutation to alanine is congruent with a histidine-null phenotype, it is likely that π-bond stacking of the polyHIS is interacting with the β-subunit binding site in the C-terminal of Snf1. Deletion of the polyHIS tract results in Snf1 being inactivated faster in glucose (Figure 2F), probably because of a loss of this level of regulation. The PSM binds (in glucose) and is released from (in ethanol) the β-subunit binding site because of its deprotonation (in glucose) and protonation (in ethanol). This was determined by upregulating Pma1 activity in ethanol (PMA1-Δ901 and hsp30Δ) (Figure 4C) or by inhibiting Pma1 with ebselen in glucose (Figure 4F), with phenotypes being suppressed by Snf1ΔH or Snf1W appropriately. This establishes a direct, distinct signaling pathway from glucose to Pma1 activity resulting in deprotonation of Snf1 and thus Snf1 inhibition. Substitution of the PSM with arginine resulted in an N-terminal truncation of Snf1 (Figure 2E) and much reduced activity (Figures 2A, S2A, and S2B). In most ways Snf1R behaves like Snf1F, but isolated PKR interacts well with the KD in ethanol media – probably because of the high positive charge interacting with phosphorylated T210. Although how exactly glucose regulates T210 activation loop phosphorylation remains unclear, it is important to note that even without a PSM, phosphorylation at T210 is still required for Snf1 activity (Figure 2A). Moreover, the interaction of the PKR with the kinase domain in the absence of glucose also requires phosphorylation at T210 (Figure 3E). Furthermore, mutation of the polyHIS tract does not prevent phosphorylation at T210 on glucose withdrawal, giving rise to a situation whereby Snf1 can be phosphorylated at T210 (Figure 2E), but a mutation in the polyHIS tract prevents dissociation of the kinase and regulatory domains (Figure 3C) and thus Snf1 kinase activity is still inhibited (Figure 2A). Similarly, although modulation of Pma1 activity affects Snf1 target gene expression (Figures 4A, S4B, and S4C) and Rod1 phosphorylation (Figures 4B, 4D, and 4E); regulation of Snf1 phosphorylation at T210 is unaffected. In other words, our work shows that Snf1 activation requires two events: T210 phosphorylation andprotonation of the polyHIS tract. Of importance, a corollary from these results is that experiments showing T210 phosphorylation do not automatically demonstrate the presence of an active Snf1. The data allows a reconstruction of a mechanism of Snf1 regulation by Pma1 activity (Figure 6). In the absence of glucose:
Figure 6

Model of the mechanism by which pH regulates Snf1.

1. In glucose medium PMA1 is phosphorylated and pumps protons out of the cell resulting in polyHIS deprotonation (mimicked by Snf1F and Snf1W).

2. On glucose deprivation, PMA1 is dephosphorylated and much less active. The polyHIS of Snf1 becomes protonated and Snf1 is phosphorylated at T210. These events are independent.

3. The protonated polyHIS tract disengages from the β-subunit binding site, and Snf4 and a β-subunit (Gal83, Sip1, Sip2) can interact.

4. Snf1-Snf4-Gal83 enters the nucleus resulting in expression of respiration genes such as ADH2 and HSP30. Hsp30 further inhibits Pma1. The nucleus in this panel is enlarged for clarity only.

Pma1 activity is inhibited, resulting in acidification of the cytoplasm and polyHIS protonation. Snf1 is phosphorylated at T210 and de-SUMOylated (independently from polyHIS protonation). The protonated polyHIS tract dissociates from the β-subunit binding site (because of the higher positive charge) and interacts with T210 phosphorylated kinase domain. β-subunits can now interact with Snf1. Snf1 is active and enters the nucleus. Snf1 activity increases Hsp30 levels, further inhibiting Pma1. Model of the mechanism by which pH regulates Snf1. 1. In glucose medium PMA1 is phosphorylated and pumps protons out of the cell resulting in polyHIS deprotonation (mimicked by Snf1F and Snf1W). 2. On glucose deprivation, PMA1 is dephosphorylated and much less active. The polyHIS of Snf1 becomes protonated and Snf1 is phosphorylated at T210. These events are independent. 3. The protonated polyHIS tract disengages from the β-subunit binding site, and Snf4 and a β-subunit (Gal83, Sip1, Sip2) can interact. 4. Snf1-Snf4-Gal83 enters the nucleus resulting in expression of respiration genes such as ADH2 and HSP30. Hsp30 further inhibits Pma1. The nucleus in this panel is enlarged for clarity only. The three mechanisms of Snf1 activation on glucose deprivation thus identified: protonation of the polyHIS tract in the PSM, phosphorylation at T210 and de-SUMOylation of K549, are all required for Snf1 activity, and are all independent from each other. Thus Snf1 acts as a signal coordinator, similar to the integration of glucose, nitrogen and phosphate sensing by Rim15 (Swinnen et al., 2006).

PolyHIS sensing is progressive

An emergent theme is that histidines can act as pH sensors, such as RasGRP1 in H. sapiens at a single histidine H212 (Vercoulen et al., 2017) or utilizing multiple histidine residues in PsaF of Y. pestis (Quinn et al., 2021) or the human prion PrP (Zahn, 2003). It is possible that pH sensing of Snf1 arises from the aggregate charge changes of the histidines in the PSM. We determined that the number of histidines in the PKR inversely correlates with Snf1 activity (Figures S1, 5A, and 5B). Thus, it appears that the crucial determinant of pH sensing by polyHIS in Snf1 is the number of local histidines, which may implicate π-bond stacking as the molecular interaction with the β-subunit binding site in the regulatory domain.

PSMs in the proteome

The polyHIS motif of the PSM is very rare. In S. cerevisiae only 12 other proteins contain polyHIS stretches with more than 7 histidines (H7+ representing 0.25% of the S. cerevisiae proteome [assuming 5815 proteins (Engel et al., 2014)] In contrast, 91 (1.6%) of S. cerevisiae proteins contain polyQ7+. The human genome contains 95 polyH7+ proteins representing 0.47% of the proteome (assuming 20,352 protein-encoding genes (Pertea et al., 2018)) whereas polyQ7+ proteins number 709 (3.5% of the proteome). The ratio of proteins with polyHIS7+ to proteins with polyQ7+ is thus constant between yeast and humans. The proteins containing polyH7+ in yeast are diverse and do not seem to have a common function, although 4 of the 12 are involved in metabolism and stress response (SNF1, TAX4, ROM2, BAG7). However, in humans, an enrichment for voltage-gated channels is observed – this could permit a progressive change in voltage to accumulate in the histidines and when a threshold is reached to effect a conformational/interactional change. Indeed, multiple histidine residues have been found to be important in the human voltage-gated channel Hv1 (Zhao et al., 2015; Cherny et al., 2018) but analysis of histdines as a group has yet to be performed. Although Snf1 has been an intensively investigated protein in the last 40 years, the function of the rarely occurring polyHIS tract at the N-terminus has hitherto not been elucidated. We have found that this region constitutes a PSM, which enables a direct progressive regulation of Snf1 by pH changes arising from the presence of glucose. Roles for pH sensing modules, such as the one we have found for the yeast AMPK ortholog are likely to govern cellular processes in many organisms, as pH variations are ubiquitous, and affected by nutrient availability and cellular metabolism.

Yeast strains, plasmids and oligonucleotides used

The strains used in this study are listed in Table 1. Plasmids are shown in Table 2. Oligonucleotides are listed in Table 3.
Table 1

Yeast strains

NumberNameGenotypeFrom
1WTW303 MATa ade2-1 his3-11,15 trp1-1 leu2-3,112 ura3-1 rad5-535 bud4Lab collection
2snf1As 1 snf1::HYGE. Young (Ratnakumar et al., 2009)
3snf1 ADE2As 1 snf1::HYG ADE2This study
4reg1As 1MATα reg1::NATE. Young (Ratnakumar et al., 2009)
5snf1 reg1As 1 snf1:: HYG reg1::NATE. Young (Ratnakumar et al., 2009)
6snf1 mig1As 1 snf1:: HYG mig1::KANThis study
7snf1 hsp30As 1 snf1::HYG hsp30::NATThis study
8snf1 Rod1-6HAAs 1 snf1::HYG Rod1-6HA::NATThis study
9snf1 Rod1-6HA hsp30As 1 snf1::HYG Rod1-6HA::NAT hsp30::KANThis study
10snf1 Rod1-6HA exp1As 1 snf1::HYG Rod1-6HA::NAT exp1::KANThis study
11PJ694trp1-901 leu2-3,112, ura3-52, his3Δ200, gal4Δ, gal80Δ, prGAL2:ADE2 lys2::prGAL1::HIS3 met2::prGAL7-LacZLab collection
12Y2H hsp30As 8 hsp30::NAT

All S. cerevisiae strains are in the W303 background except PJ694 (for Yeast 2-Hybrid assays).

Table 3

Oligonucleotides for mutagenesis of Snf1

NameDescriptionSequence
Gap Repair Mutagenesis of polyHIS

C54Y-5CTACTATTATTACTATTATTACTACTATTACTATTATGGATATggcggaagcaactcgacg
C53Y-3cataatagtaatagtagtaataatagtaataatagtagtaGCTAGAATTTGCATTGGCAGGTG
C57A-5GCAGCTGCCGCAGCAGCCGCTGCCGCAGCAGCTGCAGCAGGAGCCggcggaagcaactcgacg
C58A-3ggctcctgctgcagctgctgcggcagcggctgctgcggcagctgcGCTAGAATTTGCATTGGCAGGTG
C74F-5CTTTTTCTTCTTCTTTTTCTTCTTCTTCTTCTTCTTTGGATTTggcggaagcaactcgacg
C73F-3caaagaagaagaagaagaagaaaaagaagaagaaaaagaaGCTAGAATTTGCATTGGCAGGTG
C107K-5aaaaagaaaaagaaaaagaagaaaaaaaagaaaaaaaaaggtaagGGCGGAAGCAACTCGACG
C106K-3cttacctttttttttcttttttttcttctttttctttttctttttGCTAGAATTTGCATTGGCAGGTG
C115R-5cgtagaaggagacgtcgtaggagaagacgtaggcgtaggggtaggGGCGGAAGCAACTCGACG
C114R-3cctacccctacgcctacgtcttctcctacgacgtctccttctacgGCTAGAATTTGCATTGGCAGGTG
C104W-5tggtggtggtggtggtggtggtggtggtggtggtggtggggttggGGCGGAAGCAACTCGACG
C103W-3ccaaccccaccaccaccaccaccaccaccaccaccaccaccaccaGCTAGAATTTGCATTGGCAGGTG
C117Delta-5GGCGGAAGCAACTCGACG
C116Delta-3caccatccgctaaggacgacttgggattgtttagcgtcgagttgcttccgccGCTAGAATTTGCATTGGCAG

PKR from other species

C136Cgla-5gaagtttttttttgtaacaagttttgctacactcccttaataaagtcaacATGGAGAACAAGGAGCATCATC
C137Cgla-3ctttaccaaaggacccctctcccagcgttttgacgatttggtagttcccAACACGAGAACCATCAGCCAG
C138Dhan-5gaagtttttttttgtaacaagttttgctacactcccttaataaagtcaacATGTCAGAGCAACAAGGAAGTGATAATC
C139Dhan-3ctttaccaaaggacccctctcccagcgttttgacgatttggtagttcccTATTCTGTTTGCTGGGTTAGCC
C140Klac-5gaagtttttttttgtaacaagttttgctacactcccttaataaagtcaacATGTCGCACGACCCAAATC
C141Klac3ctttaccaaaggacccctctcccagcgttttgacgatttggtagttcccGATATGTTGGCCTTGTGCGTG
C142Kwal-5gaagtttttttttgtaacaagttttgctacactcccttaataaagtcaacATGCCAAGCCACTCGGC
C143Kwal3ctttaccaaaggacccctctcccagcgttttgacgatttggtagttcccAATCCGTTGGTCTTGCCC
C150Cpar-5gtaacaagttttgctacactcccttaataaagtcaacATGTCTGAGCCCAATGATATGCAAGTG
C151Cpar-3ctttaccaaaggacccctctcccagcgttttgacgatttggtagttcccAATGTAGTTGGCAGTATGAGCATTTTGATC
C152Ctro-5gtaacaagttttgctacactcccttaataaagtcaacATGTCTGATCAAAATCAAGGACATAACGATCAAC
C153Ctro-3ctttaccaaaggacccctctcccagcgttttgacgatttggtagttcccTATTCTATTGGCAGGATTTACATTTGGATCAATC
C154Lelo-5gtaacaagttttgctacactcccttaataaagtcaacATGTCGGAATCAAACGATATAGATGTGGATC
C155Lelo-3ctttaccaaaggacccctctcccagcgttttgacgatttggtagttcccAATTCTGTTAGCGGGGTTAGCATTTG
C164Hsap-5gaagtttttttttgtaacaagttttgctacactcccttaataaagtcaacATGCGCAGACTCAGTTCCTGGAG
C165Hsap-3CTTTACCAAAGGACCCCTCTCCCAGCGTTTTGACGATTTGGTAGTTCccgatcttcacccgcccgtc
C170Smik-5gaagtttttttttgtaacaagttttgctacactcccttaataaagtcaacATGGAAATTCTCCCCCACCATCAC
C171Smik-3CTTTACCAAAGGACCCCTCTCCCAGCGTTTTGACGATTTGGTAGTTCCCGATATGCGCACCATCCG
C195Sklu-5gaagtttttttttgtaacaagttttgctacactcccttaataaagtcaacATGTCGCAGCATCCTAATCAGGG
C196Sklu-3CTTTACCAAAGGACCCCTCTCCCAGCGTTTTGACGATTTGGTAGttgccgatgcgttggccatg

Phusion Mutagenesis of Snf1

M36∗T210E-5GTAATTTCTTAAAGGAATCTTGTGGTTC
M37∗T210E-3CATCAGTCATGATGTTTGACAAAC
M68∗4H/8H-5GGCGGAAGCAACTCGACG
M69∗4H-3ATGGTGATGGTGGCTAGAATTTGCATTGGCAGGTGCTGTG
M70∗8H-3ATGGTGGTGATGGTGGTGGTGGTGGCTAGAATTTGCATTGGCAGGTGCTGTG

Pma1 trunctation

C292PMA1-5Gaattgtaatacgactcactatagggcgaattggagctccgcttcctgaaacggagaaac
C293PMA1-3Cctcactcattaggcaccccaggctttacactttatgcttccggctccgtaaaggtatttcgcggagg
C294∗Pma1-trunc5TAAtcctgttgaagtagcatttaatc
C295∗Pma1-trunc3GACACTTCTGGTAGACTTCTTTTC

All oligonucleotides were ordered from Sigma. Desalted, no purification. Upon receipt they were resuspended to 100 μM in 5 mM Tris (pH 8.5). Primers marked with an asterisk were phosphorylated using phosphonucleotide kinase before use.

Yeast strains All S. cerevisiae strains are in the W303 background except PJ694 (for Yeast 2-Hybrid assays). Plasmids Oligonucleotides for mutagenesis of Snf1 All oligonucleotides were ordered from Sigma. Desalted, no purification. Upon receipt they were resuspended to 100 μM in 5 mM Tris (pH 8.5). Primers marked with an asterisk were phosphorylated using phosphonucleotide kinase before use.

Limitations of the study

We have shown that Snf1 is inhibited by genetically upregulating Pma1 activity which causes intracellular alkalization, and also that Snf1 remains more active when Pma1 activity is lowered (causing intracellular acidification), both chemically and genetically. Future experiments should directly measure the internal pH. We acknowledge that our model of internal and external interactions of Snf1 is inferred from yeast two-hybrid data and not by other more direct techniques. Future experiments could include co-immunoprecipitations and biophysical measurements of protein dynamics.

STAR★Methods

Key resources table

Resource availability

Lead contact

Further information and requests for strains and plasmids should be directed to and will be fulfilled by the lead contact Martin Kupiec (martin.kupiec@gmail.com).

Materials availability

All materials generated in this study are available on request to the lead contact.

Experimental model and subject details

Yeasts, plasmids, and growth conditions

Experiments were conducted using S. cerevisiae. DH5a bacteria were used for plasmid propagation and standard procedures. Strains used are listed in Table 1; plasmids used are listed in Table 2. Oligonucleotide used for mutagenesis of Snf1 and Pma1 are listed in Table 3. All strains are related to W303a except for PJ694 used for yeast 2-hybrid assays. Yeasts were transformed with DNA using the frozen lithium acetate method (Knop et al., 1999). Construction of poly-amino acid substitutions in Snf1; and exchange of selective markers in plasmids was by gap repair (Ma et al., 1987) and PCR. PMA1 (from −934 of promoter to +834 of 3'UTR) was inserted into the multiple cloning site of pRS313 and pRS314. The final 18 codons were truncated using Phusion PCR. Four or eight histidines were inserted into plasmid 847 using Phusion PCR. Primers were phosphorylated before use with polynucleotide kinase. Standard sugar concentrations were 4% for glucose (to ensure complete repression of the ADH2 promoter), 3% for glycerol and 2% for ethanol, unless stated differently. For ADH2, ACS1, and JEN1 expression assays, cells were grown in 4% glucose to ensure complete repression of expression, washed 3x with water, and resuspended in media containing 2% ethanol. Cells were grown at 30°C.

Method details

β-galactosidase assays for gene expression

β-galactosidase assays were performed using log phase cells. Cells containing prADH2/ACS1/JEN1::LacZ plasmids were grown overnight in 4% glucose synthetic defined medium, diluted in the morning and grown for an additional 3 h. A sample was taken for measuring (t = 0) and cells washed 3x with 25 mL water before resuspension in indicated medium. After 3 h β-galactosidase activity was measured. Cell concentration was determined by reading 100 μL of cells at 595 nm. 20 μL of cells was added to the β-galactosidase reaction mix (40 μL YPER (Pierce 78,990), 80 μL Z-buffer (120 mM Na2HPO4, 80 mM NaH2PO4, 20 mM KCl, 2 mM MgSO4), 24 μL ONPG (4 mg/mL), 0.4 μL β-mercaptoethanol) and incubated at 30°C for 15 min. Reactions were stopped by addition of 54 μL 1 M Na2CO3. The Eppendorf tubes were centrifuged for 1 min at full speed to pellet the cell debris, and 200 μL supernatant was removed, and absorbance read at 415 nm using a microplate reader. Miller Units were calculated by the equation Miller Units = (1000∗A415)/(time∗volume of cells∗A595-0.055, where the A415 and A595 has been corrected for blanking and path length (final path length = 1 cm). The expression was calculated as a rate of Miller unit increase per hour, and shown as a percentage of the expression rate of WT cells (=100). Typically, in wild-type cells, ADH2 is expressed at 1000 Miller Units per hour in ethanol medium and barely expressed in glucose medium, ACS1 at 1500 Miller Units/hour and JEN1 at 1350 Miller Units/hour in ethanol medium. Three biological replicates were measured. Error bars are +- 1 standard deviation.

Yeast 2-hybrid experiments

S. cerevisiae strain PJ694 (trp1-901 leu2-3,112, ura3-52, his3Δ200, gal4Δ, gal80Δ, prGAL2:ADE2 lys2::prGAL1::HIS3 met2::prGAL7-LacZ) or PJ694 hsp30Δ was transformed with the indicated plasmids expressing proteins or protein fragments fused to either GAD or GBD and PMA1-Δ901 where indicated. Cells were grown overnight in indicated medium and diluted in the morning with the same medium for an additional 3 h before β-galactosidase activity was determined as above. However, for Figure 4F cells were grown in 2% ethanol overnight. After measuring t = 0 glucose to 2% and the indicated concentration of ebselen (in DMSO, from Sigma) was added and cells grown for an additional 4 h before the interaction between the NTD and CTD was determined. Three biological replicates were measured. Error bars are + - 1 standard deviation.

Western blots

For Figures 2E, 2G, and 4B, cells were grown overnight in 4% glucose, diluted in the morning and grown for an additional 3 h. A sample was taken, and the remaining cells washed 2x with 25 mL water and resuspended in 2% ethanol medium for the time indicated (minutes). For Figure 2G (a time-course) because centrifugation activates Snf1 (Orlova et al., 2008) a sample was taken in glucose grown cells (Glc), immediately after washing and resuspension in media containing 2% ethanol (0) and thereafter at indicated times. For Figures 2F and 4E glucose was added to ethanol grown cells to 2%. For Figure 4D, cells were grown in 2% ethanol medium overnight. After harvesting a sample (ethanol) glucose to 2% (+DMSO) and ebselen at the indicated concentrations (in DMSO) was added. Cells were harvested after 15 min. Cells were harvested and protein solubilized in sample buffer by the method developed by the Kuchin group to prevent activation of Snf1 by centrifugation (Orlova et al., 2008). Cells were boiled for 5 min before resuspension in 1x TE and treatment with 0.2 M NaOH. Sample buffer volume was adjusted to give equal OD for all samples (30 μL/OD) and cells boiled for 5 min. After running samples on 8% polyacrylamide gels, the proteins were wet transferred to nitrocellulose membranes. Antibodies used were mouse α-pgk1 (Abcam) for a loading control, rabbit α-GFP (Abcam), mouse α-HA (Santa Cruz) and rabbit α-phospo-T172 (AMPK) (Cell Signalling Technology), all at 1/1000 dilution. Because Snf1-GFP and Rod1-6HA run at 100 kDa, the membrane was cut between 55 and 65 kDa (to check at 72 kDa if native Snf1 is present – it was absent from all samples), the lower half used for Pgk1 and the upper for HA, GFP or phospho-T172. Secondary antibodies were conjugated to HRP. Images were minimally processed using ImageJ.

Microscopy

snf1Δ ADE2 cells carrying prSNF1::Snf1-GFP plasmids with the indicated mutations [deletion of polyHIS (ΔH), substitution to polyA (A), polyF (F), polyK (K), polyW (W), polyY (Y), T210A] were grown overnight in medium containing 4% glucose. 1 mL of cells were washed 3x with water, and resuspended in medium containing 3% glycerol for 30 min. 5 μL of log phase cells were imaged using an EVOS microscope (60x objective) with the GFP filter for GFP and the Texas Red filter for Cherry. The dimension of each panel corresponds to 20 × 20 μm. Cells were not concentrated before imaging, to prevent perturbations to Snf1 activity (Wilson et al., 1996). Images were processed using the brightness/contrast function of Image J, to give a black background. For statistics, over 300 cells were counted. Experiments were repeated at least three times on different days. All figures shown in the manuscript are of identical magnification, a size bar is provided in Figure 2C.

Serial dilutions

snf1Δ cells were transformed with centromeric plasmids expressing prSNF1::SNF1-GFP with the indicated mutations (deletion of polyHIS (ΔH), substitution to polyA (A), polyF (F), polyW (W), polyY (Y)] or empty vector, and snf1Δ mig1Δ cells were transformed with empty vector. Cells were grown in medium lacking leucine with 4% glucose overnight, washed 3x with water and resuspended in water at a concentration of 7.5OD600/mL. Cells were serially diluted 10-fold and spotted onto-Leu plates with either 2% glucose, 2% galactose, or 3% glycerol as the carbon source and incubated at 30°C for 2 days. Cells were imaged using a Samsung A50 camera.

Quantification and statistical analysis

Data were analyzed using a two-tailed unpaired t-test (GraphPad). Unless indicated, data points were compared to the strain/condition designated as 100%. ∗ = p0.05, ∗∗ = p0.01, ∗∗∗ = p0.001, NS = Not significant.
REAGENT or RESOURCESOURCEIDENTIFIER
Antibodies

Rabbit anti-GFPAbcamAB290 RRID: Ab_303395
Rabbit anti-phospho T172 (AMPK)Cell Signalling Technologies#2535 RRID: Ab_331250
Mouse anti-HASanta CruzF-7 RRID: Ab_2894930
Mouse anti-Pgk1AbcamAb113687 RRID:Ab_10861977
Goat anti-mouse HRPJackson115035003 RRID: Ab_10015289
Goat anti-rabbit HRPJackson111035144 RRID: Ab_2307391

Bacterial and virus strains

DH5αLab StockN/A

Chemicals, peptides, and recombinant proteins

2-Nitrophenyl - beta-D-galactopyranoside (ONPG)SigmaN1127
Y-PERThermo Scientific#78990
2-mercaptoethanolFlukaM6250
Polynucleotide kinaseSigmaPNK-RO
EbselenSigmaE3520
Na2CO3Sigma791768

Deposited data

N/A

Experimental models: Organisms/strains

Please see Table 1

Oligonucleotides

Please see Table 3

Recombinant DNA

Please see Table 2

Software and algorithms

Image JImage JN/A
GraphpadGraphpad: https://www.graphpad.com/quickcalcs/ttest1/?format=SDN/A
  77 in total

1.  Oxidization without substrate unfolding triggers proteolysis of the peroxide-sensor, PerR.

Authors:  Bo-Eun Ahn; Tania A Baker
Journal:  Proc Natl Acad Sci U S A       Date:  2015-12-17       Impact factor: 11.205

2.  Ptc1 protein phosphatase 2C contributes to glucose regulation of SNF1/AMP-activated protein kinase (AMPK) in Saccharomyces cerevisiae.

Authors:  Amparo Ruiz; Xinjing Xu; Marian Carlson
Journal:  J Biol Chem       Date:  2013-09-09       Impact factor: 5.157

3.  A member of the sugar transporter family, Stl1p is the glycerol/H+ symporter in Saccharomyces cerevisiae.

Authors:  Célia Ferreira; Frank van Voorst; António Martins; Luisa Neves; Rui Oliveira; Morten C Kielland-Brandt; Cândida Lucas; Anders Brandt
Journal:  Mol Biol Cell       Date:  2005-02-09       Impact factor: 4.138

4.  Structural and functional characterization of 2-oxo-histidine in oxidized PerR protein.

Authors:  Daouda A K Traoré; Abdelnasser El Ghazouani; Lilian Jacquamet; Franck Borel; Jean-Luc Ferrer; David Lascoux; Jean-Luc Ravanat; Michel Jaquinod; Geneviève Blondin; Christelle Caux-Thang; Victor Duarte; Jean-Marc Latour
Journal:  Nat Chem Biol       Date:  2008-12-14       Impact factor: 15.040

5.  Activation of yeast Snf1 and mammalian AMP-activated protein kinase by upstream kinases.

Authors:  Seung-Pyo Hong; Fiona C Leiper; Angela Woods; David Carling; Marian Carlson
Journal:  Proc Natl Acad Sci U S A       Date:  2003-07-07       Impact factor: 11.205

6.  Snf1 Phosphorylates Adenylate Cyclase and Negatively Regulates Protein Kinase A-dependent Transcription in Saccharomyces cerevisiae.

Authors:  Raffaele Nicastro; Farida Tripodi; Marco Gaggini; Andrea Castoldi; Veronica Reghellin; Simona Nonnis; Gabriella Tedeschi; Paola Coccetti
Journal:  J Biol Chem       Date:  2015-08-26       Impact factor: 5.157

7.  Detection of endogenous Snf1 and its activation state: application to Saccharomyces and Candida species.

Authors:  Marianna Orlova; LaKisha Barrett; Sergei Kuchin
Journal:  Yeast       Date:  2008-10       Impact factor: 3.239

8.  Regulation and pH-dependent expression of a bilaterally truncated yeast plasma membrane H+-ATPase.

Authors:  A B Mason; T B Kardos; B C Monk
Journal:  Biochim Biophys Acta       Date:  1998-07-17

9.  Two novel effectors of trafficking and maturation of the yeast plasma membrane H+ -ATPase.

Authors:  Yosef Geva; Jonathan Crissman; Eric C Arakel; Natalia Gómez-Navarro; Silvia G Chuartzman; Kyle R Stahmer; Blanche Schwappach; Elizabeth A Miller; Maya Schuldiner
Journal:  Traffic       Date:  2017-08-16       Impact factor: 6.215

10.  CHESS: a new human gene catalog curated from thousands of large-scale RNA sequencing experiments reveals extensive transcriptional noise.

Authors:  Mihaela Pertea; Alaina Shumate; Geo Pertea; Ales Varabyou; Florian P Breitwieser; Yu-Chi Chang; Anil K Madugundu; Akhilesh Pandey; Steven L Salzberg
Journal:  Genome Biol       Date:  2018-11-28       Impact factor: 13.583

View more

北京卡尤迪生物科技股份有限公司 © 2022-2023.