Circular permutation profiling by deep sequencing libraries created using transposon mutagenesis

Rational Design of a Circularly Permuted Flavin-Based Fluorescent Protein

Flavin-based fluorescent proteins are oxygen-independent reporters that hold great promise for imaging anaerobic and hypoxic biological systems. In this study, we explored the feasibility of applying circular permutation, a valuable method for the creation of fluorescent sensors, to flavin-based fluorescent proteins. We used rational design and structural data to identify a suitable location for circular permutation in iLOV, a flavin-based reporter derived from A. thaliana. However, relocating the N- and C-termini to this position resulted in a significant reduction in fluorescence. This loss of fluorescence was reversible, however, by fusing dimerizing coiled coils at the new N- and C-termini to compensate for the increase in local chain entropy. Additionally, by inserting protease cleavage sites in circularly permuted iLOV, we developed two protease sensors and demonstrated their application in mammalian cells. In summary, our work establishes the first approach to engineer circularly permuted FbFPs optimized for high fluorescence and further showcases the utility of circularly permuted FbFPs to serve as a scaffold for sensor engineering.

Multiparametric domain insertional profiling of adeno-associated virus VP1

Several evolved properties of adeno-associated virus (AAV), such as broad tropism and immunogenicity in humans, are barriers to AAV-based gene therapy. Most efforts to re-engineer these properties have focused on variable regions near AAV's 3-fold protrusions and capsid protein termini. To comprehensively survey AAV capsids for engineerable hotspots, we determined multiple AAV fitness phenotypes upon insertion of six structured protein domains into the entire AAV-DJ capsid protein VP1. This is the largest and most comprehensive AAV domain insertion dataset to date. Our data revealed a surprising robustness of AAV capsids to accommodate large domain insertions. Insertion permissibility depended strongly on insertion position, domain type, and measured fitness phenotype, which clustered into contiguous structural units that we could link to distinct roles in AAV assembly, stability, and infectivity. We also identified engineerable hotspots of AAV that facilitate the covalent attachment of binding scaffolds, which may represent an alternative approach to re-direct AAV tropism.

A cellular selection identifies elongated flavodoxins that support electron transfer to sulfite reductase

Flavodoxins (Flds) mediate the flux of electrons between oxidoreductases in diverse metabolic pathways. To investigate whether Flds can support electron transfer to a sulfite reductase (SIR) that evolved to couple with a ferredoxin, we evaluated the ability of Flds to transfer electrons from a ferredoxin-NADP reductase (FNR) to a ferredoxin-dependent SIR using growth complementation of an Escherichia coli strain with a sulfur metabolism defect. We show that Flds from cyanobacteria complement this growth defect when coexpressed with an FNR and an SIR that evolved to couple with a plant ferredoxin. When we evaluated the effect of peptide insertion on Fld-mediated electron transfer, we observed a sensitivity to insertions within regions predicted to be proximal to the cofactor and partner binding sites, while a high insertion tolerance was detected within loops distal from the cofactor and within regions of helices and sheets that are proximal to those loops. Bioinformatic analysis showed that natural Fld sequence variability predicts a large fraction of the motifs that tolerate insertion of the octapeptide SGRPGSLS. These results represent the first evidence that Flds can support electron transfer to assimilatory SIRs, and they suggest that the pattern of insertion tolerance is influenced by interactions with oxidoreductase partners.

Multiparametric domain insertional profiling of Adeno-Associated Virus VP1

Evolved properties of Adeno-Associated Virus (AAV), such as broad tropism and immunogenicity in humans, are barriers to AAV-based gene therapy. Previous efforts to re-engineer these properties have focused on variable regions near AAV’s 3-fold protrusions and capsid protein termini. To comprehensively survey AAV capsids for engineerable hotspots, we determined multiple AAV fitness phenotypes upon insertion of large, structured protein domains into the entire AAV-DJ capsid protein VP1. This is the largest and most comprehensive AAV domain insertion dataset to date. Our data revealed a surprising robustness of AAV capsids to accommodate large domain insertions. There was strong positional, domain-type, and fitness phenotype dependence of insertion permissibility, which clustered into correlated structural units that we could link to distinct roles in AAV assembly, stability, and infectivity. We also identified new engineerable hotspots of AAV that facilitate the covalent attachment of binding scaffolds, which may represent an alternative approach to re-direct AAV tropism.

