Close encounters with DNA

The Potential of Liquorpheresis to Treat Leptomeningeal Disease

Leptomeningeal disease (LMD) is a devastating sequela of many cancers, with an extremely poor prognosis. Barriers to improving outcomes are related to the inability of many traditional therapies to effectively reach the cerebrospinal fluid (CSF) space within the central nervous system. Liquorpheresis is an emerging treatment modality specific to CSF diseases, the primary mechanism of action of which is direct targeted filtration of CSF content by neurosurgical access. In this review, we highlight the principles of liquorpheresis and detail how LMD can be amenable to this treatment. Further, we summarize the current in vitro and in vivo evidence supporting liquorpheresis as a feasible method to treat LMD and other central nervous system diseases as well as describe its conceivable limitations.

Comparison of gender identification using exfoliated cells obtained from toothbrush and miswak: A longitudinal study

Gender identification plays a pivotal role in forensic medicine. Among the various methods used for gender identification, deoxyribose nucleic acid (DNA) based methods are considered accurate. Exfoliated oral mucosal cells that are harvested from oral hygiene aids can be potentially used for gender identification using real-time polymerase chain rection (PCR). The aim of the present longitudinal study is to assess and compare the efficacy of toothbrush and miswak as potential tools to harvest exfoliated cells for gender identification. Forty healthy volunteers were recruited and asked to clean their teeth using new toothbrush and fresh miswak each day for 4 days. Toothbrush and miswak used by the participants were subjected to DNA analysis immediately, 1st, 2nd and 6th month. The absorbance of DNA samples were quantified and gender identification was done by amplification of sex determining gene-Sex determining region Y gene (SRY) and ALT1 genes using real-time PCR. The number of correct and positive identification for samples at various time points were tabulated and subjected to statistical analysis. Post hoc power analysis showed that the study had a power of 93%. Correct and positive gender identification was observed for the samples (100%) obtained using miswak, for tooth brush it reduced to 95%, 80%, and 35% at the end of 1st, 2nd, and 6th month. The differences seen at the end of 2nd month and 6th month were statistically significant. Miswak is a better tool to harvest exfoliated cells for gender identification when compared to a toothbrush. Hence, miswak can serve as a potential tool in forensic medicine for DNA extraction and subsequently victim identification.

Rapid Escherichia coli Cloned DNA Detection in Serum Using an Electrical Double Layer-Gated Field-Effect Transistor-Based DNA Sensor

In this study, a rapid diagnosis platform was developed for the detection of Escherichia coli O157:H7. An electrical double layer (EDL)-gated field-effect transistor-based biosensor (BioFET) as a point-of-care testing device is demonstrated with its high sensitivity, portability, high selectivity, quick response, and ease of use. The specially designed ssDNA probe was immobilized on the extended gate electrode to bind the target complementary DNA segment of E. coli, resulting in a sharp drain current change within minutes. The limit of detection for target DNA is validated to a concentration of 1 fM in buffer solution and serum. Meanwhile, the results of a Kelvin probe force microscope were shown to have reduced surface potential of the DNA immobilized sensors before and after the cDNA detection, which is consistent with the decreased drain current of the BioFET. A 1.2 kb E. coli duplex DNA synthesized in plasmid was sonicated and detected in serum samples with the sensor array. Gel electrophoresis was used to confirm the efficiency of sonication by elucidating the length of DNA. Those results show that the EDL-gated BioFET system is a promising platform for rapid identification of pathogens for future clinical needs.

New Insights into the Role of Ligand-Binding Modes in GC-DNA Condensation through Thermodynamic and Spectroscopic Studies

In biological systems, the unprompted assembly of DNA molecules by cationic ligands into condensed structures is ubiquitous. The ability of ligands to provoke DNA packaging is crucial to the molecular organization and functional control of DNA, yet their underlined physical roles have remained elusive. Here, we have examined the DNA condensation mechanism of four cationic ligands, including their primary DNA-binding modes through extensive biophysical studies. We observed contrasting changes for these ligands binding to poly[dGdC]:poly[dGdC] (GC-DNA) and poly[dAdT]:poly[dAdT] (AT-DNA). Based on a CD spectroscopic study, it was confirmed that only GC-DNA undergoes B- to Ψ-type DNA transformation in the presence of ligands. In the fluorescence displacement assay (FDA), the ability of ligands to displace GC-DNA-bound EtBr follows the order: protamine²¹⁺ > cohex³⁺ > Ni²⁺ > spermine⁴⁺, which indicates that there is no direct correlation between the ligand charge and its ability to displace the drug from the DNA, indicating that GC-DNA condensation is not just influenced by electrostatic interaction but ligand-specific interactions may also have played a crucial role. Furthermore, the detailed ITC-binding studies suggested that DNA-ligand interactions are generally driven by unfavorable enthalpy and favorable entropy. The correlations from various studies insinuate that cationic ligands show major groove binding as one of the preferred binding modes during GC-DNA condensation.

