Life continuously creates genetic innovation, e.g., genes coding for proteins with novel functions. Most new genes emerge by modification of existing genes rather than de novo (Long et al., 2013). The first who connected the generation of new functions with the accumulation of mutations and duplication of pre-existing genes was Muller, 1936. Later on, different theories arose discussing the point at which divergence takes place (Ohno, 1970; Jensen, 1976; Lynch and Force, 2000; Bergthorsson et al., 2007). Modern evolutionary models assume an ancestral gene to acquire a novel function while retaining its original function (Bergthorsson et al., 2007; Hughes, 1994; Aharoni et al., 2005; Des Marais and Rausher, 2008; Näsvall et al., 2012). If a novel function is beneficial for the organism, the gene variant is subjected to positive selection. Consecutive duplication of the gene allows the establishment of the new function in the population. This model of divergence prior to gene duplication allows for the accumulation and selection of mutations with adaptive potential while the original gene function is maintained, thus creating a ‘generalist’ gene with multiple functions (Soskine and Tawfik, 2010).

When considering the evolution of novel gene functions, the substrate range of ‘classical’, soluble enzymes is arguably the most archetypal field of study, important for the rapid evolution of antibiotic resistance, degradation of artificial xenobiotics, or use of new nutrients (Aharoni et al., 2005; Copley et al., 2023). The substrate spectrum of an enzyme is often a key determinant of the fitness of an organism. Generally, enzymes are thought of as being highly specific for a certain substrate or a group of chemically similar substrates. Nonetheless, many enzymes studied in isolation show side activities toward other substrates (Copley et al., 2023; Khersonsky and Tawfik, 2010; Andorfer et al., 2017). These side activities are often obscured in the complex setting of a biological cell, and as such contribute little to the fitness of the organism. However, these side activities are thought of as a major source for genetic innovation: proteins can evolve the ability to use new substrates by improving an already existing, but low, side activity (Janssen et al., 2005; Jørgensen et al., 2017; Caligiore et al., 2022; Zheng et al., 2020).

Membrane transporters form a special class of enzymes, as they generally do not break or make covalent bonds. Hence, they experience fewer constraints of the chemistry of their substrates. In principle, any small molecule can be transported through a lipid membrane. Thus, it is not surprising that life has evolved transport systems for virtually every metabolite, ion, and xenobiotic (Saier et al., 2014; Saier et al., 2021). In many branches of the tree of life, transporter genes underwent processes of repeated duplication and divergence, leading to extended gene families; many genomes have an array of transporter paralogs (Copley et al., 2023; Saier et al., 2021; Jack et al., 2000; Perland and Fredriksson, 2017). Typically, each paralog displays a distinct substrate specificity profile, and thus a distinctly evolved function. Some notable examples of gene family expansion are found within the APC (amino acid-polyamine-organocation) superfamily of transporters, which include: the neurotransmitter transporter family in animals (20 paralogs in human, TC 2.A.22), the AAAP family in plants (58 paralogs in rice, TC 2.A.18), the AAT family in Proteobacteria (11 paralogs in Escherichia coli, TC 2.A.3.1), and the yeast amino acid transporter (YAT) family in fungi (18 paralogs in budding yeast, TC 2.A.3.10). The APC superfamily is the second-largest group of membrane transport proteins (Vastermark et al., 2014).

Most APC superfamily members are transporters for amino acids, their analogs, or related amines. Some transporters are generalists, accepting a wide range of substrates; others are considered specialists, recognizing one or a few related substrates (Jack et al., 2000; Bianchi et al., 2019). How do such patterns arise? Does the evolution of specialists move through a generalist stage? Or do new specialists evolve through switching from one specific substrate to another? For classical enzymes, many lines of evidence point to the former (Aharoni et al., 2005; Copley et al., 2023; Khersonsky and Tawfik, 2010; Risso et al., 2013; Mascotti et al., 2018). However, the evolution of new traits in the special case of transporters is understudied, partly because of a lack of appropriate experimental evolution methods (Bali et al., 2018).

Here, we propose that APC-type transporters can evolve novel functions through generalist intermediates, similar to classical enzymes. We showcase the presence of so far unknown substrates of five out of seven studied YAT and develop a growth-based selection platform. Utilizing in vivo experimental evolution, under purifying selection for the original function and simultaneous selection for a new function, we evolve the substrate spectrum of two of these transporters. We show that the acquired mutations result either in improving an existing weak activity or in establishing activity for a new substrate. At the same time, each mutation has a distinct impact on the fitness of the organism for each of the transporter’s original substrates.

The budding yeast, Saccharomyces cerevisiae, typically encodes 18 members of the YAT family, each with a different substrate profile (Bianchi et al., 2019). This array of transporters enables yeast to use 20 different l-amino acids as the sole nitrogen source in minimal media. These include 17 of the 20 proteinogenic amino acids (the exceptions are Cys, His, and Lys) and 3 non-proteinogenic ones (γ-amino butyric acid [GABA], L-ornithine [Orn], and L-citrulline [Cit]). Here, we use an engineered yeast strain (Δ10AA) that lacks 10 membrane transporters for amino acids (Besnard et al., 2016). This strain is thus severely deficient in the uptake of amino acids, preventing it from using most of them as the nitrogen source. When an amino acid transporter is genetically reintroduced, its substrate profile directly determines which amino acids can support growth. In the following text, the gene names are used in place of the protein names for consistency (e.g. AGP1 instead of Agp1). The Δ10AA yeast strain expressing different transporter genes was used in all growth and transport assays.

Substrate spectrum of seven YAT

We investigated the substrate profile of seven YAT, by measuring the growth rates of Δ10AA strains overexpressing each transporter gene from pADHXC3GH, in the presence of 2 mM of each amino acid that served as nitrogen source. We determined growth rates in order to study each transporter’s substrate range separately. For five of the seven transporters tested, we found substrates that were not previously reported in literature. These substrates were identified by significantly increased growth rates of Δ10AA expressing the respective transporter in comparison to the strain carrying the empty plasmid (vector control) (Table 1). Because of high background growth of the vector control on Arg and Orn (probably due to the activity of the VBA5 transporter Shimazu et al., 2012), these two amino acids were excluded from further analysis. Additionally, some growth of the empty vector control was observed on Gly, Trp, Leu, Met, and Gln, which complicates the interpretation of the data. The observed growth could indicate an unknown substrate specificity of an endogenous transporter. We observed that some of the newly characterized substrates support growth rates in similar ranges as the transporter’s substrates already reported in literature, while other substrates support very slow growth rates, which we call weak or promiscuous activities (Figure 1 and Figure 1—figure supplement 1).

Table 1

Transporter gene	Uniprot accession number (UniProt Consortium, 2023)	Substrates described in literature (reviewed in Bianchi et al., 2019)	Substrates found in the present study
AGP1	P25376	Ala, Asn, Asp, Cys*, Gln, Glu, Gly, His*, Ile, Leu, Met, Phe, Pro, Ser, Thr, Trp, Tyr, Val, GABA	Ala, Asn, Asp, GABA, Gln, Glu, Ile, Leu, Met, Phe, Pro, Ser, Thr, Trp, Tyr, Val, (Cit)
BAP2	P38084	Ala, Cys*, Ile, Leu, Met, Phe, Tyr, Val	Ala, Asn, Asp, Cit, Gln, Glu, Gly, Ile, Leu, Met, Phe, Pro, Ser, Thr, Tyr, Val
CAN1	P04817	Arg, His*, Lys*, Orn, Ser	–
HIP1	P06775	His*	–
LYP1	P32487	Lys*, Met	Ala, Asn, Met, Phe, Ser, Val
MMP1	Q12372	S-Methylmethionine*	Val
PUT4	P15380	Ala, Gly, Pro, GABA	Ala, GABA, Pro, Ser, Val

