Abstract
The past decades have witnessed the evolution of electronic and photonic integrated circuits, from application specific to programmable1,2. Although liquid-phase DNA circuitry holds the potential for massive parallelism in the encoding and execution of algorithms3,4, the development of general-purpose DNA integrated circuits (DICs) has yet to be explored. Here we demonstrate a DIC system by integration of multilayer DNA-based programmable gate arrays (DPGAs). We find that the use of generic single-stranded oligonucleotides as a uniform transmission signal can reliably integrate large-scale DICs with minimal leakage and high fidelity for general-purpose computing. Reconfiguration of a single DPGA with 24 addressable dual-rail gates can be programmed with wiring instructions to implement over 100 billion distinct circuits. Furthermore, to control the intrinsically random collision of molecules, we designed DNA origami registers to provide the directionality for asynchronous execution of cascaded DPGAs. We exemplify this by a quadratic equation-solving DIC assembled with three layers of cascade DPGAs comprising 30 logic gates with around 500 DNA strands. We further show that integration of a DPGA with an analog-to-digital converter can classify disease-related microRNAs. The ability to integrate large-scale DPGA networks without apparent signal attenuation marks a key step towards general-purpose DNA computing.
This is a preview of subscription content, access via your institution
Access options
Access Nature and 54 other Nature Portfolio journals
Get Nature+, our best-value online-access subscription
24,99 € / 30 days
cancel any time
Subscribe to this journal
Receive 51 print issues and online access
185,98 € per year
only 3,65 € per issue
Buy this article
- Purchase on Springer Link
- Instant access to full article PDF
Prices may be subject to local taxes which are calculated during checkout
Similar content being viewed by others
Data availability
The data that support the findings of this study are available in the manuscript or the Supplementary Information. Source data are provided with this paper. All other data is available on request.
Code availability
Source codes used in this study (Visual DSD, MATLAB, Python) are available from GitHub (https://github.com/FeiWANG-SJTU/DPGA). All other codes are available from the corresponding authors on reasonable request.
References
Burks, A. W. in A History of Computing in the Twentieth Century (ed. Metropolis, N.) 311–344 (Elsevier, 1980).
Bogaerts, W. et al. Programmable photonic circuits. Nature 586, 207–216 (2020).
Chen, X. & Ellington, A. D. Shaping up nucleic acid computation. Curr. Opin. Biotechnol. 21, 392–400 (2010).
Li, J., Green, A. A., Yan, H. & Fan, C. Engineering nucleic acid structures for programmable molecular circuitry and intracellular biocomputation. Nat. Chem. 9, 1056–1067 (2017).
Benenson, Y. et al. Programmable and autonomous computing machine made of biomolecules. Nature 414, 430–434 (2001).
Qian, L. & Winfree, E. Scaling up digital circuit computation with DNA strand displacement cascades. Science 332, 1196–1201 (2011).
Seelig, G., Soloveichik, D., Zhang, D. Y. & Winfree, E. Enzyme-free nucleic acid logic circuits. Science 314, 1585–1588 (2006).
Wang, F. et al. Implementing digital computing with DNA-based switching circuits. Nat. Commun. 11, 121 (2020).
Thubagere, A. J. et al. A cargo-sorting DNA robot. Science 357, eaan6558 (2017).
Chao, J. et al. Solving mazes with single-molecule DNA navigators. Nat. Mater. 18, 273–279 (2019).
Cherry, K. M. & Qian, L. L. Scaling up molecular pattern recognition with DNA-based winner-take-all neural networks. Nature 559, 370–376 (2018).
Kishi, J. Y., Schaus, T. E., Gopalkrishnan, N., Xuan, F. & Yin, P. Programmable autonomous synthesis of single-stranded DNA. Nat. Chem. 10, 155–164 (2018).
Zhang, Y. et al. Selective transformations between nanoparticle superlattices via the reprogramming of DNA-mediated interactions. Nat. Mater. 14, 840–847 (2015).
Douglas, S. M., Bachelet, I. & Church, G. M. A logic-gated nanorobot for targeted transport of molecular payloads. Science 335, 831–834 (2012).
Lopez, R., Wang, R. & Seelig, G. A molecular multi-gene classifier for disease diagnostics. Nat. Chem. 10, 746–754 (2018).
Zhang, C. et al. Cancer diagnosis with DNA molecular computation. Nat. Nanotechnol. 15, 709–715 (2020).
Hills, G. et al. Modern microprocessor built from complementary carbon nanotube transistors. Nature 572, 595–602 (2019).
Debnath, S. et al. Demonstration of a small programmable quantum computer with atomic qubits. Nature 536, 63–66 (2016).
Athanas, P. M. & Silverman, H. F. Processor reconfiguration through instruction-set metamorphosis. Computer 26, 11–18 (1993).
Ruiz-Rosero, J., Ramirez-Gonzalez, G. & Khanna, R. Field programmable gate array applications—a scientometric review. Computation 7, 63 (2019).
Benenson, Y. Biomolecular computing systems: principles, progress and potential. Nat. Rev. Genet. 13, 455–468 (2012).
Pei, R., Matamoros, E., Liu, M., Stefanovic, D. & Stojanovic, M. N. Training a molecular automaton to play a game. Nat. Nanotechnol. 5, 773–777 (2010).
Woods, D. et al. Diverse and robust molecular algorithms using reprogrammable DNA self-assembly. Nature 567, 366–372 (2019).
Rabaey, J. M. Digital Integrated Circuits: A Design Perspective (Prentice Hall, 1995).
Wan, N. H. et al. Large-scale integration of artificial atoms in hybrid photonic circuits. Nature 583, 226–231 (2020).
Klosin, A. et al. Phase separation provides a mechanism to reduce noise in cells. Science 367, 464–468 (2020).
Chatterjee, G., Dalchau, N., Muscat, R. A., Phillips, A. & Seelig, G. A spatially localized architecture for fast and modular DNA computing. Nat. Nanotechnol. 12, 920–927 (2017).
Bian, X., Zhang, Z., Xiong, Q., De Camilli, P. & Lin, C. A programmable DNA-origami platform for studying lipid transfer between bilayers. Nat. Chem. Biol. 15, 830–837 (2019).
Joesaar, A. et al. DNA-based communication in populations of synthetic protocells. Nat. Nanotechnol. 14, 369–378 (2019).
Zhang, D. Y. & Winfree, E. Control of DNA strand displacement kinetics using toehold exchange. J. Am. Chem. Soc. 131, 17303–17314 (2009).
Li, W., Zhang, F., Yan, H. & Liu, Y. DNA based arithmetic function: a half adder based on DNA strand displacement. Nanoscale 8, 3775–3784 (2016).
Song, T. et al. Fast and compact DNA logic circuits based on single-stranded gates using strand-displacing polymerase. Nat. Nanotechnol. 14, 1075–1081 (2019).
Lund, K. et al. Molecular robots guided by prescriptive landscapes. Nature 465, 206–210 (2010).
Ke, Y. G., Lindsay, S., Chang, Y., Liu, Y. & Yan, H. Self-assembled water-soluble nucleic acid probe tiles for label-free RNA hybridization assays. Science 319, 180–183 (2008).
Praetorius, F. et al. Biotechnological mass production of DNA origami. Nature 552, 84–87 (2017).
Lauback, S. et al. Real-time magnetic actuation of DNA nanodevices via modular integration with stiff micro-levers. Nat. Commun. 9, 1446 (2018).
Genot, A. J. et al. High-resolution mapping of bifurcations in nonlinear biochemical circuits. Nat. Chem. 8, 760��767 (2016).
Rothemund, P. W. K. Folding DNA to create nanoscale shapes and patterns. Nature 440, 297–302 (2006).
Wagenbauer, K. F. et al. How we make DNA origami. ChemBioChem 18, 1873–1885 (2017).
Ouldridge, T. E., Louis, A. A. & Doye, J. P. K. Structural, mechanical, and thermodynamic properties of a coarse-grained DNA model. J. Chem. Phys. 134, 085101 (2011).
Snodin, B. E. K. et al. Introducing improved structural properties and salt dependence into a coarse-grained model of DNA. J. Chem. Phys. 142, 234901 (2015).
Suma, A. et al. TacoxDNA: a user-friendly web server for simulations of complex DNA structures, from single strands to origami. J. Comput. Chem. 40, 2586–2595 (2019).
Sulc, P. Coarse-grained Modelling of Nucleic Acids (Univ. of Oxford, 2014).
Doye, J. P. et al. The oxDNA coarse-grained model as a tool to simulate DNA origami. Methods Mol. Biol. 2639, 93–112 (2023).
Poppleton, E. et al. Design, optimization and analysis of large DNA and RNA nanostructures through interactive visualization, editing and molecular simulation. Nucleic Acids Res. 48, e72 (2020).
Acknowledgements
We thank L. Qian (Caltech) and Y. Huang (Peking University) for helpful discussions. This work was supported by the National Key R&D Program of China (grant no. 2021YFF1200300), the National Natural Science Foundation of China (grant nos. T2188102, 21991134, 21904060, 22025404 and 22104088), the Science Foundation of Shanghai Municipal Science and Technology Commission (grant nos. 20dz1101000 and 21TQ1400222) and the New Cornerstone Investigator Program. Molecular dynamics simulations were run on the π2.0 cluster supported by the Center for High Performance Computing at Shanghai Jiao Tong University.
Author information
Authors and Affiliations
Contributions
C.F. and F.W. conceived the research. F.W. designed the circuits and wrote the sequence-generating and -compiling programmes. H.L. performed the majority of the experiments. N.X. performed AFM experiments. C.S. and Q.Z. performed magnetic field-related experiments. M.L. and F.W. performed simulations. L.Z. trained the nonlinear classification model. F.W., C.F. and H.L. analysed data and wrote the manuscript. M.D., J.L. and H.C. participated in data analysis and discussions. C.F., H.C., J.L., N.X., M.D. and X.Z. reviewed and edited the manuscript.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature thanks the anonymous reviewers for their contribution to the peer review of this work. Peer reviewer reports are available.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data figures and tables
Extended Data Fig. 1 Outline of the design strategy for general-purpose DICs.
a, Architecture for electronic chips integration. b, Hierarchical illustration of scalable DPGA integration (shown by the logical arrangement). We referred architectural properties of programmable electronic integrated circuits to design general-purpose DICs. In electronic integrated circuits, general-purpose chips can be physically integrated, with the information exchange between chips and storage realized via electrons. Analogous to electronic signal, DNA-UTS is used to transmit information. Inter-gate and inter-DPGA information transmission are all enabled by DNA-UTS. With uniform transmitted signals, integrability is permitted both at the gate and the DPGA levels. Asynchronous execution of cascaded DPGAs interdicts molecular diffusion between DPGAs, allowing DPGA integration. Hence, the scalability is enabled with the use of DNA-UTS and DNA origami register. In addition, arbitrary gate connection is allowed in a DPGA, providing rich programming space. In all, the programmability and the scalability support general-purpose computing with DICs.
Extended Data Fig. 2 Representative DNA circuits experimentally implemented via multi-level programming of DPGAs.
The DIC is composed of three levels: DNA gates; DPGAs and DPGA networks.
Extended Data Fig. 3 Demonstrated advantages of dual-rail gates.
a, A single-rail gate receives one molecule as an input and generates one molecule as output. The concentration of output increases rapidly when the output is 1 while slowly when the output is 0. b, A dual-rail gate receives a molecule representing 1 or another molecule representing 1 as an input, and generates one molecule as output = 0 or another molecule as output = 1. The output signal is represented by the difference of two output signals. When output = 1, the output signal increases; when output = 0, the output signal decreases. c, Implemention of a dual-rail XOR gate with AND-OR gates requires six gates. d, Left, heatmap showing the result values with all possible combinations of high and low signals, when a result is supposed to be 0. Triangular region above the upper red line represents the obtained dual-rail results smaller than 0.4. Right, heatmap of error flag with all possible combinations of high and low signals. Regions inside the green box have error flag values lower than 0.4. e, Five possible computing states for a result supposed to be 0. f, Left, heatmap showing the result values with all possible combinations of high and low signals, when a result is supposed to be 1. Triangular region below the lower red line represents the obtained dual-rail result larger than 0.6. Right, heatmap of error flag with all possible combinations of high and low signals. Regions inside green boxed have error flag values lower than 0.4. g, Five possible computing states for a result supposed to be 1.
Extended Data Fig. 4 Internal structures and operating mechanisms of the four dual-rail gates.
a, Schematics showing internal structures of the dual-rail gates. AND gate has two series switches to respond to in1H and in2H respectively, generating outH in the presence of both in1H and in2H. Another switch responds to both L inputs, generating ouL in the presence of either in1L or in2L. OR gate has an opposite internal structure of AND gate. NOT gate has two switches, one for L input and one for H input. XOR contains four input-controlled switches. b–e, Signal transmission paths with all possible input combinations for AND (b), OR (c), NOT (d) and XOR (e) gates. High signals are shown in red and low signals in green. Single-stranded inputs bind to logic gates and release the corresponding output strands through SDRs. For AND gate, both in1L (representing in1 = 0) and in2L hybridize with s1 to displace s2 (=outL). in2H hybridizes with S3 to displace S5, exposing toehold γ for in1H to replace S4 (=outH).
Extended Data Fig. 5 The complete simulated time trajectories (a) for Fig. 2e and distance distributions during the simulation time (b).
With 1-nt gap, the toehold domain bound reversibly to S3 (middles panel of Distance 1 trajectory, inset), which did not lead to further branch migration (high fluctuation for Distance 2 and 3). Therefore, the presence of only in1H cannot generate a fault result, permitting neglectable leakage similar to that of 0-nt gap. However, 2-nt gap allowed stable binding of in1H (left panel of Distance 1 trajectory, indicated by H-bond arrow), and the output strand was replaced via branch migration (left panels of Distance 2 and 3 trajectories, indicated by H-bond arrows). Insets in (b): Distributions of Distance 1 ranging from 0.5 nm to 2.5 nm.
Extended Data Fig. 6 Experimental optimization and performance evaluation of WIR2 and WIR3.
a, Molecular reactions for wiring (for WIR2) or reading out (for WIR3) an output signal from a gate. Low (green) and high (red) signals are transmitted independently. The dashed lines indicate the upstream binding region. b, To introduce the threshold-over-amplifier binding priority, the duplex region was shortened by 2 bp, which minimized non-specific input-threshold binding and the leakage. c, Without Threshold, output from a gate could be amplified to close to 1. However, weak signal leakage could also be amplified, leading to false result. Thus, the threshold is essential to suppress leakage before amplification. d, We used a Threshold molecule (Th) that could interact with output quickly. With 0.4× leakage signal, we found Th with a concentration higher than 0.4× can effectively suppress leakage. With 1× output signal, we found signal transmission speed decreased with Th concentration. To balance leakage suppression and computing speed, we used 0.4× to 0.6× Threshold for experiments. e, Signal wiring for High signal (left) and Low signal (right) using 0.4× threshold.
Extended Data Fig. 7 A representative compiling process of the DNA Compiler.
Statements containing different type of operations and different priorities can be compiled into wiring instructions to configure DPGA.
Extended Data Fig. 8 The signal relay function of WIR4s.
a, Possible transient binding reactions between unmatched molecules that may affect the circuit performance. b, Numerical simulation of the computing process from upstream inputs to downstream inputs considering the transient binding. With the increase of the circuit size, the exposed toeholds have a higher chance of being occupied by unmatched strands via transient binding, leading to reduced computing speed. c, Experimental computing kinetics of an OR gate when 0, 500, 1000 and 1500 nM unmatched threshold molecules were added, respectively. d, Dual-rail results of the OR gate followed by a WIR3 at different unmatched threshold levels. e, The 11-layer cascade circuit in Fig. 3g when implemented with a single DPGA. f, Simulated computing kinetics of 1- to 11-layer subcircuits implemented with a single DPGA without WIR4. g, Simulated dual-rail results of the cascade circuit showing decayed performance with the increase of the circuit depth, which is generally consistent with the experimental results. h, Circuit diagram of a cascade circuit containing 11 layers of dual-rail gates, which was divided and implemented by three configured DPGAs. i, Simulated computing kinetics of the subcircuits in the cascaded DPGAs connected by WIR4s.
Extended Data Fig. 9 Signal decay during the cascade of DPGAs.
a–b, Illustration of the circuit depth. The largest circuit depth (a) and the corresponding SDR pathway that contains 30 steps of reactions (b), beyond what can be achieved with a single reaction. c, Schematic illustrations showing tested systems with direct ideal input and that transmitted by WIR4s. d, Paired comparison showing the increase of the half-competition time (t1/2). d–g, Output, leakage and error showing insignificant difference between DPGAs using direct input and transmitted input.
Extended Data Fig .10 Computing kinetics of 20 experimentally tested samples in addition to that in Fig. 6g in test set.
The values in bracket show the normalized concentration of mir-200c, mir-204 and mir-887, respectively.
Supplementary information
Supplementary Information
Methods, texts, Figs. 1–62, Tables 1–7 and references that extend and support the data and discussion presented in the main text.
Supplementary Data 1
DNA sequences used in each figure.
Supplementary Data 2
DNA sequences for dual-rails gates and wiring instructions.
Supplementary Data 3
Information and experimental results of 103 circuits for robustness evaluation.
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Lv, H., Xie, N., Li, M. et al. DNA-based programmable gate arrays for general-purpose DNA computing. Nature 622, 292–300 (2023). https://doi.org/10.1038/s41586-023-06484-9
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41586-023-06484-9
This article is cited by
-
DNA as a universal chemical substrate for computing and data storage
Nature Reviews Chemistry (2024)
-
Chemical unclonable functions based on operable random DNA pools
Nature Communications (2024)
-
Chemical reservoir computation in a self-organizing reaction network
Nature (2024)
-
A spatially localized DNA linear classifier for cancer diagnosis
Nature Communications (2024)
-
A signal transmission strategy driven by gap-regulated exonuclease hydrolysis for hierarchical molecular networks
Communications Biology (2024)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.