Computer Aided Development of Nucleic Acid Applications in Nanotechnologies

Utilization of nucleic acids (NAs) in nanotechnologies and nanotechnology-related applications is a growing field with broad application potential, ranging from biosensing up to targeted cell delivery. Computer simulations are useful techniques that can aid design and speed up development in this field. This review focuses on computer simulations of hybrid nanomaterials composed of NAs and other components. Current state-of-the-art molecular dynamics simulations, empirical force fields (FFs), and coarse-grained approaches for the description of deoxyribonucleic acid and ribonucleic acid are critically discussed. Challenges in combining biomacromolecular and nanomaterial FFs are emphasized. Recent applications of simulations for modeling NAs and their interactions with nano- and biomaterials are overviewed in the fields of sensing applications, targeted delivery, and NA templated materials. Future perspectives of development are also highlighted.

Capillary-Assisted Molecular Pendulum Bioanalysis

The development of robust biosensing strategies that can be easily implemented in everyday life remains a challenge for the future of modern biosensor research. While several reagentless approaches have attempted to address this challenge, they often achieve user-friendliness through sacrificing sensitivity or universality. While acceptable for certain applications, these trade-offs hinder the widespread adoption of reagentless biosensing technologies. Here, we report a novel approach to reagentless biosensing that achieves high sensitivity, rapid detection, and universality using the SARS-CoV-2 virus as a model target. Universality is achieved by using nanoscale molecular pendulums, which enables reagentless electrochemical biosensing through a variable antibody recognition element. Enhanced sensitivity and rapid detection are accomplished by incorporating the coffee-ring phenomenon into the sensing scheme, allowing for target preconcentration on a ring-shaped electrode. Using this approach, we obtained limits of detection of 1 fg/mL and 20 copies/mL for the SARS-CoV-2 nucleoproteins and viral particles, respectively. In addition, clinical sample analysis showed excellent agreement with Ct values from PCR-positive SARS-CoV-2 patients.

Equilibrium fluctuations of DNA plectonemes

Plectonemes are intertwined helically looped domains which form when a DNA molecule is supercoiled, i.e., over- or underwound. They are ubiquitous in cellular DNA, and their physical properties have attracted significant interest both from the experimental side and from the modeling side. In this paper, we investigate fluctuations of the end-point distance z of supercoiled linear DNA molecules subject to external stretching forces. Our analysis is based on a two-phase model, which describes the supercoiled DNA as composed of a stretched phase and a plectonemic phase. A variety of mechanisms are found to contribute to extension fluctuations, characterized by the variance 〈Δz^{2}〉. We find the dominant contribution to 〈Δz^{2}〉 to originate from phase-exchange fluctuations, the transient shrinking and expansion of plectonemes, which is accompanied by an exchange of molecular length between the two phases. We perform Monte Carlo simulations of the twistable wormlike chain and analyze the fluctuation of various quantities, the results of which are found to agree with the two-phase model predictions. Furthermore, we show that the extension and its variance at high forces are very well captured by the two-phase model, provided that one goes beyond quadratic approximations.

Systematic mapping of DNAzymes research from 1995 to 2019

