An orthogonalized PYR1-based CID module with reprogrammable ligand-binding specificity.

Plants sense abscisic acid (ABA) using chemical-induced dimerization (CID) modules, including the receptor PYR1 and HAB1, a phosphatase inhibited by ligand-activated PYR1. This system is unique because of the relative ease with which ligand recognition can be reprogrammed. To expand the PYR1 system, we designed an orthogonal '*' module, which harbors a dimer interface salt bridge; X-ray crystallographic, biochemical and in vivo analyses confirm its orthogonality. We used this module to create PYR1*MANDI/HAB1* and PYR1*AZIN/HAB1*, which possess nanomolar sensitivities to their activating ligands mandipropamid and azinphos-ethyl. Experiments in Arabidopsis thaliana and Saccharomyces cerevisiae demonstrate the sensitive detection of banned organophosphate contaminants using living biosensors and the construction of multi-input/output genetic circuits. Our new modules enable ligand-programmable multi-channel CID systems for plant and eukaryotic synthetic biology that can empower new plant-based and microbe-based sensing modalities.

Standardizing cassette-based deep mutagenesis by Golden Gate assembly.

Protocols for the construction of large, deeply mutagenized protein encoding libraries via Golden Gate assembly of synthetic DNA cassettes employ disparate, system-specific methodology. Here we present a standardized Golden Gate method for building user-defined libraries. We demonstrate that a 25 µL reaction, using 40 fmol of input DNA, can generate a library on the order of 1 × 106 members and that reaction volume or input DNA concentration can be scaled up with no losses in transformation efficiency. Such libraries can be constructed from dsDNA cassettes generated either by degenerate oligonucleotides or oligo pools. We demonstrate its real-world effectiveness by building custom, user-defined libraries on the order of 104 -107 unique protein encoding variants for two orthogonal protein engineering systems. We include a detailed protocol and provide several general-use destination vectors.

A Closed Form Model for Molecular Ratchet-Type Chemically Induced Dimerization Modules.

Chemical-induced dimerization (CID) modules enable users to implement ligand-controlled cellular and biochemical functions for a number of problems in basic and applied biology. A special class of CID modules occur naturally in plants and involve a hormone receptor that binds a hormone, triggering a conformational change in the receptor that enables recognition by a second binding protein. Two recent reports show that such hormone receptors can be engineered to sense dozens of structurally diverse compounds. As a closed form model for molecular ratchets would be of immense utility in forward engineering of biological systems, here we have developed a closed form model for these distinct CID modules. These modules, which we call molecular ratchets, are distinct from more common CID modules called molecular glues in that they engage in saturable binding kinetics and are characterized well by a Hill equation. A defining characteristic of molecular ratchets is that the sensitivity of the response can be tuned by increasing the molar ratio of the hormone receptor to the binding protein. Thus, the same molecular ratchet can have a pico- or micromolar EC50 depending on the concentration of the different receptor and binding proteins. Closed form models are derived for a base elementary reaction rate model, for ligand-independent complexation of the receptor and binding protein, and for homodimerization of the hormone receptor. Useful governing equations for a variety of in vitro and in vivo applications are derived, including enzyme-linked immunosorbent assay-like microplate assays, transcriptional activation in prokaryotes and eukaryotes, and ligand-induced split protein complementation.

A method for generating user-defined circular single-stranded DNA from plasmid DNA using Golden Gate intramolecular ligation.

Construction of user-defined long circular single stranded DNA (cssDNA) and linear single stranded DNA (lssDNA) is important for various biotechnological applications. Many current methods for synthesis of these ssDNA molecules do not scale to multikilobase constructs. Here we present a robust methodology for generating user-defined cssDNA employing Golden Gate assembly, a nickase, and exonuclease degradation. Our technique is demonstrated for three plasmids with insert sizes ranging from 2.1 to 3.4 kb, requires no specialized equipment, and can be accomplished in 5 h with a yield of 33%-43% of the theoretical. To produce lssDNA, we evaluated different CRISPR-Cas9 cleavage conditions and reported a 52 ± 8% cleavage efficiency of cssDNA. Thus, our current method does not compete with existing protocols for lssDNA generation. Nevertheless, our protocol can make long, user-defined cssDNA readily available to biotechnology researchers.

