Search

We present a machine learning-based interpretive framework (whatprot) for analyzing single molecule protein sequencing data produced by fluorosequencing, a recently developed proteomics technology that determines sparse amino acid sequences for many individual peptide molecules in a highly parallelized fashion. Whatprot uses Hidden Markov Models (HMMs) to represent the states of each peptide undergoing the various chemical processes during fluorosequencing, and applies these in a Bayesian classifier, in combination with pre-filtering by a k-Nearest Neighbors (kNN) classifier trained on large volumes of simulated fluorosequencing data. We have found that by combining the HMM based Bayesian classifier with the kNN pre-filter, we are able to retain the benefits of both, achieving both tractable runtimes and acceptable precision and recall for identifying peptides and their parent proteins from complex mixtures, outperforming the capabilities of either classifier on its own. Whatprot's hybrid kNN-HMM approach enables the efficient interpretation of fluorosequencing data using a full proteome reference database and should now also enable improved sequencing error rate estimates., Competing Interests: I have read the journal’s policy and the authors of this manuscript have the following competing interests: E.M.M and Z.B.S. are co-founders and shareholders of Erisyon, Inc., and are co-inventors on granted patents or pending patent applications related to single-molecule protein sequencing. Z.B.S. is an employee of Erisyon, Inc. E.M.M. serves on the scientific advisory board., (Copyright: © 2023 Smith et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Molecular encoding in sequence-defined polymers shows promise as a new paradigm for data storage. Here, we report what is, to our knowledge, the first use of self-immolative oligourethanes for storing and reading encoded information. As a proof of principle, we describe how a text passage from Jane Austen's Mansfield Park was encoded in sequence-defined oligourethanes and reconstructed via self-immolative sequencing. We develop Mol.E-coder, a software tool that uses a Huffman encoding scheme to convert the character table to hexadecimal. The oligourethanes are then generated by a high-throughput parallel synthesis. Sequencing of the oligourethanes by self-immolation is done concurrently in a parallel fashion, and the liquid chromatography-mass spectrometry (LC-MS) information decoded by our Mol.E-decoder software. The passage is capable of being reproduced wholly intact by a third-party, without any purifications or the use of tandem MS (MS/MS), despite multiple rounds of compression, encoding, and synthesis., Competing Interests: DECLARATION OF INTERESTS The authors declare no competing interests.

Psychiatric disorders are characterized by major fluctuations in psychological function over the course of weeks and months, but the dynamic characteristics of brain function over this timescale in healthy individuals are unknown. Here, as a proof of concept to address this question, we present the MyConnectome project. An intensive phenome-wide assessment of a single human was performed over a period of 18 months, including functional and structural brain connectivity using magnetic resonance imaging, psychological function and physical health, gene expression and metabolomics. A reproducible analysis workflow is provided, along with open access to the data and an online browser for results. We demonstrate dynamic changes in brain connectivity over the timescales of days to months, and relations between brain connectivity, gene expression and metabolites. This resource can serve as a testbed to study the joint dynamics of human brain and metabolic function over time, an approach that is critical for the development of precision medicine strategies for brain disorders.

The economic theory of comparative advantage postulates that beneficial trading relationships can be arrived at by two self-interested entities producing the same goods as long as they have opposing relative efficiencies in producing those goods. The theory predicts that upon entering trade, in order to maximize consumption both entities will specialize in producing the good they can produce at higher efficiency, that the weaker entity will specialize more completely than the stronger entity, and that both will be able to consume more goods as a result of trade than either would be able to alone. We extend this theory to the realm of unicellular organisms by developing mathematical models of genetic circuits that allow trading of a common good (specifically, signaling molecules) required for growth in bacteria in order to demonstrate comparative advantage interactions. In Conception 1, the experimenter controls production rates via exogenous inducers, allowing exploration of the parameter space of specialization. In Conception 2, the circuits self-regulate via feedback mechanisms. Our models indicate that these genetic circuits can demonstrate comparative advantage, and that cooperation in such a manner is particularly favored under stringent external conditions and when the cost of production is not overly high. Further work could involve implementing the models in living bacteria and searching for naturally occurring cooperative relationships between bacteria that conform to the principles of comparative advantage., (Copyright © 2014 The Authors. Published by Elsevier Ltd.. All rights reserved.)