DNAzymes (catalytic DNA) have gained significant diagnostic and therapeutic applications with increasing research output over the years. Functional oligonucleotides are used as molecular recognition elements within biosensors for detection of analytes and viral infections such as SARS-CoV-2. DNAzymes are also applied for silencing and regulating cancer specific genes. However, there has not been any report on systematic analysis to track research status, reveal hotspots, and map knowledge in this field. Therefore, in the present study, research articles on DNAzymes from 1995 to 2019 were extracted from Web of Science (SCI-Expanded) after which, 1037 articles were imported into Rstudio (version 3.6.2) and analysed accordingly. The highest number of articles was published in 2019 (n = 138), while the least was in 1995 (n = 1). The articles were published across 216 journals by 2344 authors with 2337 multi-author and 7 single authors. The most prolific authors were Li Y (n = 47), Liu J (n = 46), Wang L (n = 33), Willner I (n = 33) and Zhang L (n = 33). The top three most productive countries were China (n = 2018), USA (n = 447) and Canada (n = 251). The most productive institutions were Hunan University, China (n = 141), University of Illinois, USA (n = 139) and Fuzhou University, China (n = 101). Despite the increasing interest in this field, international collaborations between institutions were very low which requires immediate attention to mitigate challenges such as limited funding, access to facilities, and existing knowledge gap.

Toothbrushes as a Source of DNA for Gender and Human Identification-A Systematic Review

Background: Few studies have reported the use of toothbrushes as a reliable source of DNA for human or gender identification. The present systematic review with the available information was conducted to answer the focus question “Is a toothbrush a reliable source of DNA for human or gender identification?”. Methods: The keyword combination “Toothbrush” and “DNA” was used to search databases including MEDLINE, Scopus, and Web of Science along with a manual search of reference lists of relevant articles. Duplicates and irrelevant articles were excluded, and the remaining articles were fully read for the final selection of articles. The risk of bias of the included studies was evaluated using the Appraisal tool for Cross-Sectional Studies (AXIS tool). Results: Of the 130 articles obtained, 122 duplicates or irrelevant articles were eliminated. Following the full-text reading of eight articles, five articles were selected based on eligibility criteria. The five studies reported that a toothbrush is a good source of DNA irrespective of the time interval. In a few studies some samples were not sufficient for complete DNA profiling due to factors such as the method of DNA extraction. Conclusion: Although a toothbrush is an excellent source of DNA for human and gender identification, future studies with a larger sample size, appropriate control group, and standardized technique of DNA extraction need to be conducted. Additionally, factors influencing the quantity and quality of DNA in toothbrushes need to be determined with standardized techniques.

Ionic Atmosphere Effect on the Absorption Spectrum of a Flavoprotein: A Reminder to Consider Solution Ions

This study utilizes the FMN-dependent NADH:quinone oxidoreductase from Pseudomonas aeruginosa PAO1 to investigate the effect of introducing an active site negative charge on the flavin absorption spectrum both in the absence and presence of a long-range electrostatic potential coming from solution ions. There were no observed changes in the flavin UV-visible spectrum when an active site tyrosine (Y277) becomes deprotonated in vitro. These results could only be reproduced computationally using average solvent electrostatic configuration (ASEC) QM/MM simulations that include both positive and negative solution ions. The same calculations performed with minimal ions to neutralize the total protein charge predicted that deprotonating Y277 would significantly alter the flavin absorption spectrum. Analyzing the distribution of solution ions indicated that the ions reorganize around the protein surface upon Y277 deprotonation to cancel the effect of the tyrosinate on the flavin absorption spectrum. Additional biochemical experiments were performed to test this hypothesis.

Reagentless biomolecular analysis using a molecular pendulum

The development of reagentless sensors that can detect molecular analytes in biological fluids could enable a broad range of applications in personalized health monitoring. However, only a limited set of molecular inputs can currently be detected using reagentless sensors. Here, we report a sensing mechanism that is compatible with the analysis of proteins that are important physiological markers of stress, allergy, cardiovascular health, inflammation and cancer. The sensing method is based on the motion of an inverted molecular pendulum that exhibits field-induced transport modulated by the presence of a bound analyte. We measure the sensor's electric field-mediated transport using the electron-transfer kinetics of an attached reporter molecule. Using time-resolved electrochemical measurements that enable unidirectional motion of our sensor, the presence of an analyte bound to our sensor complex can be tracked continuously in real time. We show that this sensing approach is compatible with making measurements in blood, saliva, urine, tears and sweat and that the sensors can collect data in situ in living animals.

Binding Mode Multiplicity and Multiscale Chirality in the Supramolecular Assembly of DNA and a π-Conjugated Polymer