Figure 1 with 2 supplements see all

Download asset Open asset

Figure 1—source data 1

Growth rates of yeast amino acid transporter (YAT) expressing yeast.

Growth rate values of Δ10AA expressing either one of the seven different wild-type YAT genes (AGP1, BAP2, CAN1, HIP1, LYP1, MMP1, PUT4) from pADHXC3GH and the empty vector control on 2 mM of each amino acid. The growth rates were calculated based on Figure 1—figure supplement 1.

https://cdn.elifesciences.org/articles/93971/elife-93971-fig1-data1-v1.xlsx

Download elife-93971-fig1-data1-v1.xlsx

For example, BAP2, previously reported to transport seven amino acids, was found to support growth on 16 amino acids (specific growth rate µ=0.07–0.31 hr^–1), making it a broad-range transporter. The differences on the observed growth rates on each amino acid can be explained by differences in the activity of the overexpressed transporter for these substrates or changes in the expression level of the endogenous transporter. Specifically, the twofold higher growth rate on Leu compared to Pro could be caused by a higher transcription rate of the endogenous transporter in the presence of Leu (Didion et al., 1996; Magasanik and Kaiser, 2002). Another transporter which showed a broader substrate range than previously reported is PUT4. The PUT4 transporter was found to support fast growth on Ser (µ=0.21 hr^–1) and slow growth on Val (µ=0.03 hr^–1), in addition to the already reported transport of Ala, Gly, Pro, and GABA. AGP1, an already known broad-range transporter, additionally supports very slow growth (µ=0.03 hr^–1) on Cit (l-citrulline), which is lower than the growth rate on any of the previously known substrates (µ=0.11–0.31 hr^–1). As a comparison, growth on ammonium, which does not depend on amino acid transporters, yielded the highest measured growth rate of µ=0.36 hr^–1 (vector control shown in Figure 3—figure supplement 1).

Experimental evolution of membrane transporter specificity

We tested whether the substrate range of amino acid transporters could be changed by experimental evolution, either by improving promiscuous activities or by evolving toward completely new substrates. For the former, we chose to study AGP1 under a selective pressure for Cit uptake. For the latter, we studied PUT4 under selective pressures for uptake of Asp and Glu. For the in vivo evolution, the OrthoRep system was used, which allows for random mutagenesis of a gene while it is actively expressed in yeast (Ravikumar et al., 2014; Ravikumar et al., 2018). Each of the two transporter genes was encoded on this system and introduced into the transport-deficient Δ10 strain, with additional Δura3 and Δhis3 deletions needed for selection.

To select for mutants with changed substrate specificity, we employed a dual-selection scheme: The cultures were grown in a minimal medium containing a nitrogen-limiting mixture of natively transported amino acids (1 mM final concentration) to ensure low-level growth of the culture. The mixture of amino acids was composed of equimolar amounts of 17 original substrates for AGP1 and 2 original substrates for PUT4. A concentration of 3 mM of the target substrate was added in the selection for evolved mutants. We expect that in this setup nonfunctional transporters are purged from the evolving population, as these do not support uptake, neither of the original substrates nor of the target amino acid. Transporters with the original substrate range are kept in low numbers, and transporter variants with novel substrates can increase to high numbers as they support growth on the target substrate present at threefold higher concentration than the original substrates.

The cultures of AGP1 and PUT4 evolution strains were passaged for multiple generations in the selective media. Colonies that supported growth on the target amino acid were isolated, and the transporter gene was sequenced. The gene was cloned to the pADHXC3GH expression vector and reintroduced into Δ10AA cells. Three AGP1 variants from Cit selection and three PUT4 variants each from Asp and Glu selection (Table 2) conferred the ability to grow on the respective amino acids (Figure 2—figure supplement 1 and Figure 2—figure supplement 2).

Table 2

Evolved variant	Nucleotide substitutions (nonsynonymous in bold)	Amino acid substitutions	Abbreviated names of site-directed mutants
AGP1-Cit1	T117C T1001A T1191A A1234G T1329C C1461G T1899C	I334N (=N) I412V (=V)	AGP1-N AGP1-V AGP1-NV AGP1-G AGP1-T
AGP1-Cit3	T117C A308G T915C T1001A C1451G C1461G T1684C	D103G I334N (=N) A484G (=G) F562L
AGP1-Cit11	T111C T121C G1450A A1609G	S41P A484T (=T) I537V
PUT4-Asp1	T620C T1474C	L207S (=S) F492L	PUT4-S
PUT4-Asp2	T620C T1415C T1474C	L207S (=S) V472A F492L
PUT4-Asp3	T191C T620C C734T G956A	I64T L207S (=S) S245F G319D
PUT4-Glu1	G225A T620C	L207S (=S)
PUT4-Glu2	T191C T620C G956A	I64T L207S (=S) G319D
PUT4-Glu3	T191C T620C	I64T L207S (=S)

Each evolved AGP1 variant had multiple nonsynonymous mutations, most of which occurred in the predicted transmembrane part of the protein (Figure 2 and Figure 2—figure supplement 3). In two of the three variants, the I334N (transmembrane helix TM6) mutation was found. The I334 position is located at the permeation pathway of the protein as predicted from AlphaFold (Jumper et al., 2021). A mutation in position A484 (TM10) was found in two variants (A484G and A484T). Additionally, mutations were found at the intracellular N-terminus (S41P, D103G) and the predicted transmembrane helices further away from the permeation pathway (I537V in TM11, F562L in TM12).

Figure 2 with 4 supplements see all

Download asset Open asset

Most PUT4 variants had multiple nonsynonymous mutations, again mostly in the predicted transmembrane part of the protein (Figure 2 and Figure 2—figure supplement 3). Interestingly, each variant from both Asp and Glu selections had the same L207S mutation (TM3). Overlay of homologous and structurally similar proteins indicates that this position is part of the substrate-binding site of the transporter (Figure 2—figure supplement 4; Gibrat et al., 1996). A change from a large hydrophobic residue to a small hydrophilic residue is thus expected to considerably change the properties of the substrate-binding site. Furthermore, in two of the three variants from Asp selection, the F492L mutation (TM10) was found. Mutations were also found at the intracellular N-terminus (I64T), between TM4 and TM5 (S245F), in TM6 (G319D), and in TM9 (V472A).

Single mutations are sufficient for altered substrate specificity

We wondered whether the altered substrate specificities could be conferred by single mutations. We made site-directed mutations in the wild-type transporter genes and tested the variants for uptake of Cit (AGP1) and Asp/Glu (PUT4). Generally, only mutations in the transmembrane domains were considered. The shorthand names for each site-directed mutant can be found in Table 2.

