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.

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.

DIMPLE: deep insertion, deletion, and missense mutation libraries for exploring protein variation in evolution, disease, and biology.

Insertions and deletions (indels) enable evolution and cause disease. Due to technical challenges, indels are left out of most mutational scans, limiting our understanding of them in disease, biology, and evolution. We develop a low cost and bias method, DIMPLE, for systematically generating deletions, insertions, and missense mutations in genes, which we test on a range of targets, including Kir2.1. We use DIMPLE to study how indels impact potassium channel structure, disease, and evolution. We find deletions are most disruptive overall, beta sheets are most sensitive to indels, and flexible loops are sensitive to deletions yet tolerate insertions.

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.

A large-scale survey of pairwise epistasis reveals a mechanism for evolutionary expansion and specialization of PDZ domains.

Deep mutational scanning (DMS) facilitates data-driven models of protein structure and function. Here, we adapted Saturated Programmable Insertion Engineering (SPINE) as a programmable DMS technique. We validate SPINE with a reference single mutant dataset in the PSD95 PDZ3 domain and then characterize most pairwise double mutants to study epistasis. We observe wide-spread proximal negative epistasis, which we attribute to mutations affecting thermodynamic stability, and strong long-range positive epistasis, which is enriched in an evolutionarily conserved and function-defining network of "sector" and clade-specifying residues. Conditional neutrality of mutations in clade-specifying residues compensates for deleterious mutations in sector positions. This suggests that epistatic interactions between these position pairs facilitated the evolutionary expansion and specialization of PDZ domains. We propose that SPINE provides easy experimental access to reveal epistasis signatures in proteins that will improve our understanding of the structural basis for protein function and adaptation.

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.

Combinatorial Assembly of Lumitoxins.

Ion channels are among the most important proteins in neuroscience and serve as drug targets for many brain disorders. During development, learning, disease progression, and other processes, the activity levels of specific ion channels are tuned in a cell-type specific manner. However, it is difficult to assess how cell-specific changes in ion channel activity alter emergent brain functions. We have developed a protein architecture for fully genetically encoded light-activated modulation of endogenous ion channel activity. Fusing a genetically encoded photoswitch and an ion channel-modulating peptide toxin in a computationally designed fashion, this reagent, which we call Lumitoxins, can mediate light-modulation of specific endogenous ion channel activities in targeted cells. The modular lumitoxin architecture may be useful in a diversity of neuroscience tools. Here, we delineate how to construct lumitoxin genes from synthesized components, and provide a general outline for how to test their function in mammalian cell culture.

Uncovering divergent evolution of α/ß-hydrolases: a surprising residue substitution needed to convert Hevea brasiliensis hydroxynitrile lyase into an esterase.

Hevea brasiliensis hydroxynitrile lyase (HbHNL) and salicylic acid binding protein 2 (SABP2, an esterase) share 45% amino acid sequence identity, the same protein fold, and even the same catalytic triad of Ser-His-Asp. However, they catalyze different reactions: cleavage of hydroxynitriles and hydrolysis of esters, respectively. To understand how other active site differences in the two enzymes enable the same catalytic triad to catalyze different reactions, we substituted amino acid residues in HbHNL with the corresponding residues from SABP2, expecting hydroxynitrile lyase activity to decrease and esterase activity to increase. Previous mechanistic studies and x-ray crystallography suggested that esterase activity requires removal of an active site lysine and threonine from the hydroxynitrile lyase. The Thr11Gly Lys236Gly substitutions in HbHNL reduced hydroxynitrile lyase activity for cleavage of mandelonitrile 100-fold, but increased esterase activity only threefold to kcat ~ 0.1 min-1 for hydrolysis of p-nitrophenyl acetate. Adding a third substitution - Glu79His - increased esterase activity more than tenfold to kcat ~ 1.6 min-1. The specificity constant (kcat/KM) for this triple substitution variant versus wild type HbHNL shifted more than one million-fold from hydroxynitrile lyase activity (acetone cyanohydrin substrate) to esterase activity (p-nitrophenyl acetate substrate). The contribution of Glu79His to esterase activity was surprising since esterases and lipases contain many different amino acids at this position, including glutamate. Saturation mutagenesis at position 79 showed that 13 of 19 possible amino acid substitutions increased esterase activity, suggesting that removal of glutamate, not addition of histidine, increased esterase activity. Molecular modeling indicates that Glu79 disrupts esterase activity in HbHNL when its negatively charged side chain distorts the orientation of the catalytic histidine. Naturally occurring glutamate at the corresponding location of Candida lipases is uncharged due to other active site differences and does not cause the same distortion. This example of the fine tuning of the same catalytic triad for different types of catalysis by subtle interactions with other active site residues shows how difficult it is to design new catalytic reactions of enzymes.

Increasing the reaction rate of hydroxynitrile lyase from Hevea brasiliensis toward mandelonitrile by copying active site residues from an esterase that accepts aromatic esters.

The natural substrate of hydroxynitrile lyase from rubber tree (HbHNL, Hevea brasiliensis) is acetone cyanohydrin, but synthetic applications usually involve aromatic cyanohydrins such as mandelonitrile. To increase the activity of HbHNL toward this unnatural substrate, we replaced active site residues in HbHNL with the corresponding ones from esterase SABP2 (salicylic acid binding protein 2). Although this enzyme catalyzes a different reaction (hydrolysis of esters), its natural substrate (methyl salicylate) contains an aromatic ring. Three of the eleven single-amino-acid-substitution variants of HbHNL reacted more rapidly with mandelonitrile. The best was HbHNL-L121Y, with a kcat 4.2 times higher and high enantioselectivity. Site-saturation mutagenesis at position 121 identified three other improved variants. We hypothesize that the smaller active site orients the aromatic substrate more productively.