Water-soluble π-conjugated polymers are increasingly considered for DNA biosensing. However, the conformational rearrangement, supramolecular organization and dynamics upon interaction with DNA have been overlooked, which prevents the rational design of such detection tools. To elucidate the binding of a cationic polythiophene (CPT) to DNA with atomistic resolution, we performed molecular simulations of their supramolecular assembly. Comparison of replicated simulations show a multiplicity of CPT binding geometries that contribute to the wrapping of CPT around DNA. The different binding geometries are stabilized by both electrostatic interactions between CPT lateral cations and DNA phosphodiesters and van der Waals interactions between the CPT backbone and the DNA grooves. Simulated circular dichroism (CD) spectra show that the induced CD signal stems from a conserved geometrical feature across the replicated simulations, i. e. the presence of segments of syn configurations between thiophene units along the CPT chain. At the macromolecular scale, we inspected the different shapes related to the CPT binding modes around the DNA through symmetry metrics. Altogether, molecular dynamics (MD) simulations, model Hamiltonian calculations of the CD spectra, and symmetry indices provide insights into the origin of induced chirality from the atomic to the macromolecular scale. Our multidisciplinary approach points out the hierarchical aspect of CPT chiral organization induced by DNA.

Improving the Performance of the Amber RNA Force Field by Tuning the Hydrogen-Bonding Interactions

Molecular dynamics (MD) simulations became a leading tool for investigation of structural dynamics of nucleic acids. Despite recent efforts to improve the empirical potentials (force fields, ffs), RNA ffs have persisting deficiencies, which hamper their utilization in quantitatively accurate simulations. Previous studies have shown that at least two salient problems contribute to difficulties in the description of free-energy landscapes of small RNA motifs: (i) excessive stabilization of the unfolded single-stranded RNA ensemble by intramolecular base-phosphate and sugar-phosphate interactions and (ii) destabilization of the native folded state by underestimation of stability of base pairing. Here, we introduce a general ff term (gHBfix) that can selectively fine-tune nonbonding interaction terms in RNA ffs, in particular, the H bonds. The gHBfix potential affects the pairwise interactions between all possible pairs of the specific atom types, while all other interactions remain intact; i.e., it is not a structure-based model. In order to probe the ability of the gHBfix potential to refine the ff nonbonded terms, we performed an extensive set of folding simulations of RNA tetranucleotides and tetraloops. On the basis of these data, we propose particular gHBfix parameters to modify the AMBER RNA ff. The suggested parametrization significantly improves the agreement between experimental data and the simulation conformational ensembles, although our current ff version still remains far from being flawless. While attempts to tune the RNA ffs by conventional reparametrizations of dihedral potentials or nonbonded terms can lead to major undesired side effects, as we demonstrate for some recently published ffs, gHBfix has a clear promising potential to improve the ff performance while avoiding introduction of major new imbalances.

An Overview of Molecular Modeling for Drug Discovery with Specific Illustrative Examples of Applications

In this paper we review the current status of high-performance computing applications in the general area of drug discovery. We provide an introduction to the methodologies applied at atomic and molecular scales, followed by three specific examples of implementation of these tools. The first example describes in silico modeling of the adsorption of small molecules to organic and inorganic surfaces, which may be applied to drug delivery issues. The second example involves DNA translocation through nanopores with major significance to DNA sequencing efforts. The final example offers an overview of computer-aided drug design, with some illustrative examples of its usefulness.

Molecular Mechanism of Spontaneous Nucleosome Unraveling

Meters of DNA wrap around histone proteins to form nucleosomes and fit inside the micron-diameter nucleus. For the genetic information encoded in the DNA to become available for transcription, replication, and repair, the DNA-histone assembly must be disrupted. Experiment has indicated that the outer stretches of nucleosomal DNA "breathe" by spontaneously detaching from and reattaching to the histone core. Here, we report direct observation of spontaneous DNA breathing in atomistic molecular dynamics simulations, detailing a microscopic mechanism of the DNA breathing process. According to our simulations, the outer stretches of nucleosomal DNA detach in discrete steps involving 5 or 10 base pairs, with the detachment process being orchestrated by the motion of several conserved histone residues. The inner stretches of nucleosomal DNA are found to be more stably associated with the histone core by more abundant nonspecific DNA-protein contacts, providing a microscopic interpretation of nucleosome unraveling experiments. The CG content of nucleosomal DNA is found to anticorrelate with the extent of unwrapping, supporting the possibility that AT-rich segments may signal the start of transcription by forming less stable nucleosomes.

