Force-free activation of Notch with DNA origami.

The Notch signaling pathway is a highly conserved, fundamental process to embryogenesis and neurogenesis. While force-induced conformational change is known to activate Notch receptors, Smyrlaki et al. recently used DNA origami to reveal an additional, force-independent mode of Notch activation via soluble presentation of spatially controlled ligand nanopatterns.

Engineering Exciton Dynamics with Synthetic DNA Scaffolds.

Excitons are the molecular-scale currency of electronic energy. Control over excitons enables energy to be directed and harnessed for light harvesting, electronics, and sensing. Excitonic circuits achieve such control by arranging electronically active molecules to prescribe desired spatiotemporal dynamics. Photosynthetic solar energy conversion is a canonical example of the power of excitonic circuits, where chromophores are positioned in a protein scaffold to perform efficient light capture, energy transport, and charge separation. Synthetic systems that aim to emulate this functionality include self-assembled aggregates, molecular crystals, and chromophore-modified proteins. While the potential of this approach is clear, these systems lack the structural precision to control excitons or even test the limits of their power. In recent years, DNA origami has emerged as a designer material that exploits biological building blocks to construct nanoscale architectures. The structural precision afforded by DNA origami has enabled the pursuit of naturally inspired organizational principles in a highly precise and scalable manner. In this Account, we describe recent developments in DNA-based platforms that spatially organize chromophores to construct tunable excitonic systems. The high fidelity of DNA base pairing enables the formation of programmable nanoscale architectures, and sequence-specific placement allows for the precise positioning of chromophores within the DNA structure. The integration of a wide range of chromophores across the visible spectrum introduces spectral tunability. These excitonic DNA-chromophore assemblies not only serve as model systems for light harvesting, solar conversion, and sensing but also lay the groundwork for the integration of coupled chromophores into larger-scale nucleic acid architectures.We have used this approach to generate DNA-chromophore assemblies of strongly coupled delocalized excited states through both sequence-specific self-assembly and the covalent attachment of chromophores. These strategies have been leveraged to independently control excitonic coupling and system-bath interaction, which together control energy transfer. We then extended this framework to identify how scaffold configurations can steer the formation of symmetry-breaking charge transfer states, paving the way toward the design of dual light-harvesting and charge separation DNA machinery. In an orthogonal application, we used the programmability of DNA chromophore assemblies to change the optical emission properties of strongly coupled dimers, generating a series of fluorophore-modified constructs with separable emission properties for fluorescence assays. Upcoming advances in the chemical modification of nucleotides, design of large-scale DNA origami, and predictive computational methods will aid in constructing excitonic assemblies for optical and computing applications. Collectively, the development of DNA-chromophore assemblies as a platform for excitonic circuitry offers a pathway to identifying and applying design principles for light harvesting and molecular electronics.

Binding Site Programmable Self-Assembly of 3D Hierarchical DNA Origami Nanostructures.

DNA nanotechnology has broad applications in biomedical drug delivery and programmable materials. Characterization of the self-assembly of DNA origami and quantum dots (QDs) is necessary for the development of new DNA-based nanostructures. We use computation and experiment to show that the self-assembly of 3D hierarchical nanostructures can be controlled by programming the binding site number and their positions on DNA origami. Using biotinylated pentagonal pyramid wireframe DNA origamis and streptavidin capped QDs, we demonstrate that DNA origami with 1 binding site at the outer vertex can assemble multimeric origamis with up to 6 DNA origamis on 1 QD, and DNA origami with 1 binding site at the inner center can only assemble monomeric and dimeric origamis. Meanwhile, the yield percentages of different multimeric origamis are controlled by the QD:DNA-origami stoichiometric mixing ratio. DNA origamis with 2 binding sites at the αγ positions (of the pentagon) make larger nanostructures than those with binding sites at the αß positions. In general, increasing the number of binding sites leads to increases in the nanostructure size. At high DNA origami concentration, the QD number in each cluster becomes the limiting factor for the growth of nanostructures. We find that reducing the QD size can also affect the self-assembly because of the reduced access to the binding sites from more densely packed origamis.

