Selecting for Altered Substrate Specificity Reveals the Evolutionary Flexibility of ATP-Binding Cassette Transporters

Specialization restricts the evolutionary paths available to yeast sugar transporters

Functional innovation at the protein level is a key source of evolutionary novelties. The constraints on functional innovations are likely to be highly specific in different proteins, which are shaped by their unique histories and the extent of global epistasis that arises from their structures and biochemistries. These contextual nuances in the sequence-function relationship have implications both for a basic understanding of the evolutionary process and for engineering proteins with desirable properties. Here, we have investigated the molecular basis of novel function in a model member of an ancient, conserved, and biotechnologically relevant protein family. These Major Facilitator Superfamily sugar porters are a functionally diverse group of proteins that are thought to be highly plastic and evolvable. By dissecting a recent evolutionary innovation in an α-glucoside transporter from the yeast Saccharomyces eubayanus, we show that the ability to transport a novel substrate requires high-order interactions between many protein regions and numerous specific residues proximal to the transport channel. To reconcile the functional diversity of this family with the constrained evolution of this model protein, we generated new, state-of-the-art genome annotations for 332 Saccharomycotina yeast species spanning approximately 400 million years of evolution. By integrating phylogenetic and phenotypic analyses across these species, we show that the model yeast α-glucoside transporters likely evolved from a multifunctional ancestor and became subfunctionalized. The accumulation of additive and epistatic substitutions likely entrenched this subfunction, which made the simultaneous acquisition of multiple interacting substitutions the only reasonably accessible path to novelty.

A phage-encoded RNA-binding protein inhibits the antiviral activity of a toxin-antitoxin system

Bacteria harbor diverse mechanisms to defend themselves against their viral predators, bacteriophages. In response, phages can evolve counter-defense systems, most of which are poorly understood. In T4-like phages, the gene tifA prevents bacterial defense by the type III toxin-antitoxin (TA) system toxIN, but the mechanism by which TifA inhibits ToxIN remains unclear. Here, we show that TifA directly binds both the endoribonuclease ToxN and RNA, leading to the formation of a high molecular weight ribonucleoprotein complex in which ToxN is inhibited. The RNA binding activity of TifA is necessary for its interaction with and inhibition of ToxN. Thus, we propose that TifA inhibits ToxN during phage infection by trapping ToxN on cellular RNA, particularly the abundant 16S rRNA, thereby preventing cleavage of phage transcripts. Taken together, our results reveal a novel mechanism underlying inhibition of a phage-defensive RNase toxin by a small, phage-encoded protein.

Advances in ligand-specific biosensing for structurally similar molecules

The specificity of biological systems makes it possible to develop biosensors targeting specific metabolites, toxins, and pollutants in complex medical or environmental samples without interference from structurally similar compounds. For the last two decades, great efforts have been devoted to creating proteins or nucleic acids with novel properties through synthetic biology strategies. Beyond augmenting biocatalytic activity, expanding target substrate scopes, and enhancing enzymes' enantioselectivity and stability, an increasing research area is the enhancement of molecular specificity for genetically encoded biosensors. Here, we summarize recent advances in the development of highly specific biosensor systems and their essential applications. First, we describe the rational design principles required to create libraries containing potential mutants with less promiscuity or better specificity. Next, we review the emerging high-throughput screening techniques to engineer biosensing specificity for the desired target. Finally, we examine the computer-aided evaluation and prediction methods to facilitate the construction of ligand-specific biosensors.

Selection of a de novo gene that can promote survival of Escherichia coli by modulating protein homeostasis pathways

Cellular novelty can emerge when non-functional loci become functional genes in a process termed de novo gene birth. But how proteins with random amino acid sequences beneficially integrate into existing cellular pathways remains poorly understood. We screened ~108 genes, generated from random nucleotide sequences and devoid of homology to natural genes, for their ability to rescue growth arrest of Escherichia coli cells producing the ribonuclease toxin MazF. We identified ~2,000 genes that could promote growth, probably by reducing transcription from the promoter driving toxin expression. Additionally, one random protein, named Random antitoxin of MazF (RamF), modulated protein homeostasis by interacting with chaperones, leading to MazF proteolysis and a consequent loss of its toxicity. Finally, we demonstrate that random proteins can improve during evolution by identifying beneficial mutations that turned RamF into a more efficient inhibitor. Our work provides a mechanistic basis for how de novo gene birth can produce functional proteins that effectively benefit cells evolving under stress.

