TY - JOUR AU - Rondelez, Yannick AB - Abstract In vitro molecular circuits, based on DNA-programmable chemistries, can perform an increasing range of high-level functions, such as molecular level computation, image or chemical pattern recognition and pattern generation. Most reported demonstrations, however, can only accept nucleic acids as input signals. Real-world applications of these programmable chemistries critically depend on strategies to interface them with a variety of non-DNA inputs, in particular small biologically relevant chemicals. We introduce here a general strategy to interface DNA-based circuits with non-DNA signals, based on input-translating modules. These translating modules contain a DNA response part and an allosteric protein sensing part, and use a simple design that renders them fully tunable and modular. They can be repurposed to either transmit or invert the response associated with the presence of a given input. By combining these translating-modules with robust and leak-free amplification motifs, we build sensing circuits that provide a fluorescent quantitative time-response to the concentration of their small-molecule input, with good specificity and sensitivity. The programmability of the DNA layer can be leveraged to perform DNA based signal processing operations, which we demonstrate here with logical inversion, signal modulation and a classification task on two inputs. The DNA circuits are also compatible with standard biochemical conditions, and we show the one-pot detection of an enzyme through its native metabolic activity. We anticipate that this sensitive small-molecule-to-DNA conversion strategy will play a critical role in the future applications of molecular-level circuitry. INTRODUCTION The field of molecular programming (MP) develops systematic methods for the synthesis of complex molecular systems with desired behaviours or functions, in vitro (1). In particular, the use of synthetic DNA oligonucleotides as information-encoding molecules has allowed fast progress in our ability to construct multi-species reaction networks. These systems are able to emulate, in vitro, electronic circuits (2,3), cellular regulation networks (4,5), ecological webs (6,7) or neural architectures (8). The modularity and chemical simplicity of this approach allow rational and scalable circuit assembly. Analog molecular computation, for example, was demonstrated by Song et al. in the case of polynomial functions (9). One critical limit of most in vitro DNA-based approaches, however, is that they only accept inputs presented as nucleic acid species. In contrast, natural networks, as well as many in vivo synthetic circuits, accept a variety of chemicals, including small molecules, as inputs. For example, the historically important lac operon is a small regulatory circuit that optimises the metabolic response of a bacteria on a 2-dimensional glucose/lactose concentration space (10). Synthetic biology also builds on the rich catalogue of regulatory elements, such as allosteric regulators or riboswitches, to engineer a variety of circuits responding to small-molecule inputs (11,12). Following the biological inspiration, early studies explored two strategies, a direct or an indirect one, to connect nucleic acid reactions networks and small molecule chemistry. The direct strategy leverages specific, sequence-dependent DNA or RNA folds, termed aptamers (13,14), that display tailored binding affinity for a given target molecule (15). Inserted in the synthetic strand design, these elements can modulate the circuit's reactivity according to their ligand's concentration (16). The indirect approach uses an intermediate molecule to bridge nucleic and small molecule chemistries. A prominent example is provided by some natural or engineered (17) transcription factors (TFs), which recognise and bind specific sequences in double-stranded DNA, and whose binding affinity is strongly modulated by small effectors. Both approaches have allowed the building of simple in vitro sensors (18–20). For example, L-tryptophan has been detected via the release of TF-induced inhibition of the rolling circle amplification (RCA) on a template displaying the operator site (18). Using an aptamer Oishi et al. designed an ATP-sensitive circuit embedded in a DNA gel (21). However, little is known concerning the potential of TFs to provide a general interface between small molecule stimuli and DNA-based molecular circuits. In this work, we set out to develop a versatile approach that both leverages state-of-the-art capabilities in terms of MP, in particular neural architectures (8), and connects it to physiologically relevant concentrations of small molecule inputs (i.e. micro to millimolar). We show the benefit of the molecular signal-processing layer, which we use to invert the signal, or generate multi-input classification. In addition, we adapted the system to work in mild and bio-compatible conditions, thereby increasing the range of possible applications of these circuits. We demonstrate this last advantage by building a DNA circuit that senses the activity of enzymes, whose catalytic action is in principle completely orthogonal to the chemistry of nucleic acids. MATERIALS AND METHODS Oligonucleotides High-performance liquid chromatography-purified DNA oligonucleotides were obtained from biomers.net. The sequences are listed in Supplementary Table S2. Most templates were modified at their 5′-end by phosphorothioate bonds and at their 3′-end by a phosphate group to respectively prevent their degradation by the exonuclease and their extension by the polymerase. An exception are the protein sensing templates, which were not susceptible to the issue. Reporter template (rT) was conjugated to Atto663 and BHQ2 at 5′- and 3′-ends (instead of the phosphate group), for fluorescent reporting. Thermodynamic parameters were calculated using the DINAMelt server. Protein purification LacI and TrpR CDSs were cloned from NEB 5-alpha competent cells via Gibson assembly and expressed in NEB T7 express on pIVEX-based vector containing C-terminal His-tag kindly provided by Rémi Sieskind. Purification was carried out at 4°C on GE Healthcare Life Sciences His GraviTrap TALON columns and stored in 20 mM Tris–HCl pH 8.0, 0.5 M NaCl, 50% glycerol prior to use. TrpB variants were purified in a similar manner but using KRX cells. The PfB gene was ordered as a gBlock and NEB KLD kit was used to generate point mutations in genes from both organisms (See Supplementary Note 6 for cloning materials). Enzymatic activity of EcB L-trp synthesis was confirmed through HPLC evaluation. Briefly, 200 mM phosphate buffer pH 8.0 with 3 mM l-serine, 2.4 mM indole, 40 μM pyridoxal 5'-phosphate (PLP) was prepared. 5.6 μM EcB was added for an overnight incubation at 37°C along with negative control that contained no enzyme and was not incubated. Enzyme was filtered out using Amicon Ultra 0.5 ml centrifugal filters MWCO 3kDa. Samples were run on 2.1 × 100 mm C-18 silica column with acetonitrile/water: 0% acetonitrile for 1 min, 0% to 100% over 5 min, 100% for 4 min and then a wash with water for 10 min before the following sample. Using absorbance measurements at 290 nm, the negative control only contained a peak at 5.36 min, corresponding to the indole standard, whilst overnight incubation produced a single peak at 3.24 min, corresponding to L-trp standard. Reaction assembly Reactions were performed in a buffer containing 20 mM Tris–HCl pH8.9, 10 mM (NH4)2SO4, 40 mM KCl, 10 mM MgSO4, 10 mM NaCl, 2 μM Netropsin (Sigma-Aldrich), 200 μg ml–1 BSA (New England Biolabs, NEB), 0.1% Synperonic F108 (Sigma-Aldrich), 0.4× EvaGreen dye (Biotium) and 25 μM of each dNTP. All experiments also contained Bst DNA polymerase, large fragment, the nicking endonucleases Nb.BsmI and Nt.BstNBI (all from NEB) respectively used at 20, 200 and 0.5 U ml–1 (0.25, 2 and 0.005% final dilutions of the commercial stock solutions). The thermophilic 5′ → 3′ exonuclease ttRecJ was purified in the laboratory, stored in Diluent A (NEB)+ 0.1% Triton X-100 at a concentration of 1.53 μM and used at 23 nM. Where relevant, LacI was added at 110 nM and TrpR at 210 nM of dimer. Reactions were assembled in a total volume of 10 μl and run at 37°C in a CFX96 real-time PCR detection system (Biorad). As reported for polymerase/nickase amplification schemes (27,72), we sometimes observed the emergence of template-independent amplification products at longer incubation times. These parasitic reactions—that may occur from both the high or low state of the switches—are easily distinguished from their fluorescent signature and were not taken into account for the analysis. Concentrations of aT and rT were fixed at 50 nM and 10 nM accordingly as in previous works (22–25), sT 8 pM (minimally saturating concentration) for pskT-only circuits, pT at 7 nM, unless otherwise specified whilst for the remainder: Figure 1: F. psTLacI 5 nM; G. psTLacI 2.5 nM. Figure 2: B. psTTrpR 1 nM; D. pT 4 nM, pskTLacI 10nM; F. pT 4 nM, pskTTrpR 10 nM. Figure 3: A. pskTTrpR 10 nM, psTLacI 10 nM; B. psTTrpR 0.13 nM, psTLacI 2.5 nM; C. pT 4 nM, pskTTrpR 1.5 nM, pskTLacI 7.5 nM; D. pT 1 nM, psTTrpR 0.063 nM, pskTLacI 40 nM. Figure 4: 50 nM of rT was used, as well as 1.5 times concentration of Nb.BsmI and 1.8 times concentration of Bst DNA polymerase, 50 nM of TrpR dimer. pT 5nM, sT 2.5 nM, 10 mM of L-serine, 80 μM PLP, 5 mM indole, pskTTrpR 10 nM. Data acquisition and analysis A template conjugated to Atto663 fluorophore on the 5′-end and BHQ2 quencher on 3′-end (rT) was used to specifically monitor switch activation. The signal of the fluorophore is quenched when the template is in a hairpin conformation, which opens upon trigger binding and extension. Nonspecific reactions were monitored from the fluorescence signal produced by a double-stranded DNA intercalating dye (EvaGreen). Fluorescence signals were acquired every minute in the channels corresponding to the two dyes present in solution. ‘Start time’ was derived from Atto663 time traces through a workflow implemented in Excel. In brief, each trace is normalised so that fluorescence values lie between 1 and 0; differences between 2 points 10 min apart are then calculated (last 10 min are excluded) and the lower time-point corresponding to maximum difference value is extracted. The data are manually curated to identify flat traces with no distinctive switch-driven fluorescence increase which are assigned 1/Cq = 0. Where appropriate, Hill Function fitting was performed using GraphPadPrism7. The estimated error determination is based on the maximum single point error for the extraction of Cq values in triplicate experiments performed for Figure 1G, and applied to all other plots. Figure 1. Open in new tabDownload slide PEN-DNA toolbox elements for small molecule sensing. Polymerase, exonuclease and nickases (PEN) are used to drive a series of modules. (A) The autocatalytic module where ssDNA signal B11 binds the dual-repeat template aT, is elongated and nicked, thus releasing two identical triggers. (B) The pseudo-template module where the signal binds pT, is extended with a TT 3′-end (orange), and released. The TT-tailed signal cannot initiate an aT cycle and will eventually be degraded. The pT reaction is faster than the aT one, but has limited capacity, leading to an adjustable threshold effect (C) Reporting module where the signal is extended over the hairpin structure, therefore dequenching a fluorophore, giving rise to fluorescence. (D) The combination of aT, pT and rT provides a bistable amplification switch, whose state is controlled by the concentration of B11 signal. (E) TF (e.g. LacI) occludes the Nt.BstNBI nicking site on psT blocking its ability to produce signal via sequential nicking and polymerisation cycles. Inducer addition results in TF unbinding and signal production. (F) Time-response to a range of LacI dimer added to the switch with psTLacI. Inset – normalised fluorescent time traces. (G) Time-response to a range of IPTG (circles) and Lactose (squares) added to the switch with psTLacI and LacI mixture. Figure 1. Open in new tabDownload slide PEN-DNA toolbox elements for small molecule sensing. Polymerase, exonuclease and nickases (PEN) are used to drive a series of modules. (A) The autocatalytic module where ssDNA signal B11 binds the dual-repeat template aT, is elongated and nicked, thus releasing two identical triggers. (B) The pseudo-template module where the signal binds pT, is extended with a TT 3′-end (orange), and released. The TT-tailed signal cannot initiate an aT cycle and will eventually be degraded. The pT reaction is faster than the aT one, but has limited capacity, leading to an adjustable threshold effect (C) Reporting module where the signal is extended over the hairpin structure, therefore dequenching a fluorophore, giving rise to fluorescence. (D) The combination of aT, pT and rT provides a bistable amplification switch, whose state is controlled by the concentration of B11 signal. (E) TF (e.g. LacI) occludes the Nt.BstNBI nicking site on psT blocking its ability to produce signal via sequential nicking and polymerisation cycles. Inducer addition results in TF unbinding and signal production. (F) Time-response to a range of LacI dimer added to the switch with psTLacI. Inset – normalised fluorescent time traces. (G) Time-response to a range of IPTG (circles) and Lactose (squares) added to the switch with psTLacI and LacI mixture. RESULTS AND DISCUSSION Design of the bistable amplification switch We started with the PEN-DNA framework (22), a powerful molecular programming language which has allowed the rational building of a variety of molecular systems, including batch oscillators (23), programmable pattern formation (24), or ultrasensitive multi-input sensors (25). The PEN-DNA employs Polymerase, Exonuclease and Nickase (PEN) enzymes as the chemical machinery. Short protected single-stranded oligos (templates) with an input-output structure provide instructions. Data flows are implemented via the endogenous production, exchange and decay of smaller strands of DNA called signals. The sequence of these signal strands is designed in such a way that (i) the Tm of a single hybridised signal is similar to the reaction temperature, yielding dynamic binding and unbinding to templates and (ii) it carries recognition sequence site enabling productive nicking/polymerisation cycles on the template. Here, we build on three essential modules from the PEN-DNA toolbox: the first one is the autocatalytic module (aT, for autocatalytic template) (Figure 1A), whose encoding oligo consists of a dual repeat of the complement of an arbitrary ‘signal’ sequence. When a cognate trigger binds to an aT, it is extended by the polymerase and the resulting full duplex is nicked, releasing two copies of the signal. aTs thus provide positive feedback loops and exponential amplification of signals. However, as noted in a number of reports (26,27), this motif is intrinsically sensitive to spurious initiation. The second module, a pseudotemplate (pT), provides a solution to this instability issue. pT (Figure 1B) also has a trigger binding site, which is followed by a short poly-T tail (in this paper 2 bases). Additionally, the binding energy of the trigger for pT versus aT is biased in favor of pT by removing a few bases to the 3′ extremity of the aT. At low concentration, triggers thus bind preferentially to pT, are extended by a short tail, which deactivates them. Above a certain concentration threshold, this sink gets saturated and exponential amplification ensues. The combination aT/pT thus creates a bistable switch that can be adapted to sensitive detection of a variety of inputs (25,28,29). The dynamics of these systems can be observed with non-specific dsDNA fluorescent dyes or by a variety of sequence-specific fluorescent reporters based on DNA-dye conjugates. Here we will use a trigger-specific fluorescent probe, we called reporter template (rT) (Figure 1C). It consists of a hairpin-structure DNA with the two ends carrying a quencher and a fluorescent dye. Upon binding and extension of the trigger, the hairpin structure is disrupted, increasing the separation between the two moieties and giving rise to fluorescence. Since the trigger is short, this form is relatively unstable, but it is captured by the polymerase, which extends the trigger along the probe, thus making the unquenching reaction essentially irreversible. We shall refer to this 3-module assembly (aT, pT, rT) as a ‘switch’ throughout the text (Figure 1D). Similar switches have already been used for the sensitive detection of nucleic acids, down to the femtomolar range (25,28). However, these circuits were typically run at temperature ranging from 45 to 50°C, which is not appropriate if one wishes to use them to detect, or in combination with, protein components coming from a mesophilic organism. The first step was therefore to adapt the working temperature of the switches down to 37°C. Although the PEN enzymes are thermophilic, with optimal temperature ranging from 55°C to 70°C, we found that shortening the signal sequence and decreasing its melting temperature to around 26°C could still provide reasonable kinetics. We then optimised the concentration of switch components to achieve robust bistability, yet reasonable threshold of activation (see Supplementary Note 1, Tables S1 and S2 for details). Design of a sensing module for sensitive small biomolecule detection As a first step to connect the PEN-based networks to non-DNA inputs, we designed a new module called the protein sensing template (psT) (Figure 1E). A psT is a hairpin-shaped structure containing in the stem a nickase site followed by two variable elements: a protein binding site (operator) and an output-encoding sequence. For the nickase Nt.BstNBI, the cleavage site is located 4 bp downstream of the recognition site, so it can be designed to overlap with the operator site without sequence constraints (Supplementary Figure S3). In absence of the cognate operator-binding protein, this template acts simply as a constant source of its output sequence, through self-primed extension-and-nicking cycles. Here, the output sequence is the switch's ‘signal’ sequence, so a sufficient concentration of free psT will eventually turn the switch on. When the DNA-binding protein attaches to the psT, its blocks access for the nickase because of steric hindrance. In that configuration, the output production rate is expected to drop, and the switch start-up time to be delayed or abolished (see Supplementary Figure S4 for full circuit design). To validate this strategy we expressed and purified Escherichia coli LacI, a well-characterised transcription factor (30). LacI dimers bind DNA containing one of the three naturally occurring operator sequences with an affinity ranging from 1 pM to 100 nM (31–34). In parallel, we designed the corresponding psT, with the particular output ‘signal’ sequence chosen—‘B11’ (an 11 nucleotide sequence) and LacI O1 sequence—the operator sequence with the lowest Kd. We combined psTLacI with the switch and tested a range of LacI concentrations. As expected, we observed a direct effect on the switch response, with more TF leading to a delayed initiation (Figure 1F). The repressive effect seems to saturate close to 4 nM of LacI dimer, consistent with the concentration of psT (5 nM), and the reported tight affinity of the TF for its operator (35) (∼1 pM). We define the time of switching (Cq) as the time point where the gradient of rT fluorescent trace is steepest. Since in some conditions (strongly repressed psT) the switch does not start at all in the experimental timescale, we find it useful to use 1/Cq values. This allows to represent non-starters as 0 and to assign the highest values to the fastest reactions (see Supplementary Figure S5 for details on data processing). The affinity of many DNA binding proteins to their target site is allosterically modulated by small ligands. The corresponding psT thus also provides an opportunity to interface the circuits to these effectors. For example, the classic inducer isopropyl β-d-1-thiogalactopyranoside (IPTG) can decrease LacI affinity to DNA by over 4 orders of magnitude (35). We therefore tested whether this strategy could be used to connect PEN-DNA reactions to small molecule inputs. A range of IPTG going from 0 to 1000 μM was added to the DNA circuit described above. We observed indeed that the DNA system responds, within a couple of hours, to the various IPTG concentrations (Figure 1G). The behavior is reminiscent of the typical sigmoidal bacterial induction curve for LacI: more IPTG led to smaller Cq, that is, higher switch activity. The apparent Kd, close to 100 μM, appears reasonable, provided the reported one is around 23 μM (35) at pH 9.2, close to our pH 8.9. IPTG is an analogue of naturally occurring allolactose, the true effector of this system, and we achieved similar results with the latter (Supplementary Figure S6). As a control, we used the same concentration range of lactose, also a close chemical analogue – a metabolic precursor of allolactose with no significant allosteric effect on LacI. As expected, the switch circuit did not respond at all to lactose (Figure 1G). The results above thus suggest that the specificity of the natural allosteric transducer (36) has been maintained upon embedding into an artificial in vitro circuit. Using appropriate calibration curves, this isothermal DNA circuit can thus be used to specifically quantify the concentration of a DNA-binding protein, or its small molecule effector. Expanding the range of signal inputs The binding of allosteric transcription factor is either activated or repressed by their ligand. While IPTG, or allolactose disrupts the LacI-DNA complex, other TFs may follow the opposite logic: for example, TrpR, the transcription factor responsible for tryptophan (trp) metabolic regulation, can only bind its operator sequence in presence of L-trp (37,38). To test whether this class of sensing proteins could be used with DNA-based circuits, we purified E. coli TrpR and designed the corresponding psT presenting the TrpR operator site and the same B11 signal as output sequence (Figure 2A). When the reaction was assembled as before, but with TrpR and psTtrpR as input module, we observed that the switch indeed responded to the presence of L-trp. However, the response was inverted compared to the previous case, with more L-trp leading to larger Cq. Here as well, we could confirm the specificity (39) of the sensing circuit using another amino acid as analogue: a concentration of l-threonine (l-thr) as high as 1 mM had no detectable effect on the switch. The Hill curve fitting on the l-trp range yielded a Kd value of 24 μM, again consistent with the reported value of 20 μM (40,41) (Figure 2B). Figure 2. Open in new tabDownload slide PEN-DNA small molecule sensing modules and their dynamics. (A) TF blocks the Nt.BstNBI site on psT upon repressor addition and stops B11 production via sequential nicking and polymerisation cycles. (B) Time-response to a range of L-trp (circles) and L-thr (squares) added to the switch with TrpR and psTTrpR mixture. (C) TF occludes the Nt.BstNBI nicking site on pskT, blocking its ability to produce antiB11 (aB11) via sequential nicking and polymerisation cycles. Inducer addition results in TF unbinding and aB11 production leading to B11 consumption and switch deactivation. (D) Time–response to a range of IPTG added to the switch with pskTLacI and LacI mixture. (E) aB11 is produced from pskT via sequential nicking and polymerisation cycles and the switch is initially deactivated. Repressor addition results in TF binding and aB11 production disruption leading to switch activation. (F) Time-response to a range of l-trp added to the switch with pskTTrpR and TrpR mixture. Figure 2. Open in new tabDownload slide PEN-DNA small molecule sensing modules and their dynamics. (A) TF blocks the Nt.BstNBI site on psT upon repressor addition and stops B11 production via sequential nicking and polymerisation cycles. (B) Time-response to a range of L-trp (circles) and L-thr (squares) added to the switch with TrpR and psTTrpR mixture. (C) TF occludes the Nt.BstNBI nicking site on pskT, blocking its ability to produce antiB11 (aB11) via sequential nicking and polymerisation cycles. Inducer addition results in TF unbinding and aB11 production leading to B11 consumption and switch deactivation. (D) Time–response to a range of IPTG added to the switch with pskTLacI and LacI mixture. (E) aB11 is produced from pskT via sequential nicking and polymerisation cycles and the switch is initially deactivated. Repressor addition results in TF binding and aB11 production disruption leading to switch activation. (F) Time-response to a range of l-trp added to the switch with pskTTrpR and TrpR mixture. In practical applications, being dependent on the TF internal allosteric logic, i.e. its positive or negative response to the molecular stimulus, may seem inconvenient. Fortunately, the PEN toolbox offers straightforward mechanisms for logical inversion (42,43). Here we decided to replace the B11 output of the psT with the complementary sequence antiB11 (aB11). For clarity we named these repressing input modules protein sensing killer Template (pskT), in line with recently described killer Templates (kT) (25) (Figure 2C, E). At the same time, the mixture was complemented with an unconditional mechanism to produce the signal B11, which we call a source template (sT). An sT is similar to a psT but lacks the operator element. Triggering of the switch is then controlled by the opposing effects of sT and pskT, whose balance is shifted by the TF binding on pskT. Figure 2D, F shows the result of inverted LacI and TrpR-based sensors using this approach: the cognate small molecule now acts as inhibitor (respectively activator) of the DNA-based bistable switch. See Supplementary Note 2 and Supplementary Figures S7 and S8 for details about sT, psT and pskT function and design. In these systems, while signal inversion behaved as expected, the apparent Kd values increased: Kd were around 500 and 69 μM, for LacI and TrpR-based systems respectively. We attribute this decrease in performance to the fact that pskT product has to overcome the continuous production of triggers by the sT, in order to maintain the switch in an OFF state. Interestingly however, the competitive mechanism implemented by this additional inversion layer could result in a much sharper transition than the corresponding direct design. In the case of TrpR, we observed a Hill coefficient shifting from 1.9 to 10.5. This observation is reminiscent of similar ultrasensitive behavior in system based on inhibitor titration (44,45), and shows that the DNA circuit can also be used to reshape the functional response of the allosteric element (see also Supplementary Note 3). Circuit integration Transducing small molecule concentration into DNA reactions in principle allows us to deploy the full toolbox of molecular information processing, including logic functions, pattern recognition, signal amplification, etc, to the sensing of metabolites. To demonstrate this idea, we created two-input molecular classifiers, based on IPTG and l-trp. Using the logical inversion mechanism reported above, we designed the four possible 2D systems (Figure 3). Figure 3. Open in new tabDownload slide Combinations of protein sensing (killer) templates with ranges of concentrations of two inputs and LacI with TrpR in all samples. Circuit schematics are shown above the figures. Colour range represents relative 1/Cq variation within individual plots. (A) pskTTrpR and psTLacI for l-trp inverted sensing, IPTG direct sensing. (B) psTTrpR and psTLacI for l-trp direct sensing, IPTG direct sensing. (C) pskTTrpR and pskTLacI for l-trp direct sensing, IPTG inverted sensing. (D) psTTrpR and pskTLacI for l-trp inverted sensing, IPTG inverted sensing. Figure 3. Open in new tabDownload slide Combinations of protein sensing (killer) templates with ranges of concentrations of two inputs and LacI with TrpR in all samples. Circuit schematics are shown above the figures. Colour range represents relative 1/Cq variation within individual plots. (A) pskTTrpR and psTLacI for l-trp inverted sensing, IPTG direct sensing. (B) psTTrpR and psTLacI for l-trp direct sensing, IPTG direct sensing. (C) pskTTrpR and pskTLacI for l-trp direct sensing, IPTG inverted sensing. (D) psTTrpR and pskTLacI for l-trp inverted sensing, IPTG inverted sensing. Combining two input channels into an integrated response, whatever the logic, implies that the time scales and reaction strength are adequately matched. Fortunately, in DNA-based chemical reaction networks these parameters are controlled by the reaction kinetics, and ultimately by the encoding DNA concentrations. We could therefore balance the systems, by adjusting psT/pskT concentrations to change the inputs strength, and pT to adapt switch sensitivity. Given these adjustments, we obtained a series of plots showing rotational symmetry. This tunability is both useful for rapid rewiring of a circuit and for adaptation to different transcription factors. In the plots shown in Figure 3, we tuned the system to obtain an approximate weighted log-sum on the inputs, where the weight can be positive or negative, but other responses are possible. The weight and shape of the contribution of a particular input can be varied, which can be used to prioritise one analyte over another or binarise the response (see Supplementary Note 4, for examples of different 2D sensing functions). This again demonstrates the benefit of the DNA circuit as a versatile signal-processing layer acting to modulate the intrinsic molecular characteristics of the functional components – here, the TFs. In situ metabolic sensing One step upstream of the small molecule world stand the metabolic enzymes that create, consume or modify them. For example, l-trp is the product of Tryptophan synthase subunit B (TrpB), which synthesises it from indole and l-serine. This reaction is usually coupled with indole production in subunit A of the TrpAB dimer, and is slow in the absence of TrpA (46). Standalone TrpB, which could be used in a variety of drug synthesis schemes, has therefore recently been the focus of rich directed evolution studies (47–50) and provided a suitable target. Some of the improved variants reported are Escherichia coli TrpBM149T N171D (EcB from now on) and Pyrococcus furiosus TrpB T292S (from now on PfB T292S, or PfB wt for wildtype). We thus purified these variants and endeavoured to optimise our reaction assembly to accommodate for enzymatic synthesis as well as L-trp product sensing, over a range of enzyme concentrations (Figure 4A). Figure 4. Open in new tabDownload slide Sensing of L-trp production in situ. (A) L-serine and indole are converted to L-trp by one of the TrpB variants listed in the box (from E. coli or P. furiosus), a metabolic activity that is sensed in situ by the switch. (B) Switch response to L-trp production by various concentrations of EcB (circles), PfB T292S (squares), PfB wt (diamonds). Figure 4. Open in new tabDownload slide Sensing of L-trp production in situ. (A) L-serine and indole are converted to L-trp by one of the TrpB variants listed in the box (from E. coli or P. furiosus), a metabolic activity that is sensed in situ by the switch. (B) Switch response to L-trp production by various concentrations of EcB (circles), PfB T292S (squares), PfB wt (diamonds). As we could expect, for equivalent enzyme concentrations and at 37°C, we found that the mesophilic homologue performed best, followed by enhanced P. furiosus variant (Figure 4B). TrpB co-factor pyridoxal 5'-phosphate (PLP), as well as the reactants, do not seem to affect the circuit. Only very high enzyme concentrations (>10 μM) showed a non-specific disruptive effect to the DNA circuit. CONCLUSION This work demonstrates that DNA-based molecular programs can be interfaced with small molecule inputs using engineered modular converter elements. Compared to other small molecule sensing approaches (e.g. optical (51) or electrochemical (52)), which provide physical readouts, our approach uses chemical-level transducers. Once the molecular signal is forwarded to the DNA layer, a wide range of off-the-shelf molecular programming techniques and modules can be used for signal processing in moleculo. Here, we used signal amplification, thresholding, shaping, inversion and integration. This approach is modular and tunable, as DNA reactions can be designed with limited cross-over, between themselves or with other chemistries. It is thus possible to add a branch – or replace one – with the operator sequence and the TF of interest, without the need to redesign the other components. A limit of this design, however, may lie in its dependence on existing allosteric DNA-binding proteins. Fortunately, large libraries of natural and artificial TFs are readily available, due to their common use in a variety of applications (17,53,54). Faulon et al. for example demonstrated the modular use of various TFs for signal integration in bacterial context (55,56), and compiled a catalogue containing over 4600 entries of characterised TFs (57). This is complemented with a CAD tool that suggests enzymatic pathways to convert small molecules of interest to those detectable by a TF (58). The latter is also compatible with our circuit. Transcription factors are usually very specific with respect to their effectors and operators, and could therefore be used to create highly multi-dimensional metabolome profiling networks. Concerning non-biologically relevant molecule sensing, we note that a number of studies have reported a surprising mutational plasticity of natural allosteric TF (59). TrpR for example, can be repurposed from L-trp to specific bromotryptophan sensing with only a few mutations (60). Although we have focused here on small organic molecules, a number of transcription factors are sensitive to metal ions (61), pH (62) or physical signal (63), widening the range of potential inputs. Likewise, as we show, TFs themselves can be targeted for detection (Figure 1F and Supplementary Note 3). In principle, the modularity of the design makes it compatible with design automation, which could leverage the existing databases and computational tools (53,54). Finally, the protein mediated strategy explored here is not limited to transcription factors. Split DNA- or RNA-processing enzymes, for example RNA polymerase (64) or CRISPR–Cas9 (65), can be made sensitive to a variety of ligands, and mediate their presence to a DNA/RNA circuit. Similarly, optogenetics explores strategies to modulate genetic processes using light (66), which could be reused in vitro. In-emulsion molecular programs offer massive assay parallelisation at little extra cost, provided all operations are performed autonomously in the closed microcompartment (67). We therefore envisage the use of our sensors in applications involving high-throughput screening. For example, metabolic profile of a collection of cells distributed in emulsion could be examined on single-cell basis for millions of cells in parallel, and combined with ultrasensitive miRNA quantification (28). Libraries of enzymatic variants could be assayed for activity by distributing the cells containing mutant pathways or individual DNA fragments coupled to their protein products in emulsion. Although this paper explores the chemical interface upstream of the DNA circuit, the downstream connection will prove equally important for ‘real-world’ applications (68). For example, the actuation of enzymes with DNA has been demonstrated (69) and implemented in particular as an output of PEN-DNA modules (70). Another option is the control of genetic replication processes, which open the route to programmable and autonomous directed evolution protocols (71). As a proof of principle, by replacing the L-trp circuit's fluorescent output by cognate primers production, we created a one-pot circuit that connects TrpB activity to polymerase chain reaction (PCR) gene amplification (Supplementary Note 5). Such a circuit brings the possibility of autonomous and programmable selection of metabolic enzymes. Overall, our generic approach makes molecular programming strategies relevant to a new range of biological metabolite signals. We believe this may open innovative application areas, such as directed evolution and point-of-care medical diagnostics, through providing parallelised molecular level computation on molecular inputs. DATA AVAILABILITY The data that support the findings of this study are available from the corresponding author upon request. Sequences for synthetic oligonucleotides and plasmids are provided in the Supporting information. SUPPLEMENTARY DATA Supplementary Data are available at NAR Online. ACKNOWLEDGEMENTS We would like to thank Gabriel Ragala for his help in transcription factor production and experimentally testing them. We would also like to thank Emmanuelle Demey and Joelle Vinh for helpful scientific discussion and LC-UV analysis (SMBP, ESPCI Paris). Author contributions: Y.R. conceived this project. V.S., G.G. and Y.R. designed the experiments. V.S. and C.L. performed the experiments. V.S., G.G. and Y.R. analysed the results and wrote the paper. All authors have given approval to the final version of the manuscript. FUNDING This work was supported by the European Research Council [647275 (ProFF)]; the LC system was supported by DIM Analytics from the Conseil Régional d’Ile de France. Conflict of interest statement. None declared. REFERENCES 1. Hagiya M. From molecular computing to molecular programming . Lecture Notes in Computer Science (including subseries Lecture Notes in Artificial Intelligence and Lecture Notes in Bioinformatics) . 2001 ; 1716 : 89 – 102 . Google Scholar Crossref Search ADS Google Preview WorldCat COPAC 2. Seelig G. , Soloveichik D., Zhang D.Y., Winfree E. Enzyme-free nucleic acid logic circuits . Science . 2006 ; 314 : 1585 – 1588 . Google Scholar Crossref Search ADS PubMed WorldCat 3. Zhang D.Y. , Turberfield A.J., Yurke B., Winfree E. Engineering entropy-driven reactions and networks catalyzed by DNA . Science . 2007 ; 318 : 1121 – 1125 . Google Scholar Crossref Search ADS PubMed WorldCat 4. Genot A.J. , Fujii T., Rondelez Y. In vitro regulatory models for systems biology . Biotechnol. Adv. 2013 ; 31 : 789 – 796 . Google Scholar Crossref Search ADS PubMed WorldCat 5. Tyson J.J. , Novák B. Functional motifs in biochemical reaction networks . Annu. Rev. Phys. Chem. 2010 ; 61 : 219 – 240 . Google Scholar Crossref Search ADS PubMed WorldCat 6. Rondelez Y. Competition for catalytic resources alters biological network dynamics . Phys. Rev. Lett. 2012 ; 108 : 018102 . Google Scholar Crossref Search ADS PubMed WorldCat 7. Genot A.J. , Fujii T., Rondelez Y. Computing with competition in biochemical networks . Phys. Rev. Lett. 2012 ; 109 : 208102 . Google Scholar Crossref Search ADS PubMed WorldCat 8. Qian L. , Winfree E., Bruck J. Neural network computation with DNA strand displacement cascades . Nature . 2011 ; 475 : 368 – 372 . Google Scholar Crossref Search ADS PubMed WorldCat 9. Song T. , Garg S., Mokhtar R., Bui H., Reif J. Analog computation by DNA strand displacement circuits . ACS Synth. Biol. 2016 ; 5 : 898 – 912 . Google Scholar Crossref Search ADS PubMed WorldCat 10. Jacob F. , Monod J. Genetic regulatory mechanisms in the synthesis of proteins . J. Mol. Biol. 1961 ; 3 : 318 – 356 . Google Scholar Crossref Search ADS PubMed WorldCat 11. Westbrook A.M. , Lucks J.B. Achieving large dynamic range control of gene expression with a compact RNA transcription-translation regulator . Nucleic Acids Res. 2017 ; 45 : 5614 – 5624 . Google Scholar Crossref Search ADS PubMed WorldCat 12. Chen Y. , Kim J.K., Hirning A.J., Josić K., Bennett M.R. Emergent genetic oscillations in a synthetic microbial consortium . Science . 2015 ; 349 : 986 – 989 . Google Scholar Crossref Search ADS PubMed WorldCat 13. Tuerk C. , Gold L. Systematic evolution of ligands by exponential enrichment: RNA ligands to bacteriophage T4 DNA polymerase . Science . 1990 ; 249 : 505 – 510 . Google Scholar Crossref Search ADS PubMed WorldCat 14. Ellington A.D. , Szostak J.W. In vitro selection of RNA molecules that bind specific ligands . Nature . 1990 ; 346 : 818 – 822 . Google Scholar Crossref Search ADS PubMed WorldCat 15. Porchetta A. , Vallée-Bélisle A., Plaxco K.W., Ricci F. Using distal-site mutations and allosteric inhibition to tune, extend, and narrow the useful dynamic range of aptamer-based sensors . J. Am. Chem. Soc. 2012 ; 134 : 20601 – 20604 . Google Scholar Crossref Search ADS PubMed WorldCat 16. Fabry-Wood A. , Fetrow M.E., Oloyede A., Yang K.A., Stojanovic M.N., Stefanovic D., Graves S.W., Carroll N.J., Lakin M.R. Microcompartments for Protection and Isolation of Nanoscale DNA Computing Elements . ACS Appl. Mater. Interfaces . 2019 ; 11 : 11262 – 11269 . Google Scholar Crossref Search ADS PubMed WorldCat 17. Li J.W. , Zhang X.Y., Wu H., Bai Y.P. Transcription factor engineering for high-throughput strain evolution and organic acid bioproduction: a review . Front. Bioeng. Biotechnol. 2020 ; 8 : 98 . Google Scholar Crossref Search ADS PubMed WorldCat 18. Zhao G. , Hu T., Li J., Wei H., Shang H., Guan Y. A novel strategy to analyze l-tryptophan through allosteric Trp repressor based on rolling circle amplification . Biosens. Bioelectron. 2015 ; 71 : 103 – 107 . Google Scholar Crossref Search ADS PubMed WorldCat 19. Yao Y. , Li S., Cao J., Liu W., Fan K., Xiang W., Yang K., Kong D., Wang W. Development of small molecule biosensors by coupling the recognition of the bacterial allosteric transcription factor with isothermal strand displacement amplification . Chem. Commun. 2018 ; 54 : 4774 – 4777 . Google Scholar Crossref Search ADS WorldCat 20. Song S. , Wang L., Li J., Fan C., Zhao J. Aptamer-based biosensors . TrAC - Trends Anal. Chem. 2008 ; 27 : 108 – 117 . Google Scholar Crossref Search ADS WorldCat 21. Yang X. , Bing T., Mei H., Fang C., Cao Z., Shangguan D. Characterization and application of a DNA aptamer binding to l-tryptophan . Analyst . 2011 ; 136 : 577 – 585 . Google Scholar Crossref Search ADS PubMed WorldCat 22. Baccouche A. , Montagne K., Padirac A., Fujii T., Rondelez Y. Dynamic DNA-toolbox reaction circuits: a walkthrough . Methods . 2014 ; 67 : 234 – 249 . Google Scholar Crossref Search ADS PubMed WorldCat 23. Fujii T. , Rondelez Y. Predator - Prey molecular ecosystems . ACS Nano . 2013 ; 7 : 27 – 34 . Google Scholar Crossref Search ADS PubMed WorldCat 24. Rondelez Y. , Urtel G., Galas J.-C., Estevez-Torres A., Zadorin A.S., Gines G., Zambrano A., Dilhas V. Synthesis and materialization of a reaction–diffusion French flag pattern . Nat. Chem. 2017 ; 9 : 990 – 996 . Google Scholar PubMed OpenURL Placeholder Text WorldCat 25. Rondelez Y. , Gines G. Multiplex digital microRNA detection using cross-inhibitory DNA circuits . ACS sensors . 2020 ; 5 : 2430 – 2437 . Google Scholar Crossref Search ADS PubMed WorldCat 26. Ogata N. , Miura T. Creation of genetic information by DNA polymerase of the thermophilic bacterium Thermus thermophilus . Nucleic Acids Res. 1998 ; 26 : 4657 – 4661 . Google Scholar Crossref Search ADS PubMed WorldCat 27. Irvine B. , Voelkerding K., Tan E., Niemz A., Buechel M., Dames S., Erwin B., Ferguson T. Specific versus nonspecific isothermal DNA amplification through thermophilic polymerase and nicking enzyme activities . Biochemistry . 2008 ; 47 : 9987 – 9999 . Google Scholar PubMed OpenURL Placeholder Text WorldCat 28. Gines G. , Menezes R., Nara K., Kirstetter A.S., Taly V., Rondelez Y. Isothermal digital detection of microRNAs using background-free molecular circuit . Sci. Adv. 2020 ; 6 : eaay5952 . Google Scholar Crossref Search ADS PubMed WorldCat 29. Montagne K. , Gines G., Fujii T., Rondelez Y. Boosting functionality of synthetic DNA circuits with tailored deactivation . Nat. Commun. 2016 ; 7 : 13474 . Google Scholar Crossref Search ADS PubMed WorldCat 30. Lewis M. A tale of two repressors—a historical perspective . J. Mol. Biol. 2011 ; 409 : 14 – 27 . Google Scholar Crossref Search ADS PubMed WorldCat 31. Gilbert W. , Muller-Hill B. Isolation of the Lac Repressor . Proc. Natl. Acad. Sci. 1966 ; 56 : 1891 – 1898 . Google Scholar Crossref Search ADS WorldCat 32. Gilbert W. , Maxam A. The nucleotide sequence of the lac operator . Proc. Natl. Acad. Sci. U.S.A. 1973 ; 70 : 3581 – 3584 . Google Scholar Crossref Search ADS PubMed WorldCat 33. Platt T. , Files J.G., Weber K. Lac repressor . J. Biol. Chem. 1973 ; 248 : 110 – 121 . Google Scholar Crossref Search ADS PubMed WorldCat 34. Romanuka J. , Folkers G.E., Biris N., Tishchenko E., Wienk H., Bonvin A.M.J.J., Kaptein R., Boelens R. Specificity and affinity of lac repressor for the auxiliary operators O2 and O3 are explained by the structures of their protein – DNA complexes . J. Mol. Biol. 2009 ; 390 : 478 – 489 . Google Scholar Crossref Search ADS PubMed WorldCat 35. Falcon C.M. , Matthews K.S. Glycine insertion in the hinge region of lactose repressor protein alters DNA binding . J. Biol. Chem. 1999 ; 274 : 30849 – 30857 . Google Scholar Crossref Search ADS PubMed WorldCat 36. Matthews K.S. , Nichols J.C. Lactose repressor protein: functional properties and structure . Prog. Nucleic Acid Res. Mol. Biol. 1997 ; 58 : 127 – 150 . Google Scholar Crossref Search ADS WorldCat 37. Zhang R. , Lawson C.L., Otwinowski Z., Sigler P.B., Schevitz R.W., Joachimiak A. The crystal structure of trp aporepressor at 1.8 Å shows how binding tryptophan enhances DNA affinity . Nature . 1987 ; 327 : 591 – 597 . Google Scholar Crossref Search ADS PubMed WorldCat 38. Arvidson D.N. , Bruce C., Gunsalus R.P. Interaction of the Escherichia coli trp aporepressor with its ligand, L-tryptophan . J. Biol. Chem. 1986 ; 261 : 238 – 243 . Google Scholar Crossref Search ADS PubMed WorldCat 39. Marmorstein R.Q. , Sigler P.B. Stereochemistry effects of L-tryptophan and its analogues on trp repressor's affinity for operator-DNA . J. Biol. Chem. 1989 ; 264 : 9149 – 9154 . Google Scholar Crossref Search ADS PubMed WorldCat 40. Joachimiak A. , Kelley R.L., Gunsalus R.P., Yanofsky C., Sigler P.B. Purification and characterization of trp aporepressor . Proc. Natl. Acad. Sci. U.S.A. 1983 ; 80 : 668 – 672 . Google Scholar Crossref Search ADS PubMed WorldCat 41. Hurlburt B.K. , Yanofsky C. trp Repressor/trp operator interaction. Equilibrium and kinetic analysis of complex formation and stability . J. Biol. Chem. 1992 ; 267 : 16783 – 16789 . Google Scholar Crossref Search ADS PubMed WorldCat 42. Montagne K. , Plasson R., Sakai Y., Fujii T., Rondelez Y. Programming an in vitro DNA oscillator using a molecular networking strategy . Mol. Syst. Biol. 2011 ; 7 : 466 . Google Scholar Crossref Search ADS PubMed WorldCat 43. Baccouche A. , Montagne K., Padirac A., Fujii T., Rondelez Y. Dynamic DNA-toolbox reaction circuits: a walkthrough . Methods . 2014 ; 67 : 234 – 249 . Google Scholar Crossref Search ADS PubMed WorldCat 44. Buchler N.E. , Cross F.R. Protein sequestration generates a flexible ultrasensitive response in a genetic network . Mol. Syst. Biol. 2009 ; 5 : 272 . Google Scholar Crossref Search ADS PubMed WorldCat 45. Lu M.S. , Mauser J.F., Prehoda K.E. Ultrasensitive synthetic protein regulatory networks using mixed decoys . ACS Synth. Biol. 2012 ; 1 : 65 – 72 . Google Scholar Crossref Search ADS PubMed WorldCat 46. Dunn M.F. Allosteric regulation of substrate channeling and catalysis in the tryptophan synthase bienzyme complex . Arch. Biochem. Biophys. 2012 ; 519 : 154 – 166 . Google Scholar Crossref Search ADS PubMed WorldCat 47. Buller A.R. , Brinkmann-Chen S., Romney D.K., Herger M., Murciano-Calles J., Arnold F.H. Directed evolution of the tryptophan synthase β-subunit for stand-alone function recapitulates allosteric activation . Proc. Natl. Acad. Sci. U.S.A. 2015 ; 112 : 14599 – 14604 . Google Scholar Crossref Search ADS PubMed WorldCat 48. Herger M. , Van Roye P., Romney D.K., Brinkmann-Chen S., Buller A.R., Arnold F.H. Synthesis of β-branched tryptophan analogues using an engineered subunit of tryptophan synthase . J. Am. Chem. Soc. 2016 ; 138 : 8388 – 8391 . Google Scholar Crossref Search ADS PubMed WorldCat 49. Murciano-Calles J. , Romney D.K., Brinkmann-Chen S., Buller A.R., Arnold F.H. A panel of TrpB biocatalysts derived from tryptophan synthase through the transfer of mutations that mimic allosteric activation . Angew. Chemie - Int. Ed. 2016 ; 55 : 11577 – 11581 . Google Scholar Crossref Search ADS WorldCat 50. Buller A.R. , Van Roye P., Murciano-Calles J., Arnold F.H. Tryptophan synthase uses an atypical mechanism to achieve substrate specificity . Biochemistry . 2016 ; 55 : 7043 – 7046 . Google Scholar Crossref Search ADS PubMed WorldCat 51. Sutariya P.G. , Pandya A., Lodha A., Menon S.K. A unique fluorescence biosensor for selective detection of tryptophan and histidine . Analyst . 2014 ; 139 : 4794 – 4798 . Google Scholar Crossref Search ADS PubMed WorldCat 52. Bagheri Hashkavayi A. , Raoof J.B. Ultrasensitive and reusable electrochemical aptasensor for detection of tryptophan using of [Fe(bpy)3](p-CH3C6H4SO2)2 as an electroactive indicator . J. Pharm. Biomed. Anal. 2019 ; 163 : 180 – 187 . Google Scholar Crossref Search ADS PubMed WorldCat 53. Mahr R. , Frunzke J. Transcription factor-based biosensors in biotechnology: current state and future prospects . Appl. Microbiol. Biotechnol. 2016 ; 100 : 79 – 90 . Google Scholar Crossref Search ADS PubMed WorldCat 54. Engstrom M.D. , Pfleger B.F. Transcription control engineering and applications in synthetic biology . Synth. Syst. Biotechnol. 2017 ; 2 : 176 – 191 . Google Scholar Crossref Search ADS PubMed WorldCat 55. Pandi A. , Koch M., Voyvodic P.L., Soudier P., Bonnet J., Kushwaha M., Faulon J.L. Metabolic perceptrons for neural computing in biological systems . Nat. Commun. 2019 ; 10 : 3880 . Google Scholar Crossref Search ADS PubMed WorldCat 56. Voyvodic P.L. , Pandi A., Koch M., Conejero I., Valjent E., Courtet P., Renard E., Faulon J.L., Bonnet J. Plug-and-play metabolic transducers expand the chemical detection space of cell-free biosensors . Nat. Commun. 2019 ; 10 : 1697 . Google Scholar Crossref Search ADS PubMed WorldCat 57. Koch M. , Pandi A., Delépine B., Faulon J.L. A dataset of small molecules triggering transcriptional and translational cellular responses . Data Br. 2018 ; 17 : 1374 – 1378 . Google Scholar Crossref Search ADS WorldCat 58. Delépine B. , Libis V., Carbonell P., Faulon J.L. SensiPath: computer-aided design of sensing-enabling metabolic pathways . Nucleic Acids Res. 2016 ; 44 : W226 – W231 . Google Scholar Crossref Search ADS PubMed WorldCat 59. Snoek T. , Chaberski E.K., Ambri F., Kol S., Bjørn S.P., Pang B., Barajas J.F., Welner D.H., Jensen M.K., Keasling J.D. Evolution-guided engineering of small-molecule biosensors . Nucleic Acids Res. 2020 ; 48 : e3 . Google Scholar Crossref Search ADS PubMed WorldCat 60. Ellefson J.W. , Ledbetter M.P., Ellington A.D. Directed evolution of a synthetic phylogeny of programmable Trp repressors article . Nat. Chem. Biol. 2018 ; 14 : 361 – 367 . Google Scholar Crossref Search ADS PubMed WorldCat 61. Rutherford J.C. , Bird A.J. Metal-responsive transcription factors that regulate iron, zinc, and copper homeostasis in eukaryotic cells . Eukaryot. Cell . 2004 ; 3 :https://doi.org/10.1128/EC.3.1.1-13.2004. Google Scholar OpenURL Placeholder Text WorldCat 62. Rollins J.A. , Dickman M.B. pH signaling in Sclerotinia sclerotiorum: identification of a pacC/RIM1 homolog . Appl. Environ. Microbiol. 2001 ; 67 : 75 – 81 . Google Scholar Crossref Search ADS PubMed WorldCat 63. Janmey P.A. , Wells R.G., Assoian R.K., McCulloch C.A. From tissue mechanics to transcription factors . Differentiation . 2013 ; 86 : 112 – 120 . Google Scholar Crossref Search ADS PubMed WorldCat 64. Pu J. , Zinkus-Boltz J., Dickinson B.C. Evolution of a split RNA polymerase as a versatile biosensor platform . Nat. Chem. Biol. 2017 ; 13 : 432 – 438 . Google Scholar Crossref Search ADS PubMed WorldCat 65. Schmelas C. , Grimm D. Split Cas9, not hairs − advancing the therapeutic index of CRISPR technology . Biotechnol. J. 2018 ; 13 : e1700432 . Google Scholar Crossref Search ADS PubMed WorldCat 66. Repina N.A. , Rosenbloom A., Mukherjee A., Schaffer D. V., Kane R.S. At light speed: advances in optogenetic systems for regulating cell signaling and behavior . Annu. Rev. Chem. Biomol. Eng. 2017 ; 8 : 13 – 39 . Google Scholar Crossref Search ADS PubMed WorldCat 67. Genot A.J. , Baccouche A., Sieskind R., Aubert-Kato N., Bredeche N., Bartolo J.F., Taly V., Fujii T., Rondelez Y. High-resolution mapping of bifurcations in nonlinear biochemical circuits . Nat. Chem. 2016 ; 8 : 760 – 767 . Google Scholar Crossref Search ADS PubMed WorldCat 68. Green L.N. , Subramanian H.K.K., Mardanlou V., Kim J., Hariadi R.F., Franco E. Autonomous dynamic control of DNA nanostructure self-assembly . Nat. Chem. 2019 ; 11 : 510 – 520 . Google Scholar Crossref Search ADS PubMed WorldCat 69. Janssen B.M.G. , Engelen W., Merkx M. DNA-directed control of enzyme-inhibitor complex formation: a modular approach to reversibly switch enzyme activity . ACS Synth. Biol. 2015 ; 4 : 547 – 553 . Google Scholar Crossref Search ADS PubMed WorldCat 70. Meijer L.H.H. , Joesaar A., Steur E., Engelen W., Van Santen R.A., Merkx M., De Greef T.F.A. Hierarchical control of enzymatic actuators using DNA-based switchable memories . Nat. Commun. 2017 ; 8 : 1117 . Google Scholar Crossref Search ADS PubMed WorldCat 71. Dramé-Maigné A. , Espada R., McCallum G., Sieskind R., Rondelez Y. Directed evolution of enzymes based on in vitro programmable self-replication . 2021 ; bioRxiv doi: 22 April 2021, preprint: not peer reviewedhttps://doi.org/10.1101/2021.04.22.440993. 72. Urtel G. , Van Der Hofstadt M., Galas J.C., Estevez-Torres A. REXPAR: an isothermal amplification scheme that is robust to autocatalytic parasites . Biochemistry . 2019 ; 58 : 2675 – 2681 . Google Scholar Crossref Search ADS PubMed WorldCat © The Author(s) 2021. Published by Oxford University Press on behalf of Nucleic Acids Research. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited. TI - A small-molecule chemical interface for molecular programs JF - Nucleic Acids Research DO - 10.1093/nar/gkab470 DA - 2021-07-05 UR - https://www.deepdyve.com/lp/oxford-university-press/a-small-molecule-chemical-interface-for-molecular-programs-QtRec0hVnb SP - 1 EP - 1 VL - Advance Article IS - DP - DeepDyve ER -