Showing papers on "Molecular models of DNA published in 2014"

All-atom polarizable force field for DNA based on the classical drude oscillator model

Abstract: Presented is a first generation atomistic force field for DNA in which electronic polarization is modeled based on the classical Drude oscillator formalism. The DNA model is based on parameters for small molecules representative of nucleic acids, including alkanes, ethers, dimethylphosphate, and the nucleic acid bases and empirical adjustment of key dihedral parameters associated with the phosphodiester backbone, glycosidic linkages and sugar moiety of DNA. Our optimization strategy is based on achieving a compromise between satisfying the properties of the underlying model compounds in the gas phase targeting QM data and reproducing a number of experimental properties of DNA duplexes in the condensed phase. The resulting Drude force field yields stable DNA duplexes on the 100 ns time scale and satisfactorily reproduces (1) the equilibrium between A and B forms of DNA and (2) transitions between the BI and BII sub-states of B form DNA. Consistency with the gas phase QM data for the model compounds is significantly better for the Drude model as compared to the CHARMM36 additive force field, which is suggested to be due to the improved response of the model to changes in the environment associated with the explicit inclusion of polarizability. Analysis of dipole moments associated with the nucleic acid bases shows the Drude model to have significantly larger values than those present in CHARMM36, with the dipoles of individual bases undergoing significant variations during the MD simulations. Additionally, the dipole moment of water was observed to be perturbed in the grooves of DNA.

"Nano-oddities": unusual nucleic acid assemblies for DNA-based nanostructures and nanodevices.

Abstract: CONSPECTUS: DNA is an attractive polymer building material for nanodevices and nanostructures due to its ability for self-recognition and self-assembly. Assembly relies on the formation of base-specific interactions that allow strands to adopt structures in a controllable fashion. Most DNA-based higher order structures such as DNA cages, 2D and 3D DNA crystals, or origamis are based on DNA double helices stabilized by Watson-Crick complementarity. A number of nonclassical pairing patterns are possible between or among DNA strands; these interactions result in formation of unusual structures that include, but are not limited to, G-quadruplexes, i-motifs, triplexes, and parallel-stranded duplexes. These structures create greater diversity of DNA-based building blocks for nanomaterials and have certain advantages over conventional duplex DNA, such as enhanced thermal stability and sensitivity to chemical stimuli. In this Account, we briefly introduce these alternative DNA structures and describe in detail their utilization in a variety of nanomaterials and nanomachines. The field of DNA "nano-oddities" emerged in the late 1990s when for the first time a DNA nanomachine was designed based on equilibrium between B-DNA and noncanonical, left-handed Z-DNA. Soon after, "proof-of-principle" DNA nanomachines based on several DNA "oddities" were reported. These machines were set in motion by the addition of complementary strands (a principle used by many B-DNA-based nanodevices), by the addition of selected cations, small molecules, or proteins, or by a change in pH or temperature. Today, we have fair understanding of the mechanism of action of these devices, excellent control over their performance, and knowledge of basic principles of their design. pH sensors and pH-controlled devices occupy a central niche in the field. They are usually based on i-motifs or triplex DNA, are amazingly simple, robust, and reversible, and create no waste apart from salt and water. G-quadruplex based nanostructures have unusually high stability, resist DNase and temperature, and display high selectivity toward certain cations. The true power of using these "nano-oddities" comes from combining them with existing nanomaterials (e.g., DNA origami, gold nanoparticles, graphene oxide, or mesoporous silica) and integrating them into existing mechanical and optoelectronic devices. Creating well-structured junctions for these interfaces, finding appropriate applications for the vast numbers of reported "nano-oddities", and proving their biological innocence comprise major challenges in the field. Our Account is not meant to be an all-inclusive review of the field but should give a reader a firm grasp of the current state of DNA nanotechnology based on noncanonical DNA structures.

Developmental self-assembly of a DNA tetrahedron.

Abstract: Kinetically controlled isothermal growth is fundamental to biological development, yet it remains challenging to rationally design molecular systems that self-assemble isothermally into complex geometries via prescribed assembly and disassembly pathways. By exploiting the programmable chemistry of base pairing, sophisticated spatial and temporal control have been demonstrated in DNA self-assembly, but largely as separate pursuits. By integrating temporal with spatial control, here we demonstrate the “developmental” self-assembly of a DNA tetrahedron, where a prescriptive molecular program orchestrates the kinetic pathways by which DNA molecules isothermally self-assemble into a well-defined three-dimensional wireframe geometry. In this reaction, nine DNA reactants initially coexist metastably, but upon catalysis by a DNA initiator molecule, navigate 24 individually characterizable intermediate states via prescribed assembly pathways, organized both in series and in parallel, to arrive at the tetrahedral final product. In contrast to previous work on dynamic DNA nanotechnology, this developmental program coordinates growth of ringed substructures into a three-dimensional wireframe superstructure, taking a step toward the goal of kinetically controlled isothermal growth of complex three-dimensional geometries.

