Quantitative prediction of 3D solution shape and flexibility of nucleic acid nanostructures

DNA nanotechnology enables the programmed synthesis of intricate nanometer-scale structures for diverse applications in materials and biological science. Precise control over the 3D solution shape and mechanical flexibility of target designs is important to achieve desired functionality. Because experimental validation of designed nanostructures is time-consuming and cost-intensive, predictive physical models of nanostructure shape and flexibility have the capacity to enhance dramatically the design process. Here, we significantly extend and experimentally validate a computational modeling framework for DNA origami previously presented as CanDo [Castro,C.E., Kilchherr,F., Kim,D.-N., Shiao,E.L., Wauer,T., Wortmann,P., Bathe,M., Dietz,H. (2011) A primer to scaffolded DNA origami. Nat. Meth., 8, 221-229.]. 3D solution shape and flexibility are predicted from basepair connectivity maps now accounting for nicks in the DNA double helix, entropic elasticity of single-stranded DNA, and distant crossovers required to model wireframe structures, in addition to previous modeling (Castro,C.E., et al.) that accounted only for the canonical twist, bend and stretch stiffness of double-helical DNA domains. Systematic experimental validation of nanostructure flexibility mediated by internal crossover density probed using a 32-helix DNA bundle demonstrates for the first time that our model not only predicts the 3D solution shape of complex DNA nanostructures but also their mechanical flexibility. Thus, our model represents an important advance in the quantitative understanding of DNA-based nanostructure shape and flexibility, and we anticipate that this model will increase significantly the number and variety of synthetic nanostructures designed using nucleic acids.     

Probing tethered targets of a single biomolecular complex with atomic force microscopy

DNA origami shows tremendous promise as templates for the assembly of nano‐components and detection of molecular recognition events. So far, the method of choice for evaluating these structures has been atomic force microscopy (AFM), a powerful tool for imaging nanoscale objects. In most cases, tethered targets on DNA origami have proven to be highly effective samples for investigation. Still, while maximal assembly of the nanostructures might benefit from the greatest flexibility in the tether, AFM imaging requires a sufficient stability of the adsorbed components. The balance between the tether flexibility and sample stability is a major, poorly understood, concern in such studies. Here, we investigated the dependence of the tethering length on molecular capture events monitored by AFM. In our experiments, single biotin molecules were attached to DNA origami templates with various linker lengths of thymidine nucleotides, and their interaction with streptavidin was observed with AFM. Our results show that the streptavidin‐biotin complexes are easily detected with short tethered lengths, and that their morphological features clearly change with the tethering length. We identify the functionally useful tether lengths for these investigations, which are also expected to prove useful in the construction and further application of DNA origami in bio‐nanotechnology studies. Copyright © 2013 John Wiley & Sons, Ltd.     

Formation of non-base-pairing DNA microgels using directed phase transition of amphiphilic monomers

Programmability of DNA sequences enables the formation of synthetic DNA nanostructures and their macromolecular assemblies such as DNA hydrogels. The base pair-level interaction of DNA is a foundational and powerful mechanism to build DNA structures at the nanoscale; however, its temperature sensitivity and weak interaction force remain a barrier for the facile and scalable assembly of DNA structures toward higher-order structures. We conducted this study to provide an alternative, non-base-pairing approach to connect nanoscale DNA units to yield micrometer-sized gels based on the sequential phase transition of amphiphilic unit structures. Strong electrostatic interactions between DNA nanostructures and polyelectrolyte spermines led to the formation of giant phase-separated aggregates of monomer units. Gelation could be initiated by the addition of NaCl, which weakened the electrostatic DNA-spermine interaction while attractive interactions between cholesterols created stable networks by crosslinking DNA monomers. In contrast to the conventional DNA gelation techniques, our system used solid aggregates as a precursor for DNA microgels. Therefore, in situ gelation could be achieved by depositing aggregates on the desired substrate and subsequently initiating a phase transition. Our approach can expand the utility and functionality of DNA hydrogels by using more complex nucleic acid assemblies as unit structures and combining the technique with top-down microfabrication methods.     

Cooperativity-based modeling of heterotypic DNA nanostructure assembly

