No Arabic abstract
In this paper, we study a two-lane totally asymmetric simple exclusion process (TASEP) coupled with random attachment and detachment of particles (Langmuir kinetics) in both lanes under open boundary conditions. Our model can describe the directed motion of molecular motors, attachment and detachment of motors, and free inter-lane transition of motors between filaments. In this paper, we focus on some finite-size effects of the system because normally the sizes of most real systems are finite and small (e.g., size $leq 10,000$). A special finite-size effect of the two-lane system has been observed, which is that the density wall moves left first and then move towards the right with the increase of the lane-changing rate. We called it the jumping effect. We find that increasing attachment and detachment rates will weaken the jumping effect. We also confirmed that when the size of the two-lane system is large enough, the jumping effect disappears, and the two-lane system has a similar density profile to a single-lane TASEP coupled with Langmuir kinetics. Increasing lane-changing rates has little effect on density and current after the density reaches maximum. Also, lane-changing rate has no effect on density profiles of a two-lane TASEP coupled with Langmuir kinetics at a large attachment/detachment rate and/or a large system size. Mean-field approximation is presented and it agrees with our Monte Carlo simulations.
Protein contacts contain important information for protein structure and functional study, but contact prediction from sequence information remains very challenging. Recently evolutionary coupling (EC) analysis, which predicts contacts by detecting co-evolved residues (or columns) in a multiple sequence alignment (MSA), has made good progress due to better statistical assessment techniques and high-throughput sequencing. Existing EC analysis methods predict only a single contact map for a given protein, which may have low accuracy especially when the protein under prediction does not have a large number of sequence homologs. Analogous to ab initio folding that usually predicts a few possible 3D models for a given protein sequence, this paper presents a novel structure learning method that can predict a set of diverse contact maps for a given protein sequence, in which the best solution usually has much better accuracy than the first one. Our experimental tests show that for many test proteins, the best out of 5 solutions generated by our method has accuracy at least 0.1 better than the first one when the top L/5 or L/10 (L is the sequence length) predicted long-range contacts are evaluated, especially for protein families with a small number of sequence homologs. Our best solutions also have better quality than those generated by the two popular EC methods Evfold and PSICOV.
We consider the exclusion process on a ring with time-dependent defective bonds at which the hoping rate periodically switches between zero and one. This system models main roads in city traffics, intersecting with perpendicular streets. We explore basic properties of the system, in particular dependence of the vehicular flow on the parameters of signalization as well as the system size and the car density. We investigate various types of the spatial distribution of the vehicular density, and show existence of a shock profile. We also measure waiting time behind traffic lights, and examine its relationship with the traffic flow.
Non-coding RNA molecules fold into precise base pairing patterns to carry out critical roles in genetic regulation and protein synthesis. We show here that coupling systematic mutagenesis with high-throughput SHAPE chemical mapping enables accurate base pair inference of domains from ribosomal RNA, ribozymes, and riboswitches. For a six-RNA benchmark that challenged prior chemical/computational methods, this mutate-and-map strategy gives secondary structures in agreement with crystallographic data (2 % error rates), including a blind test on a double-glycine riboswitch. Through modeling of partially ordered RNA states, the method enables the first test of an interdomain helix-swap hypothesis for ligand-binding cooperativity in a glycine riboswitch. Finally, the mutate-and-map data report on tertiary contacts within non-coding RNAs; coupled with the Rosetta/FARFAR algorithm, these data give nucleotide-resolution three-dimensional models (5.7 {AA} helix RMSD) of an adenine riboswitch. These results highlight the promise of a two-dimensional chemical strategy for inferring the secondary and tertiary structures that underlie non-coding RNA behavior.
We investigate the dynamics of a one-dimensional asymmetric exclusion process with Langmuir kinetics and a fluctuating wall. At the left boundary, particles are injected onto the lattice; from there, the particles hop to the right. Along the lattice, particles can adsorb or desorb, and the right boundary is defined by a wall particle. The confining wall particle has intrinsic forward and backward hopping, a net leftward drift, and cannot desorb. Performing Monte Carlo simulations and using a moving-frame finite segment approach coupled to mean field theory, we find the parameter regimes in which the wall acquires a steady state position. In other regimes, the wall will either drift to the left and fall off the lattice at the injection site, or drift indefinitely to the right. Our results are discussed in the context of non-equilibrium phases of the system, fluctuating boundary layers, and particle densities in the lab frame versus the frame of the fluctuating wall.
In this paper, we investigate the dynamics of synchronous totally asymmetric exclusion processes on lattices with a multiple-input single-output (MISO) junction, which consists of m subchains for the input and one main chain for the output. A MISO junction is a type of complex geometry that is relevant to many biological processes as well as vehicular and pedestrian traffic flow. A mean-field approach is developed to deal with the junction that connects the subchains and the main chain. Theoretical results for stationary particle currents, density profiles, and a phase diagram have been obtained. It is found that the phase boundary moves toward the left in the phase diagram with an increase of the number of subchains. The nonequilibrium stationary states, stationary-state phases, and phase boundaries are determined by the boundary conditions of the system as well as by the number of subchains. The density profiles obtained from computer simulations show very good agreement with our theoretical analysis.