A grand challenge in synthetic biology is to use our current knowledge of RNA science to perform the automatic engineering of completely synthetic sequences encoding functional RNAs in living cells. We report here a fully automated design methodology and experimental validation of synthetic RNA interaction circuits working in a cellular environment. The computational algorithm, based on a physicochemical model, produces novel RNA sequences by exploring the space of possible sequences compatible with predefined structures. We tested our methodology in Escherichia coli by designing several positive riboregulators with diverse structures and interaction models, suggesting that only the energy of formation and the activation energy (free energy barrier to overcome for initiating the hybridization reaction) are sufficient criteria to engineer RNA interaction and regulation in bacteria. The designed sequences exhibit nonsignificant similarity to any known noncoding RNA sequence. Our riboregulatory devices work independently and in combination with transcription regulation to create complex logic circuits. Our results demonstrate that a computational methodology based on first-principles can be used to engineer interacting RNAs with allosteric behavior in living cells.

Motivation: The biological solution for synthesis or remediation of organic compounds using living organisms, particularly bacteria and yeast, has been promoted because of the cost reduction with respect to the non-living chemical approach. In that way, computational frameworks can profit from the previous knowledge stored in large databases of compounds, enzymes and reactions. In addition, the cell behavior can be studied by modeling the cellular context. Results: We have implemented a Monte Carlo algorithm (DESHARKY) that finds a metabolic pathway from a target compound by exploring a database of enzymatic reactions. DESHARKY outputs a biochemical route to the host metabolism together with its impact in the cellular context by using mathematical models of the cell resources and metabolism. Furthermore, we provide the sequence of amino acids for the enzymes involved in the route closest phylogenetically to the considered organism. We provide examples of designed metabolic pathways with their genetic load characterizations. Here, we have used Escherichia coli as host organism. In addition, our bioinformatic tool can be applied for biodegradation or biosynthesis and its performance scales with the database size. Availability: Software, a tutorial and examples are freely available and open source at http://soft.synth-bio.org/desharky.html Contact:alfonso.jaramillo@polytechnique.fr Supplementary information: Supplementary data are available at Bioinformatics online.

Background: Understanding the molecular mechanisms plants have evolved to adapt their biological activities to a constantly changing environment is an intriguing question and one that requires a systems biology approach. Here we present a network analysis of genome-wide expression data combined with reverse-engineering network modeling to dissect the transcriptional control of Arabidopsis thaliana. The regulatory network is inferred by using an assembly of microarray data containing steady-state RNA expression levels from several growth conditions, developmental stages, biotic and abiotic stresses, and a variety of mutant genotypes. Results: We show that the A. thaliana regulatory network has the characteristic properties of hierarchical networks. We successfully applied our quantitative network model to predict the full transcriptome of the plant for a set of microarray experiments not included in the training dataset. We also used our model to analyze the robustness in expression levels conferred by network motifs such as the coherent feed-forward loop. In addition, the meta-analysis presented here has allowed us to identify regulatory and robust genetic structures. Conclusions: These data suggest that A. thaliana has evolved high connectivity in terms of transcriptional regulation among cellular functions involved in response and adaptation to changing environments, while gene networks constitutively expressed or less related to stress response are characterized by a lower connectivity. Taken together, these findings suggest conserved regulatory strategies that have been selected during the evolutionary history of this eukaryote.

Motivation: The rational design of biological networks with prescribed functions is limited to gene circuits of a few genes. Larger networks involve complex interactions with many parameters and the use of automated computational tools can be very valuable. We propose a new tool to design transcriptional networks with targeted behavior that could be used to better understand the design principles of genetic circuits. Results: We have implemented a Simulated Annealing optimization algorithm that explores throughout the space of transcription networks to obtain a specific behavior. The software outputs a transcriptional network with all the corresponding kinetic parameters in SBML format. We provide examples of transcriptional circuits with logical and oscillatory behaviors. Our tool can also be applied to design networks with multiple external input and output genes. Availability: The software, a tutorial manual, parameter sets and examples are freely available at http://synth-bio.yi.org/genetdes.html Contact:Alfonso.Jaramillo@polytechnique.edu

Synthetic biology has seen an explosive growth in the capability of engineering artificial gene circuits from transcription factors (TFs), particularly in bacteria. However, most artificial networks still employ the same core set of TFs (for example LacI, TetR and cI). The TFs mostly function via repression and it is difficult to integrate multiple inputs in promoter logic. Here we present to our knowledge the first set of dual activator-repressor switches for orthogonal logic gates, based on bacteriophage λ cI variants and multi-input promoter architectures. Our toolkit contains 12 TFs, flexibly operating as activators, repressors, dual activator–repressors or dual repressor–repressors, on up to 270 synthetic promoters. To engineer non cross-reacting cI variants, we design a new M13 phagemid-based system for the directed evolution of biomolecules. Because cI is used in so many synthetic biology projects, the new set of variants will easily slot into the existing projects of other groups, greatly expanding current engineering capacities.

An automatic protein design procedure is used to select amino acid sequences that optimize the folding free energy function for a given protein. The only information used in designing the sequences is a set of known backbone structures for each protein, a rotamer library, and a well established classical empirical force field, which relies on basic physical chemical principles that underlie molecular interactions and protein stability, and has not been adjusted to yield native-like sequences. Applying the procedure to 7 different known protein folds, representing a total of 45 different native protein structures, yields ensembles of designed sequences displaying remarkable similarity to their natural counterparts in the protein core, but which are distinctly non-native on the protein surface. We show that natural and designed sequences for a given fold score significantly higher than random sequences against profiles derived from both, designed and natural sequence ensembles. Furthermore, we find that designed sequence profiles can be used to retrieve the native sequences for many of the analyzed proteins using standard PSI-BLAST searches in sequence databases. These findings may have important implications for our understanding the selection pressures operating on natural protein sequences and hold promise for improving fold recognition.