Extending the reach of homology by using successive computational filters to find yeast pheromone genes

The mating of fungi depends on pheromones that mediate communication between two mating types. Most species use short peptides as pheromones, which are either unmodified (e.g., α-factor in Saccharomyces cerevisiae) or C-terminally farnesylated (e.g., a-factor in S. cerevisiae). Peptide pheromones have been found by genetics or biochemistry in a small number of fungi, but their short sequences and modest conservation make it impossible to detect homologous sequences in most species. To overcome this problem, we used a four-step computational pipeline to identify candidate a-factor genes in sequenced genomes of the Saccharomycotina, the fungal clade that contains most of the yeasts: we require that candidate genes have a C-terminal prenylation motif, are shorter than 100 amino acids long, and contain a proteolytic-processing motif upstream of the potential mature pheromone sequence and that closely related species contain highly conserved homologs of the potential mature pheromone sequence. Additional manual curation exploits the observation that many species carry more than one a-factor gene, encoding identical or nearly identical pheromones. From 332 Saccharomycotina genomes, we identified strong candidate pheromone genes in 241 genomes, covering 13 clades that are each separated from each other by at least 100 million years, the time required for evolution to remove detectable sequence homology among small pheromone genes. For one small clade, the Yarrowia, we demonstrated that our algorithm found the a-factor genes: deleting all four related genes in the a-mating type of Yarrowia lipolytica prevents mating.

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

Drug efflux is a common resistance mechanism found in bacteria and cancer cells, but studies providing comprehensive functional insights are scarce. In this study, we performed deep mutational scanning (DMS) on the bacterial ABC transporter EfrCD to determine the drug efflux activity profile of more than 1,430 single variants. These systematic measurements revealed that the introduction of negative charges at different locations within the large substrate binding pocket results in strongly increased efflux activity toward positively charged ethidium, whereas additional aromatic residues did not display the same effect. Data analysis in the context of an inward-facing cryogenic electron microscopy structure of EfrCD uncovered a high-affinity binding site, which releases bound drugs through a peristaltic transport mechanism as the transporter transits to its outward-facing conformation. Finally, we identified substitutions resulting in rapid Hoechst influx without affecting the efflux activity for ethidium and daunorubicin. Hence, single mutations can convert EfrCD into a drug-specific ABC importer.

On-target inhibition of Cryptosporidium parvum by nitazoxanide (NTZ) and paclitaxel (PTX) validated using a novel MDR1-transgenic host cell model and algorithms to quantify the effect on the parasite target

Cryptosporidium parvum is a globally distributed zoonotic protozoan parasite that causes moderate to severe, sometime deadly, watery diarrhea in humans and animals, for which fully effective treatments are yet unavailable. In studying the mechanism of action of drugs against intracellular pathogens, it is important to validate whether the observed anti-infective activity is attributed to the drug action on the pathogen or host target. For the epicellular parasite Cryptosporidium, we have previously developed a concept that the host cells with significantly increased drug tolerance by transient overexpression of the multidrug resistance protein-1 (MDR1) could be utilized to evaluate whether and how much the observed anti-cryptosporidial activity of an inhibitor was attributed to the inhibitor’s action on the parasite target. However, the transient transfection model was only applicable to evaluating native MDR1 substrates. Here we report an advanced model using stable MDR1-transgenic HCT-8 cells that allows rapid development of novel resistance to non-MDR1 substrates by multiple rounds of drug selection. Using the new model, we successfully validated that nitazoxanide, a non-MDR1 substrate and the only FDA-approved drug to treat human cryptosporidiosis, killed C. parvum by fully (100%) acting on the parasite target. We also confirmed that paclitaxel acted fully on the parasite target, while several other inhibitors including mitoxantrone, doxorubicin, vincristine and ivermectin acted partially on the parasite targets. Additionally, we developed mathematical models to quantify the proportional contribution of the on-parasite-target effect to the observed anti-cryptosporidial activity and to evaluate the relationships between several in vitro parameters, including antiparasitic efficacy (ECi), cytotoxicity (TCi), selectivity index (SI) and Hill slope (h). Owning to the promiscuity of the MDR1 efflux pump, the MDR1-transgenic host cell model could be applied to assess the on-parasite-target effects of newly identified hits/leads, either substrates or non-substrates of MDR1, against Cryptosporidium or other epicellular pathogens.

