Saturation mutagenesis of α-synuclein reveals monomer fold that modulates aggregation.

α-Synuclein (aSyn) aggregation underlies neurodegenerative synucleinopathies. aSyn seeds are proposed to replicate and propagate neuronal pathology like prions. Seeding of aSyn can be recapitulated in cellular systems of aSyn aggregation; however, the mechanism of aSyn seeding and its regulation are not well understood. We developed an mEos-based aSyn seeding assay and performed saturation mutagenesis to identify with single-residue resolution positive and negative regulators of aSyn aggregation. We not only found the core regions that govern aSyn aggregation but also identified mutants outside of the core that enhance aggregation. We identified local structure within the N terminus of aSyn that hinders the fibrillization propensity of its aggregation-prone core. Based on the screen, we designed a minimal aSyn fragment that shows a ~4-fold enhancement in seeding activity and enabled discrimination of synucleinopathies. Our study expands the basic knowledge of aSyn aggregation and advances the design of cellular systems of aSyn aggregation to diagnose synucleinopathies based on protein conformation.

Integrative taxonomy reveals two new narrowly-endemic crayfish species (Decapoda: Cambaridae) from the Yadkin River Basin in western North Carolina, USA.

Two new species of stream-dwelling crayfish, Cambarus lapidosus, the Stony Fork Crayfish, and Cambarus burchfielae, the Falls Crayfish, are described from the Yadkin River basin in western North Carolina, USA, using an integrative taxonomic approach consisting of morphological, genetic, and biogeographic data. Both species were previously considered to be members of the widely distributed Cambarus species C complex, which occurs throughout mid-Atlantic Slope river basins; however, they are in fact morphologically and genetically more similar to the Cambarus robustus species complex from interior basins in the south-central Appalachians, indicating Atlantic basin stream capture of an Interior basin faunal group has occurred in this region. Both new species described herein can be differentiated from these two complexes, and each other, by several morphological characteristics. Cambarus lapidosus and C. burchfielae are differentiated from C. species C by the absence of cervical spines and the presence of 1-2 subpalmar tubercles on the chelae in most specimens; both species are less-punctate across the areola than C. aff. robustus. Cambarus burchfielae is further differentiated from C. lapidosus, C. species C, and C. aff. robustus by the presence of a narrower and sparsely punctate areola, a single weak row of tubercles on the mesial margin of the palm, and the absence of tubercles on the dorsal surface of the dactyl. The newly described species are genetically sister taxa and together are most similar by genetic distance to undescribed members of the C. robustus species complex from the nearby (~30 km) Watauga River basin and most phylogenetically similar to C. aff. robustus from the adjacent New River basin; both are interior basin drainages. The newly-described species are endemic to the Blue Ridge Mountains of western North Carolina and have restricted distributions (<100 km2-150 km2) in small (<4th order) tributaries to the Yadkin River. Cambarus lapidosus is known only from the upper Stony Fork watershed and C. burchfielae is known only from the upper reaches of the adjacent Lewis Fork watershed; both species likely qualify for conservation status protections under narrow geographic range criteria.

Undersampling and the inference of coevolution in proteins.

Protein structure, function, and evolution depend on local and collective epistatic interactions between amino acids. A powerful approach to defining these interactions is to construct models of couplings between amino acids that reproduce the empirical statistics (frequencies and correlations) observed in sequences comprising a protein family. The top couplings are then interpreted. Here, we show that as currently implemented, this inference unequally represents epistatic interactions, a problem that fundamentally arises from limited sampling of sequences in the context of distinct scales at which epistasis occurs in proteins. We show that these issues explain the ability of current approaches to predict tertiary contacts between amino acids and the inability to obviously expose larger networks of functionally relevant, collectively evolving residues called sectors. This work provides a necessary foundation for more deeply understanding and improving evolution-based models of proteins.

The Antibiotic Efflux Protein TolC Is a Highly Evolvable Target under Colicin E1 or TLS Phage Selection.