For the AGP1 gene, mutations in positions that occurred multiple times were investigated (i.e. I334N, A484G, and A484T; respective abbreviations AGP1-N, AGP1-G, and AGP1-T). Additionally, the I412V mutation was included in the study (shorthand AGP1-V), along with a double mutant combining I334N and I412V (abbreviation AGP1-NV, recreating the genotype observed in AGP1-Cit1). To compare the growth of the mutants on different media, the growth rate of Δ10AA expressing each variant from the plasmid pADHXC3GH was measured. The surface expression of the transporters was estimated by analyzing the fluorescence at the cell envelope (Figure 3—figure supplement 1A) and was similar for the wild-type and the evolved proteins (Figure 3—figure supplement 1B). Furthermore, there was no statistically significant difference between the growth rate of the strains under conditions where the nitrogen source was not dependent on amino acid uptake (Figure 3—figure supplement 1C). Strikingly, when grown in media with Cit as the sole nitrogen source, each of the variants showed a large increase in growth rate relative to the wild-type (Figure 3A). Specifically, the growth rate increased 1.9-fold for AGP1-N, 2.4-fold for AGP1-V, 3.2-fold for AGP1-G, and 2.4-fold for AGP1-T. The double mutant AGP1-NV showed a growth rate increase that is higher than each of its single mutants (2.8-fold). In separate uptake assays with radiolabeled Cit, we found a similar general trend of higher uptake rate of Cit by the mutants (Figure 3B and Figure 3—figure supplement 3A). We thus conclude that each of the tested single mutations is adaptive for growth on Cit as the nitrogen source.

Figure 3 with 3 supplements see all

Download asset Open asset

Regarding the PUT4 gene, we focused on the L207S single mutation since it occurred in all evolved variants (mutant abbreviation PUT4-S). The growth rates of Δ10AA expressing either the wild-type or the PUT4-S variant from pADHXC3GH were measured. The growth rate of the strains under conditions where the nitrogen source was not dependent on amino acid uptake was not significantly different (Figure 3—figure supplement 2C). While the wild-type transporter conferred no growth in the presence of Asp and Glu, the PUT4-S variant allowed growth on both of these substrates, showing that the single mutation is sufficient to change the specificity range of the transporter (Figure 3C). To confirm that the observed growth is due to transport of a novel substrate by PUT4-S, we conducted a whole-cell uptake assay with radiolabeled Glu. Indeed, Glu was only taken up by the mutant, whereas the wild-type PUT4 behaved like the vector control (Figure 3D and Figure 3—figure supplement 3B). Even though the surface expression of PUT4-S was higher than that of the wild-type (Figure 3—figure supplement 2B), this difference is unlikely to reason for the observed gain of function.

Single mutations have distinct impacts on the fitness for the original substrates

Above, we established that amino acid transporters can be evolved toward new specificity or increased activity of a substrate by just one mutation. Next, we wanted to know how these adaptive mutations affect the original substrate range and specificity of the transporters. Are the original substrates still transported by the variants? And if so, how does the transport compare to the wild-type? For that, we compared the fitness effects of the AGP1 and PUT4 variants for their respective substrates. The relative fitness of each mutant in the studied media was calculated from the growth rate of yeast expressing that mutant, normalized to the growth rate of yeast expressing the wild-type transporter. Additionally, we measured the uptake rates of a set of original substrates in whole cells to investigate if the growth fitness effects were due to changed transport activity.

For the AGP1 variants, we measured growth rates for 17 substrates (Figure 4 and Figure 4—figure supplement 1) and found that none of the single and double mutants had lost the capacity to grow on any of the original substrates. As a general trend, the AGP1 variants had a lower fitness than the wild-type, although some had a higher fitness for a given substrate.

Figure 4 with 2 supplements see all

Download asset Open asset

The AGP1 variants had different growth rates on the original substrates (Table 3 and Figure 4B). The amino acid substitutions I334N and I412V had a different effect on the relative fitness depending on whether they were analyzed individually or together in the double mutant. Specifically, the AGP1-V variant showed significantly increased relative fitness on GABA as nitrogen source, while showing no significant differences in any of the other substrates. The AGP1-N variant showed decreased relative fitness on the negatively charged amino acids Asp and Glu, as well as on Ala. The double mutant AGP1-NV showed decreased relative fitness on the majority of the tested amino acids (a total of 12 out of 17), including the ones in which the single mutant AGP1-N also exhibited decreased relative fitness. Interestingly, the AGP1-NV did not show significantly increased relative fitness in GABA, in contrast to AGP1-V. Despite the higher relative fitness of the double mutant on Cit (2.8) compared to that of the single mutants, the variant’s low relative fitness on the rest of the amino acids has a greater cost on the overall fitness of that strain (0.78). Also, the observed mistargeting of the transporters to internal membranes in the cases of AGP1-N and AGP1-NV could indicate misfolding of a fraction of the proteins, which could impact the relative fitness (Figure 5A).

Table 3

AGP1 variant	Substrates with significantly increased relative fitness	Substrates with significantly decreased relative fitness	Overall relative fitness (mean of 17 original substrates)	Relative fitness on Cit
AGP1-N		Ala, Asp, Glu	0.92	1.9
AGP1-V	GABA		0.98	2.4
AGP1-NV		Ala, Asn, Asp, Glu, Ile, Leu, Met, Phe, Pro, Ser, Thr, Tyr	0.78	2.8
AGP1-G	GABA		1.0	3.2
AGP1-T		Asp, Gln, Glu, Ile, Leu, Met, Phe, Thr, Val	0.86	2.4

Figure 5 with 2 supplements see all

Download asset Open asset

The fitness effects were drastically different between AGP1-G and -T, both having A484 substituted. The AGP1-G variant showed significantly increased relative fitness on GABA and a trend of overall increased relative fitness on small aliphatic amino acids of intermediate hydrophilicity (Pro, Gly). Simultaneously, its relative fitness on Cit was higher than any of the other mutants. AGP1-T showed mainly negative relative fitness effects, with decreased relative fitness for charged amino acids and the majority of non-polar amino acids (a total of 9 out of 17).

The idea that the observed fitness effects are due to changed transport activity and/or affinity constant for the substrate was investigated by measuring the uptake rates of radiolabeled Phe and Glu. The uptake was measured at two different amino acid concentrations, namely 0.1 mM and 2 mM, which are respectively lower and higher than the K_m of the wild-type protein for amino acids (K_m = 0.2–0.79 mM [reviewed in Bianchi et al., 2019]; Figure 5 and Figure 5—figure supplement 1). For Glu it was found that the uptake rates at both concentrations (Figure 5C and D) correlate well with the relative fitness. An almost identical pattern was observed when the uptake and the growth rates of the variants were compared at 2 mM Glu (Figure 5B and D). The fact that AGP1-N, AGP1-NV, and AGP1-T showed a pattern of both reduced growth rate and reduced uptake rate compared to the wild-type indicates that the relative fitness indeed is linked to the activity of the transporter (Figure 5—figure supplement 2B). Equivalently, in the case of Phe as nitrogen source, the relative fitness and uptake rate of the AGP1 variants followed similar trends. At 2 mM Phe, the effects on uptake were almost congruent to the growth rates (Figure 5E and G, Figure 5—figure supplement 2A). The reduced uptake and growth rate of AGP1-NV compared to the wild-type indicates that the lower fitness of the AGP1-NV on Phe possibly stems from its reduced transport rate, but an effect of the internalization of the transporter cannot be excluded (Figure 5A).

