Engineering orthogonal signalling pathways reveals the sparse occupancy of sequence space.

Gene duplication is a common and powerful mechanism by which cells create new signalling pathways1,2, but recently duplicated proteins typically must become insulated from each other and from other paralogues to prevent unwanted crosstalk3. A similar challenge arises when new sensors or synthetic signalling pathways are engineered within cells or transferred between genomes. How easily new pathways can be introduced into cells depends on the density and distribution of paralogous pathways in the sequence space that is defined by their specificity-determining residues4,5. Here we directly investigate how crowded this sequence space is, by generating novel two-component signalling proteins in Escherichia coli using cell sorting coupled to deep sequencing to analyse large libraries designed on the basis of coevolutionary patterns. We produce 58 insulated pathways comprising functional kinase-substrate pairs that have different specificities than their parent proteins, and demonstrate that several of these new pairs are orthogonal to all 27 paralogous pathways in E. coli. Additionally, from the kinase-substrate pairs generated, we identify sets consisting of six pairs that are mutually orthogonal to each other, which considerably increases the two-component signalling capacity of E. coli. These results indicate that sequence space is not densely occupied. The relative sparsity of paralogues in sequence space suggests that new insulated pathways can arise easily during evolution, or be designed de novo. We demonstrate the latter by engineering a signalling pathway in E. coli that responds to a plant cytokinin, without crosstalk to extant pathways. Our work also demonstrates how coevolution-guided mutagenesis and the mapping of sequence space can be used to design large sets of orthogonal protein-protein interactions.

Definition of a saxitoxin (STX) binding code enables discovery and characterization of the anuran saxiphilin family.

American bullfrog (Rana castesbeiana) saxiphilin (RcSxph) is a high-affinity "toxin sponge" protein thought to prevent intoxication by saxitoxin (STX), a lethal bis-guanidinium neurotoxin that causes paralytic shellfish poisoning (PSP) by blocking voltage-gated sodium channels (NaVs). How specific RcSxph interactions contribute to STX binding has not been defined and whether other organisms have similar proteins is unclear. Here, we use mutagenesis, ligand binding, and structural studies to define the energetic basis of Sxph:STX recognition. The resultant STX "recognition code" enabled engineering of RcSxph to improve its ability to rescue NaVs from STX and facilitated discovery of 10 new frog and toad Sxphs. Definition of the STX binding code and Sxph family expansion among diverse anurans separated by ∼140 My of evolution provides a molecular basis for understanding the roles of toxin sponge proteins in toxin resistance and for developing novel proteins to sense or neutralize STX and related PSP toxins.

High-throughput dissection of the thermodynamic and conformational properties of a ubiquitous class of RNA tertiary contact motifs.

Despite RNA's diverse secondary and tertiary structures and its complex conformational changes, nature utilizes a limited set of structural "motifs"-helices, junctions, and tertiary contact modules-to build diverse functional RNAs. Thus, in-depth descriptions of a relatively small universe of RNA motifs may lead to predictive models of RNA tertiary conformational landscapes. Motifs may have different properties depending on sequence and secondary structure, giving rise to subclasses that expand the universe of RNA building blocks. Yet we know very little about motif subclasses, given the challenges in mapping conformational properties in high throughput. Previously, we used "RNA on a massively parallel array" (RNA-MaP), a quantitative, high-throughput technique, to study thousands of helices and two-way junctions. Here, we adapt RNA-MaP to study the thermodynamic and conformational properties of tetraloop/tetraloop receptor (TL/TLR) tertiary contact motifs, analyzing 1,493 TLR sequences from different classes. Clustering analyses revealed variability in TL specificity, stability, and conformational behavior. Nevertheless, natural GAAA/11ntR TL/TLRs, while varying in tertiary stability by ∼2.5 kcal/mol, exhibited conserved TL specificity and conformational properties. Thus, RNAs may tune stability without altering the overall structure of these TL/TLRs. Furthermore, their stability correlated with natural frequency, suggesting thermodynamics as the dominant selection pressure. In contrast, other TL/TLRs displayed heterogenous conformational behavior and appear to not be under strong thermodynamic selection. Our results build toward a generalizable model of RNA-folding thermodynamics based on the properties of isolated motifs, and our characterized TL/TLR library can be used to engineer RNAs with predictable thermodynamic and conformational behavior.

Molecular physiology of chemical defenses in a poison frog.

Poison frogs sequester small molecule lipophilic alkaloids from their diet of leaf litter arthropods for use as chemical defenses against predation. Although the dietary acquisition of chemical defenses in poison frogs is well documented, the physiological mechanisms of alkaloid sequestration has not been investigated. Here, we used RNA sequencing and proteomics to determine how alkaloids impact mRNA or protein abundance in the little devil frog (Oophaga sylvatica), and compared wild-caught chemically defended frogs with laboratory frogs raised on an alkaloid-free diet. To understand how poison frogs move alkaloids from their diet to their skin granular glands, we focused on measuring gene expression in the intestines, skin and liver. Across these tissues, we found many differentially expressed transcripts involved in small molecule transport and metabolism, as well as sodium channels and other ion pumps. We then used proteomic approaches to quantify plasma proteins, where we found several protein abundance differences between wild and laboratory frogs, including the amphibian neurotoxin binding protein saxiphilin. Finally, because many blood proteins are synthesized in the liver, we used thermal proteome profiling as an untargeted screen for soluble proteins that bind the alkaloid decahydroquinoline. Using this approach, we identified several candidate proteins that interact with this alkaloid, including saxiphilin. These transcript and protein abundance patterns suggest that the presence of alkaloids influences frog physiology and that small molecule transport proteins may be involved in toxin bioaccumulation in dendrobatid poison frogs.

