BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Yunxin Zhang

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

ResearcherUnclaimed source record
Biological PhysicsSubcellular Processescond-mat.stat-mechmath-phmath.MPcond-mat.softphysics.chem-phBiomoleculesmath.APmath.DGcond-mat.mes-hallGenomicsmath.CAmath.DSmath.PRMolecular Networksphysics.ao-phphysics.med-ph

Catalog footprint

What is connected

39works
18topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

39 published item(s)

preprint2022arXiv

Dual-domain Attention-based Deep Network for Sparse-view CT Artifact Reduction

Due to the wide applications of X-ray computed tomography (CT) in medical imaging activities, radiation exposure has become a major concern for public health. Sparse-view CT is a promising approach to reduce the radiation dose by down-sampling the total number of acquired projections. However, the CT images reconstructed by this sparse-view imaging approach suffer from severe streaking artifacts and structural information loss. In this work, an end-to-end dual-domain attention-based deep network (DDANet) is proposed to solve such an ill-posed CT image reconstruction problem. The image-domain CT image and the projection-domain sinogram are put into the two parallel sub-networks of the DDANet to independently extract the distinct high-level feature maps. In addition, a specified attention module is introduced to fuse the aforementioned dual-domain feature maps to allow complementary optimizations of removing the streaking artifacts and mitigating the loss of structure. Numerical simulations, anthropomorphic thorax phantom and in vivo pre-clinical experiments are conducted to verify the sparse-view CT imaging performance of the DDANet. Results demonstrate that this newly developed approach is able to robustly remove the streaking artifacts while maintaining the fine structures. As a result, the DDANet provides a promising solution in achieving high quality sparse-view CT imaging.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Existence and uniqueness of solution of the differential equation describing the TASEP-LK coupled transport process

In this paper, the existence and uniqueness of solution of a specific differential equation is studied. This equation originates from the description of a coupled process by totally asymmetric simple exclusion process (TASEP) and Langmuir kinetics (LK). In the fields of physics and biology, the properties of the TASEP-LK coupled process have been extensively studied by Monte Carlo simulations and numerical calculations, as well as detailed experiments. However, so far, no rigorous mathematical analysis has been given to the corresponding differential equations, especially their existence and uniqueness of solution. In this paper, using the upper and lower solution method, the existence of solution of the steady state equation is obtained. Then using a generalized maximum principle, we show that the solution constructed from the upper and lower solution method is actually the unique solution in C∞ space. Moreover, the existence and uniqueness of solution of the time dependent differential equation are also obtained in one specific space X\b{eta}. Our results imply that the previous results obtained by numerical calculations and Monte Carlo simulations are theoretically correct, especially the most important phase diagram of particle density along the travel track under different model parameters. The study in this paper provides theoretical foundations for the analysis of TASEP-LK coupled process. At the same time, the methods used in this paper may be instructive for studies about the more general cases of the TASEP-LK process, such as the one with multiple travel tracks or the one with multiple particle species.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Existence and uniqueness of solution of the differential equation describing the TASEP-LK coupled transport process

We study the existence and uniqueness of solution of a evolutionary partial differential equation originating from the continuum limit of a coupled process of totally asymmetric simple exclusion process (TASEP) and Langmuir kinetics (LK). In the fields of physics and biology, the TASEP-LK coupled process has been extensively studied by Monte Carlo simulations, numerical computations, and detailed experiments. However, no rigorous mathematical analysis so far has been given for the corresponding differential equations, especially the existence and uniqueness of their solutions. In this paper, we prove the existence of the $C^\infty[0,1]$ steady-state solution by the method of upper and lower solution, and the uniqueness in both $W^{1,2}(0,1)$ and $L^\infty(0,1)$ by a generalized maximum principle. We further prove the global existence and uniqueness of the time-dependent solution in $C([0,1]\times [0,+\infty))\cap C^{2,1}([0,1]\times (0,+\infty))$, which, for any continuous initial value, converges to the steady-state solution in $C[0,1]$ (global attractivity). Our results support the numerical calculations and Monte Carlo simulations, and provide theoretical foundations for the TASEP-LK coupled process, especially the most important phase diagram of particle density along the travel track under different model parameters, which is difficult because the boundary layers (at one or both boundaries) and domain wall (separating high and low particle densities) may appear as the length of the travel track tends to infinity. The methods used in this paper may be instructive for studies of the more general cases of the TASEP-LK process, such as the one with multiple travel tracks and/or multiple particle species.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Optimization of a class of heat engines with explicit solution