The growth rates of Δ10AA expressing PUT4 and PUT4-S showed that the mutation considerably affects the relative fitness of the strain in three of the original substrates (Figure 6B). Specifically, PUT4-S was associated with a significant fitness loss during growth on 2 mM Ala and GABA, and a fitness gain in Val. To investigate the phenotype further, we followed the uptake of radiolabeled Ala and GABA at a concentration of 10 μM, well below the reported K_m (reviewed in Bianchi et al., 2019; Figure 6 and Figure 6—figure supplement 1). As with AGP1 variants, the relative fitness and uptake rate trends coincided, indicating that the slower amino acid uptake may cause the slower growth of PUT4-S on both Ala and GABA (Figure 6C and D). Additionally, we investigated the uptake of Gly, another known substrate of PUT4. Due to high background growth of the Δ10AA strain on Gly, there was no significant difference in the growth rates between PUT4-expressing and vector control cultures. However, the uptake assays showed a clear transport activity of the wild-type PUT4 strain, which decreased dramatically in the case of PUT4-S (Figure 6E). There is thus a clear cost of the mutation on the uptake of Gly.

Figure 6 with 1 supplement see all

Download asset Open asset

The wild-type transporter supported very slow growth on Val (µ=0.03 hr^–1) and we therefore consider its transport a promiscuous activity. This coincides with the larger size and hydrophobicity of the Val molecule relative to Ala, GABA, Gly, Pro, and Ser. For Val, the fitness of PUT4-S was dramatically increased compared to wild-type (2.4-fold), which could be explained with the predicted size increase of the substrate-binding site of the transporter (see above). Thus, the gain-of-function mutation L207S affected the transport of at least four original substrates, three of which negatively, and one positively.

Single mutations can significantly broaden the substrate range of amino acid transporters

All variants with increased fitness for one substrate can still transport the wild-type’s original substrates. For AGP1, each of the four tested single mutants as well as the double mutant retained the ability to confer growth on all amino acids that we tested.

Similarly, for PUT4, the L207S mutation added transport of Glu and Asp, while still retaining the ability to transport the original substrates Ala, GABA, Gly, Pro, Ser, and Val. Since we established that one single evolutionary mutation extended the strain’s growth by two additional amino acids, we wondered whether this mutation also impacted the overall substrate range of the transporter. Therefore, growth assays for the PUT4-S strain were conducted with 18 amino acids (Figure 7A). Remarkably, the transporter variant conferred growth on 15 of the studied substrates. Specifically, the evolved variant carrying the L207S mutation transports Asp and Glu (the amino acids used for the initial selection), as well as Asn, Cit, Gln, Thr, Ile, Leu, and Met in addition to the six substrates of the wild-type transporter PUT4. Regarding the growth on Gln, we note that although the calculated growth rates of PUT4-S and PUT4 were not significantly different, the shape of the growth curve implies that Gln is indeed taken up (Figure 4—figure supplement 1). In the case of growth on Phe and Trp, the results were not conclusive due to the high background growth of the vector control strain.

Figure 7

Download asset Open asset

Upon grouping the amino acids according to their size and hydropathy index (Kyte and Doolittle, 1982; Zamyatnin, 1972), a clear trend emerged for the novel substrates of PUT4-S. In addition to the small substrates of low and intermediate hydrophilicity, the evolved variant facilitates the uptake of larger hydrophobic amino acids (Ile, Leu, Met) as well as hydrophilic amino acids (Thr, Asp, Asn, Gln, Glu, Cit) (Figure 7B, colored areas 2 and 1, respectively). Thus, we showcase that a single adaptive mutation is sufficient to evolve a narrow-range transporter (PUT4: 6 substrates) to a broad-range transporter (PUT4-S: 15 substrates).

Different theoretical models have been developed to explain the origin of new functions in proteins. Many favor the idea that new functions are acquired prior to gene duplication while the original functions are maintained (Soskine and Tawfik, 2010). While these models have been used to describe the evolution of functions in soluble enzymes, we showcase that the theory also applies to membrane transporters. To simulate natural transporter evolution, we leveraged in vivo evolution of two YAT, AGP1 and PUT4, with selective pressure for transport of novel substrates.

We show that the substrate spectrum of five wild-type YAT is substantially broader than previously described. The use of a growth-based screening process enabled us to identify so far unknown substrates (up to nine in the case of BAP2). The transporter genes were expressed from a multicopy plasmid (2µ origin of replication), where the number of plasmid copies that each daughter cell receives is stochastic (Chan et al., 2013). Cells with different copy numbers could potentially have an advantage in our growth rate-based measurements of transporter activity. Thus, the mean copy number of the culture is expected to change during the 3-day cultivation period. The growth rate is therefore a measure of the fitness of the strain when the transporter gene copy number is optimized, which is probably the reason why the growth-based assay is so sensitive for promiscuous transporter activities that have not been described before.

It is apparent that the relative fitness of each strain expressing a transporter depends on the amino acid provided. For the strain expressing AGP1, Cit supports very low growth rates indicative of a promiscuous uptake activity by AGP1 that has not been reported previously. We show that upon the presence of selective pressure for the uptake of Cit, this weak transport activity of AGP1 is increased by adaptive mutations. Each of the five studied mutants resulted in a gain of relative fitness for growth on Cit. The differential Cit uptake rates of the single and double mutants showcase how in vivo experimentally generated single mutations can alter the substrate spectrum of the YAT without the loss of ancestral functions. Our observations are in line with studies of substrate promiscuity being the result of adaptive mutations in different membrane transporters (Regenberg and Kielland-Brandt, 2001; Tutulan-Cunita et al., 2005; Adler and Bibi, 2005; Gros and Schuldiner, 2010; Ghaddar et al., 2014; Young et al., 2014; Meier et al., 2023). The altered specificity for a substrate, which is mirrored in our case in the altered fitness, could be a result of a different affinity constant for the amino acid and/or maximal transport rate.

We also established that through applying evolutionary pressure upon the PUT4 membrane transporter, single mutations can emerge in the population which provide a gain of function. The newly obtained function is essential for the survival of the organism in the new environment. Interestingly, the same mutation (L207S) in the PUT4 gene was found when evolutionary pressure was applied for the growth on media supplemented with either Asp or Glu as the sole nitrogen source. The L207 residue is predicted to be part of the substrate-binding site of the transporter. Structurally equivalent sites on other APC superfamily transporters have been shown to affect the function of the studied proteins (Edwards et al., 2018; Ghaddar et al., 2014). Specifically, the non-conserved V104 in Aquifex aeolicus LeuT shapes the binding pocket and influences the substrate specificity (Singh et al., 2008; Yamashita et al., 2005), while T180 in CAN1 is expected to contribute to the transporter’s interaction with its arginine substrate (Ghaddar et al., 2014; Gournas et al., 2015). In the case of PUT4-S, the replacement of the large hydrophobic leucine residue with a small hydrophilic serine residue is likely to increase both the size and the hydrophilicity of the binding pocket. This could be the reason that the evolved transporter now facilitates the uptake of hydrophilic as well as large hydrophobic amino acids.

The evolved mutations in AGP1 are also positioned in the transmembrane helices that create the permeation pathway of LeuT-fold transporters (Kazmier et al., 2017). I334 in TM6 of AGP1 is conserved in some members of the YAT family (DIP5, CAN1, LYP1, ALP1, TAT3). I412 in TM8 is part of the conserved sequence [IV]-L-[ITAVL] and thus in AGP1-V and AGP1-NV, Ile is changed to the only other common residue, Val, which is found in DIP5, LYP1, and ALP1 (Gournas et al., 2015). Reportedly, single site substitutions in conserved regions of CAN1 and GNP1 have been shown to broaden the substrate range of these YAT and accommodate for the uptake of Cit (Regenberg and Kielland-Brandt, 2001). Substitutions in non-conserved regions can also have an effect on a transporter’s substrate spectrum. The non-conserved T456 in CAN1 is important for the substrate specificity of the protein. Single site substitution of T456 can broaden the substrate spectrum of CAN1 by one amino acid, while S176N/T456S abolishes the uptake of arginine, a high affinity substrate of CAN1, and supports high affinity for the uptake of lysine similar to that of LYP1 (Ghaddar et al., 2014). Interestingly, T456 lies two residues away from N454 in CAN1, whose equivalent position is A484 in AGP1 (Gournas et al., 2015).