We have developed a theoretical framework for developing patterns in multiple dimensions using controllable diffusion and designed reactions implemented in DNA. This includes so-called strand displacement reactions in which one single-stranded DNA hybridizes to a hemi-duplex DNA and displaces another single-stranded DNA, reversibly or irreversibly. These reactions can be designed to proceed with designed rate and molecular specificity. By also controlling diffusion by partial complementarity to a stationary, cross-linked DNA, we can generate predictable patterns. We demonstrate this with several simulations showing deterministic, predictable shapes in space.

DEAD-box helicase proteins accelerate folding and rearrangements of highly structured RNAs and RNA-protein complexes (RNPs) in many essential cellular processes. Although DEAD-box proteins have been shown to use ATP to unwind short RNA helices, it is not known how they disrupt RNA tertiary structure. Here, we use single molecule fluorescence to show that the DEAD-box protein CYT-19 disrupts tertiary structure in a group I intron using a helix capture mechanism. CYT-19 binds to a helix within the structured RNA only after the helix spontaneously loses its tertiary contacts, and then CYT-19 uses ATP to unwind the helix, liberating the product strands. Ded1, a multifunctional yeast DEAD-box protein, gives analogous results with small but reproducible differences that may reflect its in vivo roles. The requirement for spontaneous dynamics likely targets DEAD-box proteins toward less stable RNA structures, which are likely to experience greater dynamic fluctuations, and provides a satisfying explanation for previous correlations between RNA stability and CYT-19 unfolding efficiency. Biologically, the ability to sense RNA stability probably biases DEAD-box proteins to act preferentially on less stable misfolded structures and thereby to promote native folding while minimizing spurious interactions with stable, natively folded RNAs. In addition, this straightforward mechanism for RNA remodeling does not require any specific structural environment of the helicase core and is likely to be relevant for DEAD-box proteins that promote RNA rearrangements of RNP complexes including the spliceosome and ribosome., Competing Interests: The authors have declared that no competing interests exist.

Readily programmable chemical networks are important tools as the scope of chemistry expands from individual molecules to larger molecular systems. Although many complex systems are constructed using conventional organic and inorganic chemistry, the programmability of biological molecules such as nucleic acids allows for precise, high-throughput and automated design, as well as simple, rapid and robust implementation. Here we show that systematic and quantitative control over the diffusivity and reactivity of DNA molecules yields highly programmable chemical reaction networks (CRNs) that execute at the macroscale. In particular, we designed and implemented non-enzymatic DNA circuits capable of performing pattern-transformation algorithms such as edge detection. We also showed that it is possible to fine-tune and multiplex such circuits. We believe these strategies will provide programmable platforms on which to prototype CRNs, discover bottom-up construction principles and generate patterns in materials.

The power of electronic computation is due in part to the development of modular gate structures that can be coupled to carry out sophisticated logical operations and whose performance can be readily modelled. However, the equivalences between electronic and biochemical operations are far from obvious. In order to help cross between these disciplines, we develop an analogy between complementary metal oxide semiconductor and transcriptional logic gates. We surmise that these transcriptional logic gates might prove to be useful in amorphous computations and model the abilities of immobilized gates to form patterns. Finally, to begin to implement these computations, we design unique hairpin transcriptional gates and then characterize these gates in a binary latch similar to that already demonstrated by Kim et al. (Kim, White & Winfree 2006 Mol. Syst. Biol. 2, 68 (doi:10.1038/msb4100099)). The hairpin transcriptional gates are uniquely suited to the design of a complementary NAND gate that can serve as an underlying basis of molecular computing that can output matter rather than electronic information.