A specific class of stochastic heat engines driven cyclically by time-dependent potential, which is defined in the half-line ($0<x<+\infty$), is analysed. For such engines, most of their physical quantities can be obtained explicitly, including the entropy and internal energy of the heat engine, as well as output work (power) and heat exchange with the environment during a finite time interval. The optimisation method based on the external potential to reduce {\it irreversible} work and increase energy efficiency is presented. With this optimised potential, efficiency $η^*$ and its particular value at maximum power $η^*_{\rm EMP}$ are calculated and discussed briefly.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Optimization of stochastic thermodynamic machines

The study of stochastic thermodynamic machines is one of the main topics in nonequilibrium thermodynamics. In this study, within the framework of Fokker-Planck equation, and using the method of characteristics of partial differential equation as well as the variational method, performance of stochastic thermodynamic machines is optimized according to the external potential, with the irreversible work $W_{irr}$, or the total entropy production $ΔS_{\rm tot}$ equivalently, reaching its lower bound. Properties of the optimal thermodynamic machines are discussed, with explicit expressions of upper bounds of work output $W$, power $P$, and energy efficiency $η$ are presented. To illustrate the results obtained, typical examples with optimal protocols (external potentials) are also presented.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Work needed to drive a thermodynamic system between two distributions

In this study, the minimum amount of work needed to drive a thermodynamic system from one initial distribution to another in a given time duration is discussed. Equivalently, for given amount of work, the minimum time duration required to complete such a transition is obtained. Results show that the minimum amount of work is used to achieve the following three objectives, to increase the internal energy of the system, to decrease the system entropy, to change the mean position of the system, and with other nonzero part dissipated into environment. To illustrate the results, an example with explicit solutions is presented.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Mean velocity and effective diffusion constant for translocation of biopolymer chains across membrane

Chaperone-assisted translocation through a nanopore embedded in membrane holds a prominent role in the transport of biopolymers. Inspired by classical Brownian ratchet, we develop a theoretical framework characterizing such translocation process through a master equation approach. In this framework, the polymer chain, provided with reversible binding of chaperones, undergoes forward/backward diffusion, which is rectified by chaperones. We drop the assumption of timescale separation and keep the length of a polymer chain finite, both of which happen to be the key points in most of the previous studies. Our framework makes it accessible to derive analytical expressions for mean translocation velocity and effective diffusion constant in stationary state, which is the basis of a comprehensive understanding towards the dynamics of such process. Generally, the translocation of polymer chain across membrane consists of three subprocesses: initiation, termination, and translocation of the main body part of a polymer chain, where the translocation of the main body part depends on the binding/unbinding kinetics of chaperones. That is the main concern of this study. Our results show that the increase of forward/backward diffusion rate of a polymer chain and the binding/unbinding ratio of chaperones both raise the mean translocation velocity of a polymer chain, and roughly speaking, the dependence of effective diffusion constant on these two factors achieves similar behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Theoretical model of transcription based on torsional mechanics of DNA template

Transcription is the first step of gene expression, in which a particular segment of DNA is copied to RNA by the enzyme RNA polymerase (RNAP). Despite many details of the complex interactions between DNA and RNA synthesis disclosed experimentally, much of physical behavior of transcription remains largely unknown. Understanding torsional mechanics of DNA and RNAP together with its nascent RNA and RNA-bound proteins in transcription maybe the first step towards deciphering the mechanism of gene expression. In this study, based on the balance between viscous drag on RNA synthesis and torque resulted from untranscribed supercoiled DNA template, a simple model is presented to describe mechanical properties of transcription. With this model, the rotation and supercoiling density of the untranscribed DNA template are discussed in detail. Two particular cases of transcription are considered, transcription with constant velocity and transcription with torque dependent velocity. Our results show that, during the initial stage of transcription, rotation originated from the transcribed part of DNA template is mainly released by the rotation of RNAP synthesis. During the intermediate stage, the rotation is usually released by both the supercoiling of the untranscribed part of DNA template and the rotation of RNAP synthesis, with proportion depending on the friction coefficient in environment and the length of nascent RNA. However, with the approaching to the upper limit of twisting of the untranscribed DNA template, the rotation resulted from transcription will then be mainly released by the rotation of RNAP synthesis.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

How does a protein reach its binding locus: sliding along DNA chain or not?

