Browse
Submit Data
Databases
API
Help

Dataset Information

0 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

A Monte Carlo Study of Knots in Long Double-Stranded DNA Chains.

ABSTRACT: We determine knotting probabilities and typical sizes of knots in double-stranded DNA for chains of up to half a million base pairs with computer simulations of a coarse-grained bead-stick model: Single trefoil knots and composite knots which include at least one trefoil as a prime factor are shown to be common in DNA chains exceeding 250,000 base pairs, assuming physiologically relevant salt conditions. The analysis is motivated by the emergence of DNA nanopore sequencing technology, as knots are a potential cause of erroneous nucleotide reads in nanopore sequencing devices and may severely limit read lengths in the foreseeable future. Even though our coarse-grained model is only based on experimental knotting probabilities of short DNA strands, it reproduces the correct persistence length of DNA. This indicates that knots are not only a fine gauge for structural properties, but a promising tool for the design of polymer models.

SUBMITTER: Rieger FC

PROVIDER: S-EPMC5025000 | biostudies-literature | 2016 Sep

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Publications

A Monte Carlo Study of Knots in Long Double-Stranded DNA Chains.

Rieger Florian C FC Virnau Peter P

PLoS computational biology 20160915 9

We determine knotting probabilities and typical sizes of knots in double-stranded DNA for chains of up to half a million base pairs with computer simulations of a coarse-grained bead-stick model: Single trefoil knots and composite knots which include at least one trefoil as a prime factor are shown to be common in DNA chains exceeding 250,000 base pairs, assuming physiologically relevant salt conditions. The analysis is motivated by the emergence of DNA nanopore sequencing technology, as knots a ...[more]

PMID: 27631891

Similar Datasets

Long-lived metastable knots in polyampholyte chains.

Project description:Knots in proteins and DNA are known to have significant effect on their equilibrium and dynamic properties as well as on their function. While knot dynamics and thermodynamics in electrically neutral and uniformly charged polymer chains are relatively well understood, proteins are generally polyampholytes, with varied charge distributions along their backbones. Here we use simulations of knotted polymer chains to show that variation in the charge distribution on a polyampholyte chain with zero net charge leads to significant variation in the resulting knot dynamics, with some charge distributions resulting in long-lived metastable knots that escape the (open-ended) chain on a timescale that is much longer than that for knots in electrically neutral chains. The knot dynamics in such systems can be described, quantitatively, using a simple one-dimensional model where the knot undergoes biased Brownian motion along a "reaction coordinate", equal to the knot size, in the presence of a potential of mean force. In this picture, long-lived knots result from charge sequences that create large electrostatic barriers to knot escape. This model allows us to predict knot lifetimes even when those times are not directly accessible by simulations.

| S-EPMC10266668 | biostudies-literature

Towards long double-stranded chains and robust DNA-based data storage using the random code system

Project description:DNA has become a popular choice for next-generation storage media due to its high storage density and stability. As the storage medium of life’s information, DNA has significant storage capacity and low-cost, low-power replication and transcription capabilities. However, utilizing long double-stranded DNA for storage can introduce unstable factors that make it difficult to meet the constraints of biological systems. To address this challenge, we have designed a highly robust coding scheme called the “random code system,” inspired by the idea of fountain codes. The random code system includes the establishment of a random matrix, Gaussian preprocessing, and random equilibrium. Compared to Luby transform codes (LT codes), random code (RC) has better robustness and recovery ability of lost information. In biological experiments, we successfully stored 29,390 bits of data in 25,700 bp chains, achieving a storage density of 1.78 bits per nucleotide. These results demonstrate the potential for using long double-stranded DNA and the random code system for robust DNA-based data storage.

| S-EPMC10294226 | biostudies-literature

Molecular renormalization group coarse-graining of polymer chains: application to double-stranded DNA.