Nondestructive Chemical Sensing within Bulk Soil Using 1000 Biosensors Per Gram of Matrix

Gene expression can be monitored in hard-to-image environmental materials using gas-reporting biosensors, but these outputs have only been applied in autoclaved matrices that are hydrated with rich medium. To better understand the compatibility of indicator gas reporting with environmental samples, we evaluated how matrix hydration affects the gas signal of an engineered microbe added to a sieved soil. A gas-reporting microbe presented a gas signal in a forest soil (Alfisol) when hydrated to an environmentally relevant osmotic pressure. When the gas signal was concentrated prior to analysis, a biosensor titer of 10³ cells/gram of soil produced a significant signal when soil was supplemented with halides. A signal was also observed without halide amendment, but a higher cell titer (10⁶ cells/gram of soil) was required. A sugar-regulated gas biosensor was able to report with a similar level of sensitivity when added to an unsterilized soil matrix, illustrating how gas concentration enables biosensing within a soil containing environmental microbes. These results establish conditions where engineered microbes can report on gene expression in living environmental matrices with decreased perturbation of the soil environment compared to previously reported approaches, using biosensor titers that are orders of magnitude lower than the number of cells typically observed in a gram of soil.

Determinants of Multiheme Cytochrome Extracellular Electron Transfer Uncovered by Systematic Peptide Insertion

The multiheme cytochrome MtrA enables microbial respiration by transferring electrons across the outer membrane to extracellular electron acceptors. While structural studies have identified residues that mediate the binding of MtrA to hemes and to other cytochromes that facilitate extracellular electron transfer (EET), the relative importance of these interactions for EET is not known. To better understand EET, we evaluated how insertion of an octapeptide across all MtrA backbone locations affects Shewanella oneidensis MR-1 respiration on Fe(III). The EET efficiency was found to be inversely correlated with the proximity of the insertion to the heme prosthetic groups. Mutants with decreased EET efficiencies also arose from insertions in a subset of the regions that make residue-residue contacts with the porin MtrB, while all sites contacting the extracellular cytochrome MtrC presented high peptide insertion tolerance. MtrA variants having peptide insertions within the CXXCH motifs that coordinate heme cofactors retained some ability to support respiration on Fe(III), although these variants presented significantly decreased EET efficiencies. Furthermore, the fitness of cells expressing different MtrA variants under Fe(III) respiration conditions correlated with anode reduction. The peptide insertion profile, which represents the first comprehensive sequence-structure-function map for a multiheme cytochrome, implicates MtrA as a strategic protein engineering target for the regulation of EET.

High-throughput split-protein profiling by combining transposon mutagenesis and regulated protein-protein interactions with deep sequencing

Splitting a protein at a position may lead to self- or assisted-complementary fragments depending on whether two resulting fragments can reconstitute to maintain the native function spontaneously or require assistance from two interacting molecules. Assisted complementary fragments with high contrast are an important tool for probing biological interactions. However, only a small number of assisted-complementary split-variants have been identified due to manual, labour-intensive optimization of a candidate gene. Here, we introduce a technique for high-throughput split-protein profiling (HiTS) that allows fast identification of self- and assisted complementary positions by transposon mutagenesis, a rapamycin-regulated FRB-FKBP protein interaction pair, and deep sequencing. We test this technique by profiling three antibiotic-resistant genes (fosfomycin-resistant gene, fosA3, erythromycin-resistant gene, ermB, and chloramphenicol-resistant gene, catI). Self- and assisted complementary fragments discovered by the high-throughput technique were subsequently confirmed by low-throughput testing of individual split positions. Thus, the HiTS technique provides a quicker alternative for discovering the proteins with suitable self- and assisted-complementary split positions when combining with a readout such as fluorescence, bioluminescence, cell survival, gene transcription or genome editing.