In gene expression, various kinds of proteins need to bind to specific locus of DNA. It is still not clear how these proteins find their target locus. In this study, the mean first-passage time (FPT) of protein binding to its target locus on DNA chain is discussed by a chain-space coupled model. Our results show that the 1-dimensional diffusion constant has a critical value, with which the mean time spent by a protein to find its target locus is almost independent of the binding rate of protein to DNA chain and the detachment rate from DNA chain. Which implies that, the frequency of protein binding to DNA and the sliding time on DNA chain have little influence on the search efficiency, and therefore whether or not the 1-dimensional sliding on DNA chain increases the search efficiency depends on the 1-dimensional diffusion constant of the protein on DNA chain. This study also finds that only protein bindings to DNA loci which are close to the target locus help to increase the search efficiency, while bindings to those loci which are far from the target locus might delay the target binding process. As expected, the mean FPT increases with the distance between the initial position of protein in cell space and its target locus on DNA chain. The direct binding probability, which can be regarded as one index to describe if the 1-dimensional sliding along DNA chain is helpful to increase the search efficiency is calculated. Our results show that the influence of 1-dimensional sliding along DNA chain on the search process depends on both diffusion constants of protein in cell space and on the 1-dimensional DNA chain.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Mechanism of helicase translocation along nucleic acid

In cells, helicase translocation along nucleic acid is essential for many biological processes. However, so far, the mechanism of this translocation is not fully understood. Recent studies show that helicase might translocate through two processes, active process and passive process, with different translocation rate. In this study, a model including such two processes is presented. In which, each of these two processes consists of two sub-processes, chemical sub-process in which needed translocation factors are attached, and mechanochemical sub-process in which helicase makes a forward translocation step. Helicase can switch stochastically between these two processes with external force dependent rates. By this model, ribosome translocation along message RNA is detailed discussed. We found that, with the increase of external force, the mean translocation rate of ribosome increases from one lower limit to one upper limit, and both of these two limits increase with concentrations of the translocation factors. Under high external force, ribosome translocates mainly through active process. At saturating concentrations of translocation factors, the translocation is governed by mechanochemical sub-processes.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A weaker rigidity theorem for pairs of hyperquadrics and its application

In this short article, we establish a rigidity theorem for pairs of hyperquadrics in a weaker sense, i.e., we impose a condition that minimal rational curves are preserved, which is stronger than inheriting a sub-VMRT structure, a notion raised by Mok & Zhang (2014). This problem has its source in a theorem of Tsai (1993), and the main result of this article can be applied back to give a more intrinsic proof of Tsai's theorem.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Comment on "Optimal Reaction Time for Surface-Mediated Diffusion"

In recent letter [Phys. Rev. Lett {\bf 105}, 150606 (2010)], the surface-mediated diffusion problem is theoretically discussed, and interesting results have been obtained. However, for more general cases, the ansatz of solutions of the diffusion equation, which is the starting point of their analysis, might not be appropriate. In this comment, suggested ansatz and corresponding methods will be presented.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Correlations between promoter activity and its nucleotide positions in spacing region

Transcription is one of the essential processes for cells to read genetic information encoded in genes, which is initiated by the binding of RNA polymerase to related promoter. Experiments have found that the nucleotide sequence of promoter has great influence on gene expression strength, or promoter activity. In synthetic biology, one interesting question is how we can synthesize a promoter with given activity, and which positions of promoter sequence are important for determining its activity. In this study, based on recent experimental data, correlations between promoter activity and its sequence positions are analyzed by various methods. Our results show that, except nucleotides in the two highly conserved regions, $-35$ box and $-10$ box, influences of nucleotides in other positions are also not neglectable. For example, modifications of nucleotides around position $-19$ in spacing region may change promoter activity in a large scale. The results of this study might be helpful to our understanding of biophysical mechanism of gene transcription, and may also be helpful to the design of synthetic cell factory.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Mean-field analysis of two-species TASEP with attachment and detachment

In cells, most of cargos are transported by motor proteins along microtubule. Biophysically, unidirectional motion of large number of motor proteins along a single track can be described by totally asymmetric simple exclusion process (TASEP). From which many meaningful properties, such as the appearance of domain wall (defined as the borderline of high density and low density of motor protein along motion track) and boundary layers, can be obtained. However, it is biologically obvious that a single track may be occupied by different motor species. So previous studies based on TASEP of one particle species are not reasonable enough to find more detailed properties of the motion of motors along a single track. To address this problem, TASEP with two particle species is discussed in this study. Theoretical methods to get densities of each particle species are provided. Using these methods, phase transition related properties of particle densities are obtained. Our analysis show that domain wall and boundary layer of single species densities always appear simultaneously with those of the total particle density. The height of domain wall of total particle density is equal to the summation of those of single species. Phase diagrams for typical model parameters are also presented. The methods presented in this study can be generalized to analyze TASEP with more particle species.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Theoretical analysis of transcription process with polymerase stalling