Feline Interleukin-31 Shares Overlapping Epitopes with the Oncostatin M Receptor and IL-31RA.

Interleukin-31 (IL-31) is a major protein involved in severe inflammatory skin disorders. Its signaling pathway is mediated through two type I cytokine receptors, IL-31RA (also known as the gp130-like receptor) and the oncostatin M receptor (OSMR). Understanding molecular details in these interactions would be helpful for developing antagonist anti-IL-31 monoclonal antibodies (mAbs) as potential therapies. Previous studies suggest that human IL-31 binds to IL-31RA and then recruits OSMR to form a ternary complex. In this model, OSMR cannot interact with IL-31 in the absence of IL-31RA. In this work, we show that feline IL-31 (fIL-31) binds independently with feline OSMR using surface plasmon resonance, an enzyme-linked immunosorbent assay, and yeast surface display. Moreover, competition experiments suggest that OSMR shares a partially overlapping epitope with IL-31RA. We then used deep mutational scanning to map the binding sites of both receptors on fIL-31. In agreement with previous studies of the human homologue, the binding site for IL31-RA contains fIL-31 positions E20 and K82, while the binding site for OSMR comprises the "PADNFERK" motif (P103-K110) and position G38. However, our results also revealed a new overlapping site, composed of positions R69, R72, P73, D76, D81, and E97, between both receptors that we called the "shared site". The conformational epitope of an anti-feline IL-31 mAb that inhibits both OSMR and IL-31RA also mapped to this shared site. Combined, our results show that fIL-31 binds IL-31RA and OSMR independently through a partially shared epitope. These results suggest reexamination of the putative canonical mechanisms for IL-31 signaling in higher animals.

Impact of In Vivo Protein Folding Probability on Local Fitness Landscapes.

It is incompletely understood how biophysical properties like protein stability impact molecular evolution and epistasis. Epistasis is defined as specific when a mutation exclusively influences the phenotypic effect of another mutation, often at physically interacting residues. In contrast, nonspecific epistasis results when a mutation is influenced by a large number of nonlocal mutations. As most mutations are pleiotropic, the in vivo folding probability-governed by basal protein stability-is thought to determine activity-enhancing mutational tolerance, implying that nonspecific epistasis is dominant. However, evidence exists for both specific and nonspecific epistasis as the prevalent factor, with limited comprehensive data sets to support either claim. Here, we use deep mutational scanning to probe how in vivo enzyme folding probability impacts local fitness landscapes. We computationally designed two different variants of the amidase AmiE with statistically indistinguishable catalytic efficiencies but lower probabilities of folding in vivo compared with wild-type. Local fitness landscapes show slight alterations among variants, with essentially the same global distribution of fitness effects. However, specific epistasis was predominant for the subset of mutations exhibiting positive sign epistasis. These mutations mapped to spatially distinct locations on AmiE near the initial mutation or proximal to the active site. Intriguingly, the majority of specific epistatic mutations were codon dependent, with different synonymous codons resulting in fitness sign reversals. Together, these results offer a nuanced view of how protein folding probability impacts local fitness landscapes and suggest that transcriptional-translational effects are as important as stability in determining evolutionary outcomes.

The importance and future of biochemical engineering.

Today's Biochemical Engineer may contribute to advances in a wide range of technical areas. The recent Biochemical and Molecular Engineering XXI conference focused on "The Next Generation of Biochemical and Molecular Engineering: The role of emerging technologies in tomorrow's products and processes". On the basis of topical discussions at this conference, this perspective synthesizes one vision on where investment in research areas is needed for biotechnology to continue contributing to some of the world's grand challenges.

Trade-offs between enzyme fitness and solubility illuminated by deep mutational scanning.