Project description:Coarse-graining of atomistic force fields allows us to investigate complex biological problems, occurring at long timescales and large length scales. In this work, we have developed an accurate coarse-grained model for double-stranded DNA chain, derived systematically from atomistic simulations. Our approach is based on matching correlators obtained from atomistic and coarse-grained simulations, for observables that explicitly enter the coarse-grained Hamiltonian. We show that this requirement leads to equivalency of the corresponding partition functions, resulting in a one-step renormalization. Compared to prior works exploiting similar ideas, the main novelty of this work is the introduction of a highly compact set of Hamiltonian basis functions, based on molecular interaction potentials. We demonstrate that such compactification allows us to reproduce many-body effects, generated by one-step renormalization, at low computational cost. In addition, compact Hamiltonians greatly increase the likelihood of finding unique solutions for the coarse-grained force-field parameter values. By successfully applying our molecular renormalization group coarse-graining technique to double-stranded DNA, we solved, for the first time, a long-standing problem in coarse-graining polymer systems, namely, how to accurately capture the correlations among various polymeric degrees of freedom. Excellent agreement is found among atomistic and coarse-grained distribution functions for various structural observables, including those not included in the Hamiltonian. We also suggest higher-order generalization of this method, which may allow capturing more subtle correlations in biopolymer dynamics.

| S-EPMC2712212 | biostudies-literature

Double Excitation Energies from Quantum Monte Carlo Using State-Specific Energy Optimization.

Project description:We show that recently developed quantum Monte Carlo methods, which provide accurate vertical transition energies for single excitations, also successfully treat double excitations. We study the double excitations in medium-sized molecules, some of which are challenging for high-level coupled-cluster calculations to model accurately. Our fixed-node diffusion Monte Carlo excitation energies are in very good agreement with reliable benchmarks, when available, and provide accurate predictions for excitation energies of difficult systems where reference values are lacking.

| S-EPMC9648196 | biostudies-literature

MERIT: Controlling Monte-Carlo error rate in large-scale Monte-Carlo hypothesis testing.

Project description:The use of Monte-Carlo (MC) p$$ p $$ -values when testing the significance of a large number of hypotheses is now commonplace. In large-scale hypothesis testing, we will typically encounter at least some p$$ p $$ -values near the threshold of significance, which require a larger number of MC replicates than p$$ p $$ -values that are far from the threshold. As a result, some incorrect conclusions can be reached due to MC error alone; for hypotheses near the threshold, even a very large number (eg, 106$$ 1{0}^6 $$ ) of MC replicates may not be enough to guarantee conclusions reached using MC p$$ p $$ -values. Gandy and Hahn (GH)6-8 have developed the only method that directly addresses this problem. They defined a Monte-Carlo error rate (MCER) to be the probability that any decisions on accepting or rejecting a hypothesis based on MC p$$ p $$ -values are different from decisions based on ideal p$$ p $$ -values; their method then makes decisions by controlling the MCER. Unfortunately, the GH method is frequently very conservative, often making no rejections at all and leaving a large number of hypotheses "undecided". In this article, we propose MERIT, a method for large-scale MC hypothesis testing that also controls the MCER but is more statistically efficient than the GH method. Through extensive simulation studies, we demonstrate that MERIT controls the MCER while making more decisions that agree with the ideal p$$ p $$ -values than GH does. We also illustrate our method by an analysis of gene expression data from a prostate cancer study.

| S-EPMC10909586 | biostudies-literature

Streaming Variational Monte Carlo.

Project description:Nonlinear state-space models are powerful tools to describe dynamical structures in complex time series. In a streaming setting where data are processed one sample at a time, simultaneous inference of the state and its nonlinear dynamics has posed significant challenges in practice. We develop a novel online learning framework, leveraging variational inference and sequential Monte Carlo, which enables flexible and accurate Bayesian joint filtering. Our method provides an approximation of the filtering posterior which can be made arbitrarily close to the true filtering distribution for a wide class of dynamics models and observation models. Specifically, the proposed framework can efficiently approximate a posterior over the dynamics using sparse Gaussian processes, allowing for an interpretable model of the latent dynamics. Constant time complexity per sample makes our approach amenable to online learning scenarios and suitable for real-time applications.

| S-EPMC10082974 | biostudies-literature

Long-distance lateral diffusion of human Rad51 on double-stranded DNA.