Charge Neutralization Drives the Shape Reconfiguration of DNA Nanotubes

Reconfiguration of membrane protein channels for gated transport is highly regulated under physiological conditions. However, a mechanistic understanding of such channels remains challenging owing to the difficulty in probing subtle gating-associated structural changes. Herein, we show that charge neutralization can drive the shape reconfiguration of a biomimetic 6-helix bundle DNA nanotube (6HB). Specifically, 6HB adopts a compact state when its charge is neutralized by Mg2+ ; whereas Na+ switches it to the expanded state, as revealed by MD simulations, small-angle X-ray scattering (SAXS), and FRET characterization. Furthermore, partial neutralization of the DNA backbone charges by chemical modification renders 6HB compact and insensitive to ions, suggesting an interplay between electrostatic and hydrophobic forces in the channels. This system provides a platform for understanding the structure-function relationship of biological channels and designing rules for the shape control of DNA nanostructures in biomedical applications.

A pH-induced conformational switch in a tyrosine kinase inhibitor identified by electronic spectroscopy and quantum chemical calculations

Tyrosine kinase inhibitors (TKIs) are a major class of drug utilised in the clinic. During transit to their cognate kinases, TKIs will encounter different pH environments that could have a major influence on TKI structure. To address this, we report UV-Vis spectroscopic and computational studies of the TKI, AG1478, as a function of pH. The electronic absorption spectrum of AG1478 shifted by 10 nm (from 342 nm to 332 nm) from acid to neutral pH and split into two peaks (at 334 nm and 345 nm) in highly alkaline conditions. From these transitions, the pKa value was calculated as 5.58 ± 0.01. To compute structures and spectra, time-dependent density functional theory (TD-DFT) calculations were performed along with conductor-like polarizable continuum model (CPCM) to account for implicit solvent effect. On the basis of the theoretical spectra, we could assign the AG1478 experimental spectrum at acidic pH to a mixture of two twisted conformers (71% AG1478 protonated at quinazolyl nitrogen N(1) and 29% AG1478 protonated at quinazolyl nitrogen N(3)) and at neutral pH to the neutral planar conformer. The AG1478 absorption spectrum (pH 13.3) was fitted to a mixture of neutral (70%) and NH-deprotonated species (30%). These studies reveal a pH-induced conformational transition in a TKI.

Order and interactions in DNA arrays: Multiscale molecular dynamics simulation

While densely packed DNA arrays are known to exhibit hexagonal and orthorhombic local packings, the detailed mechanism governing the associated phase transition remains rather elusive. Furthermore, at high densities the atomistic resolution is paramount to properly account for fine details, encompassing the DNA molecular order, the contingent ordering of counterions and the induced molecular ordering of the bathing solvent, bringing together electrostatic, steric, thermal and direct hydrogen-bonding interactions, resulting in the observed osmotic equation of state. We perform a multiscale simulation of dense DNA arrays by enclosing a set of 16 atomistically resolved DNA molecules within a semi-permeable membrane, allowing the passage of water and salt ions, and thus mimicking the behavior of DNA arrays subjected to external osmotic stress in a bathing solution of monovalent salt and multivalent counterions. By varying the DNA density, local packing symmetry, and counterion type, we obtain osmotic equation of state together with the hexagonal-orthorhombic phase transition, and full structural characterization of the DNA subphase in terms of its positional and angular orientational fluctuations, counterion distributions, and the solvent local dielectric response profile with its order parameters that allow us to identify the hydration force as the primary interaction mechanism at high DNA densities.

Single-molecule dissection of stacking forces in DNA

We directly measured at the single-molecule level the forces and lifetimes of DNA base-pair stacking interactions for all stack sequence combinations. Our experimental approach combined dual-beam optical tweezers with DNA origami components to allow positioning of blunt-end DNA helices so that the weak stacking force could be isolated. Base-pair stack arrays that lacked a covalent backbone connection spontaneously dissociated at average rates ranging from 0.02 to 500 per second, depending on the sequence combination and stack array size. Forces in the range from 2 to 8 piconewtons that act along the helical direction only mildly accelerated the stochastic unstacking process. The free-energy increments per stack that we estimate from the measured forward and backward kinetic rates ranged from -0.8 to -3.4 kilocalories per mole, depending on the sequence combination. Our data contributes to understanding the mechanics of DNA processing in biology, and it is helpful for designing the kinetics of DNA-based nanoscale devices according to user specifications.

