Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userThe looping probability of random heteropolymers helps to understand the scaling properties of biopolymers
60
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Random heteropolymers are a minimal description of biopolymers and can provide a theoretical framework to the investigate the formation of loops in biophysical experiments. A two--state model provides a consistent and robust way to study the scaling properties of loop formation in polymers of the size of typical biological systems. Combining it with self--adjusting simulated--tempering simulations, we can calculate numerically the looping properties of several realizations of the random interactions within the chain. Differently from homopolymers, random heteropolymers display at different temperatures a continuous set of scaling exponents. The necessity of using self--averaging quantities makes finite--size effects dominant at low temperatures even for long polymers, shadowing the length--independent character of looping probability expected in analogy with homopolymeric globules. This could provide a simple explanation for the small scaling exponents found in experiments, for example in chromosome folding.
rate research
Read More
Polarons, introduced by Davydov to explain energy transport in $alpha$-helices, correspond to electrons localised on a few lattice sites because of their interaction with phonons. While the static polaron field configurations have been extensively studied, their displacement is more difficult to explain. In this paper we show that, when the next to nearest neighbour interactions are included, for physical values of the parameters, polarons can spontaneously move, at T=0, on bent chains that exhibit a positive gradient in their curvature. At room temperature polarons perform a random walk but a curvature gradient can induce a non-zero average speed similar to the one observed at zero temperature. We also show that at zero temperature a polaron bounces on sharply kinked junctions. We interpret these results in light of the energy transport by transmembrane proteins.
We propose a model for the formation of chromatin loops based on the diffusive sliding of a DNA-bound factor which can dimerise to form a molecular slip-link. Our slip-links mimic the behaviour of cohesin-like molecules, which, along with the CTCF protein, stabilize loops which organize the genome. By combining 3D Brownian dynamics simulations and 1D exactly solvable non-equilibrium models, we show that diffusive sliding is sufficient to account for the strong bias in favour of convergent CTCF-mediated chromosome loops observed experimentally. Importantly, our model does not require any underlying, and energetically costly, motor activity of cohesin. We also find that the diffusive motion of multiple slip-links along chromatin may be rectified by an intriguing ratchet effect that arises if slip-links bind to the chromatin at a preferred loading site. This emergent collective behaviour is driven by a 1D osmotic pressure which is set up near the loading point, and favours the extrusion of loops which are much larger than the ones formed by single slip-links.
In single molecule laser optical tweezer (LOT) pulling experiments a protein or RNA is juxtaposed between DNA handles that are attached to beads in optical traps. The LOT generates folding trajectories under force in terms of time-dependent changes in the distance between the beads. How to construct the full intrinsic folding landscape (without the handles and the beads) from the measured time series is a major unsolved problem. By using rigorous theoretical methods---which account for fluctuations of the DNA handles, rotation of the optical beads, variations in applied tension due to finite trap stiffness, as well as environmental noise and the limited bandwidth of the apparatus---we provide a tractable method to derive intrinsic free energy profiles. We validate the method by showing that the exactly calculable intrinsic free energy profile for a Generalized Rouse Model, which mimics the two-state behavior in nucleic acid hairpins, can be accurately extracted from simulated time series in a LOT setup regardless of the stiffness of the handles. We next apply the approach to trajectories from coarse grained LOT molecular simulations of a coiled-coil protein based on the GCN4 leucine zipper, and obtain a free energy landscape that is in quantitative agreement with simulations performed without the beads and handles. Finally, we extract the intrinsic free energy landscape from experimental LOT measurements for the leucine zipper, which is independent of the trap parameters.
We develop a theoretical approach to the protein folding problem based on out-of-equilibrium stochastic dynamics. Within this framework, the computational difficulties related to the existence of large time scale gaps in the protein folding problem are removed and simulating the entire reaction in atomistic details using existing computers becomes feasible. In addition, this formalism provides a natural framework to investigate the relationships between thermodynamical and kinetic aspects of the folding. For example, it is possible to show that, in order to have a large probability to remain unchanged under Langevin diffusion, the native state has to be characterized by a small conformational entropy. We discuss how to determine the most probable folding pathway, to identify configurations representative of the transition state and to compute the most probable transition time. We perform an illustrative application of these ideas, studying the conformational evolution of alanine di-peptide, within an all-atom model based on the empiric GROMOS96 force field.
Motivation: DNA data is transcribed into single-stranded RNA, which folds into specific molecular structures. In this paper we pose the question to what extent sequence- and structure-information correlate. We view this correlation as structural semantics of sequence data that allows for a different interpretation than conventional sequence alignment. Structural semantics could enable us to identify more general embedded patterns in DNA and RNA sequences. Results: We compute the partition function of sequences with respect to a fixed structure and connect this computation to the mutual information of a sequence-structure pair for RNA secondary structures. We present a Boltzmann sampler and obtain the a priori probability of specific sequence patterns. We present a detailed analysis for the three PDB-structures, 2JXV (hairpin), 2N3R (3-branch multi-loop) and 1EHZ (tRNA). We localize specific sequence patterns, contrast the energy spectrum of the Boltzmann sampled sequences versus those sequences that refold into the same structure and derive a criterion to identify native structures. We illustrate that there are multiple sequences in the partition function of a fixed structure, each having nearly the same mutual information, that are nevertheless poorly aligned. This indicates the possibility of the existence of relevant patterns embedded in the sequences that are not discoverable using alignments.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