Determinants of trafficking, conduction, and disease within a K+ channel revealed through multiparametric deep mutational scanning

A long-standing goal in protein science and clinical genetics is to develop quantitative models of sequence, structure, and function relationships to understand how mutations cause disease. Deep mutational scanning (DMS) is a promising strategy to map how amino acids contribute to protein structure and function and to advance clinical variant interpretation. Here, we introduce 7429 single-residue missense mutations into the inward rectifier K⁺ channel Kir2.1 and determine how this affects folding, assembly, and trafficking, as well as regulation by allosteric ligands and ion conduction. Our data provide high-resolution information on a cotranslationally folded biogenic unit, trafficking and quality control signals, and segregated roles of different structural elements in fold stability and function. We show that Kir2.1 surface trafficking mutants are underrepresented in variant effect databases, which has implications for clinical practice. By comparing fitness scores with expert-reviewed variant effects, we can predict the pathogenicity of 'variants of unknown significance' and disease mechanisms of known pathogenic mutations. Our study in Kir2.1 provides a blueprint for how multiparametric DMS can help us understand the mechanistic basis of genetic disorders and the structure-function relationships of proteins.

Probing ion channel functional architecture and domain recombination compatibility by massively parallel domain insertion profiling

Protein domains are the basic units of protein structure and function. Comparative analysis of genomes and proteomes showed that domain recombination is a main driver of multidomain protein functional diversification and some of the constraining genomic mechanisms are known. Much less is known about biophysical mechanisms that determine whether protein domains can be combined into viable protein folds. Here, we use massively parallel insertional mutagenesis to determine compatibility of over 300,000 domain recombination variants of the Inward Rectifier K⁺ channel Kir2.1 with channel surface expression. Our data suggest that genomic and biophysical mechanisms acted in concert to favor gain of large, structured domain at protein termini during ion channel evolution. We use machine learning to build a quantitative biophysical model of domain compatibility in Kir2.1 that allows us to derive rudimentary rules for designing domain insertion variants that fold and traffic to the cell surface. Positional Kir2.1 responses to motif insertion clusters into distinct groups that correspond to contiguous structural regions of the channel with distinct biophysical properties tuned towards providing either folding stability or gating transitions. This suggests that insertional profiling is a high-throughput method to annotate function of ion channel structural regions.

Comprehensive deletion landscape of CRISPR-Cas9 identifies minimal RNA-guided DNA-binding modules

Proteins evolve through the modular rearrangement of elements known as domains. Extant, multidomain proteins are hypothesized to be the result of domain accretion, but there has been limited experimental validation of this idea. Here, we introduce a technique for genetic minimization by iterative size-exclusion and recombination (MISER) for comprehensively making all possible deletions of a protein. Using MISER, we generate a deletion landscape for the CRISPR protein Cas9. We find that the catalytically-dead Streptococcus pyogenes Cas9 can tolerate large single deletions in the REC2, REC3, HNH, and RuvC domains, while still functioning in vitro and in vivo, and that these deletions can be stacked together to engineer minimal, DNA-binding effector proteins. In total, our results demonstrate that extant proteins retain significant modularity from the accretion process and, as genetic size is a major limitation for viral delivery systems, establish a general technique to improve genome editing and gene therapy-based therapeutics.

A systematic approach to inserting split inteins for Boolean logic gate engineering and basal activity reduction