Proteins are marginally stable, and an understanding of the sequence determinants for improved protein solubility is highly desired. For enzymes, it is well known that many mutations that increase protein solubility decrease catalytic activity. These competing effects frustrate efforts to design and engineer stable, active enzymes without laborious high-throughput activity screens. To address the trade-off between enzyme solubility and activity, we performed deep mutational scanning using two different screens/selections that purport to gauge protein solubility for two full-length enzymes. We assayed a TEM-1 beta-lactamase variant and levoglucosan kinase (LGK) using yeast surface display (YSD) screening and a twin-arginine translocation pathway selection. We then compared these scans with published experimental fitness landscapes. Results from the YSD screen could explain 37% of the variance in the fitness landscapes for one enzyme. Five percent to 10% of all single missense mutations improve solubility, matching theoretical predictions of global protein stability. For a given solubility-enhancing mutation, the probability that it would retain wild-type fitness was correlated with evolutionary conservation and distance to active site, and anticorrelated with contact number. Hybrid classification models were developed that could predict solubility-enhancing mutations that maintain wild-type fitness with an accuracy of 90%. The downside of using such classification models is the removal of rare mutations that improve both fitness and solubility. To reveal the biophysical basis of enhanced protein solubility and function, we determined the crystallographic structure of one such LGK mutant. Beyond fundamental insights into trade-offs between stability and activity, these results have potential biotechnological applications.

Plasmid-based one-pot saturation mutagenesis.

Deep mutational scanning is a foundational tool for addressing the functional consequences of large numbers of mutants, but a more efficient and accessible method for construction of user-defined mutagenesis libraries is needed. Here we present nicking mutagenesis, a robust, single-day, one-pot saturation mutagenesis method performed on routinely prepped plasmid dsDNA. The method can be used to produce comprehensive or single- or multi-site saturation mutagenesis libraries.

Characterization of Individual Human Antibodies That Bind Pertussis Toxin Stimulated by Acellular Immunization.

Despite high vaccination rates, the incidence of whooping cough has steadily been increasing in developing countries for several decades. The current acellular pertussis (aP) vaccines all include the major protective antigen pertussis toxin (PTx) and are safer, but they appear to be less protective than infection or older, whole-cell vaccines. To better understand the attributes of individual antibodies stimulated by aP, we isolated plasmablast clones recognizing PTx after booster immunization of two donors. Five unique antibody sequences recognizing native PTx were recovered and expressed as recombinant human IgG1 antibodies. The antibodies all bind different epitopes on the PTx S1 subunit, B oligomer, or S1-B subunit interface, and just one clone neutralized PTx in an in vitro assay. To better understand the epitopes bound by the nonneutralizing S1-subunit antibodies, comprehensive mutagenesis with yeast display provided a detailed map of the epitope recognized by antibodies A8 and E12. Residue R76 is required for antibody A8 binding and is present on the S1 surface but is only partially exposed in the holotoxin, providing a structural explanation for A8's inability to neutralize holotoxin. The B-subunit-specific antibody D8 inhibited PTx binding to a model receptor and neutralized PTx in vitro as well as in an in vivo leukocytosis assay. This is the first study, to our knowledge, to identify individual human antibodies stimulated by the acellular pertussis vaccine and demonstrates the feasibility of using these approaches to address outstanding issues in pertussis vaccinology, including mechanisms of accelerated waning of protective immunity despite repeated aP immunization.

Pro region engineering of nerve growth factor by deep mutational scanning enables a yeast platform for conformational epitope mapping of anti-NGF monoclonal antibodies.

Nerve growth factor (NGF) plays a central role in multiple chronic pain conditions. As such, anti-NGF monoclonal antibodies (mAbs) that function by antagonizing NGF downstream signaling are leading drug candidates for non-opioid pain relief. To evaluate anti-canine NGF (cNGF) mAbs we sought a yeast surface display platform of cNGF. Both mature cNGF and pro-cNGF displayed on the yeast surface but bound conformationally sensitive mAbs at most 2.5-fold in mean fluorescence intensity above background, suggesting that cNGF was mostly misfolded. To improve the amount of folded, displayed cNGF, we used comprehensive mutagenesis, FACS, and deep sequencing to identify point mutants in the pro-region of canine NGF that properly enhance the folded protein displayed on the yeast surface. Out of 1,737 tested single point mutants in the pro region, 49 increased the amount of NGF recognized by conformationally sensitive mAbs. These gain-of-function mutations cluster around residues A-61-P-26. Gain-of-function mutants were additive, and a construct containing three mutations increased amount of folded cNGF to 23-fold above background. Using this new cNGF construct, fine conformational epitopes for tanezumab and three anti-cNGF mAbs were evaluated. The epitope revealed by the yeast experiments largely overlapped with the tanezumab epitope previously determined by X-ray crystallography. The other mAbs showed site-specific differences with tanezumab. As the number of binding epitopes of functionally neutralizing anti-NGF mAbs on NGF are limited, subtle differences in the individual interacting residues on NGF that bind each mAb contribute to the understanding of each antibody and variations in its neutralizing activity. These results demonstrate the potential of deep sequencing-guided protein engineering to improve the production of folded surface-displayed protein, and the resulting cNGF construct provides a platform to map conformational epitopes for other anti-neurotrophin mAbs.