Experimental evidences show that in gene transcription, RNA polymerase has the possibility to be stalled at certain position of the transcription template. This may be due to the template damage, or protein barriers. Once stalled, polymerase may backtrack along the template to the previous nucleotide to wait for the repair of the damaged site, or simply bypass the barrier or damaged site and consequently synthesize an incorrect messenger RNA, or degrade and detach from the template. Thus, the {\it effective} transcription rate (the rate to synthesize correct product mRNA) and the transcription {\it effectiveness} (the ratio of the {\it effective} transcription rate to the {\it effective} transcription initiation rate) are both influenced by polymerase stalling events. This study shows that, Without backtracking, detachment of stalled polymerase can also help to increase the {\it effective} transcription rate and transcription {\it effectiveness}. Generally, the increase of bypass rate of the stalled polymerase will lead to the decrease of the {\it effective} transcription rate and transcription {\it effectiveness}. But when both detachment rate and backtracking rate of the stalled polymerase vanish, the {\it effective} transcription rate may also be increased by the bypass mechanism.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Translation with frameshifting of ribosome along mRNA transcript

Translation is an important process for prokaryotic and eukaryotic cells to produce necessary proteins for cell growth. Numerious experiments have been performed to explore the translational properties. Diverse models have also been developed to determine the biochemical mechanism of translation. However, to simplify the majority of the existing models, the frameshifting of ribosome along the mRNA transcript is neglected, which actually occurs in real cells and has been extensively experimentally studied. The frameshifting of ribosome evidently influences the efficiency and speed of translation, considering that the peptide chains synthesized by shifted ribosomes will not fold into functional proteins and will degrade rapidly. In this study, a theoretical model is presented to describe the translational process based on the model for totally asymmetric simple exclusion process. In this model, the frameshifting of the ribosome along the mRNA transcript and the attachment/detachment of the ribosome to/from the main body of mRNA codons during translation elongation process, are explicitly included. The results show that, with ribosome frameshifing, the speed of correctly synthesized peptide chains may increase first and then decrease with both the translation initiation rate $α$ and the ribosome detachment rate $ω_d$. This results indicates that regulating the translation process to reach maximal synthesized speed of proteins is theoretically feasible. Traffic-related problems of ribosome motion along the mRNA transcript are also addressed theoretically. Depending on parameter values, shock wave (or domain wall) may exist for ribosome probabilities along the mRNA.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Admissible pairs of Hermitian symmetric spaces in the perspective of the theory of varieties of minimal rational tangents

We study a pair (\mathcal{S}_0,\mathcal{S}) of irreducible Hermitian Symmetric Spaces of compact type (cHSS) in this paper, with the first aim being classifying all the admissible pairs (\mathcal{S}_0,\mathcal{S})). This notion is a natural generalization of the pairs of sub-diagram type originated by Jaehyun Hong and Ngaiming Mok ([HoM 10]). Based on this classification, we partially solve the rigidity problem for the admissible pairs (\mathcal{S}_0,\mathcal{S}) which was raised by Mok and Zhang (2014) ([MoZ 14]), culminating in determining a sufficient condition for the pairs being non-rigid and proving that special pairs, which show up in the classification procedure, are algebraic, as a weaker result than being rigid. However, whether special pairs are rigid or not remains unknown and needs further investigation in the framework of VMRT theory.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Relationship between promoter sequence and its strength in gene expression

In this study, through various tests one theoretical model is presented to describe the relationship between promoter strength and its nucleotide sequence. Our analysis shows that, promoter strength is greatly influenced by nucleotide groups with three adjacent nucleotides in its sequence. Meanwhile, nucleotides in different regions of promoter sequence have different effects on promoter strength. Based on experimental data for {\it E. coli} promoters, our calculations indicate, nucleotides in -10 region, -35 region, and the discriminator region of promoter sequence are more important than those in spacing region for determining promoter strength. With model parameter values obtained by fitting to experimental data, four promoter libraries are theoretically built for the corresponding experimental environments under which data for promoter strength in gene expression has been measured previously.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Improvement of risk estimate on wind turbine tower buckled by hurricane