Split inteins are powerful tools for seamless ligation of synthetic split proteins. Yet, their use remains limited because the already intricate split site identification problem is often complicated by the requirement of extein junction sequences. To address this, we augment a mini-Mu transposon-based screening approach and devise the intein-assisted bisection mapping (IBM) method. IBM robustly reveals clusters of split sites on five proteins, converting them into AND or NAND logic gates. We further show that the use of inteins expands functional sequence space for splitting a protein. We also demonstrate the utility of our approach over rational inference of split sites from secondary structure alignment of homologous proteins, and that basal activities of highly active proteins can be mitigated by splitting them. Our work offers a generalizable and systematic route towards creating split protein-intein fusions for synthetic biology.

A mini-Tn5-derived transposon with reportable and selectable markers enables rapid generation and screening of insertional mutants in Gram-negative bacteria

We re-engineered a classic tool for mutagenesis and gene expression studies in Gram-negative bacteria. Our modified Tn5-based transposon contains multiple features that allow rapid selection for mutants, direct quantification of gene expression and straightforward cloning of the inactivated gene. The promoter-less gfp-km cassette provides selection and reporter assay depending on the activity of the promoter upstream of the transposon insertion site. The cat gene facilitates positive antibiotic selection for mutants, while the narrow R6Kγ replication origin forces transposition in recipient strains lacking the pir gene and enables cloning of the transposon flanked with the disrupted gene from the chromosome. The suicide vector pCKD100, a plasmid that could be delivered into recipient cells through biparental mating or electroporation, harbours the modified transposon. We used the transposon to mutagenize Pectobacterium versatile KD100, Pseudumonas coronafaciens PC27R and Escherichia coli 35150N. The fluorescence intensities of mutants expressing high GFP could be quantified and detected qualitatively. Transformation efficiency from conjugation ranged from 1600 to 1900 CFU per ml. We sequenced the upstream flanking regions, identified the putative truncated genes and demonstrated the restoration of the GFP phenotype through marker exchange. The mini-Tn5 transposon was also utilized to construct mutant a library of P. versatile for forward genetic screens.

Targeted insertional mutagenesis libraries for deep domain insertion profiling

Domain recombination is a key principle in protein evolution and protein engineering, but inserting a donor domain into every position of a target protein is not easily experimentally accessible. Most contemporary domain insertion profiling approaches rely on DNA transposons, which are constrained by sequence bias. Here, we establish Saturated Programmable Insertion Engineering (SPINE), an unbiased, comprehensive, and targeted domain insertion library generation technique using oligo library synthesis and multi-step Golden Gate cloning. Through benchmarking to MuA transposon-mediated library generation on four ion channel genes, we demonstrate that SPINE-generated libraries are enriched for in-frame insertions, have drastically reduced sequence bias as well as near-complete and highly-redundant coverage. Unlike transposon-mediated domain insertion that was severely biased and sparse for some genes, SPINE generated high-quality libraries for all genes tested. Using the Inward Rectifier K⁺ channel Kir2.1, we validate the practical utility of SPINE by constructing and comparing domain insertion permissibility maps. SPINE is the first technology to enable saturated domain insertion profiling. SPINE could help explore the relationship between domain insertions and protein function, and how this relationship is shaped by evolutionary forces and can be engineered for biomedical applications.

Protein tolerance to random circular permutation correlates with thermostability and local energetics of residue-residue contacts

Adenylate kinase (AK) orthologs with a range of thermostabilities were subjected to random circular permutation, and deep mutational scanning was used to evaluate where new protein termini were nondisruptive to activity. The fraction of circularly permuted variants that retained function in each library correlated with AK thermostability. In addition, analysis of the positional tolerance to new termini, which increase local conformational flexibility, showed that bonds were either functionally sensitive to cleavage across all homologs, differentially sensitive, or uniformly tolerant. The mobile AMP-binding domain, which displays the highest calculated contact energies, presented the greatest tolerance to new termini across all AKs. In contrast, retention of function in the lid and core domains was more dependent upon AK melting temperature. These results show that family permutation profiling identifies primary structure that has been selected by evolution for dynamics that are critical to activity within an enzyme family. These findings also illustrate how deep mutational scanning can be applied to protein homologs in parallel to differentiate how topology, stability, and local energetics govern mutational tolerance.