Prevalence of pvmrp1 Polymorphisms and Its Contribution to Antimalarial Response

As more sporadic cases of chloroquine resistance occur (CQR) in Plasmodium vivax (P. vivax) malaria, molecular markers have become an important tool to monitor the introduction and spread of drug resistance. P. vivax multidrug resistance-associated protein 1 (PvMRP1), as one of the members of the ATP-binding cassette (ABC) transporters, may modulate this phenotype. In this study, we investigated the gene mutations and copy number variations (CNVs) in the pvmrp1 in 102 P. vivax isolates from China, the Republic of Korea (ROK), Myanmar, Papua New Guinea (PNG), Pakistan, the Democratic People’s Republic of Korea (PRK), and Cambodia. And we also obtained 72 available global pvmrp1 sequences deposited in the PlasmoDB database to investigate the genetic diversity, haplotype diversity, natural selection, and population structure of pvmrp1. In total, 29 single nucleotide polymorphisms reflected in 23 non-synonymous, five synonymous mutations and one gene deletion were identified, and CNVs were found in 2.9% of the isolates. Combined with the antimalarial drug susceptibility observed in the previous in vitro assays, except the prevalence of S354N between the two CQ sensitivity categories revealed a significant difference, no genetic mutations or CNVs associated with drug sensitivity were found. The genetic polymorphism analysis of 166 isolates worldwide found that the overall nucleotide diversity (π) of pvmrp1 was 0.0011, with 46 haplotypes identified (Hd = 0.9290). The ratio of non-synonymous to synonymous mutations (dn/ds = 0.5536) and the neutrality tests statistic Fu and Li’s D* test (Fu and Li’s D* = −3.9871, p < 0.02) suggests that pvmrp1 had evolved under a purifying selection. Due to geographical differences, genetic differentiation levels of pvmrp1 in different regions were different to some extent. Overall, this study provides a new idea for finding CQR molecular monitoring of P. vivax and provides more sequences of pvmrp1 in Asia for subsequent research. However, further validation is still needed through laboratory and epidemiological field studies of P. vivax samples from more regions.

An Identification and Expression Analysis of the ABCG Genes Related to Benzaldehyde Transportation among Three Prunus Species

The plants of Prunus mostly bloom in early spring, and the flowers of various species possess their individual floral scent characteristics; Prunus mume, especially, can volatilize a large amount of benzenoid compounds into the air during the flowering phase. In order to elucidate the molecular basis of the differences in the volatile capacity of aromatic substances among Prunus flowers, the endogenous and the headspace volatile components and the expression of ABCG genes were studied among P. mume, P. armeniaca, and P. persica. We detected the floral components in the three species by gas chromatography-mass spectrometry (GC-MS), and we found that benzaldehyde was the key component. Meanwhile, the volatilization efficiency of benzaldehyde in P. mume and P. armeniaca were much higher than that in P. persica. Furthermore, 130, 135, and 133 ABC family members from P. mume, P. armeniaca, and P. persica were identified, respectively. WGCNA analysis demonstrated that candidate ABCG genes were positively correlated with benzaldehyde volatilization efficiency. Moreover, quantitative Real-time PCR indicated that ABCG17 was more likely to be involved in the transmembrane transport of benzaldehyde. This study aimed to provide a theoretical basis for elucidating the transmembrane transport of benzaldehyde and to further the valuable information for fragrant flower breeding in Prunus.

Characterization and tissue localization of zebrafish homologs of the human ABCB1 multidrug transporter