Wind is one of the important reasonable resources. However, wind turbine towers are sure to be threatened by hurricanes. In this paper, method to estimate the number of wind turbine towers that would be buckled by hurricanes is discussed. Monte Carlo simulations show that our method is much better than the previous one. Since in our method, the probability density function of the buckling probability of a single turbine tower in a single hurricane is obtained accurately but not from one approximated expression. The result in this paper may be useful to the design and maintenance of wind farms.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Properties of sodium-driven bacterial flagellar motor: A two-state model approach

Bacterial flagellar motor (BFM) is one of the ion-driven molecular machines, which drives the rotation of flagellar filaments and enable bacteria to swim in viscous solutions. Understanding its mechanism is one challenge in biophysics. Based on previous models and inspired by the idea used in description of motor proteins, in this study one two-state model is provided. Meanwhile, according to corresponding experimental data, mathematical relationship between BFM membrane voltage and pH value of the environment, and relationship between internal and external sodium concentrations are given. Therefore, with model parameter values obtained by fitting theoretical results of torque-speed relation to recent experimental data, many biophysical properties of bacterial flagellar motor can be obtained for any pH values and any external sodium concentrations. Including the rotation speed, stall torque (i.e. the torque generated by BFM), rotation dispersion, and rotation randomness. In this study, the single-stator BFM will be firstly analyzed, and then properties of multiple-stator BFM are addressed briefly.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

ATP binding to a multisubunit enzyme: statistical thermodynamics analysis

Due to inter-subunit communication, multisubunit enzymes usually hydrolyze ATP in a concerted fashion. However, so far the principle of this process remains poorly understood. In this study, from the viewpoint of statistical thermodynamics, a simple model is presented. In this model, we assume that the binding of ATP will change the potential of the corresponding enzyme subunit, and the degree of this change depends on the state of its adjacent subunits. The probability of enzyme in a given state satisfies the Boltzmann's distribution. Although it looks much simple, this model can fit the recent experimental data of chaperonin TRiC/CCT well. From this model, the dominant state of TRiC/CCT can be obtained. This study provided a new way to understand biophysical processes by statistical thermodynamics analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Cargo transportation by two species of motor protein

The cargo motion in living cells transported by two species of motor protein with different intrinsic directionality is discussed in this study. Similar to single motor movement, cargo steps forward and backward along microtubule stochastically. Recent experiments found that, cargo transportation by two motor species has a memory, it does not change its direction as frequently as expected, which means that its forward and backward step rates depends on its previous motion trajectory. By assuming cargo has only the least memory, i.e. its step direction depends only on the direction of its last step, two cases of cargo motion are detailed analyzed in this study: {\bf (I)} cargo motion under constant external load; and {\bf (II)} cargo motion in one fixed optical trap. Due to the existence of memory, for the first case, cargo can keep moving in the same direction for a long distance. For the second case, the cargo will oscillate in the trap. The oscillation period decreases and the oscillation amplitude increases with the motor forward step rates, but both of them decrease with the trap stiffness. The most likely location of cargo, where the probability of finding the oscillated cargo is maximum, may be the same as or may be different with the trap center, which depends on the step rates of the two motor species. Meanwhile, if motors are robust, i.e. their forward to backward step rate ratios are high, there may be two such most likely locations, located on the two sides of the trap center respectively. The probability of finding cargo in given location, the probability of cargo in forward/backward motion state, and various mean first passage times of cargo to give location or given state are also analyzed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Microtubule length dependence of motor traffic in cells

In living cells, motor proteins, such as kinesin and dynein can move processively along microtubule (MT), and also detach from or attach to MT stochastically. Experiments have found that, the traffic of motor might be jammed, and various theoretical models have been designed to understand this traffic jam phenomenon. But previous studies mainly focus on motor attachment/detachment rate dependent properties. Recent experiment of Leduc {\it et al.} found that the traffic jam formation of motor protein kinesin depends also on the length of MT [Proc. Natl. Acad. Sci. U.S.A. {\bf 109}, 6100-6105 (2012)]. In this study, the MT length dependent properties of motor traffic will be analyzed. We found that MT length has one {\it critical value} $N_c$, traffic jam occurs only when MT length $N>N_c$. The jammed length of MT increases with total MT length, while the non-jammed MT length might not change monotonically with the total MT length. The critical value $N_c$ increases with motor detachment rate from MT, but decreases with motor attachment rate to MT.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

One-Dimensional Birth-Death Process and Delbrück-Gillespie Theory of Mesoscopic Nonlinear Chemical Reactions