Computational investigation of the impact of core sequence on immobile DNA four-way junction structure and dynamics.

Immobile four-way junctions (4WJs) are core structural motifs employed in the design of programmed DNA assemblies. Understanding the impact of sequence on their equilibrium structure and flexibility is important to informing the design of complex DNA architectures. While core junction sequence is known to impact the preferences for the two possible isomeric states that junctions reside in, previous investigations have not quantified these preferences based on molecular-level interactions. Here, we use all-atom molecular dynamics simulations to investigate base-pair level structure and dynamics of four-way junctions, using the canonical Seeman J1 junction as a reference. Comparison of J1 with equivalent single-crossover topologies and isolated nicked duplexes reveal conformational impact of the double-crossover motif. We additionally contrast J1 with a second junction core sequence termed J24, with equal thermodynamic preference for each isomeric configuration. Analyses of the base-pair degrees of freedom for each system, free energy calculations, and reduced-coordinate sampling of the 4WJ isomers reveal the significant impact base sequence has on local structure, isomer bias, and global junction dynamics.

Versatile Dehydration-Assisted Functionalization of Quantum Dots and Rods.

Functionalization of quantum dots (QDs) and quantum rods (QRs) with ligands is essential for their further practical application across various domains. Dehydration-assisted functionalization (DAF) is a versatile method applicable to a wide range of hydrophilic ligands with an affinity to the surface of QDs and QRs. This approach facilitates rapid one-pot ligand exchange and dense modification by efficiently transferring these ligands onto the surface of QDs and QRs. This study demonstrates the efficacy of DAF in preparing chiral QRs, engineering the surface charge of QDs, utilizing QR aggregates, and conjugating dense DNA onto cadmium-free InP/ZnS QDs. DAF therefore offers a versatile solution for hydrophilic ligand functionalization of QDs and QRs applicable to diverse applications.

Rapid prototyping of arbitrary 2D and 3D wireframe DNA origami.

Wireframe DNA origami assemblies can now be programmed automatically from the top-down using simple wireframe target geometries, or meshes, in 2D and 3D, using either rigid, six-helix bundle (6HB) or more compliant, two-helix bundle (DX) edges. While these assemblies have numerous applications in nanoscale materials fabrication due to their nanoscale spatial addressability and high degree of customization, no easy-to-use graphical user interface software yet exists to deploy these algorithmic approaches within a single, standalone interface. Further, top-down sequence design of 3D DX-based objects previously enabled by DAEDALUS was limited to discrete edge lengths and uniform vertex angles, limiting the scope of objects that can be designed. Here, we introduce the open-source software package ATHENA with a graphical user interface that automatically renders single-stranded DNA scaffold routing and staple strand sequences for any target wireframe DNA origami using DX or 6HB edges, including irregular, asymmetric DX-based polyhedra with variable edge lengths and vertices demonstrated experimentally, which significantly expands the set of possible 3D DNA-based assemblies that can be designed. ATHENA also enables external editing of sequences using caDNAno, demonstrated using asymmetric nanoscale positioning of gold nanoparticles, as well as providing atomic-level models for molecular dynamics, coarse-grained dynamics with oxDNA, and other computational chemistry simulation approaches.

Isometric Hamming embeddings of weighted graphs.

A mapping α:V(G)→V(H) from the vertex set of one graph G to another graph H is an isometric embedding if the shortest path distance between any two vertices in G equals the distance between their images in H. Here, we consider isometric embeddings of a weighted graph G into unweighted Hamming graphs, called Hamming embeddings, when G satisfies the property that every edge is a shortest path between its endpoints. Using a Cartesian product decomposition of G called its canonical isometric representation, we show that every Hamming embedding of G may be partitioned into a canonical partition, whose parts provide Hamming embeddings for each factor of the canonical isometric representation of G. This implies that G permits a Hamming embedding if and only if each factor of its canonical isometric representation is Hamming embeddable. This result extends prior work on unweighted graphs that showed that an unweighted graph permits a Hamming embedding if and only if each factor is a complete graph. When a graph G has nontrivial isometric representation, determining whether G has a Hamming embedding can be simplified to checking embeddability of two or more smaller graphs.