Understanding the origin of liquid crystal ordering of ultrashort double-stranded DNA

Recent experiments have shown that short double-stranded DNA (dsDNA) fragments having six- to 20-base pairs exhibit various liquid crystalline phases. This violates the condition of minimum molecular shape anisotropy that analytical theories demand for liquid crystalline ordering. It has been hypothesized that the liquid crystalline ordering is the result of end-to-end stacking of dsDNA to form long supramolecular columns which satisfy the shape anisotropy criterion necessary for ordering. To probe the thermodynamic feasibility of this process, we perform molecular dynamics simulations on ultrashort (four base pair long) dsDNA fragments, quantify the strong end-to-end attraction between them, and demonstrate that the nematic ordering of the self-assembled stacked columns is retained for a large range of temperature and salt concentration.

Sequence-Dependent Mechanism of DNA Oligonucleotide Dehybridization Resolved through Infrared Spectroscopy

Despite its important role in biology and nanotechnology, many questions remain regarding the molecular mechanism and dynamics by which oligonucleotides recognize and hybridize to their complementary sequence. The thermodynamics and kinetics of DNA oligonucleotide hybridization and dehybridization are often assumed to involve an all-or-nothing two-state dissociation pathway, but deviations from this behavior can be considerable even for short sequences. We introduce a new strategy to characterize the base-pair-specific thermal dissociation mechanism of DNA oligonucleotides through steady-state and time-resolved infrared spectroscopy. Experiments are interpreted with a lattice model to provide a structure-specific interpretation. This method is applied to a model set of self-complementary 10-base-pair sequences in which the placement of GC base pairs is varied in an otherwise AT strand. Through a combination of Fourier transform infrared and two-dimensional infrared spectroscopy, experiments reveal varying degrees of deviation from simple two-state behavior. As the temperature is increased, duplexes dissociate through a path in which the terminal bases fray, without any significant contribution from loop configurations. Transient temperature jump experiments reveal time scales of 70-100 ns for fraying and 10-30 μs for complete dissociation near the melting temperature. Whether or not frayed states are metastable intermediates or short-lived configurations during the full dissociation of the duplex is dictated by the nucleobase sequence.

Multiscale coarse-grained modelling of chromatin components: DNA and the nucleosome

To model large biomolecular systems, such as cell and organelles an atomistic description is not currently achievable and is not generally practical. Therefore, simplified coarse-grained (CG) modelling becomes a necessity. One of the most important cellular components is chromatin, a large DNA-protein complex where DNA is highly compacted. Recent progress in coarse graining modelling of the major chromatin components, double helical DNA and the nucleosome core particle (NCP) is presented. First, general principles and approaches allowing rigorous bottom-to-top generation of interaction potentials in the CG models are presented. Then, recent CG models of DNA are reviewed and their adequacy is benchmarked against experimental data on the salt dependence of DNA flexibility (persistence length). Furthermore, a few recent CG models of the NCP are described and their application for studying salt-dependent NCP-NCP interaction is discussed. An example of a multiscale approach to CG modelling of chromatin is presented where interactions and self-assembly of thousands of NCPs in solution are observed.

Introducing improved structural properties and salt dependence into a coarse-grained model of DNA

We introduce an extended version of oxDNA, a coarse-grained model of deoxyribonucleic acid (DNA) designed to capture the thermodynamic, structural, and mechanical properties of single- and double-stranded DNA. By including explicit major and minor grooves and by slightly modifying the coaxial stacking and backbone-backbone interactions, we improve the ability of the model to treat large (kilobase-pair) structures, such as DNA origami, which are sensitive to these geometric features. Further, we extend the model, which was previously parameterised to just one salt concentration ([Na(+)] = 0.5M), so that it can be used for a range of salt concentrations including those corresponding to physiological conditions. Finally, we use new experimental data to parameterise the oxDNA potential so that consecutive adenine bases stack with a different strength to consecutive thymine bases, a feature which allows a more accurate treatment of systems where the flexibility of single-stranded regions is important. We illustrate the new possibilities opened up by the updated model, oxDNA2, by presenting results from simulations of the structure of large DNA objects and by using the model to investigate some salt-dependent properties of DNA.