As a mathematical theory for the stochasstic, nonlinear dynamics of individuals within a population, Delbrück-Gillespie process (DGP) $n(t)\in\mathbb{Z}^N$, is a birth-death system with state-dependent rates which contain the system size $V$ as a natural parameter. For large $V$, it is intimately related to an autonomous, nonlinear ordinary differential equation as well as a diffusion process. For nonlinear dynamical systems with multiple attractors, the quasi-stationary and stationary behavior of such a birth-death process can be underestood in terms of a separation of time scales by a $T^*\sim e^{αV}$ $(α>0)$: a relatively fast, intra-basin diffusion for $t\ll T^*$ and a much slower inter-basin Markov jump process for $t\gg T^*$. In the present paper for one-dimensional systems, we study both stationary behavior ($t=\infty$) in terms of invariant distribution $p_n^{ss}(V)$, and finite time dynamics in terms of the mean first passsage time (MFPT) $T_{n_1\rightarrow n_2}(V)$. We obtain an asymptotic expression of MFPT in terms of the "stochastic potential" $Φ(x,V)=-(1/V)\ln p^{ss}_{xV}(V)$. We show in general no continuous diffusion process can provide asymptotically accurate representations for both the MFPT and the $p_n^{ss}(V)$ for a DGP. When $n_1$ and $n_2$ belong to two different basins of attraction, the MFPT yields the $T^*(V)$ in terms of $Φ(x,V)\approx ϕ_0(x)+(1/V)ϕ_1(x)$. For systems with a saddle-node bifurcation and catastrophe, discontinuous "phase transition" emerges, which can be characterized by $Φ(x,V)$ in the limit of $V\rightarrow\infty$. In terms of time scale separation, the relation between deterministic, local nonlinear bifurcations and stochastic global phase transition is discussed. The one-dimensional theory is a pedagogic first step toward a general theory of DGP.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Comment on "Length-dependent translation of messenger RNA by ribosomes"

In the recent paper of Valleriani {\it et al} [Phys. Rev. E {\bf 83}, 042903 (2011)], a simple model for describing the translation of messenger RNA (mRNA) by ribosomes is presented, and an expression of the translational ratio $r$, defined as the ratio of translation rate $ω_{\rm tl}$ of protein from mRNA to degradation rate $ω_p$ of protein, is obtained. The key point to get this ratio $r$ is to get the translation rate $ω_{\rm tl}$. In the study of Valleriani {\it et al}, $ω_{\rm tl}$ is assumed to be the mean value of measured translation rate, i.e. the mean value of ratio of the translation number of protein to the lifetime of mRNA. However, in experiments different methods might be used to get $ω_{\rm tl}$. Therefore, for the sake of future application of their model to more experimental data analysis, in this comment three methods to get the translation rate $ω_{\rm tl}$, and consequently the translational ratio $r$, are provided. Based on one of the methods which might be employed in most of the experiments, we find that the translational ratio $r$ decays exponentially with the length of mRNA in prokaryotic cells, and decays reciprocally with the length of mRNA in eukaryotic cells. This result is slight different from that obtained in Valleriani's study.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Growth and shortening of microtubules: a two-state model approach

In this study, a two-state mechanochemical model is presented to describe the dynamic instability of microtubules (MTs) in cells. The MTs switches between two states, assembly state and disassembly state. In assembly state, the growth of MTs includes two processes: free GTP-tubulin binding to the tip of protofilament (PF) and conformation change of PF, during which the first tubulin unit which curls outwards is rearranged into MT surface using the energy released from the hydrolysis of GTP in the penultimate tubulin unit. In disassembly state, the shortening of MTs includes also two processes, the release of GDP-tibulin from the tip of PF and one new tubulin unit curls out of the MT surface. Switches between these two states, which are usually called rescue and catastrophe, happen stochastically with external force dependent rates. Using this two-state model with parameters obtained by fitting the recent experimental data, detailed properties of MT growth are obtained, we find that MT is mainly in assembly state, its mean growth velocity increases with external force and GTP-tubulin concentration, MT will shorten in average without external force. To know more about the external force and GTP-tubulin concentration dependent properties of MT growth, and for the sake of the future experimental verification of this two-state model, eleven {\it critical forces} are defined and numerically discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Loose mechanochemical coupling of molecular motors

In living cells, molecular motors convert chemical energy into mechanical work. Its thermodynamic energy efficiency, i.e. the ratio of output mechanical work to input chemical energy, is usually high. However, using two-state models, we found the motion of molecular motors is loosely coupled to the chemical cycle. Only part of the input energy can be converted into mechanical work. Others is dissipated into environment during substeps without contributions to the macro scale unidirectional movement.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Phenomenological analysis of ATP dependence of motor protein

