Source author record

Petr Šulc

Petr Šulc appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2020arXiv

Coarse-Grained Nucleic Acid-Protein Model for Hybrid Nanotechnology

The emerging field of hybrid DNA - protein nanotechnology brings with it the potential for many novel materials which combine the addressability of DNA nanotechnology with versatility of protein interactions. However, the design and computational study of these hybrid structures is difficult due to the system sizes involved. To aid in the design and in silico analysis process, we introduce here a coarse-grained DNA/RNA-protein model that extends the oxDNA/oxRNA models of DNA/RNA with a coarse-grained model of proteins based on an anisotropic network model representation. Fully equipped with analysis scripts and visualization, our model aims to facilitate hybrid nanomaterial design towards eventual experimental realization, as well as enabling study of biological complexes. We further demonstrate its usage by simulating DNA-protein nanocage, DNA wrapped around histones, and a nascent RNA in polymerase.

preprint2020arXiv

Designing patchy interactions to self-assemble arbitrary structures

One of the fundamental goals of nanotechnology is to exploit selective and directional interactions between molecules to design particles that self-assemble into desired structures, from capsids, to nano-clusters, to fully formed crystals with target properties (e.g. optical, mechanical, etc.). Here we provide a general framework which transforms the inverse problem of self-assembly of colloidal crystals into a Boolean satisfiability problem for which solutions can be found numerically. Given a reference structure and the desired number of components, our approach produces designs for which the target structure is an energy minimum, and also allows to exclude solutions that correspond to competing structures. We demonstrate the effectiveness of our approach by designing model particles that spontaneously nucleate milestone structures such as the cubic diamond, the pyrochlore and the clathrate lattices.

preprint2015arXiv

Coarse-grained modelling of supercoiled RNA

We study the behaviour of double-stranded RNA under twist and tension using oxRNA, a recently developed coarse-grained model of RNA. Introducing explicit salt-dependence into the model allows us to directly compare our results to data from recent single-molecule experiments. The model reproduces extension curves as a function of twist and stretching force, including the buckling transition and the behaviour of plectoneme structures. For negative supercoiling, we predict denaturation bubble formation in plectoneme end-loops, suggesting preferential plectoneme localisation in weak base sequences. OxRNA exhibits a positive twist-stretch coupling constant, in agreement with recent experimental observations.

preprint2014arXiv

A comparison between parallelization approaches in molecular dynamics simulations on GPUs

We test the performances of two different approaches to the computation of forces for molecular dynamics simulations on Graphics Processing Units. A "vertex-based" approach, where a computing thread is started per particle, is compared to a newly proposed "edge-based" approach, where a thread is started per each potentially non-zero interaction. We find that the former is more efficient for systems with many simple interactions per particle, while the latter is more efficient if the system has more complicated interactions or fewer of them. By comparing computation times on more and less recent GPU technology, we predict that, if the current trend of increasing the number of processing cores - as opposed to their computing power - remains, the "edge-based" approach will gradually become the most efficient choice in an increasing number of cases.

preprint2014arXiv

A nucleotide-level coarse-grained model of RNA

We present a new, nucleotide-level model for RNA, oxRNA, based on the coarse-graining methodology recently developed for the oxDNA model of DNA. The model is designed to reproduce structural, mechanical and thermodynamic properties of RNA, and the coarse-graining level aims to retain the relevant physics for RNA hybridization and the structure of single- and double-stranded RNA. In order to explore its strengths and weaknesses, we test the model in a range of nanotechnological and biological settings. Applications explored include the folding thermodynamics of a pseudoknot, the formation of a kissing loop complex, the structure of a hexagonal RNA nanoring, and the unzipping of a hairpin motif. We argue that the model can be used for efficient simulations of the structure of systems with thousands of base pairs, and for the assembly of systems of up to hundreds of base pairs. The source code implementing the model is released for public use.

preprint2014arXiv

Modelling toehold-mediated RNA strand displacement

We study the thermodynamics and kinetics of an RNA toehold-mediated strand displacement reaction with a recently developed coarse-grained model of RNA. Strand displacement, during which a single strand displaces a different strand previously bound to a complementary substrate strand, is an essential mechanism in active nucleic acid nanotechnology and has also been hypothesized to occur in vivo. We study the rate of displacement reactions as a function of the length of the toehold and temperature and make two experimentally testable predictions: that the displacement is faster if the toehold is placed at the 5' end of the substrate and that the displacement slows down with increasing temperature for longer toeholds.

