DeCoDe: degenerate codon design for complete protein-coding DNA libraries

MOTIVATION

High-throughput protein screening is a critical technique for dissecting and designing protein function. Libraries for these assays can be created through a number of means, including targeted or random mutagenesis of a template protein sequence or direct DNA synthesis. However, mutagenic library construction methods often yield vastly more nonfunctional than functional variants and, despite advances in large-scale DNA synthesis, individual synthesis of each desired DNA template is often prohibitively expensive. Consequently, many protein-screening libraries rely on the use of degenerate codons (DCs), mixtures of DNA bases incorporated at specific positions during DNA synthesis, to generate highly diverse protein-variant pools from only a few low-cost synthesis reactions. However, selecting DCs for sets of sequences that covary at multiple positions dramatically increases the difficulty of designing a DC library and leads to the creation of many undesired variants that can quickly outstrip screening capacity.

RESULTS

We introduce a novel algorithm for total DC library optimization, degenerate codon design (DeCoDe), based on integer linear programming. DeCoDe significantly outperforms state-of-the-art DC optimization algorithms and scales well to more than a hundred proteins sharing complex patterns of covariation (e.g. the lab-derived avGFP lineage). Moreover, DeCoDe is, to our knowledge, the first DC design algorithm with the capability to encode mixed-length protein libraries. We anticipate DeCoDe to be broadly useful for a variety of library generation problems, ranging from protein engineering attempts that leverage mutual information to the reconstruction of ancestral protein states.

AVAILABILITY AND IMPLEMENTATION

github.com/OrensteinLab/DeCoDe.

CONTACT

yaronore@bgu.ac.il.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

An Open-Source, Programmable Pneumatic Setup for Operation and Automated Control of Single- and Multi-Layer Microfluidic Devices

Microfluidic technologies have been used across diverse disciplines (e.g. high-throughput biological measurement, fluid physics, laboratory fluid manipulation) but widespread adoption has been limited in part due to the lack of openly disseminated resources that enable non-specialist labs to make and operate their own devices. Here, we report the open-source build of a pneumatic setup capable of operating both single and multilayer (Quake-style) microfluidic devices with programmable scripting automation. This setup can operate both simple and complex devices with 48 device valve control inputs and 18 sample inputs, with modular design for easy expansion, at a fraction of the cost of similar commercial solutions. We present a detailed step-by-step guide to building the pneumatic instrumentation, as well as instructions for custom device operation using our software, Geppetto, through an easy-to-use GUI for live on-chip valve actuation and a scripting system for experiment automation. We show robust valve actuation with near real-time software feedback and demonstrate use of the setup for high-throughput biochemical measurements on-chip. This open-source setup will enable specialists and novices alike to run microfluidic devices easily in their own laboratories.

BET-seq: Binding energy topographies revealed by microfluidics and high-throughput sequencing

Biophysical models of transcriptional regulation rely on energetic measurements of the binding affinities between transcription factors (TFs) and target DNA binding sites. Historically, assays capable of measuring TF-DNA binding affinities have been relatively low-throughput (measuring ~10³ sequences in parallel) and have required significant specialized equipment, limiting their use to a handful of laboratories. Recently, we developed an experimental assay and analysis pipeline that allows measurement of binding energies between a single TF and up to 10⁶ DNA species in a single experiment (Binding Energy Topography by sequencing, or BET-seq) (Le et al., 2018). BET-seq employs the Mechanically Induced Trapping of Molecular Interactions (MITOMI) platform to purify DNA bound to a TF at equilibrium followed by high coverage sequencing to reveal relative differences in binding energy for each sequence. While we have previously used BET-seq to refine the binding affinity landscapes surrounding high-affinity DNA consensus target sites, we anticipate this technique will be applied in future work toward measuring a wide variety of TF-DNA landscapes. Here, we provide detailed instructions and general considerations for DNA library design, performing BET-seq assays, and analyzing the resulting data.