Direct Simulation of the Self-Assembly of a Small DNA Origami

By using oxDNA, a coarse-grained nucleotide-level model of DNA, we are able to directly simulate the self-assembly of a small 384-base-pair origami from single-stranded scaffold and staple strands in solution. In general, we see attachment of new staple strands occurring in parallel, but with cooperativity evident for the binding of the second domain of a staple if the adjacent junction is already partially formed. For a system with exactly one copy of each staple strand, we observe a complete assembly pathway in an intermediate temperature window; at low temperatures successful assembly is prevented by misbonding while at higher temperature the free-energy barriers to assembly become too large for assembly on our simulation time scales. For high-concentration systems involving a large staple strand excess, we never see complete assembly because there are invariably instances where two copies of the same staple both bind to the scaffold, creating a kinetic trap that prevents the complete binding of either staple. This mutual staple blocking could also lead to aggregates of partially formed origamis in real systems, and helps to rationalize certain successful origami design strategies.

Free-energy landscape and characteristic forces for the initiation of DNA unzipping

DNA unzipping, the separation of its double helix into single strands, is crucial in modulating a host of genetic processes. Although the large-scale separation of double-stranded DNA has been studied with a variety of theoretical and experimental techniques, the minute details of the very first steps of unzipping are still unclear. Here, we use atomistic molecular-dynamics simulations, coarse-grained simulations, and a statistical-mechanical model to study the initiation of DNA unzipping by an external force. Calculation of the potential of mean force profiles for the initial separation of the first few terminal basepairs in a DNA oligomer revealed that forces ranging between 130 and 230 pN are needed to disrupt the first basepair, and these values are an order of magnitude larger than those needed to disrupt basepairs in partially unzipped DNA. The force peak has an echo of ∼50 pN at the distance that unzips the second basepair. We show that the high peak needed to initiate unzipping derives from a free-energy basin that is distinct from the basins of subsequent basepairs because of entropic contributions, and we highlight the microscopic origin of the peak. To our knowledge, our results suggest a new window of exploration for single-molecule experiments.

Computing Nonequilibrium Conformational Dynamics of Structured Nucleic Acid Assemblies

Synthetic nucleic acids can be programmed to form precise three-dimensional structures on the nanometer-scale. These thermodynamically stable complexes can serve as structural scaffolds to spatially organize functional molecules including multiple enzymes, chromophores, and force-sensing elements with internal dynamics that include substrate reaction-diffusion, excitonic energy transfer, and force-displacement response that often depend critically on both the local and global conformational dynamics of the nucleic acid assembly. However, high molecular weight assemblies exhibit long time-scale and large length-scale motions that cannot easily be sampled using all-atom computational procedures such as molecular dynamics. As an alternative, here we present a computational framework to compute the overdamped conformational dynamics of structured nucleic acid assemblies and apply it to a DNA-based tweezer, a nine-layer DNA origami ring, and a pointer-shaped DNA origami object, which consist of 204, 3,600, and over 7,000 basepairs, respectively. The framework employs a mechanical finite element model for the DNA nanostructure combined with an implicit solvent model to either simulate the Brownian dynamics of the assembly or alternatively compute its Brownian modes. Computational results are compared with an all-atom molecular dynamics simulation of the DNA-based tweezer. Several hundred microseconds of Brownian dynamics are simulated for the nine-layer ring origami object to reveal its long time-scale conformational dynamics, and the first ten Brownian modes of the pointer-shaped structure are predicted.

Refinement of the Sugar-Phosphate Backbone Torsion Beta for AMBER Force Fields Improves the Description of Z- and B-DNA

Z-DNA duplexes are a particularly complicated test case for current force fields. We performed a set of explicit solvent molecular dynamics (MD) simulations with various AMBER force field parametrizations including our recent refinements of the ε/ζ and glycosidic torsions. None of these force fields described the ZI/ZII and other backbone substates correctly, and all of them underpredicted the population of the important ZI substate. We show that this underprediction can be attributed to an inaccurate potential for the sugar-phosphate backbone torsion angle β. We suggest a refinement of this potential, β(OL1), which was derived using our recently introduced methodology that includes conformation-dependent solvation effects. The new potential significantly increases the stability of the dominant ZI backbone substate and improves the overall description of the Z-DNA backbone. It also has a positive (albeit small) impact on another important DNA form, the antiparallel guanine quadruplex (G-DNA), and improves the description of the canonical B-DNA backbone by increasing the population of BII backbone substates, providing a better agreement with experiment. We recommend using β(OL1) in combination with our previously introduced corrections, εζ(OL1) and χ(OL4), (the combination being named OL15) as a possible alternative to the current β torsion potential for more accurate modeling of nucleic acids.

