No Arabic abstract
We show how coarse-grained modelling combined with umbrella sampling using distance-based order parameters can be applied to compute the free-energy landscapes associated with mechanical deformations of large DNA nanostructures. We illustrate this approach for the strong bending of DNA nanotubes and the potentially bistable landscape of twisted DNA origami sheets. The homogeneous bending of the DNA nanotubes is well described by the worm-like chain model; for more extreme bending the nanotubes reversibly buckle with the bending deformations localized at one or two kinks. For a twisted one-layer DNA origami, the twist is coupled to the bending of the sheet giving rise to a free-energy landscape that has two nearly-degenerate minima that have opposite curvatures. By contrast, for a two-layer origami, the increased stiffness with respect to bending leads to a landscape with a single free-energy minimum that has a saddle-like geometry. The ability to compute such landscapes is likely to be particularly useful for DNA mechanotechnology and for understanding stress accumulation during the self-assembly of origamis into higher-order structures.
The fractal properties of the total potential energy V as a function of time t are studied for a number of systems, including realistic models of proteins (PPT, BPTI and myoglobin). The fractal dimension of V(t), characterized by the exponent gamma, is almost independent of temperature and increases with time, more slowly the larger the protein. Perhaps the most striking observation of this study is the apparent universality of the fractal dimension, which depends only weakly on the type of molecular system. We explain this behavior by assuming that fractality is caused by a self-generated dynamical noise, a consequence of intermode coupling due to anharmonicity. Global topological features of the potential energy landscape are found to have little effect on the observed fractal behavior.
In this paper we first analyzed the inductive bias underlying the data scattered across complex free energy landscapes (FEL), and exploited it to train deep neural networks which yield reduced and clustered representation for the FEL. Our parametric method, called Information Distilling of Metastability (IDM), is end-to-end differentiable thus scalable to ultra-large dataset. IDM is also a clustering algorithm and is able to cluster the samples in the meantime of reducing the dimensions. Besides, as an unsupervised learning method, IDM differs from many existing dimensionality reduction and clustering methods in that it neither requires a cherry-picked distance metric nor the ground-true number of clusters, and that it can be used to unroll and zoom-in the hierarchical FEL with respect to different timescales. Through multiple experiments, we show that IDM can achieve physically meaningful representations which partition the FEL into well-defined metastable states hence are amenable for downstream tasks such as mechanism analysis and kinetic modeling.
Cooperative events requiring anomalously large fluctuations are a defining characteristic for the onset of glassy relaxation across many materials. The importance of such intermittent events has been noted in systems as diverse as superconductors, metallic glasses, gels, colloids, and granular piles. Here, we show that prohibiting the attainment of new record-high energy fluctuations -- by explicitly imposing a ``lid on the fluctuation spectrum -- impedes further relaxation in the glassy phase. This lid allows us to directly measure the impact of record events on the evolving system in extensive simulations of aging in such vastly distinct glass formers as spin glasses and tapped granular piles. Interpreting our results in terms of a dynamics of records succeeds in explaining the ubiquity of both, the logarithmic decay of the energy and the memory effects encoded in the scaling of two-time correlation functions of aging systems.
Enhanced sampling techniques have become an essential tool in computational chemistry and physics, where they are applied to sample activated processes that occur on a time scale that is inaccessible to conventional simulations. Despite their popularity, it is well known that they have constraints that hinder their applications to complex problems. The core issue lies in the need to describe the system using a small number of collective variables (CVs). Any slow degree of freedom that is not properly described by the chosen CVs will hinder sampling efficiency. However, exploration of configuration space is also hampered by including variables that are not relevant to describe the activated process under study. This paper presents the Adaptive Topography of Landscape for Accelerated Sampling (ATLAS), a new biasing method capable of working with many CVs. The root idea of ATLAS is to apply a divide-and-conquer strategy where the high-dimensional CVs space is divided into basins, each of which is described by an automatically-determined, low-dimensional set of variables. A well-tempered metadynamics-like bias is constructed as a function of these local variables. Indicator functions associated with the basins switch on and off the local biases, so that the sampling is performed on a collection of low-dimensional CV spaces, that are smoothly combined to generate an effectively high-dimensional bias. The unbiased Boltzmann distribution is recovered through reweigting, making the evaluation of conformational and thermodynamic properties straightforward. The decomposition of the free-energy landscape in local basins can be updated iteratively as the simulation discovers new (meta)stable states.
A key objective in DNA-based material science is understanding and precisely controlling the mechanical properties of DNA hydrogels. We perform microrheology measurements using diffusing-wave spectroscopy (DWS) to investigate the viscoelastic behavior of a hydrogel made of Y-shaped DNA nano-stars over a wide range of frequencies and temperatures. Results show a clear liquid-to-equilibrium-gel transition as the temperature cycles up and down across the melting-temperature region for which the Y-DNA bind to each other. These first measurements reveal the crossover of the elastic G({omega}) and loss modulus G({omega}) when the DNA-hydrogel formed at low temperatures is heated to a fluid phase of DNA nano-stars well above the melt temperature Tm. We show that the crossover relates to the life-time of the DNA-bond and also that percolation coincides with the systems Tm. The approach demonstrated here can be easily extended to more complicated DNA hydrogel systems and provides guidance for the future design of such transient, semi-flexible networks that can be adapted to the application of molecular sensing and controlled release.