Project description:Rad51 is the primary eukaryotic recombinase responsible for initiating DNA strand exchange during homologous recombination. Although the subject of intense study for over a decade, many molecular details of the reactions promoted by Rad51 and related recombinases remain unknown. Using total internal reflection fluorescence microscopy, we directly visualized the behavior of individual Rad51 complexes on double-stranded DNA (dsDNA) molecules suspended in an extended configuration above a lipid bilayer. Here we show that complexes of Rad51 can bind to and slide freely along the helical axis of dsDNA. Sliding is bidirectional, does not require ATP hydrolysis, and displays properties consistent with a 1D random walk driven solely by thermal diffusion. The proteins move freely on the DNA for long periods of time; however, sliding terminates and the proteins become immobile upon encountering the free end of a linear dsDNA molecule. This study provides previously uncharacterized insights into the behaviors of human Rad51, which may apply to other members of the RecA-like family of recombinases.

| S-EPMC1345706 | biostudies-literature

Thermodynamics of long supercoiled molecules: insights from highly efficient Monte Carlo simulations.

Project description:Supercoiled DNA polymer models for which the torsional energy depends on the total twist of molecules (Tw) are a priori well suited for thermodynamic analysis of long molecules. So far, nevertheless, the exact determination of Tw in these models has been based on a computation of the writhe of the molecules (Wr) by exploiting the conservation of the linking number, Lk=Tw+Wr, which reflects topological constraints coming from the helical nature of DNA. Because Wr is equal to the number of times the main axis of a DNA molecule winds around itself, current Monte Carlo algorithms have a quadratic time complexity, O(L(2)), with respect to the contour length (L) of the molecules. Here, we present an efficient method to compute Tw exactly, leading in principle to algorithms with a linear complexity, which in practice is O(L(1.2)). Specifically, we use a discrete wormlike chain that includes the explicit double-helix structure of DNA and where the linking number is conserved by continuously preventing the generation of twist between any two consecutive cylinders of the discretized chain. As an application, we show that long (up to 21 kbp) linear molecules stretched by mechanical forces akin to magnetic tweezers contain, in the buckling regime, multiple and branched plectonemes that often coexist with curls and helices, and whose length and number are in good agreement with experiments. By attaching the ends of the molecules to a reservoir of twists with which these can exchange helix turns, we also show how to compute the torques in these models. As an example, we report values that are in good agreement with experiments and that concern the longest molecules that have been studied so far (16 kbp).

| S-EPMC4571024 | biostudies-literature

Addressing nonlinearities in Monte Carlo.

Project description:Monte Carlo is famous for accepting model extensions and model refinements up to infinite dimension. However, this powerful incremental design is based on a premise which has severely limited its application so far: a state-variable can only be recursively defined as a function of underlying state-variables if this function is linear. Here we show that this premise can be alleviated by projecting nonlinearities onto a polynomial basis and increasing the configuration space dimension. Considering phytoplankton growth in light-limited environments, radiative transfer in planetary atmospheres, electromagnetic scattering by particles, and concentrated solar power plant production, we prove the real-world usability of this advance in four test cases which were previously regarded as impracticable using Monte Carlo approaches. We also illustrate an outstanding feature of our method when applied to acute problems with interacting particles: handling rare events is now straightforward. Overall, our extension preserves the features that made the method popular: addressing nonlinearities does not compromise on model refinement or system complexity, and convergence rates remain independent of dimension.

| S-EPMC6125312 | biostudies-literature

Sequential Monte Carlo multiple testing.

Project description:MotivationIn molecular biology, as in many other scientific fields, the scale of analyses is ever increasing. Often, complex Monte Carlo simulation is required, sometimes within a large-scale multiple testing setting. The resulting computational costs may be prohibitively high.ResultsWe here present MCFDR, a simple, novel algorithm for false discovery rate (FDR) modulated sequential Monte Carlo (MC) multiple hypothesis testing. The algorithm iterates between adding MC samples across tests and calculating intermediate FDR values for the collection of tests. MC sampling is stopped either by sequential MC or based on a threshold on FDR. An essential property of the algorithm is that it limits the total number of MC samples whatever the number of true null hypotheses. We show on both real and simulated data that the proposed algorithm provides large gains in computational efficiency.AvailabilityMCFDR is implemented in the Genomic HyperBrowser (http://hyperbrowser.uio.no/mcfdr), a web-based system for genome analysis. All input data and results are available and can be reproduced through a Galaxy Pages document at: http://hyperbrowser.uio.no/mcfdr/u/sandve/p/mcfdr.Contactgeirksa@ifi.uio.no.

| S-EPMC3223366 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

Tweets

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data