A toxic environment selects for specialist microbiome in poison frogs.

Shifts in microbiome community composition can have large effects on host health. It is therefore important to understand how perturbations, like those caused by the introduction of exogenous chemicals, modulate microbiome community composition. In poison frogs within the family Dendrobatidae, the skin microbiome is exposed to the alkaloids that the frogs sequester from their diet and use for defense. Given the demonstrated antimicrobial effects of these poison frog alkaloids, these compounds may be structuring the skin microbial community. To test this, we first characterized microbial communities from chemically defended and closely related non-defended frogs from Ecuador. Then we conducted a laboratory experiment to monitor the effect of the alkaloid decahydroquinoline (DHQ) on the microbiome of a single frog species. In both the field and lab experiments, we found that alkaloid-exposed microbiomes are more species rich and phylogenetically diverse, with an increase in rare taxa. To better understand the strain-specific behavior in response to alkaloids, we cultured microbial strains from poison frog skin and found the majority of strains exhibited either enhanced growth or were not impacted by the addition of DHQ. Additionally, stable isotope tracing coupled to nanoSIMS suggests that some of these strains are able to metabolize DHQ. Taken together, these data suggest that poison frog chemical defenses open new niches for skin-associated microbes with specific adaptations, including the likely metabolism of alkaloids, that enable their survival in this toxic environment. This work helps expand our understanding of how exposure to exogenous compounds like alkaloids can impact host microbiomes.

Binding and sequestration of poison frog alkaloids by a plasma globulin.

Alkaloids are important bioactive molecules throughout the natural world, and in many animals they serve as a source of chemical defense against predation. Dendrobatid poison frogs bioaccumulate alkaloids from their diet to make themselves toxic or unpalatable to predators. Despite the proposed roles of plasma proteins as mediators of alkaloid trafficking and bioavailability, the responsible proteins have not been identified. We use chemical approaches to show that a ~50 kDa plasma protein is the principal alkaloid-binding molecule in blood of poison frogs. Proteomic and biochemical studies establish this plasma protein to be a liver-derived alkaloid-binding globulin (ABG) that is a member of the serine-protease inhibitor (serpin) family. In addition to alkaloid-binding activity, ABG sequesters and regulates the bioavailability of 'free' plasma alkaloids in vitro. Unexpectedly, ABG is not related to saxiphilin, albumin, or other known vitamin carriers, but instead exhibits sequence and structural homology to mammalian hormone carriers and amphibian biliverdin-binding proteins. ABG represents a new small molecule binding functionality in serpin proteins, a novel mechanism of plasma alkaloid transport in poison frogs, and more broadly points toward serpins acting as tunable scaffolds for small molecule binding and transport across different organisms.

Molecular physiology of pumiliotoxin sequestration in a poison frog.

Poison frogs bioaccumulate alkaloids for chemical defense from their arthropod diet. Although many alkaloids are accumulated without modification, some poison frog species can metabolize pumiliotoxin (PTX 251D) into the more potent allopumiliotoxin (aPTX 267A). Despite extensive research characterizing the chemical arsenal of poison frogs, the physiological mechanisms involved in the sequestration and metabolism of individual alkaloids remain unclear. We first performed a feeding experiment with the Dyeing poison frog (Dendrobates tinctorius) to ask if this species can metabolize PTX 251D into aPTX 267A and what gene expression changes are associated with PTX 251D exposure in the intestines, liver, and skin. We found that D. tinctorius can metabolize PTX 251D into aPTX 267A, and that PTX 251D exposure changed the expression level of genes involved in immune system function and small molecule metabolism and transport. To better understand the functional significance of these changes in gene expression, we then conducted a series of high-throughput screens to determine the molecular targets of PTX 251D and identify potential proteins responsible for metabolism of PTX 251D into aPTX 267A. Although screens of PTX 251D binding human voltage-gated ion channels and G-protein coupled receptors were inconclusive, we identified human CYP2D6 as a rapid metabolizer of PTX 251D in a cytochrome P450 screen. Furthermore, a CYP2D6-like gene had increased expression in the intestines of animals fed PTX, suggesting this protein may be involved in PTX metabolism. These results show that individual alkaloids can modify gene expression across tissues, including genes involved in alkaloid metabolism. More broadly, this work suggests that specific alkaloid classes in wild diets may induce physiological changes for targeted accumulation and metabolism.

Poison frog dietary preference depends on prey type and alkaloid load.

The ability to acquire chemical defenses through the diet has evolved across several major taxa. Chemically defended organisms may need to balance chemical defense acquisition and nutritional quality of prey items. However, these dietary preferences and potential trade-offs are rarely considered in the framework of diet-derived defenses. Poison frogs (Family Dendrobatidae) acquire defensive alkaloids from their arthropod diet of ants and mites, although their dietary preferences have never been investigated. We conducted prey preference assays with the Dyeing Poison frog (Dendrobates tinctorius) to test the hypothesis that alkaloid load and prey traits influence frog dietary preferences. We tested size preferences (big versus small) within each of four prey groups (ants, beetles, flies, and fly larvae) and found that frogs preferred interacting with smaller prey items of the fly and beetle groups. Frog taxonomic prey preferences were also tested as we experimentally increased their chemical defense load by feeding frogs decahydroquinoline, an alkaloid compound similar to those naturally found in their diet. Contrary to our expectations, overall preferences did not change during alkaloid consumption, as frogs across groups preferred fly larvae over other prey. Finally, we assessed the protein and lipid content of prey items and found that small ants have the highest lipid content while large fly larvae have the highest protein content. Our results suggest that consideration of toxicity and prey nutritional value are important factors in understanding the evolution of acquired chemical defenses and niche partitioning.