DNA origami is a robust method for the fabrication of nanoscale 2D and 3D objects with complex features and geometries. The process of DNA origami folding has been recently studied, however quantitative understanding of it is still elusive. Here, we describe a systematic quantification of the assembly process of DNA nanostructures, focusing on the heterotypic DNA junction—in which arms are unequal—as their basic building block. Using bulk fluorescence studies we tracked this process and identified multiple levels of cooperativity from the arms in a single junction to neighboring junctions in a large DNA origami object, demonstrating that cooperativity is a central underlying mechanism in the process of DNA nanostructure assembly. We show that the assembly of junctions in which the arms are consecutively ordered is more efficient than junctions with randomly-ordered components, with the latter showing assembly through several alternative trajectories as a potential mechanism explaining the lower efficiency. This highlights consecutiveness as a new design consideration that could be implemented in DNA nanotechnology CAD tools to produce more efficient and high-yield designs. Altogether, our experimental findings allowed us to devise a quantitative, cooperativity-based heuristic model for the assembly of DNA nanostructures, which is highly consistent with experimental observations.     

Rapid prototyping of 3D DNA-origami shapes with caDNAno

DNA nanotechnology exploits the programmable specificity afforded by base-pairing to produce self-assembling macromolecular objects of custom shape. For building megadalton-scale DNA nanostructures, a long ‘scaffold’ strand can be employed to template the assembly of hundreds of oligonucleotide ‘staple’ strands into a planar antiparallel array of cross-linked helices. We recently adapted this ‘scaffolded DNA origami’ method to producing 3D shapes formed as pleated layers of double helices constrained to a honeycomb lattice. However, completing the required design steps can be cumbersome and time-consuming. Here we present caDNAno, an open-source software package with a graphical user interface that aids in the design of DNA sequences for folding 3D honeycomb-pleated shapes A series of rectangular-block motifs were designed, assembled, and analyzed to identify a well-behaved motif that could serve as a building block for future studies. The use of caDNAno significantly reduces the effort required to design 3D DNA-origami structures. The software is available at http://cadnano.org/, along with example designs and video tutorials demonstrating their construction. The source code is released under the MIT license.     

Automated macromolecular model building for X-ray crystallography using ARP/wARP version 7

ARP/wARP is a software suite to build macromolecular models in X-ray crystallography electron density maps. Structural genomics initiatives and the study of complex macromolecular assemblies and membrane proteins all rely on advanced methods for 3D structure determination. ARP/wARP meets these needs by providing the tools to obtain a macromolecular model automatically, with a reproducible computational procedure. ARP/wARP 7.0 tackles several tasks: iterative protein model building including a high-level decision-making control module; fast construction of the secondary structure of a protein; building flexible loops in alternate conformations; fully automated placement of ligands, including a choice of the best-fitting ligand from a 'cocktail'; and finding ordered water molecules. All protocols are easy to handle by a nonexpert user through a graphical user interface or a command line. The time required is typically a few minutes although iterative model building may take a few hours.     

TOPALi: software for automatic identification of recombinant sequences within DNA multiple alignments

SUMMARY: TOPALi is a new Java graphical analysis application that allows the user to identify recombinant sequences within a DNA multiple alignment (either automatically or via manual investigation). TOPALi allows a choice of three statistical methods to predict the positions of breakpoints due to past recombination. The breakpoint predictions are then used to identify putative recombinant sequences and their relationships to other sequences. In addition to its sophisticated interface, TOPALi can import many sequence formats, estimate and display phylogenetic trees and allow interactive analysis and/or automatic HTML report generation. AVAILABILITY: TOPALi is freely available from http://www.bioss.ac.uk/software.html     

Molecular scale architecture: engineered three- and four-way junctions

Biomolecular self-assembly provides a basis for the bottom-up construction of useful and diverse nanoscale architectures. DNA is commonly used to create these assemblies and is often exploited as a lattice or an array. Although geometrically rigid and highly predictable, these sheets of repetitive constructs often lack the ability to be enzymatically manipulated or elongated by standard biochemical techniques. Here, we describe two approaches for the construction of position-controlled, molecular-scale, discrete, three- and four-way DNA junctions. The first approach for constructing these junctions relies on the use of nonmigrating cruciforms generated from synthetic oligonucleotides to which large, biologically generated, double-stranded DNA segments are enzymatically ligated. The second approach utilitizes the DNA methyltransferase-based SMILing (sequence-specific methyltransferase-induced labeling of DNA) method to site-specifically incorporate a biotin within biologically derived DNA. Streptavidin is then used to form junctions between unique DNA strands. The resultant assemblies have precise and predetermined connections with lengths that can be varied by enzymatic or hybridization techniques, or geometrically controlled with standard DNA functionalization methods. These junctions are positioned with single nucleotide resolution on large, micrometer-length templates. Both approaches generate DNA assemblies which are fully compatible with standard recombinant methods and thus provide a novel basis for nanoengineering applications.     