Edge detection is a signal processing algorithm common in artificial intelligence and image recognition programs. We have constructed a genetically encoded edge detection algorithm that programs an isogenic community of E. coli to sense an image of light, communicate to identify the light-dark edges, and visually present the result of the computation. The algorithm is implemented using multiple genetic circuits. An engineered light sensor enables cells to distinguish between light and dark regions. In the dark, cells produce a diffusible chemical signal that diffuses into light regions. Genetic logic gates are used so that only cells that sense light and the diffusible signal produce a positive output. A mathematical model constructed from first principles and parameterized with experimental measurements of the component circuits predicts the performance of the complete program. Quantitatively accurate models will facilitate the engineering of more complex biological behaviors and inform bottom-up studies of natural genetic regulatory networks.

We describe a new dynamic kinetic simulation program that allows multiple data sets to be fit simultaneously to a single model based on numerical integration of the rate equations describing the reaction mechanism. Unlike other programs that allow fitting based on numerical integration of rate equations, in the dynamic simulation rate constants, output factors, and starting concentrations of reactants can be scrolled while observing the change in the shape of the simulated reaction curves. Fast dynamic simulation facilitates the exploration of initial parameters that serve as the starting point for nonlinear regression in fitting data and facilitates exploration of the relationships between individual constants and observable reactions. The exploration of parameter space by dynamic simulation provides a powerful tool for learning kinetics and for evaluating the extent to which parameters are constrained by the data. This feature is critical to avoid overly complex models that are not supported by the data.

Fitting several sets of kinetic data directly to a model based on numerical integration provides the best method to extract kinetic parameters without relying on the simplifying assumptions required to achieve analytical solutions of rate equations. However, modern computer programs make it too easy to enter an overly complex model, and standard error analysis grossly underestimates errors when a system is underconstrained and fails to reveal the full degree to which multiple parameters are linked through the complex relationships common in kinetic data. Here we describe the application of confidence contour analysis obtained by measuring the dependence of the sum square error on each pair of parameters while allowing all remaining parameters to be adjusted in seeking the best fit. The confidence contours reveal complex relationships between parameters and clearly outline the space over which parameters can vary (the "FitSpace"). The utility of the method is illustrated by examples of well-constrained fits to published data on tryptophan synthase and the kinetics of oligonucleotide binding to a ribozyme. In contrast, analysis of alanine racemase clearly refutes claims that global analysis of progress curves can be used to extract the free energy profiles of enzyme-catalyzed reactions.

Stochastic fluctuations (noise) in gene expression can cause members of otherwise genetically identical populations to display drastically different phenotypes. An understanding of the sources of noise and the strategies cells employ to function reliably despite noise is proving to be increasingly important in describing the behavior of natural organisms and will be essential for the engineering of synthetic biological systems. Here we describe the design of synthetic constructs, termed ribosome competing RNAs (rcRNAs), as a means to rationally perturb noise in cellular gene expression. We find that noise in gene expression increases in a manner proportional to the ability of an rcRNA to compete for the cellular ribosome pool. We then demonstrate that operons significantly buffer noise between coexpressed genes in a natural cellular background and can even reduce the level of rcRNA enhanced noise. These results demonstrate that synthetic genetic constructs can significantly affect the noise profile of a living cell and, importantly, that operons are a facile genetic strategy for buffering against noise.

We have designed a bacterial system that is switched between different states by red light. The system consists of a synthetic sensor kinase that allows a lawn of bacteria to function as a biological film, such that the projection of a pattern of light on to the bacteria produces a high-definition (about 100 megapixels per square inch), two-dimensional chemical image. This spatial control of bacterial gene expression could be used to 'print' complex biological materials, for example, and to investigate signalling pathways through precise spatial and temporal control of their phosphorylation steps.