Coarse-Grained Ions for Nucleic Acid Modeling

We present a general coarse-grained model of sodium, magnesium, spermidine, and chlorine in implicit solvent. The effective potentials between ions are systematically parametrized using a relative entropy coarse-graining approach [Carmichael, S. P. and M. S. Shell, J. Phys. Chem. B, 116, 8383-93 (2012)] that maximizes the information retained in a coarse-grained model. We describe the local distribution of ions in the vicinity of a recently published coarse-grained DNA model and demonstrate a dependence of persistence length on ionic strength that differs from that predicted by Odijk-Skolnick-Fixman theory. Consistent with experimental observations, we show that spermidine induces DNA condensation whereas magnesium and sodium do not. This model can be used alongside any coarse-grained DNA model that has explicit charges and an accurate reproduction of the excluded volume of dsDNA.

Computational Approaches to Nucleic Acid Origami

Recent advances in experimental DNA origami have dramatically expanded the horizon of DNA nanotechnology. Complex 3D suprastructures have been designed and developed using DNA origami with applications in biomaterial science, nanomedicine, nanorobotics, and molecular computation. Ribonucleic acid (RNA) origami has recently been realized as a new approach. Similar to DNA, RNA molecules can be designed to form complex 3D structures through complementary base pairings. RNA origami structures are, however, more compact and more thermodynamically stable due to RNA's non-canonical base pairing and tertiary interactions. With all these advantages, the development of RNA origami lags behind DNA origami by a large gap. Furthermore, although computational methods have proven to be effective in designing DNA and RNA origami structures and in their evaluation, advances in computational nucleic acid origami is even more limited. In this paper, we review major milestones in experimental and computational DNA and RNA origami and present current challenges in these fields. We believe collaboration between experimental nanotechnologists and computer scientists are critical for advancing these new research paradigms.

Extended molecular dynamics of a c-kit promoter quadruplex

The 22-mer c-kit promoter sequence folds into a parallel-stranded quadruplex with a unique structure, which has been elucidated by crystallographic and NMR methods and shows a high degree of structural conservation. We have carried out a series of extended (up to 10 μs long, ∼50 μs in total) molecular dynamics simulations to explore conformational stability and loop dynamics of this quadruplex. Unfolding no-salt simulations are consistent with a multi-pathway model of quadruplex folding and identify the single-nucleotide propeller loops as the most fragile part of the quadruplex. Thus, formation of propeller loops represents a peculiar atomistic aspect of quadruplex folding. Unbiased simulations reveal μs-scale transitions in the loops, which emphasizes the need for extended simulations in studies of quadruplex loops. We identify ion binding in the loops which may contribute to quadruplex stability. The long lateral-propeller loop is internally very stable but extensively fluctuates as a rigid entity. It creates a size-adaptable cleft between the loop and the stem, which can facilitate ligand binding. The stability gain by forming the internal network of GA base pairs and stacks of this loop may be dictating which of the many possible quadruplex topologies is observed in the ground state by this promoter quadruplex.

Martini Coarse-Grained Force Field: Extension to DNA

We systematically parameterized a coarse-grained (CG) model for DNA that is compatible with the Martini force field. The model maps each nucleotide into six to seven CG beads and is parameterized following the Martini philosophy. The CG nonbonded interactions are based on partitioning of the nucleobases between polar and nonpolar solvents as well as base-base potential of mean force calculations. The bonded interactions are fit to single-stranded DNA (ssDNA) atomistic simulations and an elastic network is used to retain double-stranded DNA (dsDNA) and other specific DNA conformations. We present the implementation of the Martini DNA model and demonstrate the properties of individual bases, ssDNA as well as dsDNA, and DNA-protein complexes. The model opens up large-scale simulations of DNA interacting with a wide range of other (bio)molecules that are available within the Martini framework.