Computational redesign of the lipid-facing surface of the outer membrane protein OmpA.

Advances in computational design methods have made possible extensive engineering of soluble proteins, but designed ß-barrel membrane proteins await improvements in our understanding of the sequence determinants of folding and stability. A subset of the amino acid residues of membrane proteins interact with the cell membrane, and the design rules that govern this lipid-facing surface are poorly understood. We applied a residue-level depth potential for ß-barrel membrane proteins to the complete redesign of the lipid-facing surface of Escherichia coli OmpA. Initial designs failed to fold correctly, but reversion of a small number of mutations indicated by backcross experiments yielded designs with substitutions to up to 60% of the surface that did support folding and membrane insertion.

Fine Epitope Mapping of Two Antibodies Neutralizing the Bordetella Adenylate Cyclase Toxin.

Adenylate cyclase toxin (ACT) is an important Bordetella pertussis virulence factor that is not included in current acellular pertussis vaccines. We previously demonstrated that immunization with the repeat-in-toxin (RTX) domain of ACT elicits neutralizing antibodies in mice and discovered the first two antibodies to neutralize ACT activities by occluding the receptor-binding site. Here, we fully characterize these antibodies and their epitopes. Both antibodies bind ACT with low nanomolar affinity and cross-react with ACT homologues produced by B. parapertussis and B. bronchiseptica. Antibody M1H5 binds B. pertussis RTX751 ∼100-fold tighter than RTX751 from the other two species, while antibody M2B10 has similar affinity for all three variants. To initially map the antibody epitopes, we generated a series of ACT chimeras and truncation variants, which implicated the repeat blocks II-III. To identify individual epitope residues, we displayed randomly mutated RTX751 libraries on yeast and isolated clones with decreased antibody binding by flow cytometry. Next-generation sequencing identified candidate epitope residues on the basis of enrichment of clones with mutations at specific positions. These epitopes form two adjacent surface patches on a predicted structural model of the RTX751 domain, one for each antibody. Notably, the cellular receptor also binds within blocks II-III and shares at least one residue with the M1H5 epitope. The RTX751 model supports the notion that the antibody and receptor epitopes overlap. These data provide insight into mechanisms of ACT neutralization and guidance for engineering more stable RTX variants that may be more appropriate vaccine antigens.

Insights into cellulase-lignin non-specific binding revealed by computational redesign of the surface of green fluorescent protein.

Biological-mediated conversion of pretreated lignocellulosic biomass to biofuels and biochemicals is a promising avenue toward energy sustainability. However, a critical impediment to the commercialization of cellulosic biofuel production is the high cost of cellulase enzymes needed to deconstruct biomass into fermentable sugars. One major factor driving cost is cellulase adsorption and inactivation in the presence of lignin, yet we currently have a poor understanding of the protein structure-function relationships driving this adsorption. In this work, we have systematically investigated the role of protein surface potential on lignin adsorption using a model monomeric fluorescent protein. We have designed and experimentally characterized 16 model protein variants spanning the physiological range of net charge (-24 to +16 total charges) and total charge density (0.28-0.40 charges per sequence length) typical for natural proteins. Protein designs were expressed, purified, and subjected to in silico and in vitro biophysical measurements to evaluate the relationship between protein surface potential and lignin adsorption properties. The designs were comparable to model fluorescent protein in terms of thermostability and heterologous expression yield, although the majority of the designs unexpectedly formed homodimers. Protein adsorption to lignin was studied at two different temperatures using Quartz Crystal Microbalance with Dissipation Monitoring and a subtractive mass balance assay. We found a weak correlation between protein net charge and protein-binding capacity to lignin. No other single characteristic, including apparent melting temperature and 2nd virial coefficient, showed correlation with lignin binding. Analysis of an unrelated cellulase dataset with mutations localized to a family I carbohydrate-binding module showed a similar correlation between net charge and lignin binding capacity. Overall, our study provides strategies to identify highly active, low lignin-binding cellulases by either rational design or by computational screening genomic databases. Biotechnol. Bioeng. 2017;114: 740-750. © 2016 Wiley Periodicals, Inc.