A Coarse-Grained DNA Model Parameterized from Atomistic Simulations by Inverse Monte Carlo

Abstract: Computer modeling of very large biomolecular systems, such as long DNA polyelectrolytes or protein-DNA complex-like chromatin cannot reach all-atom resolution in a foreseeable future and this necessitates the development of coarse-grained (CG) approximations. DNA is both highly charged and mechanically rigid semi-flexible polymer and adequate DNA modeling requires a correct description of both its structural stiffness and salt-dependent electrostatic forces. Here, we present a novel CG model of DNA that approximates the DNA polymer as a chain of 5-bead units. Each unit represents two DNA base pairs with one central bead for bases and pentose moieties and four others for phosphate groups. Charges, intra- and inter-molecular force field potentials for the CG DNA model were calculated using the inverse Monte Carlo method from all atom molecular dynamic (MD) simulations of 22 bp DNA oligonucleotides. The CG model was tested by performing dielectric continuum Langevin MD simulations of a 200 bp double helix DNA in solutions of monovalent salt with explicit ions. Excellent agreement with experimental data was obtained for the dependence of the DNA persistent length on salt concentration in the range 0.1-100 mM. The new CG DNA model is suitable for modeling various biomolecular systems with adequate description of electrostatic and mechanical properties.

Molecular processes studied at a single-molecule level using DNA origami nanostructures and atomic force microscopy

Abstract: DNA origami nanostructures allow for the arrangement of different functionalities such as proteins, specific DNA structures, nanoparticles, and various chemical modifications with unprecedented precision. The arranged functional entities can be visualized by atomic force microscopy (AFM) which enables the study of molecular processes at a single-molecular level. Examples comprise the investigation of chemical reactions, electron-induced bond breaking, enzymatic binding and cleavage events, and conformational transitions in DNA. In this paper, we provide an overview of the advances achieved in the field of single-molecule investigations by applying atomic force microscopy to functionalized DNA origami substrates.

On the consistency of Monte Carlo track structure DNA damage simulations

Abstract: Purpose: Monte Carlo track structures (MCTS) simulations have been recognized as useful tools for radiobiological modeling. However, the authors noticed several issues regarding the consistency of reported data. Therefore, in this work, they analyze the impact of various user defined parameters on simulated direct DNA damage yields. In addition, they draw attention to discrepancies in published literature in DNA strand break (SB) yields and selected methodologies. Methods: The MCTS code Geant4-DNA was used to compare radial dose profiles in a nanometer-scale region of interest (ROI) for photon sources of varying sizes and energies. Then, electron tracks of 0.28 keV–220 keV were superimposed on a geometric DNA model composed of 2.7 × 106 nucleosomes, and SBs were simulated according to four definitions based on energy deposits or energy transfers in DNA strand targets compared to a threshold energy E TH. The SB frequencies and complexities in nucleosomes as a function of incident electron energies were obtained. SBs were classified into higher order clusters such as single and double strand breaks (SSBs and DSBs) based on inter-SB distances and on the number of affected strands. Results: Comparisons of different nonuniform dose distributions lacking charged particle equilibrium may lead to erroneous conclusions regarding the effect of energy on relative biological effectiveness. The energy transfer-based SB definitions give similar SB yields as the one based on energy deposit when E TH ≈ 10.79 eV, but deviate significantly for higher E TH values. Between 30 and 40 nucleosomes/Gy show at least one SB in the ROI. The number of nucleosomes that present a complex damage pattern of more than 2 SBs and the degree of complexity of the damage in these nucleosomes diminish as the incident electron energy increases. DNA damage classification into SSB and DSB is highly dependent on the definitions of these higher order structures and their implementations. The authors’ show that, for the four studied models, different yields are expected by up to 54% for SSBs and by up to 32% for DSBs, as a function of the incident electrons energy and of the models being compared. Conclusions: MCTS simulations allow to compare direct DNA damage types and complexities induced by ionizing radiation. However, simulation results depend to a large degree on user-defined parameters, definitions, and algorithms such as: DNA model, dose distribution, SB definition, and the DNA damage clustering algorithm. These interdependencies should be well controlled during the simulations and explicitly reported when comparing results to experiments or calculations.