In this study, through phenomenological comparison of the velocity-force data of processive motor proteins, including conventional kinesin, cytoplasmic dynein and myosin V, we found that, the ratio between motor velocities of two different ATP concentrations is almost invariant for any substall, superstall or negative external loads. Therefore, the velocity of motor can be well approximated by a Michaelis-Menten like formula $V=\atp k(F)L/(\atp +K_M)$, with $L$ the step size, and $k(F)$ the external load $F$ dependent rate of one mechanochemical cycle of motor motion in saturated ATP solution. The difference of Michaelis-Menten constant $K_M$ for substall, superstall and negative external load indicates, the ATP molecule affinity of motor head for these three cases are different, though the expression of $k(F)$ as a function of $F$ might be unchanged for any external load $F$. Verifications of this Michaelis-Menten like formula has also been done by fitting to the recent experimental data.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Shock of three-state model for intracellular transport of kinesin KIF1A

Recently, a three-state model is presented to describe the intracellular traffic of unconventional (single-headed) kinesin KIF1A [Phys. Rev. Lett. {\bf 95}, 118101 (2005)], in which each motor can bind strongly or weakly to its microtubule track, and each binding site of the track might be empty or occupied by one motor. As the usual two-state model, i.e. the totally asymmetric simple exclusion process (TASEP) with motor detachment and attachment, in steady state of the system, this three-state model also exhibits shock (or domain wall separating the high-density and low density phases) and boundary layers. In this study, using mean-field analysis, the conditions of existence of shock and boundary layers are obtained theoretically. Combined with numerical calculations, the properties of shock are also studied. This study will be helpful to understand the biophysical properties of the collective transport of kinesin KIF1A.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Cargo transport by several motors

In cells, organelles and vesicles are usually transported by cooperation of several motor proteins, including plus-end directed motor kinesin and minus-end directed motor dynein. Many biophysical models have been constructed to understand the mechanism of this motion. However, so far, the basic principle about it remains unclosed. In this paper, based on the recent experimental results and existing theoretical models, a spider-like model is provided. In this model, each motor is regarded as a bead-spring system. The bead can bind to or unbind from the track stochastically, and step forward or backward with fixed step size L and force dependent transition rates. The spring connects the bead to cargo tightly. At any time, the position of cargo is determined by force balance condition. The obvious characteristics of our model are that, the motors interact with each other and they do not share the external load equally. Our results indicate, the stall force of cargo, under which the mean velocity of cargo vanishes, usually decreases with the interactions between motors. If the cargo is pulled by several motors from same motor species, the stall force of cargo is bigger than that of the single motor, but usually smaller than their sum. However, if the cargo is pulled by two motors with different directionality, the stall force of cargo might bigger or smaller than the difference of the two single motor stall forces. The results imply, cooperation of several motors is helpful to pull big cargoes, though it is not so good as might be expected and no obvious help to improve the velocity.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Properties of modified periodic one-dimensional hopping model

One-dimensional hopping model is useful to describe the motion of microscopic particle in thermal noise environment, such as motor proteins. Recent experiments about the new generation of light-driven rotary molecular motors found that, the motor in state i can jump forward to state i+1 or i+2, or backward to state i-1 or i-2 directly. In this paper, such modified periodic one-dimensional hopping model of arbitrary period N is studied mathematically. The mean velocity, effective diffusion constant, and mean dwell time in one single cycle are obtained. Corresponding results are illustrated and verified by being applied to a type of synthetic rotary molecular motors.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

van't Hoff-Arrhenius Analysis of Mesoscopic and Macroscopic Dynamics of Simple Biochemical Systems: Stochastic vs. Nonlinear Bistabilities