Capillary endothelial cells of the human blood–brain barrier (BBB) express high levels of P-glycoprotein (P-gp, encoded by ABCB1) and ABCG2 (encoded by ABCG2). However, little information is available regarding ATP-binding cassette transporters expressed at the zebrafish BBB, which has emerged as a potential model system. We report the characterization and tissue localization of two genes that are similar to ABCB1, zebrafish abcb4 and abcb5. When stably expressed in HEK293 cells, both Abcb4 and Abcb5 conferred resistance to P-gp substrates; however, Abcb5 poorly transported doxorubicin and mitoxantrone compared to zebrafish Abcb4. Additionally, Abcb5 did not transport the fluorescent P-gp probes BODIPY-ethylenediamine or LDS 751, while they were transported by Abcb4. High-throughput screening of 90 human P-gp substrates confirmed that Abcb4 has an overlapping substrate specificity profile with P-gp. In the brain vasculature, RNAscope probes for abcb4 colocalized with staining by the P-gp antibody C219, while abcb5 was not detected. The abcb4 probe also colocalized with claudin-5 in brain endothelial cells. Abcb4 and Abcb5 had different tissue localizations in multiple zebrafish tissues, potentially indicating different functions. The data suggest that zebrafish Abcb4 functionally phenocopies P-gp and that the zebrafish may serve as a model to study the role of P-gp at the BBB.

The molecular evolution of function in the CFTR chloride channel

The ATP-binding cassette (ABC) transporter superfamily includes many proteins of clinical relevance, with genes expressed in all domains of life. Although most members use the energy of ATP binding and hydrolysis to accomplish the active import or export of various substrates across membranes, the cystic fibrosis transmembrane conductance regulator (CFTR) is the only known animal ABC transporter that functions primarily as an ion channel. Defects in CFTR, which is closely related to ABCC subfamily members that bear function as bona fide transporters, underlie the lethal genetic disease cystic fibrosis. This article seeks to integrate structural, functional, and genomic data to begin to answer the critical question of how the function of CFTR evolved to exhibit regulated channel activity. We highlight several examples wherein preexisting features in ABCC transporters were functionally leveraged as is, or altered by molecular evolution, to ultimately support channel function. This includes features that may underlie (1) construction of an anionic channel pore from an anionic substrate transport pathway, (2) establishment and tuning of phosphoregulation, and (3) optimization of channel function by specialized ligand–channel interactions. We also discuss how divergence and conservation may help elucidate the pharmacology of important CFTR modulators.

CRISPR/Cas9 mutagenesis reveals a role for ABCB1 in gut immune responses to Vibrio diazotrophicus in sea urchin larvae

The ABC transporter ABCB1 plays an important role in the disposition of xenobiotics. Embryos of most species express high levels of this transporter in early development as a protective mechanism, but its native substrates are not known. Here, we used larvae of the sea urchin Strongylocentrotus purpuratus to characterize the early life expression and role of Sp-ABCB1a, a homolog of ABCB1. The results indicate that while Sp-ABCB1a is initially expressed ubiquitously, it becomes enriched in the developing gut. Using optimized CRISPR/Cas9 gene editing methods to achieve high editing efficiency in the F0 generation, we generated ABCB1a crispant embryos with significantly reduced transporter efflux activity. When infected with the opportunistic pathogen Vibrio diazotrophicus, Sp-ABCB1a crispant larvae demonstrated significantly stronger gut inflammation, immunocyte migration and cytokine Sp-IL-17 induction, as compared with infected control larvae. The results suggest an ancestral function of ABCB1 in host-microbial interactions, with implications for the survival of invertebrate larvae in the marine microbial environment.

Evolutionary history of ATP-binding cassette proteins

ATP-binding cassette (ABC) proteins are found in every sequenced genome and evolved deep in the phylogenetic tree of life. ABC proteins form one of the largest homologous protein families, with most being involved in substrate transport across biological membranes, and a few cytoplasmic members regulating in essential processes like translation. The predominant ABC protein classification scheme is derived from human members, but the increasing number of fully sequenced genomes permits to reevaluate this paradigm in the light of the evolutionary history the ABC-protein superfamily. As we study the diversity of substrates, mechanisms, and physiological roles of ABC proteins, knowledge of the evolutionary relationships highlights similarities and differences that can be attributed to specific branches in protein divergence. While alignments and trees built on natural sequence variation account for the evolutionary divergence of ABC proteins, high-throughput experiments and next-generation sequencing creating experimental sequence variation are instrumental in identifying functional constraints. The combination of natural and experimentally produced sequence variation allows a broader and more rational study of the function and physiological roles of ABC proteins.