Measuring the burden of hundreds of BioBricks defines an evolutionary limit on constructability in synthetic biology.

Engineered DNA will slow the growth of a host cell if it redirects limiting resources or otherwise interferes with homeostasis. Populations of engineered cells can rapidly become dominated by "escape mutants" that evolve to alleviate this burden by inactivating the intended function. Synthetic biologists working with bacteria rely on genetic parts and devices encoded on plasmids, but the burden of different engineered DNA sequences is rarely characterized. We measured how 301 BioBricks on high-copy plasmids affected the growth rate of Escherichia coli. Of these, 59 (19.6%) negatively impacted growth. The burden imposed by engineered DNA is commonly associated with diverting ribosomes or other gene expression factors away from producing endogenous genes that are essential for cellular replication. In line with this expectation, BioBricks exhibiting burden were more likely to contain highly active constitutive promoters and strong ribosome binding sites. By monitoring how much each BioBrick reduced expression of a chromosomal GFP reporter, we found that the burden of most, but not all, BioBricks could be wholly explained by diversion of gene expression resources. Overall, no BioBricks reduced the growth rate of E. coli by >45%, which agreed with a population genetic model that predicts such plasmids should be "unclonable" because escape mutants will take over during growth of a bacterial colony or small laboratory culture from a transformed cell. We made this model available as an interactive web tool for synthetic biology education and added our burden measurements to the iGEM Registry descriptions of each BioBrick.

Bioassay for Determining the Concentrations of Caffeine and Individual Methylxanthines in Complex Samples.

Caffeine and other methylxanthines are stimulant molecules found in formulated beverages, including sodas and energy drinks, and in brewed beverages, such as coffee and teas. Previously, we developed a bioassay for caffeine that involves monitoring the growth of a ΔguaB mutant of Escherichia coli defective in de novo guanine biosynthesis. When supplemented with a plasmid expressing the genes for an N-demethylation pathway from Pseudomonas putida CBB5, these bacteria demethylate caffeine (1,3,7-trimethylxanthine) and other methylxanthines into xanthine, which is then converted into guanine to support cell growth. A major limitation of this bioassay was that it could only measure the total concentration of all methylxanthines in a mixture. Therefore, it could not be used to measure the caffeine content of beverages like teas, which contain substantial quantities of multiple methylxanthines. To overcome this limitation, we created seven new plasmids containing all subsets of the three demethylase genes (ndmA, ndmB, and ndmC). We show that strains of ΔguaBE. coli containing each plasmid are able to demethylate specific subsets of methylxanthines and that they can be used to determine the concentrations of individual methylxanthines in complex mixtures containing multiple methylxanthines, including coffee doped with an additional methylxanthine. While validating this assay, we also discovered an unexpected demethylation event at the 1-methyl position when NdmB and NdmC were expressed in the absence of NdmA. The improved cell-based bioassay is inexpensive, is easy to use, and gives results comparable to standard high-performance liquid chromatography methods for measuring methylxanthine concentrations.IMPORTANCE Caffeine (1,3,7-trimethylxanthine) is the dominant neurostimulant found in coffee, teas, sodas, and energy drinks. Measuring the amount of caffeine and other methylxanthines in these beverages is important for quality assurance and safety in food science. Methylxanthines are also used in medicine and as performance-enhancing drugs, two contexts in which accurately determining their concentrations in bodily fluids is important. Liquid chromatography is the standard method for measuring methylxanthine concentrations in a sample, but it requires specialized equipment and expertise. We improved a previous bioassay that links E. coli growth to methylxanthine demethylation so that it can now be used to determine the amounts of individual methylxanthines in complex mixtures or beverages, such as coffee.

Synthetic Genome Defenses against Selfish DNA Elements Stabilize Engineered Bacteria against Evolutionary Failure.

Mobile genetic elements drive evolution by disrupting genes and rearranging genomes. Eukaryotes have evolved epigenetic mechanisms, including DNA methylation and RNA interference, that silence mobile elements and thereby preserve the integrity of their genomes. We created an artificial reprogrammable epigenetic system based on CRISPR interference to give engineered bacteria a similar line of defense against transposons and other selfish elements in their genomes. We demonstrate that this CRISPR interference against mobile elements (CRISPRi-ME) approach can be used to simultaneously repress two different transposon families in Escherichia coli, thereby increasing the evolutionary stability of costly protein expression. We further show that silencing a transposon in Acinetobacter baylyi ADP1 reduces mutation rates by a factor of 5, nearly as much as deleting all copies of this element from its genome. By deploying CRISPRi-ME on a broad-host-range vector, we have created a generalizable platform for stabilizing the genomes of engineered bacterial cells for applications in metabolic engineering and synthetic biology.

Rapid and Inexpensive Evaluation of Nonstandard Amino Acid Incorporation in Escherichia coli.