Multistability of mesoscopic, driven biochemical reaction systems has implications to a wide range of cellular processes. Using several simple models, we show that one class of bistable chemical systems has a deterministic counterpart in the nonlinear dynamics based on the Law of Mass Action, while another class, widely known as noise-induced stochastic bistability, does not. Observing the system's volume ($V$) playing a similar role as the inverse temperature ($β$) in classical rate theory, an van't Hoff-Arrhenius like analysis is introduced. In one-dimensional systems, a transition rate between two states, represented in terms of a barrier in the landscape for the dynamics $Φ(x,V)$, $k\propto\exp\{-VΔΦ^‡(V)\}$, can be understood from a decomposition $ΔΦ^‡(V) \approxΔϕ_0^‡ Δϕ_1^‡/V$. Nonlinear bistability means $Δϕ_0^‡>0$ while stochastic bistability has $Δϕ_0^‡<0$ but $Δϕ_1^‡>0$. Stochastic bistabilities can be viewed as remants (or "ghosts) of nonlinear bifurcations or extinction phenomenon, and $Δϕ_0^‡$ and $Δϕ_1^‡$ as "enthalpic" and "entropic" barriers to a transition.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A general two-cycle network model of molecular motors

Molecular motors are single macromolecules that generate forces at the piconewton range and nanometer scale. They convert chemical energy into mechanical work by moving along filamentous structures. In this paper, we study the velocity of two-head molecular motors in the framework of a mechanochemical network theory. The network model, a generalization of the recently work of Liepelt and Lipowsky (PRL 98, 258102 (2007)), is based on the discrete mechanochemical states of a molecular motor with multiple cycles. By generalizing the mathematical method developed by Fisher and Kolomeisky for single cycle motor (PNAS(2001) 98(14) P7748-7753), we are able to obtain an explicit formula for the velocity of a molecular motor.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Limit properties of periodic one dimensional hopping model

Periodic one dimensional hopping model is useful to study the motion of microscopic particles, which lie in thermal noise environment. The mean velocity $V_N$ and diffusion constant $D_N$ of this model have been obtained by Bernard Derrida [J. Stat. Phys. 31 (1983) 433]. In this research, we will give the limits $V_D$ and $D_D$ of $V_N$ and $D_N$ as the number $N$ of mechanochemical sates in one period tends to infinity by formal calculation. It is well known that the stochastic motion of microscopic particles also can be described by overdamped Langevin dynamics and Fokker-Planck equation. Up to now, the corresponding formulations of mean velocity and effective diffusion coefficient, $V_L$ and $D_L$ in the framework of Langevin dynamics and $V_P, D_P$ in the framework of Fokker-Planck equation, have also been known. In this research, we will find that the formulations $V_D$ and $V_L, V_P$ are theoretically equivalent, and numerical comparison indicates that $D_D, D_L$, and $D_P$ are almost the same. Through the discussion in this research, we also can know more about the relationship between the one dimensional hopping model and Fokker-Planck equation.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Properties of tug-of-war model for cargo transport by molecular motors

Molecular motors are essential components for the biophysical functions of the cell. Our current quantitative understanding of how multiple motors move along a single track is not complete; even though models and theories for single motor chemomechanics abound. Recently, M.J.I. M$\ddot{\rm u}$ller {\em et al.} have developed a tug-of-war model to describe the bidirectional movement of the cargo (PNAS(2008) 105(12) P4609-4614). Through Monte Carlo simulations, they discovered that the tug-of-war model exhibits several qualitative different motility regimes, which depend on the precise value of single motor parameters, and they suggested the sensitivity can be used by a cell to regulate its cargo traffic. In the present paper, we carry out a thorough analysis of the tug-of-war model. All the stable, i.e., biophysically observable, steady states are obtained. Depending on several parameters, the system exhibits either uni-, bi- or tristability. Based on the separating boundary of the different stable states and the initial numbers of the different motor species that are bound to the track, the steady state of the cargo movement can be predicted, and consequently the steady state velocity can be obtained. It is found that, the velocity, even the direction, of the cargo movement change with the initial numbers of the motors which are bound to the track and several other parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The efficiency of the molecular motors

Molecular motors convert chemical energy into mechanical work while operating in an environment dominated by Brownian motion. The aim of this paper is to explore the flow of energy between the molecular motors and its surroundings, in particular, its efficiency. Based on the Fokker-Planck equation with either $N$ or infinite chemical states, we find that the energy efficiency of the molecular motors, whether the Stokes efficiency or the usual thermodynamic efficiency, is strictly bounded by 1, because of the dissipation of the energy in both the overdamped surroundings and in the process of the chemical reaction.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Domain wall of the totally asymmetric exclusion process without particle number conservation

In this research, the totally asymmetric exclusion process without particle number conservation is discussed. Based on the mean field approximation and the Rankine-Hugoniot condition, the necessary and sufficient conditions of the existence of the domain wall have been obtained. Moreover, the properties of the domain wall, including the location and height, have been studied theoretically. All the theoretical results are demonstrated by the computer simulations.

0 citations0 reviewsNo code linkedOpen access