DNAlive: a tool for the physical analysis of DNA at the genomic scale

SUMMARY: DNAlive is a tool for the analysis and graphical display of structural and physical characteristics of genomic DNA. The web server implements a wide repertoire of metrics to derive physical information from DNA sequences with a powerful interface to derive 3D information on large sequences of both naked and protein-bound DNAs. Furthermore, it implements a mesoscopic Metropolis code which allows the inexpensive study of the dynamic properties of chromatin fibers. In addition, our server also surveys other protein and genomic databases allowing the user to combine and explore the physical properties of selected DNA in the context of functional features annotated on those regions. AVAILABILITY: http://mmb.pcb.ub.es/DNAlive/ ; http://www.inab.org/     

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.     

Evidence for long poly(dA).poly(dT) tracts in D. discoideum DNA at high frequencies and their preferential avoidance of nucleosomal DNA core regions

The eukaryote, Dictyostelium discoideum, has one of the most (A+T) rich genomes studied to date. Isolated nuclear D. discoideum DNA (AX3 strain) was used to qualitatively determine the frequency and length distribution of long (dA).(dT) homopolymer tracts in this genome, in comparison to the less (A+T) rich calf thymus and Schistosoma mansoni DNAs that had few observable long tracts. These experimental data accurately reflect the significantly elevated frequencies of long tracts found computationally within the D. discoideum intron and flanking sequences, but not exons. PCR amplification of long (dA).(dT) homopolymer tract containing sequences was carried out. Then experimental biotinylated (dT)18 probe hybridization to the PCR amplified DNA showed that the long (dA).(dT) homopolymer tracts were enriched in D. discoideum sequences only hundreds of base pair in length, under conditions where no equivalent hybridization was observed to S. mansoni DNA or calf DNA sequences. Similar probe hybridization to DNA isolated following micrococcal nuclease digestion of D. discoideum chromatin demonstrated that long (dA).(dT) homopolymer tracts were more highly enriched in nucleosomal DNA lengths that included the internucleosomal linker as compared to shorter linker free mononucleosomal lengths. This observation is in agreement with the frequency of tract spacing results calculated from GenBank sequence data. These frequency data indicate that adjacent long tracts plus the intervening spacer DNA are found at peak lengths (average 42 bp), exactly characteristic of the internucleosomal spacer region of D. discoideum chromatin and are in sufficient number to be found in nearly half of all nucleosomes. Compared to shuffled tract sequence controls, these lengths of adjacent long tracts plus the intervening spacer DNA were found to be significantly enriched. Lesser enrichments are observed at lengths corresponding to adjacent tracts being separated by nucleosomal core length DNA sequences (145-185 bp). These data strongly suggest that adjacent long tracts occur spaced at selected lengths so as to avoid the central core regions of nucleosomes and instead are found localized within internucleosomal DNA linker and core edge regions in D. discoideum chromatin.     

Automated 3D tree-ring detection and measurement from X-ray computed tomography

Tree ring analysis is essential to reveal the environmental information encoded in the wood structure. It provides quantitative data on the anatomical structure which can be used, for example, to measure the impact of the fluctuating environment on the tree growth, to support global vegetation models and for the dendrochronological analysis of archaeological wooden artefacts. Currently, several imaging-based methods for tree-ring detection and tree-ring feature estimation exist. However, despite advances in computer vision and edge recognition algorithms, detection of tree-rings is mostly limited to two-dimensional (2D) datasets and performed manually in some cases. This paper describes a new approach to estimate the three-dimensional (3D) structure of tree rings and their width automatically from X-ray computed tomography data. This approach relies on a modified Canny edge detection algorithm, which is capable of detecting fully connected tree-ring edges throughout the image stack. Our results show that this approach performs well on six tree species having conifer, ring-porous and diffuse-porous ring boundary structures. In our study, image denoising proved to be a critical step to achieve accurate results. A major advantage of this procedure is that it requires very little to no user interaction rendering it a reproducible procedure for tree-ring width measurements. As it also provides 3D representations of the ring edges, it also may be used in the future for the inspection of anatomical features.     

A new DNA sequence assembly program.