In both AGP1 and PUT4, the evolved variants retain their original substrate spectrum albeit with altered specificities. Our results showcase that specific adaptive mutations for a novel substrate often decrease the fitness for the original substrates, creating a situation of trade-off between new and old function. However, in some cases, the fitness for single original substrates actually increases (e.g. AGP1-G and AGP1-V for GABA, or PUT4-S for Val). Depending on the environment, a mutation can thus be adaptive for multiple substrates at the same time, including ones that were not selected for. For evolution on Cit, the mutant with the highest fitness gain (AGP1-G) showed no loss of relative fitness for the original substrates. Thus, this mutation could potentially occur in the population even without the presence of Cit as selective pressure or even outcompete others in a natural setting where the selective pressure is applied. On the contrary, the combination of two adaptive mutations (AGP1-NV) had a higher fitness cost for original substrates than the separate mutations, but also a higher fitness gain for the novel substrate, which can be viewed as the first step toward specialization. Thus, in a natural setting, such a variant would only be viable if the novel substrate gives a very high selective advantage as compared to the original ones. Alternatively, gene duplication could lead to two distinct transporter genes, one specializing in the novel substrate, while the other (re-)specializes for the original substrates (Soskine and Tawfik, 2010). To generalize our findings to natural evolution, it is important to keep in mind that our selection process was intentionally designed to emulate a newly encountered substrate, for which no suitable transporter is yet available. This led to evolution toward substrate range expansion and produced broad-spectrum transporter variants. In wild-type organisms, broad specificities are not always favorable. Each organism typically possesses multiple nutrient transporters with partial redundancy in their substrates, making it possible to finely tune-up and down-regulate specific transport activities according to the nutrient availability (Bianchi et al., 2019). Therefore, our experimentally derived transporter variants would probably not fare well in natural settings, as they come from an organism that already has evolved a full suite of amino acid transporters. They do, however, give information about evolvability of transporters when faced with a new solute, whose transport leads to increased fitness of the organism. An example of a transporter with a very special substrate can be found in the cystine transporter CYN1, which is distantly related to LYP1. Only a fraction of all yeast species has orthologs of this transporter and the respective cystine uptake activity (Yadav and Bachhawat, 2011). Seemingly, not all yeasts profit from the ability to assimilate cystine from the medium, and therefore do not possess a suitable transporter. At one point in time, however, an ancestor of LYP1 evolved toward transport of cystine and ultimately turned into the highly selective CYN1 found today. These processes take place over long periods of time and are hard to observe in the field. Our experiments emulate the first steps in such an evolutionary scenario and make them observable within laboratory time scale.

Taken together, our findings show that transporters can have promiscuous weak activities, which are often hidden in the natural context of the cell. Furthermore, transporters can evolve novel substrates through generalist intermediates, either by increasing a weak activity or by establishing a new one. This coincides with what has been shown for soluble enzymes, and thus we show that membrane transporters can follow the same evolutionary logic (Bergthorsson et al., 2007; Hughes, 1994; Aharoni et al., 2005; Des Marais and Rausher, 2008; Näsvall et al., 2012). Looking at single adaptive mutations, we quantify the fitness trade-off between novel and original substrates, indicating that each mutation impacts the original substrates differently. Some of the identified adaptive mutations are almost neutral toward the original substrates, indicating that these mutations could become present in a large population, even before a new substrate is encountered by this population. This ‘standing genetic variation’ is a hallmark of contemporary theories of genetic innovation, and can explain how new functionalities arise while the protein is still under purifying selection for its original function (Soskine and Tawfik, 2010). Furthermore, such standing variation primes organisms for rapid adaptation (Barrett and Schluter, 2008). Along the same vein, we show that five wild-type amino acid transporters have inherent side reactivities. In the light of evolution, this is a feature rather than a bug, as it can prime these transporters to enhance the side reactivities either by mutation or duplication, as observed for soluble enzymes (Copley et al., 2023; Khersonsky and Tawfik, 2010). Overall, we show how transporter evolution could help organisms to explore ecological niches where novel transport activities give a fitness benefit, while still retaining the ability to thrive in the old environment.