Bacteriophages and bacterial toxins are promising antibacterial agents to treat infections caused by multidrug-resistant (MDR) bacteria. In fact, bacteriophages have recently been successfully used to treat life-threatening infections caused by MDR bacteria (Schooley RT, Biswas B, Gill JJ, Hernandez-Morales A, Lancaster J, Lessor L, Barr JJ, Reed SL, Rohwer F, Benler S, et al. 2017. Development and use of personalized bacteriophage-based therapeutic cocktails to treat a patient with a disseminated resistant Acinetobacter baumannii infection. Antimicrob Agents Chemother. 61(10); Chan BK, Turner PE, Kim S, Mojibian HR, Elefteriades JA, Narayan D. 2018. Phage treatment of an aortic graft infected with Pseudomonas aeruginosa. Evol Med Public Health. 2018(1):60-66; Petrovic Fabijan A, Lin RCY, Ho J, Maddocks S, Ben Zakour NL, Iredell JR, Westmead Bacteriophage Therapy Team. 2020. Safety of bacteriophage therapy in severe Staphylococcus aureus infection. Nat Microbiol. 5(3):465-472). One potential problem with using these antibacterial agents is the evolution of resistance against them in the long term. Here, we studied the fitness landscape of the Escherichia coli TolC protein, an outer membrane efflux protein that is exploited by a pore forming toxin called colicin E1 and by TLS phage (Pagie L, Hogeweg P. 1999. Colicin diversity: a result of eco-evolutionary dynamics. J Theor Biol. 196(2):251-261; Andersen C, Hughes C, Koronakis V. 2000. Chunnel vision. Export and efflux through bacterial channel-tunnels. EMBO Rep. 1(4):313-318; Koronakis V, Andersen C, Hughes C. 2001. Channel-tunnels. Curr Opin Struct Biol. 11(4):403-407; Czaran TL, Hoekstra RF, Pagie L. 2002. Chemical warfare between microbes promotes biodiversity. Proc Natl Acad Sci U S A. 99(2):786-790; Cascales E, Buchanan SK, Duché D, Kleanthous C, Lloubès R, Postle K, Riley M, Slatin S, Cavard D. 2007. Colicin biology. Microbiol Mol Biol Rev. 71(1):158-229). By systematically assessing the distribution of fitness effects of ∼9,000 single amino acid replacements in TolC using either positive (antibiotics and bile salts) or negative (colicin E1 and TLS phage) selection pressures, we quantified evolvability of the TolC. We demonstrated that the TolC is highly optimized for the efflux of antibiotics and bile salts. In contrast, under colicin E1 and TLS phage selection, TolC sequence is very sensitive to mutations. Finally, we have identified a large set of mutations in TolC that increase resistance of E. coli against colicin E1 or TLS phage without changing antibiotic susceptibility of bacterial cells. Our findings suggest that TolC is a highly evolvable target under negative selection which may limit the potential clinical use of bacteriophages and bacterial toxins if evolutionary aspects are not taken into account.

RNA sectors and allosteric function within the ribosome.

The ribosome translates the genetic code into proteins in all domains of life. Its size and complexity demand long-range interactions that regulate ribosome function. These interactions are largely unknown. Here, we apply a global coevolution method, statistical coupling analysis (SCA), to identify coevolving residue networks (sectors) within the 23S ribosomal RNA (rRNA) of the large ribosomal subunit. As in proteins, SCA reveals a hierarchical organization of evolutionary constraints with near-independent groups of nucleotides forming physically contiguous networks within the three-dimensional structure. Using a quantitative, continuous-culture-with-deep-sequencing assay, we confirm that the top two SCA-predicted sectors contribute to ribosome function. These sectors map to distinct ribosome activities, and their origins trace to phylogenetic divergences across all domains of life. These findings provide a foundation to map ribosome allostery, explore ribosome biogenesis, and engineer ribosomes for new functions. Despite differences in chemical structure, protein and RNA enzymes appear to share a common internal logic of interaction and assembly.

An evolution-based model for designing chorismate mutase enzymes.

The rational design of enzymes is an important goal for both fundamental and practical reasons. Here, we describe a process to learn the constraints for specifying proteins purely from evolutionary sequence data, design and build libraries of synthetic genes, and test them for activity in vivo using a quantitative complementation assay. For chorismate mutase, a key enzyme in the biosynthesis of aromatic amino acids, we demonstrate the design of natural-like catalytic function with substantial sequence diversity. Further optimization focuses the generative model toward function in a specific genomic context. The data show that sequence-based statistical models suffice to specify proteins and provide access to an enormous space of functional sequences. This result provides a foundation for a general process for evolution-based design of artificial proteins.

Cambarus franklini, a new crayfish (Decapoda: Cambaridae) from the Catawba River Basin in western North Carolina, USA.

