A mutate-and-map protocol for inferring base pairs in structured RNA

Chemical mapping is a widespread technique for structural analysis of nucleic acids in which a molecule's reactivity to different probes is quantified at single-nucleotide resolution and used to constrain structural modeling. This experimental framework has been extensively revisited in the past decade with new strategies for high-throughput read-outs, chemical modification, and rapid data analysis. Recently, we have coupled the technique to high-throughput mutagenesis. Point mutations of a base-paired nucleotide can lead to exposure of not only that nucleotide but also its interaction partner. Carrying out the mutation and mapping for the entire system gives an experimental approximation of the molecules contact map. Here, we give our in-house protocol for this mutate-and-map strategy, based on 96-well capillary electrophoresis, and we provide practical tips on interpreting the data to infer nucleic acid structure.Comment: 22 pages, 5 figure

Genome landscapes and bacteriophage codon usage

Across all kingdoms of biological life, protein-coding genes exhibit unequal usage of synonmous codons. Although alternative theories abound, translational selection has been accepted as an important mechanism that shapes the patterns of codon usage in prokaryotes and simple eukaryotes. Here we analyze patterns of codon usage across 74 diverse bacteriophages that infect E. coli, P. aeruginosa and L. lactis as their primary host. We introduce the concept of a `genome landscape,' which helps reveal non-trivial, long-range patterns in codon usage across a genome. We develop a series of randomization tests that allow us to interrogate the significance of one aspect of codon usage, such a GC content, while controlling for another aspect, such as adaptation to host-preferred codons. We find that 33 phage genomes exhibit highly non-random patterns in their GC3-content, use of host-preferred codons, or both. We show that the head and tail proteins of these phages exhibit significant bias towards host-preferred codons, relative to the non-structural phage proteins. Our results support the hypothesis of translational selection on viral genes for host-preferred codons, over a broad range of bacteriophages.Comment: 9 Color Figures, 5 Tables, 53 Reference

An Electronic Analog of Synthetic Genetic Networks

An electronic analog of a synthetic genetic network known as the repressilator is proposed. The repressilator is a synthetic biological clock consisting of a cyclic inhibitory network of three negative regulatory genes which produces oscillations in the expressed protein concentrations. Compared to previous circuit analogs of the repressilator, the circuit here takes into account more accurately the kinetics of gene expression, inhibition, and protein degradation. A good agreement between circuit measurements and numerical prediction is observed. The circuit allows for easy control of the kinetic parameters thereby aiding investigations of large varieties of potential dynamics

Reduced stability of mRNA secondary structure near the translation-initiation site in dsDNA viruses

<p>Abstract</p> <p>Background</p> <p>Recent studies have demonstrated a selection pressure for reduced mRNA secondary-structure stability near the start codon of coding sequences. This selection pressure can be observed in bacteria, archaea, and eukaryotes, and is likely caused by the requirement of efficient translation initiation in cellular organism.</p> <p>Results</p> <p>Here, we surveyed the complete genomes of 650 dsDNA virus strains for signals of reduced stability of mRNA secondary structure near the start codon. Our analysis included viruses infecting eukaryotic, prokaryotic, and archaeic hosts. We found that many viruses showed evidence for reduced mRNA secondary-structure stability near the start codon. The effect was most pronounced in viruses infecting prokaryotes, but was also observed in viruses infecting eukaryotes and archaea. The reduction in stability generally increased with increasing genomic GC content. For bacteriophage, the reduction was correlated with a corresponding reduction of stability in the phage hosts.</p> <p>Conclusions</p> <p>We conclude that reduced stability of the mRNA secondary structure near the start codon is a common feature for dsDNA viruses, likely driven by the same selective pressures that cause it in cellular organisms.</p

Differential Trends in the Codon Usage Patterns in HIV-1 Genes

Host-pathogen interactions underlie one of the most complex evolutionary phenomena resulting in continual adaptive genetic changes, where pathogens exploit the host's molecular resources for growth and survival, while hosts try to eliminate the pathogen. Deciphering the molecular basis of host–pathogen interactions is useful in understanding the factors governing pathogen evolution and disease propagation. In host-pathogen context, a balance between mutation, selection, and genetic drift is known to maintain codon bias in both organisms. Studies revealing determinants of the bias and its dynamics are central to the understanding of host-pathogen evolution. We considered the Human Immunodeficiency Virus (HIV) type 1 and its human host to search for evolutionary signatures in the viral genome. Positive selection is known to dominate intra-host evolution of HIV-1, whereas high genetic variability underlies the belief that neutral processes drive inter-host differences. In this study, we analyze the codon usage patterns of HIV-1 genomes across all subtypes and clades sequenced over a period of 23 years. We show presence of unique temporal correlations in the codon bias of three HIV-1 genes illustrating differential adaptation of the HIV-1 genes towards the host preferred codons. Our results point towards gene-specific translational selection to be an important force driving the evolution of HIV-1 at the population level

Massively Parallel RNA Chemical Mapping with a Reduced Bias MAP-seq Protocol

Chemical mapping methods probe RNA structure by revealing and leveraging correlations of a nucleotide's structural accessibility or flexibility with its reactivity to various chemical probes. Pioneering work by Lucks and colleagues has expanded this method to probe hundreds of molecules at once on an Illumina sequencing platform, obviating the use of slab gels or capillary electrophoresis on one molecule at a time. Here, we describe optimizations to this method from our lab, resulting in the MAP-seq protocol (Multiplexed Accessibility Probing read out through sequencing), version 1.0. The protocol permits the quantitative probing of thousands of RNAs at once, by several chemical modification reagents, on the time scale of a day using a table-top Illumina machine. This method and a software package MAPseeker (http://simtk.org/home/map_seeker) address several potential sources of bias, by eliminating PCR steps, improving ligation efficiencies of ssDNA adapters, and avoiding problematic heuristics in prior algorithms. We hope that the step-by-step description of MAP-seq 1.0 will help other RNA mapping laboratories to transition from electrophoretic to next-generation sequencing methods and to further reduce the turnaround time and any remaining biases of the protocol.Comment: 22 pages, 5 figure

Eugene – A Domain Specific Language for Specifying and Constraining Synthetic Biological Parts, Devices, and Systems

BACKGROUND: Synthetic biological systems are currently created by an ad-hoc, iterative process of specification, design, and assembly. These systems would greatly benefit from a more formalized and rigorous specification of the desired system components as well as constraints on their composition. Therefore, the creation of robust and efficient design flows and tools is imperative. We present a human readable language (Eugene) that allows for the specification of synthetic biological designs based on biological parts, as well as provides a very expressive constraint system to drive the automatic creation of composite Parts (Devices) from a collection of individual Parts. RESULTS: We illustrate Eugene's capabilities in three different areas: Device specification, design space exploration, and assembly and simulation integration. These results highlight Eugene's ability to create combinatorial design spaces and prune these spaces for simulation or physical assembly. Eugene creates functional designs quickly and cost-effectively. CONCLUSIONS: Eugene is intended for forward engineering of DNA-based devices, and through its data types and execution semantics, reflects the desired abstraction hierarchy in synthetic biology. Eugene provides a powerful constraint system which can be used to drive the creation of new devices at runtime. It accomplishes all of this while being part of a larger tool chain which includes support for design, simulation, and physical device assembly

Synthetic biology: ethical ramifications 2009

During 2007 and 2008 synthetic biology moved from the manifesto stage to research programs. As of 2009, synthetic biology is ramifying; to ramify means to produce differentiated trajectories from previous determinations. From its inception, most of the players in synthetic biology agreed on the need for (a) rationalized design and construction of new biological parts, devices, and systems as well as (b) the re-design of natural biological systems for specified purposes, and that (c) the versatility of designed biological systems makes them suitable to address such challenges as renewable energy, the production of inexpensive drugs, and environmental remediation, as well as providing a catalyst for further growth of biotechnology. What is understood by these goals, however, is diverse. Those assorted understandings are currently contributing to different ramifications of synthetic biology. The Berkeley Human Practices Lab, led by Paul Rabinow, is currently devoting its efforts to documenting and analyzing these ramifications as they emerge

Adjusting Phenotypes by Noise Control

Genetically identical cells can show phenotypic variability. This is often caused by stochastic events that originate from randomness in biochemical processes involving in gene expression and other extrinsic cellular processes. From an engineering perspective, there have been efforts focused on theory and experiments to control noise levels by perturbing and replacing gene network components. However, systematic methods for noise control are lacking mainly due to the intractable mathematical structure of noise propagation through reaction networks. Here, we provide a numerical analysis method by quantifying the parametric sensitivity of noise characteristics at the level of the linear noise approximation. Our analysis is readily applicable to various types of noise control and to different types of system; for example, we can orthogonally control the mean and noise levels and can control system dynamics such as noisy oscillations. As an illustration we applied our method to HIV and yeast gene expression systems and metabolic networks. The oscillatory signal control was applied to p53 oscillations from DNA damage. Furthermore, we showed that the efficiency of orthogonal control can be enhanced by applying extrinsic noise and feedback. Our noise control analysis can be applied to any stochastic model belonging to continuous time Markovian systems such as biological and chemical reaction systems, and even computer and social networks. We anticipate the proposed analysis to be a useful tool for designing and controlling synthetic gene networks

Full design automation of multi-state RNA devices to program gene expression using energy-based optimization

[EN] Small RNAs (sRNAs) can operate as regulatory agents to control protein expression by interaction with the 59 untranslated region of the mRNA. We have developed a physicochemical framework, relying on base pair interaction energies, to design multi-state sRNA devices by solving an optimization problem with an objective function accounting for the stability of the transition and final intermolecular states. Contrary to the analysis of the reaction kinetics of an ensemble of sRNAs, we solve the inverse problem of finding sequences satisfying targeted reactions. We show here that our objective function correlates well with measured riboregulatory activity of a set of mutants. This has enabled the application of the methodology for an extended design of RNA devices with specified behavior, assuming different molecular interaction models based on Watson-Crick interaction. We designed several YES, NOT, AND, and OR logic gates, including the design of combinatorial riboregulators. In sum, our de novo approach provides a new paradigm in synthetic biology to design molecular interaction mechanisms facilitating future high-throughput functional sRNA design.Work supported by the grants FP7-ICT-043338 (BACTOCOM) to AJ, and BIO2011-26741 (Ministerio de Economia y Competitividad, Spain) to JAD. GR is supported by an EMBO long-term fellowship co-funded by Marie Curie actions (ALTF-1177-2011), and TEL by a PhD fellowship from the AXA Research Fund. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.Rodrigo Tarrega, G.; Landrain, TE.; Majer, E.; Daros Arnau, JA.; Jaramillo, A. (2013). Full design automation of multi-state RNA devices to program gene expression using energy-based optimization. PLoS Computational Biology. 9(8):1003172-1003172. https://doi.org/10.1371/journal.pcbi.1003172S1003172100317298Isaacs, F. J., Dwyer, D. J., & Collins, J. J. (2006). RNA synthetic biology. Nature Biotechnology, 24(5), 545-554. doi:10.1038/nbt1208Isaacs, F. J., Dwyer, D. J., Ding, C., Pervouchine, D. D., Cantor, C. R., & Collins, J. J. (2004). Engineered riboregulators enable post-transcriptional control of gene expression. Nature Biotechnology, 22(7), 841-847. doi:10.1038/nbt986Lucks, J. B., Qi, L., Mutalik, V. K., Wang, D., & Arkin, A. P. (2011). Versatile RNA-sensing transcriptional regulators for engineering genetic networks. Proceedings of the National Academy of Sciences, 108(21), 8617-8622. doi:10.1073/pnas.1015741108Mutalik, V. K., Qi, L., Guimaraes, J. C., Lucks, J. B., & Arkin, A. P. (2012). Rationally designed families of orthogonal RNA regulators of translation. Nature Chemical Biology, 8(5), 447-454. doi:10.1038/nchembio.919Bayer, T. S., & Smolke, C. D. (2005). Programmable ligand-controlled riboregulators of eukaryotic gene expression. Nature Biotechnology, 23(3), 337-343. doi:10.1038/nbt1069Nakashima, N., & Tamura, T. (2009). Conditional gene silencing of multiple genes with antisense RNAs and generation of a mutator strain of Escherichia coli. Nucleic Acids Research, 37(15), e103-e103. doi:10.1093/nar/gkp498Callura, J. M., Cantor, C. R., & Collins, J. J. (2012). Genetic switchboard for synthetic biology applications. Proceedings of the National Academy of Sciences, 109(15), 5850-5855. doi:10.1073/pnas.1203808109Beisel, C. L., Bayer, T. S., Hoff, K. G., & Smolke, C. D. (2008). Model‐guided design of ligand‐regulated RNAi for programmable control of gene expression. Molecular Systems Biology, 4(1), 224. doi:10.1038/msb.2008.62Qi, L., Lucks, J. B., Liu, C. C., Mutalik, V. K., & Arkin, A. P. (2012). Engineering naturally occurring trans -acting non-coding RNAs to sense molecular signals. Nucleic Acids Research, 40(12), 5775-5786. doi:10.1093/nar/gks168Carothers, J. M., Goler, J. A., Juminaga, D., & Keasling, J. D. (2011). Model-Driven Engineering of RNA Devices to Quantitatively Program Gene Expression. Science, 334(6063), 1716-1719. doi:10.1126/science.1212209Rodrigo, G., Landrain, T. E., & Jaramillo, A. (2012). De novo automated design of small RNA circuits for engineering synthetic riboregulation in living cells. Proceedings of the National Academy of Sciences, 109(38), 15271-15276. doi:10.1073/pnas.1203831109Brantl, S. (2002). Antisense-RNA regulation and RNA interference. Biochimica et Biophysica Acta (BBA) - Gene Structure and Expression, 1575(1-3), 15-25. doi:10.1016/s0167-4781(02)00280-4Majdalani, N., Vanderpool, C. K., & Gottesman, S. (2005). Bacterial Small RNA Regulators. Critical Reviews in Biochemistry and Molecular Biology, 40(2), 93-113. doi:10.1080/10409230590918702Selinger, D. W., Cheung, K. J., Mei, R., Johansson, E. M., Richmond, C. S., Blattner, F. R., … Church, G. M. (2000). RNA expression analysis using a 30 base pair resolution Escherichia coli genome array. Nature Biotechnology, 18(12), 1262-1268. doi:10.1038/82367Yelin, R., Dahary, D., Sorek, R., Levanon, E. Y., Goldstein, O., Shoshan, A., … Rotman, G. (2003). Widespread occurrence of antisense transcription in the human genome. Nature Biotechnology, 21(4), 379-386. doi:10.1038/nbt808Wang, X.-J., Gaasterland, T., & Chua, N.-H. (2005). Genome Biology, 6(4), R30. doi:10.1186/gb-2005-6-4-r30Stojanovic, M. N., & Stefanovic, D. (2003). A deoxyribozyme-based molecular automaton. Nature Biotechnology, 21(9), 1069-1074. doi:10.1038/nbt862Seelig, G., Soloveichik, D., Zhang, D. Y., & Winfree, E. (2006). Enzyme-Free Nucleic Acid Logic Circuits. Science, 314(5805), 1585-1588. doi:10.1126/science.1132493Yin, P., Choi, H. M. T., Calvert, C. R., & Pierce, N. A. (2008). Programming biomolecular self-assembly pathways. Nature, 451(7176), 318-322. doi:10.1038/nature06451Ran, T., Kaplan, S., & Shapiro, E. (2009). Molecular implementation of simple logic programs. Nature Nanotechnology, 4(10), 642-648. doi:10.1038/nnano.2009.203Penchovsky, R., & Breaker, R. R. (2005). Computational design and experimental validation of oligonucleotide-sensing allosteric ribozymes. Nature Biotechnology, 23(11), 1424-1433. doi:10.1038/nbt1155Salis, H. M., Mirsky, E. A., & Voigt, C. A. (2009). Automated design of synthetic ribosome binding sites to control protein expression. Nature Biotechnology, 27(10), 946-950. doi:10.1038/nbt.1568Laidler, K. J., & King, M. C. (1983). Development of transition-state theory. The Journal of Physical Chemistry, 87(15), 2657-2664. doi:10.1021/j100238a002Sosnick, T. R., & Pan, T. (2003). RNA folding: models and perspectives. Current Opinion in Structural Biology, 13(3), 309-316. doi:10.1016/s0959-440x(03)00066-6Yurke, B. (2003). Genetic Programming and Evolvable Machines, 4(2), 111-122. doi:10.1023/a:1023928811651Bandyra, K. J., Said, N., Pfeiffer, V., Górna, M. W., Vogel, J., & Luisi, B. F. (2012). The Seed Region of a Small RNA Drives the Controlled Destruction of the Target mRNA by the Endoribonuclease RNase E. Molecular Cell, 47(6), 943-953. doi:10.1016/j.molcel.2012.07.015Dawid, A., Cayrol, B., & Isambert, H. (2009). RNA synthetic biology inspired from bacteria: construction of transcription attenuators under antisense regulation. Physical Biology, 6(2), 025007. doi:10.1088/1478-3975/6/2/025007Lioliou, E., Romilly, C., Romby, P., & Fechter, P. (2010). RNA-mediated regulation in bacteria: from natural to artificial systems. New Biotechnology, 27(3), 222-235. doi:10.1016/j.nbt.2010.03.002Dirks, R. M., Bois, J. S., Schaeffer, J. M., Winfree, E., & Pierce, N. A. (2007). Thermodynamic Analysis of Interacting Nucleic Acid Strands. SIAM Review, 49(1), 65-88. doi:10.1137/060651100Das, R., Karanicolas, J., & Baker, D. (2010). Atomic accuracy in predicting and designing noncanonical RNA structure. Nature Methods, 7(4), 291-294. doi:10.1038/nmeth.1433Vogel, J., & Luisi, B. F. (2011). Hfq and its constellation of RNA. Nature Reviews Microbiology, 9(8), 578-589. doi:10.1038/nrmicro2615Friedland, A. E., Lu, T. K., Wang, X., Shi, D., Church, G., & Collins, J. J. (2009). Synthetic Gene Networks That Count. Science, 324(5931), 1199-1202. doi:10.1126/science.1172005Rodrigo, G., Carrera, J., Landrain, T. E., & Jaramillo, A. (2012). Perspectives on the automatic design of regulatory systems for synthetic biology. FEBS Letters, 586(15), 2037-2042. doi:10.1016/j.febslet.2012.02.031Chin, J. W. (2006). Modular approaches to expanding the functions of living matter. Nature Chemical Biology, 2(6), 304-311. doi:10.1038/nchembio789McCaskill, J. S. (1990). The equilibrium partition function and base pair binding probabilities for RNA secondary structure. Biopolymers, 29(6-7), 1105-1119. doi:10.1002/bip.360290621Chitsaz, H., Salari, R., Sahinalp, S. C., & Backofen, R. (2009). A partition function algorithm for interacting nucleic acid strands. Bioinformatics, 25(12), i365-i373. doi:10.1093/bioinformatics/btp212Hofacker, I. L., Fontana, W., Stadler, P. F., Bonhoeffer, L. S., Tacker, M., & Schuster, P. (1994). Fast folding and comparison of RNA secondary structures. Monatshefte f�r Chemie Chemical Monthly, 125(2), 167-188. doi:10.1007/bf00818163Andronescu, M., Zhang, Z. C., & Condon, A. (2005). Secondary Structure Prediction of Interacting RNA Molecules. Journal of Molecular Biology, 345(5), 987-1001. doi:10.1016/j.jmb.2004.10.082Mathews, D. H., Sabina, J., Zuker, M., & Turner, D. H. (1999). Expanded sequence dependence of thermodynamic parameters improves prediction of RNA secondary structure. Journal of Molecular Biology, 288(5), 911-940. doi:10.1006/jmbi.1999.2700Kirkpatrick, S., Gelatt, C. D., & Vecchi, M. P. (1983). Optimization by Simulated Annealing. Science, 220(4598), 671-680. doi:10.1126/science.220.4598.67