Simulating a burnt-bridges DNA motor with a coarse-grained DNA model

Abstract: We apply a recently-developed coarse-grained model of DNA, designed to capture the basic physics of nanotechnological DNA systems, to the study of a `burnt-bridges' DNA motor consisting of a single-stranded cargo that steps processively along a track of single-stranded stators. We demonstrate that the model is able to simulate such a system, and investigate the sensitivity of the stepping process to the spatial separation of stators, finding that an increased distance can suppress successful steps due to the build up of unfavourable tension. The mechanism of suppression suggests that varying the distance between stators could be used as a method for improving signal-to-noise ratios for motors that are required to make a decision at a junction of stators.

Insights into DNA-mediated interparticle interactions from a coarse-grained model

Abstract: DNA-functionalized particles have great potential for the design of complex self-assembled materials. The major hurdle in realizing crystal structures from DNA-functionalized particles is expected to be kinetic barriers that trap the system in metastable amorphous states. Therefore, it is vital to explore the molecular details of particle assembly processes in order to understand the underlying mechanisms. Molecular simulations based on coarse-grained models can provide a convenient route to explore these details. Most of the currently available coarse-grained models of DNA-functionalized particles ignore key chemical and structural details of DNA behavior. These models therefore are limited in scope for studying experimental phenomena. In this paper, we present a new coarse-grained model of DNA-functionalized particles which incorporates some of the desired features of DNA behavior. The coarse-grained DNA model used here provides explicit DNA representation (at the nucleotide level) and complementary interactions between Watson-Crick base pairs, which lead to the formation of single-stranded hairpin and double-stranded DNA. Aggregation between multiple complementary strands is also prevented in our model. We study interactions between two DNA-functionalized particles as a function of DNA grafting density, lengths of the hybridizing and non-hybridizing parts of DNA, and temperature. The calculated free energies as a function of pair distance between particles qualitatively resemble experimental measurements of DNA-mediated pair interactions.

A Nanochannel Platform for Single DNA Studies: From Crowding, Protein DNA Interaction, to Sequencing of Genomic Information

Abstract: The study of nanochannel-confined DNA is important from biotechnological and biophysical points of view. We produce nanochannels in elastomer with soft lithography and proton beam writing. Issues concerning DNA confined in such quasi one-dimensional channels are discussed. We describe DNA stretching via the control of channel diameter and buffer conditions and how the extension can be interpreted with theory and computer simulation. We then discuss the conformation of nano-confined DNA crowded by neutral polymers and like-charged proteins. As an example of a protein that has an affinity to DNA, the effect of heat-stable nucleoid-structuring protein, H-NS, on the folding and compaction of DNA is reviewed. Compaction of DNA by eukaryotic protamine and unpacking of pre-compacted DNA through an increase in salt concentration are discussed. We review results obtained with a novel, cross-channel device that allows the monitoring of the dynamic, conformational response of DNA after exposure to a ligand or protein and/or a change in buffer conditions in situ. As a biotechnological application, linearization of DNA by bottlebrush coating with a polypeptide copolymer is discussed. It is demonstrated that large-scale genomic organization can be sequenced using single DNA molecules on an array of elastomeric nanochannels. Overall, our results show that the effects of ligands and proteins on the conformation, folding, and condensation of DNA are not only related to classical controlling factors, such as osmotic pressure, charge, and binding, but that the interplay with confinement in a nanospace is of paramount importance.

Probing transient protein-mediated DNA linkages using nanoconfinement

Abstract: We present an analytic technique for probing protein-catalyzed transient DNA loops that is based on nanofluidic channels. In these nanochannels, DNA is forced in a linear configuration that makes loops appear as folds whose size can easily be quantified. Using this technique, we study the interaction between T4 DNA ligase and DNA. We find that T4 DNA ligase binding changes the physical characteristics of the DNA polymer, in particular persistence length and effective width. We find that the rate of DNA fold unrolling is significantly reduced when T4 DNA ligase and ATP are applied to bare DNA. Together with evidence of T4 DNA ligase bridging two different segments of DNA based on AFM imaging, we thus conclude that ligase can transiently stabilize folded DNA configurations by coordinating genetically distant DNA stretches.

Progress in molecular modelling of DNA materials

Abstract: The unique molecular recognition properties of DNA molecule, which store genetic information in cells, are responsible for the rise of DNA nanotechnology. In this article, we review the recent adva...

Equilibration of complexes of DNA and H-NS proteins on charged surfaces: A coarse-grained model point of view

Abstract: The Histone-like Nucleoid Structuring protein (H-NS) is a nucleoid-associated protein, which is involved in both gene regulation and DNA compaction. Although it is a key player in genome organization by forming bridges between DNA duplexes, the precise structure of complexes of DNA and H-NS proteins is still not well understood. In particular, it is not clear whether the structure of DNA/H-NS complexes in the living cell is similar to that of complexes deposited on mica surfaces, which may be observed by AFM microscopy. A coarse-grained model, which helps getting more insight into this question, is described and analyzed in the present paper. This model is able of describing both the bridging of bacterial DNA by H-NS in the bulk and the deposition and equilibration of the complex on a charged surface. Simulations performed with the model reveal that a slight attraction between DNA and the charged surface is sufficient to let DNA/H-NS complexes reorganize from 3D coils to planar plasmids bridged by H-NS proteins similar to those observed by AFM microscopy. They furthermore highlight the antagonistic effects of the interactions between DNA and the surface. Indeed, increasing these interactions slows down the equilibration of naked plasmids on the surface but, on the other hand, enables a faster equilibration of DNA/H-NS complexes. Based on the distribution of the lifetimes of H-NS bridges and the time evolution of the number of trans-binding protein dimers during equilibration of the complexes on the surface, it is argued that the decrease of the equilibration time of the complex upon increase of the interaction strength between DNA and the surface is ascribable to the associated decrease of the probability to form new bridges between DNA and the proteins.

A new building block for DNA network formation by self-assembly and polymerase chain reaction.

Abstract: The predictability of DNA self-assembly is exploited in many nanotechnological approaches. Inspired by naturally existing selfassembled DNA architectures, branched DNA has been developed that allows self-assembly to predesigned architectures with dimensions on the nanometer scale. DNA is an attractive material for generation of nanostructures due to a plethora of enzymes which modify DNA with high accuracy, providing a toolbox for many different manipulations to construct nanometer scaled objects. We present a straightforward synthesis of a rigid DNA branching building block successfully used for the generation of DNA networks by self-assembly and network formation by enzymatic DNA synthesis. The Y-shaped 3-armed DNA construct, bearing 3 primer strands is accepted by Taq DNA polymerase. The enzyme uses each arm as primer strand and incorporates the branched construct into large assemblies during PCR. The networks were investigated by agarose gel electrophoresis, atomic force microscopy, dynamic light scattering, and electron paramagnetic resonance spectroscopy. The findings indicate that rather rigid DNA networks were formed. This presents a new bottom-up approach for DNA material formation and might find applications like in the generation of functional hydrogels.

Role of microscopic flexibility in tightly curved DNA

Abstract: The genetic material in living cells is organized into complex structures in which DNA is subjected to substantial contortions. Here we investigate the difference in structure, dynamics, and flexibility between two topological states of a short (107 base pair) DNA sequence in a linear form and a covalently closed, tightly curved circular DNA form. By employing a combination of all-atom molecular dynamics (MD) simulations and elastic rod modeling of DNA, which allows capturing microscopic details while monitoring the global dynamics, we demonstrate that in the highly curved regime the microscopic flexibility of the DNA drastically increases due to the local mobility of the duplex. By analyzing vibrational entropy and Lipari-Szabo NMR order parameters from the simulation data, we propose a novel model for the thermodynamic stability of high-curvature DNA states based on vibrational untightening of the duplex. This novel view of DNA bending provides a fundamental explanation that bridges the gap between classical models of DNA and experimental studies on DNA cyclization, which so far have been in substantial disagreement.

The "sugar" coarse-grained DNA model

Abstract: More than twenty coarse-grained (CG) DNA models have been developed for simulating the behavior of this molecule under various conditions, including those required for nanotechnology. However, none of these models reproduces the DNA polymorphism associated with conformational changes in the ribose rings of the DNA backbone. These changes make an essential contribution to the DNA local deformability and provide the possibility of the transition of the DNA double helix from the B-form to the A-form during interactions with biological molecules. We propose a CG representation of the ribose conformational flexibility. We substantiate the choice of the CG sites (6 per nucleotide) needed for the "sugar" GC DNA model, and obtain the potentials of the CG interactions between the sites by the "bottom-up" approach using the all-atom AMBER force field. We show that the representation of the ribose flexibility requires one non-harmonic and one three-particle potential, the forms of both the potentials being different from the ones generally used. The model also includes (i) explicit representation of ions (in an implicit solvent) and (ii) sequence dependence. With these features, the sugar CG DNA model reproduces (with the same parameters) both the B- and A- stable forms under corresponding conditions and demonstrates both the A to B and the B to A phase transitions.

DNA topology in chromosomes: a quantitative survey and its physiological implications

Abstract: Using a simple geometric model, we propose a general method for computing the linking number of the DNA embedded in chromatin fibers. The relevance of the method is reviewed through the single molecule experiments that have been performed in vitro with magnetic tweezers. We compute the linking number of the DNA in the manifold conformational states of the nucleosome which have been evidenced in these experiments and discuss the functional dynamics of chromosomes in the light of these manifold states.

Parametric Resonance in a Mesoscopic Discrete DNA Model

Abstract: In this paper we investigate from the numerical point of view the discrete DNA model proposed in Lacitignola and Saccomandi (Bull. Math. Biol., 2014) in order to test the robustness of the parametric resonance condition found in the limit of the continuum approximation. To describe more realistically the binding of RNA polymerase to the DNA macromolecule during the first stage of the transcription process, we here consider a localized DNA-RNA polymerase interaction and a relatively high number of base-pairs. Even with these more realistic assumptions, our findings confirm the ones found in the continuum limit and indicate that the parametric resonance phenomenon can be an intrinsic property of the discrete DNA model.

Comparison of DNA strand-break simulated with different DNA models

Abstract: In Monte Carlo simulation of DNA damage, the geometric model of DNA is of great importance. To study the influence of DNA model on the simulation of DNA damage, three DNA models were created in this paper. They were a volume model and two atomic models with different parameters. Direct DNA strand-break induced by low-energy electrons were simulated respectively with the three models. The results show that most of the energy depositions in the DNA segments do not lead to strand-breaks. The simple single strand-break (SSB) tends to be the predominant damage type, and the contribution of complex double strand-break (DSB) to the total DSB cannot be neglected. Among the yields of all the three DNA target models applied here, the yields of the volume model are the highest, the yields of the atomic model with double van der Waals radii (r ) take the second place, whereas the yields of the atomic model with single r come last. On average, the ratios of SSB yields are approximately equivalent to the corresponding ratios of the models’ volume. However, there seems to be no clear relationship between the DSB yields and the models’ volume.

A DNA model for solving the hierarchical clustering problem

Abstract: The hierarchical clustering algorithm is one of the clustering algorithms used widely. The typical character of the hierarchical clustering algorithm is to find the “shortest” patterns. But it is not to ensure that to find the global optimal result by adding the shortest pattern each time. In this paper, we use DNA computing to solve this problem. Because during the ligation of the DNA strands, the bio-reaction's character is to combine all the possible single DNA strands into the shortest double DNA strand. We propose a DNA model with the coding strategy and DNA computing algorithm for the hierarchical clustering problem. we use this algorithm to solve the real problem — the retailer's region partition problem. Although we use the computer to simulate the bio-chemical reactions, we will execute this algorithm in the biological lab in the near future.

Coarse-grained DNA model capable of simulating ribose flexibility

Abstract: More than twenty coarse-grained (CG) DNA models have been developed for simulating the behavior of this molecule under various conditions, including those required for nanotechnology. However, none of these models reproduces the DNA polymorphism associated with conformational changes in the ribose rings of the DNA backbone. These changes make an essential contribution to the DNA local deformability and provide the possibility of the transition of the DNA double helix from the B-form to the A-form during interactions with biological molecules. We propose a CG representation of the ribose conformational flexibility. We substantiate the choice of the CG sites (6 per nucleotide) needed for the "sugar" GC DNA model, and obtain the potentials of the CG interactions between the sites by the "bottom-up" approach using the all-atom AMBER force field. We show that the representation of the ribose flexibility requires one non-harmonic and one three-particle potential, the forms of both the potentials being different from the ones generally used. The model also includes (i) explicit representation of ions (in an implicit solvent) and (ii) sequence dependence. With these features, the sugar CG DNA model reproduces (with the same parameters) both the B- and A- stable forms under corresponding conditions and demonstrates both the A to B and the B to A phase transitions.