Kinetics and thermodynamics of first-order Markov chain copolymerization

We report a theoretical study of stochastic processes modeling the growth of first-order Markov copolymers, as well as the reversed reaction of depolymerization. These processes are ruled by kinetic equations describing both the attachment and detachment of monomers. Exact solutions are obtained for these kinetic equations in the steady regimes of multicomponent copolymerization and depolymerization. Thermodynamic equilibrium is identified as the state at which the growth velocity is vanishing on average and where detailed balance is satisfied. Away from equilibrium, the analytical expression of the thermodynamic entropy production is deduced in terms of the Shannon disorder per monomer in the copolymer sequence. The Mayo-Lewis equation is recovered in the fully irreversible growth regime. The theory also applies to Bernoullian chains in the case where the attachment and detachment rates only depend on the reacting monomer

Noncooperative algorithms in self-assembly

We show the first non-trivial positive algorithmic results (i.e. programs whose output is larger than their size), in a model of self-assembly that has so far resisted many attempts of formal analysis or programming: the planar non-cooperative variant of Winfree's abstract Tile Assembly Model. This model has been the center of several open problems and conjectures in the last fifteen years, and the first fully general results on its computational power were only proven recently (SODA 2014). These results, as well as ours, exemplify the intricate connections between computation and geometry that can occur in self-assembly. In this model, tiles can stick to an existing assembly as soon as one of their sides matches the existing assembly. This feature contrasts with the general cooperative model, where it can be required that tiles match on \emph{several} of their sides in order to bind. In order to describe our algorithms, we also introduce a generalization of regular expressions called Baggins expressions. Finally, we compare this model to other automata-theoretic models.Comment: A few bug fixes and typo correction

Theoretical description of a DNA-linked nanoparticle self-assembly

Nanoparticles tethered with DNA strands are promising building blocks for bottom-up nanotechnology, and a theoretical understanding is important for future development. Here we build on approaches developed in polymer physics to provide theoretical descriptions for the equilibrium clustering and dynamics, as well as the self-assembly kinetics of DNA-linked nanoparticles. Striking agreement is observed between the theory and molecular modeling of DNA tethered nanoparticles.Comment: Accepted for publication in Physical Review Letter

Quantum Quenches in a Holographic Kondo Model

We study non-equilibrium dynamics and quantum quenches in a recent gauge/gravity duality model for a strongly coupled system interacting with a magnetic impurity with $SU(N)$ spin. At large $N$, it is convenient to write the impurity spin as a bilinear in Abrikosov fermions. The model describes an RG flow triggered by the marginally relevant Kondo operator. There is a phase transition at a critical temperature, below which an operator condenses which involves both an electron and an Abrikosov fermion field. This corresponds to a holographic superconductor in AdS$_2$ and models the impurity screening. We study the time dependence of the condensate induced by quenches of the Kondo coupling. The timescale for equilibration is generically given by the lowest-lying quasinormal mode of the dual gravity model. This mode also governs the formation of the screening cloud, which is obtained as the decrease of impurity degrees of freedom with time. In the condensed phase, the leading quasinormal mode is imaginary and the relaxation of the condensate is over-damped. For quenches whose final state is close to the critical point of the large $N$ phase transition, we study the critical slowing down and obtain the combination of critical exponents $z\nu=1$. When the final state is exactly at the phase transition, we find that the exponential ringing of the quasinormal modes is replaced by a power-law behaviour of the form $\sim t^{-a}\sin(b\log t)$. This indicates the emergence of a discrete scale invariance.Comment: 23 pages + appendices, 11 figure

Effects of Kinks on DNA Elasticity

We study the elastic response of a worm-like polymer chain with reversible kink-like structural defects. This is a generic model for (a) the double-stranded DNA with sharp bends induced by binding of certain proteins, and (b) effects of trans-gauche rotations in the backbone of the single-stranded DNA. The problem is solved both analytically and numerically by generalizing the well-known analogy to the Quantum Rotator. In the small stretching force regime, we find that the persistence length is renormalized due to the presence of the kinks. In the opposite regime, the response to the strong stretching is determined solely by the bare persistence length with exponential corrections due to the ``ideal gas of kinks''. This high-force behavior changes significantly in the limit of high bending rigidity of the chain. In that case, the leading corrections to the mechanical response are likely to be due to the formation of multi-kink structures, such as kink pairs.Comment: v1: 16 pages, 7 figures, LaTeX; submitted to Physical Review E; v2: a new subsection on soft kinks added to section Theory, sections Introduction and Conclusions expanded, references added, other minor changes; v3: a reference adde

Force-induced misfolding in RNA

RNA folding is a kinetic process governed by the competition of a large number of structures stabilized by the transient formation of base pairs that may induce complex folding pathways and the formation of misfolded structures. Despite of its importance in modern biophysics, the current understanding of RNA folding kinetics is limited by the complex interplay between the weak base-pair interactions that stabilize the native structure and the disordering effect of thermal forces. The possibility of mechanically pulling individual molecules offers a new perspective to understand the folding of nucleic acids. Here we investigate the folding and misfolding mechanism in RNA secondary structures pulled by mechanical forces. We introduce a model based on the identification of the minimal set of structures that reproduce the patterns of force-extension curves obtained in single molecule experiments. The model requires only two fitting parameters: the attempt frequency at the level of individual base pairs and a parameter associated to a free energy correction that accounts for the configurational entropy of an exponentially large number of neglected secondary structures. We apply the model to interpret results recently obtained in pulling experiments in the three-helix junction S15 RNA molecule (RNAS15). We show that RNAS15 undergoes force-induced misfolding where force favors the formation of a stable non-native hairpin. The model reproduces the pattern of unfolding and refolding force-extension curves, the distribution of breakage forces and the misfolding probability obtained in the experiments.Comment: 28 pages, 11 figure

Universal Formulae for Percolation Thresholds

A power law is postulated for both site and bond percolation thresholds. The formula writes $p_c=p_0[(d-1)(q-1)]^{-a}d^{\ b}$, where $d$ is the space dimension and $q$ the coordination number. All thresholds up to $d\rightarrow \infty$ are found to belong to only three universality classes. For first two classes $b=0$ for site dilution while $b=a$ for bond dilution. The last one associated to high dimensions is characterized by $b=2a-1$ for both sites and bonds. Classes are defined by a set of value for $\{p_0; \ a\}$. Deviations from available numerical estimates at $d \leq 7$ are within $\pm 0.008$ and $\pm 0.0004$ for high dimensional hypercubic expansions at $d \geq 8$. The formula is found to be also valid for Ising critical temperatures.Comment: 11 pages, latex, 3 figures not include