We describe the Genome Assembly Program (GAP), a new program for DNA sequence assembly. The program is suitable for large and small projects, a variety of strategies and can handle data from a range of sequencing instruments. It retains the useful components of our previous work, but includes many novel ideas and methods. Many of these methods have been made possible by the program's completely new, and highly interactive, graphical user interface. The program provides many visual clues to the current state of a sequencing project and allows users to interact in intuitive and graphical ways with their data. The program has tools to display and manipulate the various types of data that help to solve and check difficult assemblies, particularly those in repetitive genomes. We have introduced the following new displays: the Contig Selector, the Contig Comparator, the Template Display, the Restriction Enzyme Map and the Stop Codon Map. We have also made it possible to have any number of Contig Editors and Contig Joining Editors running simultaneously even on the same contig. The program also includes a new 'Directed Assembly' algorithm and routines for automatically detecting unfinished segments of sequence, to which it suggests experimental solutions.     

ProfDist: a tool for the construction of large phylogenetic trees based on profile distances

SUMMARY: ProfDist is a user-friendly software package using the profile-neighbor-joining method (PNJ) in inferring phylogenies based on profile distances on DNA or RNA sequences. It is a tool for reconstructing and visualizing large phylogenetic trees providing new and standard features with a special focus on time efficency, robustness and accuracy. AVAILABILITY: A Windows version of ProfDist comes with a graphical user interface and is freely available at http://profdist.bioapps.biozentrum.uni-wuerzburg.de     

C-curve: A novel 3D graphical representation of DNA sequence based on codons

In this paper, a novel 3D graphical representation of DNA sequence based on codons is proposed. Since there is not loss of information due to overlapping and containing loops, this representation will be useful for comparison of different DNA sequences. This 3D curve will be convenient for DNA mutations comparison specially. In continues we give a numerical characterization of DNA sequences based on the new 3D curve. This characterization facilitates quantitative comparisons of similarities/dissimilarities analysis of DNA sequences based on codons.     

Primer design for large scale sequencing.

We have developed PRIDE, a primer design program that automatically designs primers in single contigs or whole sequencing projects to extend the already known sequence and to double strand single-stranded regions. The program is fully integrated into the Staden package (GAP4) and accessible with a graphical user interface. PRIDE uses a fuzzy logic-based system to calculate primer qualities. The computational performance of PRIDE is enhanced by using suffix trees to store the huge amount of data being produced. A test set of 110 sequencing primers and 11 PCR primer pairs has been designed on genomic templates, cDNAs and sequences containing repetitive elements to analyze PRIDE's success rate. The high performance of PRIDE, combined with its minimal requirement of user interaction and its fast algorithm, make this program useful for the large scale design of primers, especially in large sequencing projects.     

基于Web的基因组序列数据库管理系统的设计与实现

设计一种基于网络的可用来存储和注释海量DNA数据的数据库模型。整个过程分为三部分：首先是构建数据库框架,然后对原始基因组序列数据进行批量注释并输出有效格式导入数据库,最后通过一个友好的用户交互界面,实现对基因组数据的在线读取,查询,注释等操作。设计的数据库用于解决大量产生并有待分析的基因组序列的有效存储和管理问题。     

Geno3D: automatic comparative molecular modelling of protein

Geno3D (http://geno3d-pbil.ibcp.fr) is an automatic web server for protein molecular modelling. Starting with a query protein sequence, the server performs the homology modelling in six successive steps: (i) identify homologous proteins with known 3D structures by using PSI-BLAST; (ii) provide the user all potential templates through a very convenient user interface for target selection; (iii) perform the alignment of both query and subject sequences; (iv) extract geometrical restraints (dihedral angles and distances) for corresponding atoms between the query and the template; (v) perform the 3D construction of the protein by using a distance geometry approach and (vi) finally send the results by e-mail to the user.     

Development of a simple coarse-grained DNA model for analysis of oligonucleotide complex formation

In this paper, we develop a coarse-grained nucleotide model for the purpose of simulating large-scale aptamer-based hydrogel network formation in future research. In the model, each nucleotide is represented by a single interaction site containing sugar, phosphate and base. Discontinuous molecular dynamics (DMD) simulations are performed to simulate formation and denaturation of oligonucleotide duplexes as a function of temperature. The simulated melting temperatures of oligonucleotide duplexes are calculated in simulations of systems with different sequences, lengths and concentrations of oligonucleotides, and compared to data from the OligoAnalyzer tool. The denaturation of oligonucleotide triplexes containing a hybridised structure of three different oligonucleotides is analysed using both simulations and experiments. The nucleotide model is found to be a good predictor of the oligonucleotide’s hybridised state for both duplexes and triplexes. This coarse-grained model has wide ranging applications in the development or optimisation of DNA-based technologies including DNA origami, DNA-enabled hydrogels and DNA-based biosensors.