Factorization and pseudofactorization of weighted graphs.

For unweighted graphs, finding isometric embeddings of a graph G is closely related to decompositions of G into Cartesian products of smaller graphs. When G is isomorphic to a Cartesian graph product, we call the factors of this product a factorization of G. When G is isomorphic to an isometric subgraph of a Cartesian graph product, we call those factors a pseudofactorization of G. Prior work has shown that an unweighted graph's pseudofactorization can be used to generate a canonical isometric embedding into a product of the smallest possible pseudofactors. However, for arbitrary weighted graphs, which represent a richer variety of metric spaces, methods for finding isometric embeddings or determining their existence remain elusive, and indeed pseudofactorization and factorization have not previously been extended to this context. In this work, we address the problem of finding the factorization and pseudofactorization of a weighted graph G, where G satisfies the property that every edge constitutes a shortest path between its endpoints. We term such graphs minimal graphs, noting that every graph can be made minimal by removing edges not affecting its path metric. We generalize pseudofactorization and factorization to minimal graphs and develop new proof techniques that extend the previously proposed algorithms due to Graham and Winkler [Graham and Winkler, '85] and Feder [Feder, '92] for pseudofactorization and factorization of unweighted graphs. We show that any n-vertex, m-edge graph with positive integer edge weights can be factored in O(m2) time, plus the time to find all pairs shortest paths (APSP) distances in a weighted graph, resulting in an overall running time of O(m2+n2 log log n) time. We also show that a pseudofactorization for such a graph can be computed in O(mn) time, plus the time to solve APSP, resulting in an O(mn+n2 log log n) running time.

Programming 2D Supramolecular Assemblies with Wireframe DNA Origami.

Wireframe DNA origami offers the ability to program nearly arbitrary 2D and 3D nanoscale geometries, with six-helix bundle (6HB) edge designs providing both geometric versatility and fidelity with respect to the target origami shape. Because individual DNA origami objects are limited in size by the length of the DNA scaffold, here, we introduce a hierarchical self-assembly strategy to overcome this limitation by programming supramolecular assemblies and periodic arrays using wireframe DNA origami objects as building blocks. Parallel half-crossovers are used together with lateral cohesive interactions between staples and the scaffold to introduce symmetry into supramolecular assemblies constructed from single DNA origami units that cannot be self-assembled directly using base-stacking or conventional antiparallel crossover designs. This hierarchical design approach can be applied readily to 2D wireframe DNA origami designed using the top-down sequence design strategy METIS without any prerequisites on scaffold and staple routing. We demonstrate the utility of our strategy by fabricating dimers and self-limiting hexameric superstructures using both triangular and hexagonal wireframe origami building blocks. We generalize our self-assembly approach to fabricate close-packed and non-close-packed periodic 2D arrays. Visualization using atomic force microscopy and transmission electron microscopy demonstrates that superstructures exhibit similar structural integrity to that of the individual origami building blocks designed using METIS. Our results offer a general platform for the design and fabrication of 2D materials for a variety of applications.

Random access DNA memory using Boolean search in an archival file storage system.

DNA is an ultrahigh-density storage medium that could meet exponentially growing worldwide demand for archival data storage if DNA synthesis costs declined sufficiently and if random access of files within exabyte-to-yottabyte-scale DNA data pools were feasible. Here, we demonstrate a path to overcome the second barrier by encapsulating data-encoding DNA file sequences within impervious silica capsules that are surface labelled with single-stranded DNA barcodes. Barcodes are chosen to represent file metadata, enabling selection of sets of files with Boolean logic directly, without use of amplification. We demonstrate random access of image files from a prototypical 2-kilobyte image database using fluorescence sorting with selection sensitivity of one in 106 files, which thereby enables one in 106N selection capability using N optical channels. Our strategy thereby offers a scalable concept for random access of archival files in large-scale molecular datasets.

A Bayesian framework for high-throughput T cell receptor pairing.

MOTIVATION: The study of T cell receptor (TCR) repertoires has generated new insights into immune system recognition. However, the ability to robustly characterize these populations has been limited by technical barriers and an inability to reliably infer heterodimeric chain pairings for TCRs. RESULTS: Here, we describe a novel analytical approach to an emerging immune repertoire sequencing method, improving the resolving power of this low-cost technology. This method relies upon the distribution of a T cell population across a 96-well plate, followed by barcoding and sequencing of the relevant transcripts from each T cell. Multicell Analytical Deconvolution for High Yield Paired-chain Evaluation (MAD-HYPE) uses Bayesian inference to more accurately extract TCR information, improving our ability to study and characterize T cell populations for immunology and immunotherapy applications. AVAILABILITY AND IMPLEMENTATION: The MAD-HYPE algorithm is released as an open-source project under the Apache License and is available from https://github.com/birnbaumlab/MAD-HYPE. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

DoGNet: A deep architecture for synapse detection in multiplexed fluorescence images.

Neuronal synapses transmit electrochemical signals between cells through the coordinated action of presynaptic vesicles, ion channels, scaffolding and adapter proteins, and membrane receptors. In situ structural characterization of numerous synaptic proteins simultaneously through multiplexed imaging facilitates a bottom-up approach to synapse classification and phenotypic description. Objective automation of efficient and reliable synapse detection within these datasets is essential for the high-throughput investigation of synaptic features. Convolutional neural networks can solve this generalized problem of synapse detection, however, these architectures require large numbers of training examples to optimize their thousands of parameters. We propose DoGNet, a neural network architecture that closes the gap between classical computer vision blob detectors, such as Difference of Gaussians (DoG) filters, and modern convolutional networks. DoGNet is optimized to analyze highly multiplexed microscopy data. Its small number of training parameters allows DoGNet to be trained with few examples, which facilitates its application to new datasets without overfitting. We evaluate the method on multiplexed fluorescence imaging data from both primary mouse neuronal cultures and mouse cortex tissue slices. We show that DoGNet outperforms convolutional networks with a low-to-moderate number of training examples, and DoGNet is efficiently transferred between datasets collected from separate research groups. DoGNet synapse localizations can then be used to guide the segmentation of individual synaptic protein locations and spatial extents, revealing their spatial organization and relative abundances within individual synapses. The source code is publicly available: https://github.com/kulikovv/dognet.

Programmed coherent coupling in a synthetic DNA-based excitonic circuit.

Natural light-harvesting systems spatially organize densely packed chromophore aggregates using rigid protein scaffolds to achieve highly efficient, directed energy transfer. Here, we report a synthetic strategy using rigid DNA scaffolds to similarly program the spatial organization of densely packed, discrete clusters of cyanine dye aggregates with tunable absorption spectra and strongly coupled exciton dynamics present in natural light-harvesting systems. We first characterize the range of dye-aggregate sizes that can be templated spatially by A-tracts of B-form DNA while retaining coherent energy transfer. We then use structure-based modelling and quantum dynamics to guide the rational design of higher-order synthetic circuits consisting of multiple discrete dye aggregates within a DX-tile. These programmed circuits exhibit excitonic transport properties with prominent circular dichroism, superradiance, and fast delocalized exciton transfer, consistent with our quantum dynamics predictions. This bottom-up strategy offers a versatile approach to the rational design of strongly coupled excitonic circuits using spatially organized dye aggregates for use in coherent nanoscale energy transport, artificial light-harvesting, and nanophotonics.

Structure and conformational dynamics of scaffolded DNA origami nanoparticles.

Synthetic DNA is a highly programmable nanoscale material that can be designed to self-assemble into 3D structures that are fully determined by underlying Watson-Crick base pairing. The double crossover (DX) design motif has demonstrated versatility in synthesizing arbitrary DNA nanoparticles on the 5-100 nm scale for diverse applications in biotechnology. Prior computational investigations of these assemblies include all-atom and coarse-grained modeling, but modeling their conformational dynamics remains challenging due to their long relaxation times and associated computational cost. We apply all-atom molecular dynamics and coarse-grained finite element modeling to DX-based nanoparticles to elucidate their fine-scale and global conformational structure and dynamics. We use our coarse-grained model with a set of secondary structural motifs to predict the equilibrium solution structures of 45 DX-based DNA origami nanoparticles including a tetrahedron, octahedron, icosahedron, cuboctahedron and reinforced cube. Coarse-grained models are compared with 3D cryo-electron microscopy density maps for these five DNA nanoparticles and with all-atom molecular dynamics simulations for the tetrahedron and octahedron. Our results elucidate non-intuitive atomic-level structural details of DX-based DNA nanoparticles, and offer a general framework for efficient computational prediction of global and local structural and mechanical properties of DX-based assemblies that are inaccessible to all-atom based models alone.

Inferring transient particle transport dynamics in live cells.

Live-cell imaging and particle tracking provide rich information on mechanisms of intracellular transport. However, trajectory analysis procedures to infer complex transport dynamics involving stochastic switching between active transport and diffusive motion are lacking. We applied Bayesian model selection to hidden Markov modeling to infer transient transport states from trajectories of mRNA-protein complexes in live mouse hippocampal neurons and metaphase kinetochores in dividing human cells. The software is available at http://hmm-bayes.org/.

Molecular model of J-aggregated pseudoisocyanine fibers.

Aggregated cyanines form ordered supramolecular structures with the potential to transport energy efficiently over long distances, a hallmark of photosynthetic light-harvesting complexes. In concentrated aqueous solution, pseudoisocyanine (PIC) spontaneously forms fibers with a chiral J-band red-shifted 1600 cm-1 from the monomeric 0-0 transition. A cryogenic transmission electron microscopy analysis of these fibers show an average fiber width of 2.89 nm, although the molecular-level structure of the aggregate is currently unknown. To determine a molecular model for these PIC fibers, the calculated spectra and dynamics using a Frenkel exciton model are compared to experiment. A chiral aggregate model in which the PIC monomers are neither parallel nor orthogonal to the long axis of the fiber is shown to replicate the experimental spectra most closely. This model can be physically realized by the sequential binding of PIC dimers and monomers to the ends of the fiber. These insights into the molecular aggregation model for aqueous PIC can also be applied to other similar cyanine-based supramolecular complexes with the potential for long-range energy transport, a key building block for the rational design of novel excitonic systems.

Optical determination of the electronic coupling and intercalation geometry of thiazole orange homodimer in DNA.

Sequence-selective bis-intercalating dyes exhibit large increases in fluorescence in the presence of specific DNA sequences. This property makes this class of fluorophore of particular importance to biosensing and super-resolution imaging. Here we report ultrafast transient anisotropy measurements of resonance energy transfer (RET) between thiazole orange (TO) molecules in a complex formed between the homodimer TOTO and double-stranded (ds) DNA. Biexponential homo-RET dynamics suggest two subpopulations within the ensemble: 80% intercalated and 20% non-intercalated. Based on the application of the transition density cube method to describe the electronic coupling and Monte Carlo simulations of the TOTO/dsDNA geometry, the dihedral angle between intercalated TO molecules is estimated to be 81° ± 5°, corresponding to a coupling strength of 45 ± 22 cm-1. Dye intercalation with this geometry is found to occur independently of the underlying DNA sequence, despite the known preference of TOTO for the nucleobase sequence CTAG. The non-intercalated subpopulation is inferred to have a mean inter-dye separation distance of 19 Å, corresponding to coupling strengths between 0 and 25 cm-1. This information is important to enable the rational design of energy transfer systems that utilize TOTO as a relay dye. The approach used here is generally applicable to determining the electronic coupling strength and intercalation configuration of other dimeric bis-intercalators.

Programming Self-Assembly of DNA Origami Honeycomb Two-Dimensional Lattices and Plasmonic Metamaterials.

Scaffolded DNA origami has proven to be a versatile method for generating functional nanostructures with prescribed sub-100 nm shapes. Programming DNA-origami tiles to form large-scale 2D lattices that span hundreds of nanometers to the micrometer scale could provide an enabling platform for diverse applications ranging from metamaterials to surface-based biophysical assays. Toward this end, here we design a family of hexagonal DNA-origami tiles using computer-aided design and demonstrate successful self-assembly of micrometer-scale 2D honeycomb lattices and tubes by controlling their geometric and mechanical properties including their interconnecting strands. Our results offer insight into programmed self-assembly of low-defect supra-molecular DNA-origami 2D lattices and tubes. In addition, we demonstrate that these DNA-origami hexagon tiles and honeycomb lattices are versatile platforms for assembling optical metamaterials via programmable spatial arrangement of gold nanoparticles (AuNPs) into cluster and superlattice geometries.

Structure-based model for light-harvesting properties of nucleic acid nanostructures.

Programmed self-assembly of DNA enables the rational design of megadalton-scale macromolecular assemblies with sub-nanometer scale precision. These assemblies can be programmed to serve as structural scaffolds for secondary chromophore molecules with light-harvesting properties. Like in natural systems, the local and global spatial organization of these synthetic scaffolded chromophore systems plays a crucial role in their emergent excitonic and optical properties. Previously, we introduced a computational model to predict the large-scale 3D solution structure and flexibility of nucleic acid nanostructures programmed using the principle of scaffolded DNA origami. Here, we use Förster resonance energy transfer theory to simulate the temporal dynamics of dye excitation and energy transfer accounting both for overall DNA nanostructure architecture as well as atomic-level DNA and dye chemical structure and composition. Results are used to calculate emergent optical properties including effective absorption cross-section, absorption and emission spectra and total power transferred to a biomimetic reaction center in an existing seven-helix double stranded DNA-based antenna. This structure-based computational framework enables the efficient in silico evaluation of nucleic acid nanostructures for diverse light-harvesting and photonic applications.

H2A.Z acidic patch couples chromatin dynamics to regulation of gene expression programs during ESC differentiation.

The histone H2A variant H2A.Z is essential for embryonic development and for proper control of developmental gene expression programs in embryonic stem cells (ESCs). Divergent regions of amino acid sequence of H2A.Z likely determine its functional specialization compared to core histone H2A. For example, H2A.Z contains three divergent residues in the essential C-terminal acidic patch that reside on the surface of the histone octamer as an uninterrupted acidic patch domain; however, we know little about how these residues contribute to chromatin structure and function. Here, we show that the divergent amino acids Gly92, Asp97, and Ser98 in the H2A.Z C-terminal acidic patch (H2A.Z(AP3)) are critical for lineage commitment during ESC differentiation. H2A.Z is enriched at most H3K4me3 promoters in ESCs including poised, bivalent promoters that harbor both activating and repressive marks, H3K4me3 and H3K27me3 respectively. We found that while H2A.Z(AP3) interacted with its deposition complex and displayed a highly similar distribution pattern compared to wild-type H2A.Z, its enrichment levels were reduced at target promoters. Further analysis revealed that H2A.Z(AP3) was less tightly associated with chromatin, suggesting that the mutant is more dynamic. Notably, bivalent genes in H2A.Z(AP3) ESCs displayed significant changes in expression compared to active genes. Moreover, bivalent genes in H2A.Z(AP3) ESCs gained H3.3, a variant associated with higher nucleosome turnover, compared to wild-type H2A.Z. We next performed single cell imaging to measure H2A.Z dynamics. We found that H2A.Z(AP3) displayed higher mobility in chromatin compared to wild-type H2A.Z by fluorescent recovery after photobleaching (FRAP). Moreover, ESCs treated with the transcriptional inhibitor flavopiridol resulted in a decrease in the H2A.Z(AP3) mobile fraction and an increase in its occupancy at target genes indicating that the mutant can be properly incorporated into chromatin. Collectively, our work suggests that the divergent residues in the H2A.Z acidic patch comprise a unique domain that couples control of chromatin dynamics to the regulation of developmental gene expression patterns during lineage commitment.