Code and raw data are available on Zenodo (https://zenodo.org/records/10928101).

1. 1. Karapanagioti F
   2. Atlason UA
   3. Slotboom DJ
   4. Poolman B
   5. Obermaier S

   (2024) Zenodo

   Fitness landscape of substrate-adaptive mutations in evolved APC transporters.

   https://doi.org/10.5281/zenodo.8426062

1. 1. Adler J
   2. Bibi E

   (2005) Promiscuity in the geometry of electrostatic interactions between the Escherichia coli multidrug resistance transporter MdfA and cationic substrates

   The Journal of Biological Chemistry 280:2721–2729.

   https://doi.org/10.1074/jbc.M412332200

   • PubMed
   • Google Scholar
2. 1. Aharoni A
   2. Gaidukov L
   3. Khersonsky O
   4. McQ Gould S
   5. Roodveldt C
   6. Tawfik DS

   (2005) The “evolvability” of promiscuous protein functions

   Nature Genetics 37:73–76.

   https://doi.org/10.1038/ng1482

   • PubMed
   • Google Scholar
3. 1. Andorfer MC
   2. Grob JE
   3. Hajdin CE
   4. Chael JR
   5. Siuti P
   6. Lilly J
   7. Tan KL
   8. Lewis JC

   (2017) Understanding flavin-dependent halogenase reactivity via substrate activity profiling

   ACS Catalysis 7:1897–1904.

   https://doi.org/10.1021/acscatal.6b02707

   • PubMed
   • Google Scholar
4. 1. Bali AP
   2. Genee HJ
   3. Sommer MOA

   (2018) Directed evolution of membrane transport using synthetic selections

   ACS Synthetic Biology 7:789–793.

   https://doi.org/10.1021/acssynbio.7b00407

   • PubMed
   • Google Scholar
5. 1. Baranyi J
   2. Roberts TA

   (1995) Mathematics of predictive food microbiology

   International Journal of Food Microbiology 26:199–218.

   https://doi.org/10.1016/0168-1605(94)00121-l

   • PubMed
   • Google Scholar
6. 1. Barrett RDH
   2. Schluter D

   (2008) Adaptation from standing genetic variation

   Trends in Ecology & Evolution 23:38–44.

   https://doi.org/10.1016/j.tree.2007.09.008

   • PubMed
   • Google Scholar
7. 1. Bergthorsson U
   2. Andersson DI
   3. Roth JR

   (2007) Ohno’s dilemma: evolution of new genes under continuous selection

   PNAS 104:17004–17009.

   https://doi.org/10.1073/pnas.0707158104

   • PubMed
   • Google Scholar
8. 1. Besnard J
   2. Pratelli R
   3. Zhao C
   4. Sonawala U
   5. Collakova E
   6. Pilot G
   7. Okumoto S

   (2016) UMAMIT14 is an amino acid exporter involved in phloem unloading in Arabidopsis roots

   Journal of Experimental Botany 67:6385–6397.

   https://doi.org/10.1093/jxb/erw412

   • PubMed
   • Google Scholar
9. 1. Bianchi F
   2. Van’t Klooster JS
   3. Ruiz SJ
   4. Poolman B

   (2019) Regulation of amino acid transport in Saccharomyces cerevisiae

   Microbiology and Molecular Biology Reviews 83:e00024-19.

   https://doi.org/10.1128/MMBR.00024-19

   • PubMed
   • Google Scholar
10. 1. Caligiore F
    2. Zangelmi E
    3. Vetro C
    4. Kentache T
    5. Dewulf JP
    6. Veiga-da-Cunha M
    7. Van Schaftingen E
    8. Bommer G
    9. Peracchi A

    (2022) Human cytosolic transaminases: side activities and patterns of discrimination towards physiologically available alternative substrates

    Cellular and Molecular Life Sciences 79:421.

    https://doi.org/10.1007/s00018-022-04439-3

    • PubMed
    • Google Scholar
11. 1. Chan KM
    2. Liu YT
    3. Ma CH
    4. Jayaram M
    5. Sau S

    (2013) The 2 micron plasmid of Saccharomyces cerevisiae: a miniaturized selfish genome with optimized functional competence

    Plasmid 70:2–17.

    https://doi.org/10.1016/j.plasmid.2013.03.001

    • PubMed
    • Google Scholar
12. 1. Copley SD
    2. Newton MS
    3. Widney KA

    (2023) How to recruit a promiscuous enzyme to serve a new function

    Biochemistry 62:300–308.

    https://doi.org/10.1021/acs.biochem.2c00249

    • PubMed
    • Google Scholar
13. 1. Des Marais DL
    2. Rausher MD

    (2008) Escape from adaptive conflict after duplication in an anthocyanin pathway gene

    Nature 454:762–765.

    https://doi.org/10.1038/nature07092

    • PubMed
    • Google Scholar
14. 1. Didion T
    2. Grauslund M
    3. Kielland-Brandt MC
    4. Andersen HA

    (1996) Amino acids induce expression of BAP2, a branched-chain amino acid permease gene in Saccharomyces cerevisiae

    Journal of Bacteriology 178:2025–2029.

    https://doi.org/10.1128/jb.178.7.2025-2029.1996

    • PubMed
    • Google Scholar
15. 1. Dunnett CW

    (1955) A multiple comparison procedure for comparing several treatments with a control

    Journal of the American Statistical Association 50:1096–1121.

    https://doi.org/10.1080/01621459.1955.10501294

    • Google Scholar
16. 1. Edwards N
    2. Anderson CMH
    3. Conlon NJ
    4. Watson AK
    5. Hall RJ
    6. Cheek TR
    7. Embley TM
    8. Thwaites DT

    (2018) Resculpting the binding pocket of APC superfamily LeuT-fold amino acid transporters

    Cellular and Molecular Life Sciences 75:921–938.

    https://doi.org/10.1007/s00018-017-2677-8

    • PubMed
    • Google Scholar
17. 1. García-García JD
    2. Joshi J
    3. Patterson JA
    4. Trujillo-Rodriguez L
    5. Reisch CR
    6. Javanpour AA
    7. Liu CC
    8. Hanson AD

    (2020) Potential for applying continuous directed evolution to plant enzymes: an exploratory study

    Life 10:179.

    https://doi.org/10.3390/life10090179

    • PubMed
    • Google Scholar
18. 1. Geertsma ER
    2. Dutzler R

    (2011) A versatile and efficient high-throughput cloning tool for structural biology

    Biochemistry 50:3272–3278.

    https://doi.org/10.1021/bi200178z

    • PubMed
    • Google Scholar
19. 1. Ghaddar K
    2. Krammer EM
    3. Mihajlovic N
    4. Brohée S
    5. André B
    6. Prévost M

    (2014) Converting the yeast arginine can1 permease to a lysine permease

    The Journal of Biological Chemistry 289:7232–7246.

    https://doi.org/10.1074/jbc.M113.525915

    • PubMed
    • Google Scholar
20. 1. Gibrat JF
    2. Madej T
    3. Bryant SH

    (1996) Surprising similarities in structure comparison

    Current Opinion in Structural Biology 6:377–385.

    https://doi.org/10.1016/s0959-440x(96)80058-3

    • PubMed
    • Google Scholar
21. 1. Gibson DG
    2. Young L
    3. Chuang RY
    4. Venter JC
    5. Hutchison CA
    6. Smith HO

    (2009) Enzymatic assembly of DNA molecules up to several hundred kilobases

    Nature Methods 6:343–345.

    https://doi.org/10.1038/nmeth.1318

    • PubMed
    • Google Scholar
22. 1. Gietz RD
    2. Schiestl RH

    (2007) High-efficiency yeast transformation using the LiAc/SS carrier DNA/PEG method

    Nature Protocols 2:31–34.

    https://doi.org/10.1038/nprot.2007.13

    • PubMed
    • Google Scholar
23. Book

    1. Girden ER

    (1992) ANOVA: Repeated Measures

    Sage Publications, Newbury Park.

    https://doi.org/10.4135/9781412983419

    • Google Scholar
24. 1. Gournas C
    2. Evangelidis T
    3. Athanasopoulos A
    4. Mikros E
    5. Sophianopoulou V

    (2015) The Aspergillus nidulans proline permease as a model for understanding the factors determining substrate binding and specificity of fungal amino acid transporters

    The Journal of Biological Chemistry 290:6141–6155.

    https://doi.org/10.1074/jbc.M114.612069

    • PubMed
    • Google Scholar
25. 1. Gros Y
    2. Schuldiner S

    (2010) Directed evolution reveals hidden properties of VMAT, a neurotransmitter transporter

    The Journal of Biological Chemistry 285:5076–5084.

    https://doi.org/10.1074/jbc.M109.081216

    • PubMed
    • Google Scholar
26. 1. Huang L

    (2008) Growth kinetics of Listeria monocytogenes in broth and beef frankfurters--determination of lag phase duration and exponential growth rate under isothermal conditions

    Journal of Food Science 73:E235–E242.

    https://doi.org/10.1111/j.1750-3841.2008.00785.x

    • PubMed
    • Google Scholar
27. 1. Huang L

    (2011) A new mechanistic growth model for simultaneous determination of lag phase duration and exponential growth rate and a new Bĕlehdrádek-type model for evaluating the effect of temperature on growth rate

    Food Microbiology 28:770–776.

    https://doi.org/10.1016/j.fm.2010.05.019

    • PubMed
    • Google Scholar
28. 1. Hughes AL

    (1994) The evolution of functionally novel proteins after gene duplication

    Proceedings. Biological Sciences 256:119–124.

    https://doi.org/10.1098/rspb.1994.0058

    • PubMed
    • Google Scholar
29. 1. Jack DL
    2. Paulsen IT
    3. Saier MH

    (2000) The amino acid/polyamine/organocation (APC) superfamily of transporters specific for amino acids, polyamines and organocations

    Microbiology 146:1797–1814.

    https://doi.org/10.1099/00221287-146-8-1797

    • PubMed
    • Google Scholar
30. 1. Janssen DB
    2. Dinkla IJT
    3. Poelarends GJ
    4. Terpstra P

    (2005) Bacterial degradation of xenobiotic compounds: evolution and distribution of novel enzyme activities

    Environmental Microbiology 7:1868–1882.

    https://doi.org/10.1111/j.1462-2920.2005.00966.x

    • PubMed
    • Google Scholar
31. 1. Jensen RA

    (1976) Enzyme recruitment in evolution of new function

    Annual Review of Microbiology 30:409–425.

    https://doi.org/10.1146/annurev.mi.30.100176.002205

    • PubMed
    • Google Scholar
32. 1. Jørgensen ME
    2. Xu D
    3. Crocoll C
    4. Ernst HA
    5. Ramírez D
    6. Motawia MS
    7. Olsen CE
    8. Mirza O
    9. Nour-Eldin HH
    10. Halkier BA

    (2017) Origin and evolution of transporter substrate specificity within the NPF family

    eLife 6:e19466.

    https://doi.org/10.7554/eLife.19466

    • PubMed
    • Google Scholar
33. 1. Jumper J
    2. Evans R
    3. Pritzel A
    4. Green T
    5. Figurnov M
    6. Ronneberger O
    7. Tunyasuvunakool K
    8. Bates R
    9. Žídek A
    10. Potapenko A
    11. Bridgland A
    12. Meyer C
    13. Kohl SAA
    14. Ballard AJ
    15. Cowie A
    16. Romera-Paredes B
    17. Nikolov S
    18. Jain R
    19. Adler J
    20. Back T
    21. Petersen S
    22. Reiman D
    23. Clancy E
    24. Zielinski M
    25. Steinegger M
    26. Pacholska M
    27. Berghammer T
    28. Bodenstein S
    29. Silver D
    30. Vinyals O
    31. Senior AW
    32. Kavukcuoglu K
    33. Kohli P
    34. Hassabis D

    (2021) Highly accurate protein structure prediction with AlphaFold

    Nature 596:583–589.

    https://doi.org/10.1038/s41586-021-03819-2

    • PubMed
    • Google Scholar
34. 1. Kazmier K
    2. Claxton DP
    3. Mchaourab HS

    (2017) Alternating access mechanisms of LeuT-fold transporters: trailblazing towards the promised energy landscapes

    Current Opinion in Structural Biology 45:100–108.

    https://doi.org/10.1016/j.sbi.2016.12.006

    • PubMed
    • Google Scholar
35. 1. Khersonsky O
    2. Tawfik DS

    (2010) Enzyme promiscuity: a mechanistic and evolutionary perspective

    Annual Review of Biochemistry 79:471–505.

    https://doi.org/10.1146/annurev-biochem-030409-143718

    • PubMed
    • Google Scholar
36. 1. Kunkel TA

    (1985) Rapid and efficient site-specific mutagenesis without phenotypic selection (M13 cloning vectors/silent mutations/in vitro mutagenesis/synthetic oligonucleotides/uracil-containing DNA)

    PNAS 82:488–492.

    https://doi.org/10.1073/pnas.82.2.488

    • Google Scholar
37. 1. Kyte J
    2. Doolittle RF

    (1982) A simple method for displaying the hydropathic character of A protein

    Journal of Molecular Biology 157:105–132.

    https://doi.org/10.1016/0022-2836(82)90515-0

    • PubMed
    • Google Scholar
38. 1. Long M
    2. VanKuren NW
    3. Chen S
    4. Vibranovski MD

    (2013) New gene evolution: little did we know

    Annual Review of Genetics 47:307–333.

    https://doi.org/10.1146/annurev-genet-111212-133301

    • PubMed
    • Google Scholar
39. 1. Lynch M
    2. Force A

    (2000) The probability of duplicate gene preservation by subfunctionalization

    Genetics 154:459–473.

    https://doi.org/10.1093/genetics/154.1.459

    • PubMed
    • Google Scholar
40. 1. Madeira F
    2. Pearce M
    3. Tivey ARN
    4. Basutkar P
    5. Lee J
    6. Edbali O
    7. Madhusoodanan N
    8. Kolesnikov A
    9. Lopez R

    (2022) Search and sequence analysis tools services from EMBL-EBI in 2022

    Nucleic Acids Research 50:W276–W279.

    https://doi.org/10.1093/nar/gkac240

    • PubMed
    • Google Scholar
41. 1. Magasanik B
    2. Kaiser CA

    (2002) Nitrogen regulation in Saccharomyces cerevisiae

    Gene 290:1–18.

    https://doi.org/10.1016/s0378-1119(02)00558-9

    • PubMed
    • Google Scholar
42. 1. Mascotti ML
    2. Kumar H
    3. Nguyen QT
    4. Ayub MJ
    5. Fraaije MW

    (2018) Reconstructing the evolutionary history of F₄₂₀-dependent dehydrogenases

    Scientific Reports 8:17571.

    https://doi.org/10.1038/s41598-018-35590-2

    • PubMed
    • Google Scholar
43. 1. Meier G
    2. Thavarasah S
    3. Ehrenbolger K
    4. Hutter CAJ
    5. Hürlimann LM
    6. Barandun J
    7. Seeger MA

    (2023) Deep mutational scan of a drug efflux pump reveals its structure-function landscape

    Nature Chemical Biology 19:440–450.

    https://doi.org/10.1038/s41589-022-01205-1

    • PubMed
    • Google Scholar
44. 1. Muller HJ

    (1936) Bar duplication

    Science 83:528–530.

    https://doi.org/10.1126/science.83.2161.528-a

    • PubMed
    • Google Scholar
45. 1. Näsvall J
    2. Sun L
    3. Roth JR
    4. Andersson DI

    (2012) Real-time evolution of new genes by innovation, amplification, and divergence

    Science 338:384–387.

    https://doi.org/10.1126/science.1226521

    • PubMed
    • Google Scholar
46. Book

    1. Ohno S

    (1970) Evolution by Gene Duplication

    Berlin, Heidelberg: Springer.

    https://doi.org/10.1007/978-3-642-86659-3

    • Google Scholar
47. 1. Perland E
    2. Fredriksson R

    (2017) Classification systems of secondary active transporters

    Trends in Pharmacological Sciences 38:305–315.

    https://doi.org/10.1016/j.tips.2016.11.008

    • PubMed
    • Google Scholar
48. 1. Pettersen EF
    2. Goddard TD
    3. Huang CC
    4. Meng EC
    5. Couch GS
    6. Croll TI
    7. Morris JH
    8. Ferrin TE

    (2021) UCSF ChimeraX: Structure visualization for researchers, educators, and developers

    Protein Science 30:70–82.

    https://doi.org/10.1002/pro.3943

    • PubMed
    • Google Scholar
49. 1. Ravikumar A
    2. Arrieta A
    3. Liu CC

    (2014) An orthogonal DNA replication system in yeast

    Nature Chemical Biology 10:175–177.

    https://doi.org/10.1038/nchembio.1439

    • PubMed
    • Google Scholar
50. 1. Ravikumar A
    2. Arzumanyan GA
    3. Obadi MKA
    4. Javanpour AA
    5. Liu CC

    (2018) Scalable, continuous evolution of genes at mutation rates above genomic error thresholds

    Cell 175:1946–1957.

    https://doi.org/10.1016/j.cell.2018.10.021

    • PubMed
    • Google Scholar
51. 1. Regenberg B
    2. Kielland-Brandt MC

    (2001) Amino acid residues important for substrate specificity of the amino acid permeases Can1p and Gnp1p in Saccharomyces cerevisiae

    Yeast 18:1429–1440.

    https://doi.org/10.1002/yea.792

    • PubMed
    • Google Scholar
52. 1. Risso VA
    2. Gavira JA
    3. Mejia-Carmona DF
    4. Gaucher EA
    5. Sanchez-Ruiz JM

    (2013) Hyperstability and substrate promiscuity in laboratory resurrections of Precambrian β-lactamases

    Journal of the American Chemical Society 135:2899–2902.

    https://doi.org/10.1021/ja311630a

    • PubMed
    • Google Scholar
53. 1. Saier MH
    2. Reddy VS
    3. Tamang DG
    4. Västermark A

    (2014) The transporter classification database

    Nucleic Acids Research 42:D251–D258.

    https://doi.org/10.1093/nar/gkt1097

    • PubMed
    • Google Scholar
54. 1. Saier MH
    2. Reddy VS
    3. Moreno-Hagelsieb G
    4. Hendargo KJ
    5. Zhang Y
    6. Iddamsetty V
    7. Lam KJK
    8. Tian N
    9. Russum S
    10. Wang J
    11. Medrano-Soto A

    (2021) The transporter classification database (TCDB): 2021 update

    Nucleic Acids Research 49:D461–D467.

    https://doi.org/10.1093/nar/gkaa1004

    • PubMed
    • Google Scholar
55. 1. Sanger F
    2. Nicklen S
    3. Coulson AR

    (1977) DNA sequencing with chain-terminating inhibitors

    PNAS 74:5463–5467.

    https://doi.org/10.1073/pnas.74.12.5463

    • PubMed
    • Google Scholar
56. 1. Schindelin J
    2. Arganda-Carreras I
    3. Frise E
    4. Kaynig V
    5. Longair M
    6. Pietzsch T
    7. Preibisch S
    8. Rueden C
    9. Saalfeld S
    10. Schmid B
    11. Tinevez JY
    12. White DJ
    13. Hartenstein V
    14. Eliceiri K
    15. Tomancak P
    16. Cardona A

    (2012) Fiji: an open-source platform for biological-image analysis

    Nature Methods 9:676–682.

    https://doi.org/10.1038/nmeth.2019

    • PubMed
    • Google Scholar
57. 1. Shimazu M
    2. Itaya T
    3. Pongcharoen P
    4. Sekito T
    5. Kawano-Kawada M
    6. Kakinuma Y

    (2012) Vba5p, a novel plasma membrane protein involved in amino acid uptake and drug sensitivity in Saccharomyces cerevisiae

    Bioscience, Biotechnology, and Biochemistry 76:1993–1995.

    https://doi.org/10.1271/bbb.120455

    • PubMed
    • Google Scholar
58. 1. Singh SK
    2. Piscitelli CL
    3. Yamashita A
    4. Gouaux E

    (2008) A competitive inhibitor traps LeuT in an open-to-out conformation

    Science 322:1655–1661.

    https://doi.org/10.1126/science.1166777

    • PubMed
    • Google Scholar
59. 1. Soskine M
    2. Tawfik DS

    (2010) Mutational effects and the evolution of new protein functions

    Nature Reviews. Genetics 11:572–582.

    https://doi.org/10.1038/nrg2808

    • PubMed
    • Google Scholar
60. 1. Student

    (1908) The probable error of a mean

    Biometrika 6:1.

    https://doi.org/10.2307/2331554

    • Google Scholar
61. 1. Tutulan-Cunita AC
    2. Mikoshi M
    3. Mizunuma M
    4. Hirata D
    5. Miyakawa T

    (2005) Mutational analysis of the yeast multidrug resistance ABC transporter Pdr5p with altered drug specificity

    Genes to Cells 10:409–420.

    https://doi.org/10.1111/j.1365-2443.2005.00847.x

    • PubMed
    • Google Scholar
62. 1. UniProt Consortium

    (2023) UniProt: the Universal Protein Knowledgebase in 2023

    Nucleic Acids Research 51:D523–D531.

    https://doi.org/10.1093/nar/gkac1052

    • PubMed
    • Google Scholar
63. 1. Vastermark A
    2. Wollwage S
    3. Houle ME
    4. Rio R
    5. Saier MH

    (2014) Expansion of the APC superfamily of secondary carriers

    Proteins 82:2797–2811.

    https://doi.org/10.1002/prot.24643

    • PubMed
    • Google Scholar
64. 1. Yadav AK
    2. Bachhawat AK

    (2011) CgCYN1, a plasma membrane cystine-specific transporter of Candida glabrata with orthologues prevalent among pathogenic yeast and fungi

    The Journal of Biological Chemistry 286:19714–19723.

    https://doi.org/10.1074/jbc.M111.240648

    • PubMed
    • Google Scholar
65. 1. Yamashita A
    2. Singh SK
    3. Kawate T
    4. Jin Y
    5. Gouaux E

    (2005) Crystal structure of a bacterial homologue of Na+/Cl--dependent neurotransmitter transporters

    Nature 437:215–223.

    https://doi.org/10.1038/nature03978

    • PubMed
    • Google Scholar
66. 1. Young EM
    2. Tong A
    3. Bui H
    4. Spofford C
    5. Alper HS

    (2014) Rewiring yeast sugar transporter preference through modifying a conserved protein motif

    PNAS 111:131–136.

    https://doi.org/10.1073/pnas.1311970111

    • PubMed
    • Google Scholar
67. 1. Zamyatnin AA

    (1972) Protein volume in solution

    Progress in Biophysics and Molecular Biology 24:107–123.

    https://doi.org/10.1016/0079-6107(72)90005-3

    • PubMed
    • Google Scholar
68. 1. Zheng J
    2. Guo N
    3. Wagner A

    (2020) Selection enhances protein evolvability by increasing mutational robustness and foldability

    Science 370:eabb5962.

    https://doi.org/10.1126/science.abb5962

    • PubMed
    • Google Scholar

Author details

1. Foteini Karapanagioti

   Department of Biochemistry, University of Groningen, Groningen, Netherlands

   Contribution

   Conceptualization, Data curation, Formal analysis, Validation, Investigation, Visualization, Methodology, Writing – original draft, Project administration, Writing – review and editing

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0009-0007-1700-4867

2. Úlfur Águst Atlason

   Department of Biochemistry, University of Groningen, Groningen, Netherlands

   Contribution

   Investigation, Methodology

   Competing interests

   No competing interests declared

3. Dirk J Slotboom

   Department of Biochemistry, University of Groningen, Groningen, Netherlands

   Contribution

   Supervision, Writing – original draft, Writing – review and editing

   For correspondence

   d.j.slotboom@rug.nl

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0002-5804-9689

4. Bert Poolman

   Department of Biochemistry, University of Groningen, Groningen, Netherlands

   Contribution

   Resources, Supervision, Funding acquisition, Writing – original draft, Project administration, Writing – review and editing

   For correspondence

   b.poolman@rug.nl

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0002-1455-531X

5. Sebastian Obermaier

   Department of Biochemistry, University of Groningen, Groningen, Netherlands

   Contribution

   Conceptualization, Data curation, Supervision, Validation, Investigation, Visualization, Methodology, Writing – original draft, Project administration, Writing – review and editing

   For correspondence

   s.obermaier@protonmail.com

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0002-1238-4502

• 809

  views

• 58

  downloads

• 1

  citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.