A new species of stream-dwelling crayfish, Cambarus franklini, the South Mountains crayfish, is described from the upper South Fork Catawba River basin in western North Carolina, USA using morphological and genetic data. Cambarus franklini was previously considered a member of the widespread and morphologically variable Cambarus species C complex and is morphologically most similar to an undiagnosed member of the group native to the upper Catawba River basin in NC. Cambarus franklini can be differentiated from this species group by several morphological characteristics including: lacking a well-defined double row of tubercles along the mesial margin of the palm, possessing a more weakly convergent and longer acumen, and conspicuous blue-green and red coloration, particularly throughout the telson and along the distal margins of the rami. This species is phylogenetically most similar to Cambarus johni, Cooper, 2006, another former member of the Cambarus species C group. Cambarus franklini has a limited geographic range (100 km2) and is currently known only from the Henry and Jacob Fork watersheds in the South Mountains region of the Eastern Blue Ridge foothills.

A set of experimentally validated, mutually orthogonal primers for combinatorially specifying genetic components.

The design and synthesis of novel genes and deoxyribonucleic acid (DNA) sequences is a central technique in synthetic biology. Current methods of high throughput gene synthesis use pooled oligonucleotides obtained from custom-designed DNA microarray chips, and rely on orthogonal (non-interacting) polymerase chain reaction primers to specifically de-multiplex, by amplification, the precise subset of oligonucleotides necessary to assemble a full length gene. The availability of a large validated set of mutually orthogonal primers is therefore a crucial reagent for high-throughput gene synthesis. Here, we present a set of 166 20-nucleotide primers that are experimentally verified to be non-interacting, capable of specifying 13 695 unique genes. These primers represent a valuable resource to the synthetic biology community for specifying genetic components that can be assembled through a scalable and modular architecture.

An electrostatic selection mechanism controls sequential kinase signaling downstream of the T cell receptor.

The sequence of events that initiates T cell signaling is dictated by the specificities and order of activation of the tyrosine kinases that signal downstream of the T cell receptor. Using a platform that combines exhaustive point-mutagenesis of peptide substrates, bacterial surface-display, cell sorting, and deep sequencing, we have defined the specificities of the first two kinases in this pathway, Lck and ZAP-70, for the T cell receptor ζ chain and the scaffold proteins LAT and SLP-76. We find that ZAP-70 selects its substrates by utilizing an electrostatic mechanism that excludes substrates with positively-charged residues and favors LAT and SLP-76 phosphosites that are surrounded by negatively-charged residues. This mechanism prevents ZAP-70 from phosphorylating its own activation loop, thereby enforcing its strict dependence on Lck for activation. The sequence features in ZAP-70, LAT, and SLP-76 that underlie electrostatic selectivity likely contribute to the specific response of T cells to foreign antigens.

Evolution-based design of proteins.

Statistical analysis of protein sequences indicates an architecture for natural proteins in which amino acids are engaged in a sparse, hierarchical pattern of interactions in the tertiary structure. This architecture might be a key and distinguishing feature of evolved proteins-a design principle providing not only for foldability and high-performance function but also for robustness to perturbation and the capacity for rapid adaptation to new selection pressures. Here, we describe an approach for systematically testing this design principle for natural-like proteins by (1) computational design of synthetic sequences that gradually add or remove constraints along the hierarchy of interacting residues and (2) experimental testing of the designed sequences for folding and biochemical function. By this process, we hope to understand how the constraints on fold, function, and other aspects of fitness are organized within natural proteins, a first step in understanding the process of "design" by evolution.

An interdomain sector mediating allostery in Hsp70 molecular chaperones.

Allosteric coupling between protein domains is fundamental to many cellular processes. For example, Hsp70 molecular chaperones use ATP binding by their actin-like N-terminal ATPase domain to control substrate interactions in their C-terminal substrate-binding domain, a reaction that is critical for protein folding in cells. Here, we generalize the statistical coupling analysis to simultaneously evaluate co-evolution between protein residues and functional divergence between sequences in protein sub-families. Applying this method in the Hsp70/110 protein family, we identify a sparse but structurally contiguous group of co-evolving residues called a 'sector', which is an attribute of the allosteric Hsp70 sub-family that links the functional sites of the two domains across a specific interdomain interface. Mutagenesis of Escherichia coli DnaK supports the conclusion that this interdomain sector underlies the allosteric coupling in this protein family. The identification of the Hsp70 sector provides a basis for further experiments to understand the mechanism of allostery and introduces the idea that cooperativity between interacting proteins or protein domains can be mediated by shared sectors.

Surface sites for engineering allosteric control in proteins.

