Understanding DNA-Based Nanostructures using Molecular Simulation

Deoxyribonucleic acid (DNA) is arguably the most studied and most important biological molecule. Recently, it has also been established as a potential candidate for nanoconstruction. Self-assembly of DNA molecules has emerged as a simple yet elegant technique to organize matter with sub-nanometer precision. The unique base-pairing properties which helps DNA to carry genetic information, also makes it a suitable building block for creating stable and robust nanostructures. Recent decades have witnessed a major revolution in the synthesis of different topological structures made of DNA molecules at nanoscale like, two dimensional arrays, nanotubes, polyhedra, smiley faces, three dimensional crystals etc. Due to their easier design, high fidelity and automated chemistry, DNA nanostructures have proposed applications in diverse fields of bio-nanotechnology and synthetic biology. The field of structural DNA nanotechnology is just entering in adulthood and offer paramount challenge towards the journey of DNA-based nanostructures from the laboratory to their practical implementation in the real world. The aim of my dissertation is to develop a de novo computational framework to investigate the nanoscale structure and properties of DNA-based nanostructures. This will help to understand the molecular origin of interaction governing the structure and stability of DNA nanostructures. In this thesis, we have studied the in-solution behavior of self-assembled DNA nanostructures. The state of art all atom molecular dynamics (MD) simulation has been extensively implemented to understand the various thermodynamic properties of these self-assembled soft matter systems. We expect that the results presented here will lead to better design of self-assembled DNA nanostructures to address the real world challenges. In particular, we have developed algorithms to build very accurate atomistic models of various DNA nanostructures like crossover DNA molecules, DNA nanotubes (DNTs) and DNA icosahedron (IDNA). Further, we discuss a computational framework to understand the in situ structure and dynamics of these DNA nanostructures using state-of-art MD simulation. We carried out several hundred nanosecond long MD simulations on these systems which sometimes contains close to one million atoms. Following the trajectories of nanostructures in physiological conditions, we predicted numerous properties like equilibrium solution structure behavior and elastic properties which are difficult to measure in experiments. DNTs are self-assembled tubular templates where the circular double helical domains, kept at the vertices of a polygon, are connected at crossovers junctions. Ned Seeman and co-workers at New York University have synthesized different kind of DNTs using tile-based self-assembly of oligonucleotides. To investigate their microscopic structure, stability and mechanical properties, we have come up with 3d atomistic models of various DNTs which will facilitate further studies of these nanotubes towards their proposed nanotechnological and biological applications. In chapter 3 of this thesis, we discuss the analysis of several nanoseconds long all-atom MD simulation trajectories of various DNTS in the presence of explicit salt solution. We conclude that 6-helix DNT (6HB) structures are most stable and well behaved due to the better crossover designs and geometry. There has been considerable interest to investigate and enhance the mechanical strengths of DNTs to create rigid motifs. One simple way to increase the rigidity is to add further helices to the 6HB, which is known to be the most stable design of DNT, with the same tile-based crossover method. In chapter 4, we report atomistic models of 6HB flanked symmetrically with two double helical DNA pillars (6HB+2) and 6HB flanked symmetrically by three double helical DNA pillars (6HB+3). From the fluctuation analysis of the equilibrium MD simulation trajectories, we calculated the stretch modulus and persistence length of these DNTs. The measured persistence lengths of these nanotubes are ∼10 μm, which is 2 orders of magnitude larger than that of dsDNA. We also find a gradual increase of persistence length with an increasing number of pillars, in quantitative agreement with previous experimental findings. We also carried out non-equilibrium Steered-Molecular-Dynamics (SMD) to measure the stretch modulus from the force-extension behavior of these pillared DNTs. The values of the stretch modulus calculated using contour length distribution of equilibrium MD simulations are similar to those obtained from non-equilibrium SMD simulations. The addition of pillars makes these DNTs very rigid. Engineering the synthetic nanopores through lipid bilayer membrane to access the interior of a cell is a long standing challenge in biotechnology. Recently, a new class of DNA nanopores through the lipid bilayer membranes has been characterized using advanced imaging techniques and transmembrane ionic current recordings. In chapter 5 of the thesis, we present a MD simulation study of 6HB embedded in POPC lipid bilayer membranes. The analyses of 0.2 µs long equilibrium MD simulation trajectories demonstrate that structure is stable and well behaved. We observe that the head groups of the lipid molecules close to DNT cooperatively tilt towards the hydrophilic sugar-phosphate backbone of DNA to form a toroidal structure around the patch of DNT protruding in the membrane. Based on this observation, we propose a new mechanism, which has been largely overlooked so far, to explain the stability to this DNA-lipid molecular self-assembly. We further explore the effect of monovalent ionic concentrations to the in-solution structure and stability of the nanocomposite. Transmembrane ionic current measurements during the constant electric field simulation provide the I-V characteristics of the water filled DNT lumen in lipid membrane. The conductivity of the DNT lumen turns out to be several nS and increases with ionic concentration. Recently, Krishnan’s research group at NCBS Bangalore and Chicago University have characterized DNA icosahedra (IDNA) using advance imaging techniques and validated it for biological targeting and bioimaging in vivo. A high resolution structural model of such polyhedra would be critical to widening their applications in both materials and biology. In chapter 6 of this thesis, we discuss an atomistic model of this well-characterized IDNA to study the in-solution behavior using MD simulation. We provide quantitative estimate of the surface and volume of the equilibrium structure which is essential to estimate its maximal cargo carrying capacity. Importantly, our simulation of gold nanoparticles (AuNP) encapsulated within DNA icosahedra (IAuNP) revealed enhanced stability of the AuNP loaded structure as compared to the empty icosahedra. This is consistent with experimental results that show high yields of cargo-encapsulated DNA icosahedra that have led to its diverse applications for precision targeting. These studies reveal that the stabilizing interactions between the cargo and the DNA scaffold powerfully positions DNA polyhedra as targetable nanocapsules for payload delivery. The insights from our study can be further exploited for precise molecular display for diverse biological applications. Finally, in chapter 7, we give a summary of the main results presented in this dissertation. We also briefly discuss the ongoing research work and the bright future of this emerging field of DNA nanotechnology. We believe that this thesis deepens the microscopic understanding of the recent experimental observation and provides impetus in the real world application of DNA nanostructures in vitro and well as in vivo