Producing glucose 6-phosphate from cellulosic biomass: structural insights into levoglucosan bioconversion.

The most abundant carbohydrate product of cellulosic biomass pyrolysis is the anhydrosugar levoglucosan (1,6-anhydro-ß-d-glucopyranose), which can be converted to glucose 6-phosphate by levoglucosan kinase (LGK). In addition to the canonical kinase phosphotransfer reaction, the conversion requires cleavage of the 1,6-anhydro ring to allow ATP-dependent phosphorylation of the sugar O6 atom. Using x-ray crystallography, we show that LGK binds two magnesium ions in the active site that are additionally coordinated with the nucleotide and water molecules to result in ideal octahedral coordination. To further verify the metal binding sites, we co-crystallized LGK in the presence of manganese instead of magnesium and solved the structure de novo using the anomalous signal from four manganese atoms in the dimeric structure. The first metal is required for catalysis, whereas our work suggests that the second is either required or significantly promotes the catalytic rate. Although the enzyme binds its sugar substrate in a similar orientation to the structurally related 1,6-anhydro-N-acetylmuramic acid kinase (AnmK), it forms markedly fewer bonding interactions with the substrate. In this orientation, the sugar is in an optimal position to couple phosphorylation with ring cleavage. We also observed a second alternate binding orientation for levoglucosan, and in these structures, ADP was found to bind with lower affinity. These combined observations provide an explanation for the high Km of LGK for levoglucosan. Greater knowledge of the factors that contribute to the catalytic efficiency of LGK can be used to improve applications of this enzyme for levoglucosan-derived biofuel production.

Rapid fine conformational epitope mapping using comprehensive mutagenesis and deep sequencing.

Knowledge of the fine location of neutralizing and non-neutralizing epitopes on human pathogens affords a better understanding of the structural basis of antibody efficacy, which will expedite rational design of vaccines, prophylactics, and therapeutics. However, full utilization of the wealth of information from single cell techniques and antibody repertoire sequencing awaits the development of a high throughput, inexpensive method to map the conformational epitopes for antibody-antigen interactions. Here we show such an approach that combines comprehensive mutagenesis, cell surface display, and DNA deep sequencing. We develop analytical equations to identify epitope positions and show the method effectiveness by mapping the fine epitope for different antibodies targeting TNF, pertussis toxin, and the cancer target TROP2. In all three cases, the experimentally determined conformational epitope was consistent with previous experimental datasets, confirming the reliability of the experimental pipeline. Once the comprehensive library is generated, fine conformational epitope maps can be prepared at a rate of four per day.

Determination of binding affinity upon mutation for type I dockerin-cohesin complexes from Clostridium thermocellum and Clostridium cellulolyticum using deep sequencing.

The comprehensive sequence determinants of binding affinity for type I cohesin toward dockerin from Clostridium thermocellum and Clostridium cellulolyticum was evaluated using deep mutational scanning coupled to yeast surface display. We measured the relative binding affinity to dockerin for 2970 and 2778 single point mutants of C. thermocellum and C. cellulolyticum, respectively, representing over 96% of all possible single point mutants. The interface ΔΔG for each variant was reconstructed from sequencing counts and compared with the three independent experimental methods. This reconstruction results in a narrow dynamic range of -0.8-0.5 kcal/mol. The computational software packages FoldX and Rosetta were used to predict mutations that disrupt binding by more than 0.4 kcal/mol. The area under the curve of receiver operator curves was 0.82 for FoldX and 0.77 for Rosetta, showing reasonable agreements between predictions and experimental results. Destabilizing mutations to core and rim positions were predicted with higher accuracy than support positions. This benchmark dataset may be useful for developing new computational prediction tools for the prediction of the mutational effect on binding affinities for protein-protein interactions. Experimental considerations to improve precision and range of the reconstruction method are discussed. Proteins 2016; 84:1914-1928. © 2016 Wiley Periodicals, Inc.