preprint2014arXiv

Optimal Distributed Control of Reactive Power via the Alternating Direction Method of Multipliers

We formulate the control of reactive power generation by photovoltaic inverters in a power distribution circuit as a constrained optimization that aims to minimize reactive power losses subject to finite inverter capacity and upper and lower voltage limits at all nodes in the circuit. When voltage variations along the circuit are small and losses of both real and reactive powers are small compared to the respective flows, the resulting optimization problem is convex. Moreover, the cost function is separable enabling a distributed, on-line implementation with node-local computations using only local measurements augmented with limited information from the neighboring nodes communicated over cyber channels. Such an approach lies between the fully centralized and local policy approaches previously considered. We explore protocols based on the dual ascent method and on the Alternating Direction Method of Multipliers (ADMM) and find that the ADMM protocol performs significantly better.

preprint2013arXiv

DNA hybridization kinetics: zippering, internal displacement and sequence dependence

While the thermodynamics of DNA hybridization is well understood, much less is known about the kinetics of this classic system. Filling this gap in our understanding has new urgency because DNA nanotechnology often depends critically on binding rates. Here we use a coarse-grained model to explore the hybridization kinetics of DNA oligomers, finding that strand association proceeds through a complex set of intermediate states. Successful binding events start with the formation of a few metastable base-pairing interactions, followed by zippering of the remaining bonds. However, despite reasonably strong interstrand interactions, initial contacts frequently fail to lead to zippering because the typical configurations in which they form differ from typical states of similar enthalpy in the double-stranded equilibrium ensemble. Therefore, if the association process is analyzed on the base-pair (secondary structure) level, it shows non-Markovian behavior. Initial contacts must be stabilized by two or three base pairs before full zippering is likely, resulting in negative effective activation enthalpies. Non-Arrhenius behavior is observed as the number of base pairs in the effective transition state increases with temperature. In addition, we find that alternative pathways involving misbonds can increase association rates. For repetitive sequences, misaligned duplexes frequently rearrange to form fully paired duplexes by two distinct processes which we label `pseudoknot' and `inchworm' internal displacement. We show how the above processes can explain why experimentally observed association rates of GC-rich oligomers are higher than rates of AT-rich equivalents. More generally, we argue that association rates can be modulated by sequence choice.

preprint2012arXiv

Sequence-dependent thermodynamics of a coarse-grained DNA model

We introduce a sequence-dependent parametrization for a coarse-grained DNA model [T. E. Ouldridge, A. A. Louis, and J. P. K. Doye, J. Chem. Phys. 134, 085101 (2011)] originally designed to reproduce the properties of DNA molecules with average sequences. The new parametrization introduces sequence-dependent stacking and base-pairing interaction strengths chosen to reproduce the melting temperatures of short duplexes. By developing a histogram reweighting technique, we are able to fit our parameters to the melting temperatures of thousands of sequences. To demonstrate the flexibility of the model, we study the effects of sequence on: (a) the heterogeneous stacking transition of single strands, (b) the tendency of a duplex to fray at its melting point, (c) the effects of stacking strength in the loop on the melting temperature of hairpins, (d) the force-extension properties of single strands and (e) the structure of a kissing-loop complex. Where possible we compare our results with experimental data and find a good agreement. A simulation code called oxDNA, implementing our model, is available as free software.

preprint2012arXiv

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

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.

preprint2011arXiv

Learning the nonlinear interactions from particle trajectories

Nonlinear interaction of membrane proteins with cytoskeleton and membrane leads to non-Gaussian structure of their displacement probability distribution. We propose a novel statistical analysis technique for learning the characteristics of the nonlinear potential from the cumulants of the displacement distribution. The efficiency of the approach is demonstrated on the analysis of kurtosis of the displacement distribution of the particle traveling on a membrane in a cage-type potential. Results of numerical simulations are supported by analytical predictions. We show that the approach allows robust identification of the potential for the much lower temporal resolution compare with the mean square displacement analysis.