By introducing engineered tRNA and aminoacyl-tRNA synthetase pairs into an organism, its genetic code can be expanded to incorporate nonstandard amino acids (nsAAs). The performance of these orthogonal translation systems (OTSs) varies greatly, however, with respect to the efficiency and accuracy of decoding a reassigned codon as the nsAA. To enable rapid and systematic comparisons of these critical parameters, we developed a toolkit for characterizing any Escherichia coli OTS that reassigns the amber stop codon (TAG). It assesses OTS performance by comparing how the fluorescence of strains carrying plasmids encoding a fused RFP-GFP reading frame, either with or without an intervening TAG codon, depends on the presence of the nsAA. We used this kit to (1) examine nsAA incorporation by seven different OTSs, (2) optimize nsAA concentration in growth media, (3) define the polyspecificity of an OTS, and (4) characterize evolved variants of amberless E. coli with improved growth rates.

Automated physics-based design of synthetic riboswitches from diverse RNA aptamers.

Riboswitches are shape-changing regulatory RNAs that bind chemicals and regulate gene expression, directly coupling sensing to cellular actuation. However, it remains unclear how their sequence controls the physics of riboswitch switching and activation, particularly when changing the ligand-binding aptamer domain. We report the development of a statistical thermodynamic model that predicts the sequence-structure-function relationship for translation-regulating riboswitches that activate gene expression, characterized inside cells and within cell-free transcription-translation assays. Using the model, we carried out automated computational design of 62 synthetic riboswitches that used six different RNA aptamers to sense diverse chemicals (theophylline, tetramethylrosamine, fluoride, dopamine, thyroxine, 2,4-dinitrotoluene) and activated gene expression by up to 383-fold. The model explains how aptamer structure, ligand affinity, switching free energy and macromolecular crowding collectively control riboswitch activation. Our model-based approach for engineering riboswitches quantitatively confirms several physical mechanisms governing ligand-induced RNA shape-change and enables the development of cell-free and bacterial sensors for diverse applications.

Predicting the Genetic Stability of Engineered DNA Sequences with the EFM Calculator.

Unwanted evolution can rapidly degrade the performance of genetically engineered circuits and metabolic pathways installed in living organisms. We created the Evolutionary Failure Mode (EFM) Calculator to computationally detect common sources of genetic instability in an input DNA sequence. It predicts two types of mutational hotspots: deletions mediated by homologous recombination and indels caused by replication slippage on simple sequence repeats. We tested the performance of our algorithm on genetic circuits that were previously redesigned for greater evolutionary reliability and analyzed the stability of sequences in the iGEM Registry of Standard Biological Parts. More than half of the parts in the Registry are predicted to experience >100-fold elevated mutation rates due to the inclusion of unstable sequence configurations. We anticipate that the EFM Calculator will be a useful negative design tool for avoiding volatile DNA encodings, thereby increasing the evolutionary lifetimes of synthetic biology devices.

A family of synthetic riboswitches adopts a kinetic trapping mechanism.

Riboswitches are sequences of RNA that control gene expression via RNA-ligand interactions, without the need for accessory proteins. Riboswitches consist of an aptamer that recognizes the ligand and an expression platform that couples ligand binding to a change in gene expression. Using in vitro selection, it is possible to screen large (∼ 10(13) members) libraries of RNA sequences to discover new aptamers. However, limitations in bacterial transformation efficiency make screening such large libraries for riboswitch function in intact cells impractical. Here we show that synthetic riboswitches function in an E. coli S30 extract in a manner similar to how they function in intact E. coli cells. We discovered that, although this family of riboswitches regulates the initiation of protein translation, the fate of whether an RNA message is translated is determined during transcription. Thus, ligand binding does not bias a population of rapidly equilibrating RNA structures, but rather, co-transcriptional ligand binding kinetically traps the RNA in a conformation that supports efficient translation. In addition to providing new insights into the mechanisms of action of a family of synthetic riboswitches, our experiments suggest that it may be possible to perform selections for novel synthetic riboswitches in an in vitro system.

Engineering bacteria to recognize and follow small molecules.

The ability to recognize and react to specific environmental cues allows bacteria to localize to environments favorable to their survival and growth. Synthetic biologists have begun to exploit the chemosensory pathways that control cell motility to reprogram how bacteria move in response to novel signals. Reprograming is often accomplished by designing novel protein or RNA parts that respond to specific small molecules not normally recognized by the natural chemosensory pathways. Additionally, cell motility and localization can be coupled to bacterial quorum sensing, potentially allowing consortia of cells to perform complex tasks.

Flexibility in the site of exon junction complex deposition revealed by functional group and RNA secondary structure alterations in the splicing substrate.

The exon junction complex (EJC) is critical for mammalian nonsense-mediated mRNA decay and translational regulation, but the mechanism of its stable deposition on mRNA is unknown. To examine requirements for EJC deposition, we created splicing substrates containing either DNA nucleotides or RNA secondary structure in the 5' exon. Using RNase H protection, toeprinting, and coimmunoprecipitation assays, we found that EJC location shifts upstream when a stretch of DNA or RNA secondary structure appears at the canonical deposition site. These upstream shifts occur prior to exon ligation and are often accompanied by decreases in deposition efficiency. Although the EJC core protein eIF4AIII contacts four ribose 2'OH groups in crystal structures, we demonstrate that three 2'OH groups are sufficient for deposition. Thus, the site of EJC deposition is more flexible than previously appreciated and efficient deposition appears spatially limited.