Statistical analyses of protein families reveal networks of coevolving amino acids that functionally link distantly positioned functional surfaces. Such linkages suggest a concept for engineering allosteric control into proteins: The intramolecular networks of two proteins could be joined across their surface sites such that the activity of one protein might control the activity of the other. We tested this idea by creating PAS-DHFR, a designed chimeric protein that connects a light-sensing signaling domain from a plant member of the Per/Arnt/Sim (PAS) family of proteins with Escherichia coli dihydrofolate reductase (DHFR). With no optimization, PAS-DHFR exhibited light-dependent catalytic activity that depended on the site of connection and on known signaling mechanisms in both proteins. PAS-DHFR serves as a proof of concept for engineering regulatory activities into proteins through interface design at conserved allosteric sites.

Evolutionary information for specifying a protein fold.

Classical studies show that for many proteins, the information required for specifying the tertiary structure is contained in the amino acid sequence. Here, we attempt to define the sequence rules for specifying a protein fold by computationally creating artificial protein sequences using only statistical information encoded in a multiple sequence alignment and no tertiary structure information. Experimental testing of libraries of artificial WW domain sequences shows that a simple statistical energy function capturing coevolution between amino acid residues is necessary and sufficient to specify sequences that fold into native structures. The artificial proteins show thermodynamic stabilities similar to natural WW domains, and structure determination of one artificial protein shows excellent agreement with the WW fold at atomic resolution. The relative simplicity of the information used for creating sequences suggests a marked reduction to the potential complexity of the protein-folding problem.

Natural-like function in artificial WW domains.

Protein sequences evolve through random mutagenesis with selection for optimal fitness. Cooperative folding into a stable tertiary structure is one aspect of fitness, but evolutionary selection ultimately operates on function, not on structure. In the accompanying paper, we proposed a model for the evolutionary constraint on a small protein interaction module (the WW domain) through application of the SCA, a statistical analysis of multiple sequence alignments. Construction of artificial protein sequences directed only by the SCA showed that the information extracted by this analysis is sufficient to engineer the WW fold at atomic resolution. Here, we demonstrate that these artificial WW sequences function like their natural counterparts, showing class-specific recognition of proline-containing target peptides. Consistent with SCA predictions, a distributed network of residues mediates functional specificity in WW domains. The ability to recapitulate natural-like function in designed sequences shows that a relatively small quantity of sequence information is sufficient to specify the global energetics of amino acid interactions.

The erythropoietin receptor transmembrane domain mediates complex formation with viral anemic and polycythemic gp55 proteins.

Erythropoietin receptor (EpoR) activation is crucial for mature red blood cell production. The murine EpoR can also be activated by the envelope protein of the polycythemic (P) spleen focus forming virus (SFFV), gp55-P. Due to differences in the TM sequence, gp55 of the anemic (A) strain SFFV, gp55-A, cannot efficiently activate the EpoR. Using antibody-mediated immunofluorescence co-patching, we show that the majority of EpoR forms hetero-oligomers at the cell surface with gp55-P and, surprisingly, with gp55-A. The EpoR TM domain is targeted by gp55-P and -A, as only chimeric receptors containing EpoR TM sequences oligomerized with gp55 proteins. Both gp55-P and gp55-A are homodimers on the cell surface, as shown by co-patching. However, when the homomeric interactions of the isolated TM domains were assayed by TOXCAT bacterial reporter system, only the TM sequence of gp55-P was dimerized. Thus, homo-oligomerization of gp55 proteins is insufficient for full EpoR activation, and a correct conformation of the dimer in the TM region is required. This is supported by the failure of gp55-A-->P, a mutant protein whose TM domain can homo-oligomerize, to fully activate EpoR. As unliganded EpoR forms TM-dependent but inactive homodimers, we propose that the EpoR can be activated to different extents by homodimeric gp55 proteins, depending on the conformation of the gp55 protein dimer in the TM region.

Knowledge-based potential functions in protein design.

Predicting protein sequences that fold into specific native three-dimensional structures is a problem of great potential complexity. Although the complete solution is ultimately rooted in understanding the physical chemistry underlying the complex interactions between amino acid residues that determine protein stability, recent work shows that empirical information about these first principles is embedded in the statistics of protein sequence and structure databases. This review focuses on the use of 'knowledge-based' potentials derived from these databases in designing proteins. In addition, the data suggest how the study of these empirical potentials might impact our fundamental understanding of the energetic principles of protein structure.