Removal and upgrading of lignocellulosic fermentation inhibitors by in situ biocatalysis and liquid-liquid extraction.

Hydroxycinnamic acids are known to inhibit microbial growth during fermentation of lignocellulosic biomass hydrolysates, and the ability to diminish hydroxycinnamic acid toxicity would allow for more effective biological conversion of biomass to fuels and other value-added products. In this work, we provide a proof-of-concept of an in situ approach to remove these fermentation inhibitors through constituent expression of a phenolic acid decarboxylase combined with liquid-liquid extraction of the vinyl phenol products. As a first step, we confirmed using simulated fermentation conditions in two model organisms, Escherichia coli and Saccharomyces cerevisiae, that the product 4-vinyl guaiacol is more inhibitory to growth than ferulic acid. Partition coefficients of ferulic acid, p-coumaric acid, 4-vinyl guaiacol, and 4-ethyl phenol were measured for long-chain primary alcohols and alkanes, and tetradecane was identified as a co-solvent that can preferentially extract vinyl phenols relative to the acid parent and additionally had no effect on microbial growth rates or ethanol yields. Finally, E. coli expressing an active phenolic acid decarboxylase retained near maximum anaerobic growth rates in the presence of ferulic acid if and only if tetradecane was added to the fermentation broth. This work confirms the feasibility of donating catabolic pathways into fermentative microorganisms in order to ameliorate the effects of hydroxycinnamic acids on growth rates, and suggests a general strategy of detoxification by simultaneous biological conversion and extraction.

Enhanced T cell receptor specificity through framework engineering.

Development of T cell receptors (TCRs) as immunotherapeutics is hindered by inherent TCR cross-reactivity. Engineering more specific TCRs has proven challenging, as unlike antibodies, improving TCR affinity does not usually improve specificity. Although various protein design approaches have been explored to surmount this, mutations in TCR binding interfaces risk broadening specificity or introducing new reactivities. Here we explored if TCR specificity could alternatively be tuned through framework mutations distant from the interface. Studying the 868 TCR specific for the HIV SL9 epitope presented by HLA-A2, we used deep mutational scanning to identify a framework mutation above the mobile CDR3ß loop. This glycine to proline mutation had no discernable impact on binding affinity or functional avidity towards the SL9 epitope but weakened recognition of SL9 escape variants and led to fewer responses in a SL9-derived positional scanning library. In contrast, an interfacial mutation near the tip of CDR3α that also did not impact affinity or functional avidity towards SL9 weakened specificity. Simulations indicated that the specificity-enhancing mutation functions by reducing the range of loop motions, limiting the ability of the TCR to adjust to different ligands. Although our results are likely to be TCR dependent, using framework engineering to control TCR loop motions may be a viable strategy for improving the specificity of TCR-based immunotherapies.

Rationalizing diverse binding mechanisms to the same protein fold: insights for ligand recognition and biosensor design.

The engineering of novel protein-ligand binding interactions, particularly for complex drug-like molecules, is an unsolved problem which could enable many practical applications of protein biosensors. In this work, we analyzed two engineer ed biosensors, derived from the plant hormone sensor PYR1, to recognize either the agrochemical mandipropamid or the synthetic cannabinoid WIN55,212-2. Using a combination of quantitative deep mutational scanning experiments and molecular dynamics simulations, we demonstrated that mutations at common positions can promote protein-ligand shape complementarity and revealed prominent differences in the electrostatic networks needed to complement diverse ligands. MD simulations indicate that both PYR1 protein-ligand complexes bind a single conformer of their target ligand that is close to the lowest free energy conformer. Computational design using a fixed conformer and rigid body orientation led to new WIN55,212-2 sensors with nanomolar limits of detection. This work reveals mechanisms by which the versatile PYR1 biosensor scaffold can bind diverse ligands. This work also provides computational methods to sample realistic ligand conformers and rigid body alignments that simplify the computational design of biosensors for novel ligands of interest.