Adiabatic reduction of models of stochastic gene expression with bursting

This paper considers adiabatic reduction in both discrete and continuous models of stochastic gene expression. In gene expression models, the concept of bursting is a production of several molecules simultaneously and is generally represented as a compound Poisson process of random size. In a general two-dimensional birth and death discrete model, we prove that under specific assumptions and scaling (that are characteristics of the mRNA-protein system) an adiabatic reduction leads to a one-dimensional discrete-state space model with bursting production. The burst term appears through the reduction of the first variable. In a two-dimensional continuous model, we also prove that an adiabatic reduction can be performed in a stochastic slow/fast system. In this gene expression model, the production of mRNA (the fast variable) is assumed to be bursty and the production of protein (the slow variable) is linear as a function of mRNA. When the dynamics of mRNA is assumed to be faster than the protein dynamics (due to a mRNA degradation rate larger than for the protein) we prove that, with the appropriate scaling, the bursting phenomena can be transmitted to the slow variable. We show that the reduced equation is either a stochastic differential equation with a jump Markov process or a deterministic ordinary differential equation depending on the scaling that is appropriate. These results are significant because adiabatic reduction techniques seem to have not been applied to a stochastic differential system containing a jump Markov process. Last but not least, for our particular system, the adiabatic reduction allows us to understand what are the necessary conditions for the bursting production-like of protein to occur.Comment: 24 page

Boundary value for a nonlinear transport equation emerging from a stochastic coagulation-fragmentation type model

We investigate the connection between two classical models of phase transition phenomena, the (discrete size) stochastic Becker-D\"oring, a continous time Markov chain model, and the (continuous size) deterministic Lifshitz-Slyozov model, a nonlinear transport partial differential equation. For general coefficients and initial data, we introduce a scaling parameter and prove that the empirical measure associated to the stochastic Becker-D\"oring system converges in law to the weak solution of the Lifshitz-Slyozov equation when the parameter goes to 0. Contrary to previous studies, we use a weak topology that includes the boundary of the state space (\ie\ the size $x=0$) allowing us to rigorously derive a boundary value for the Lifshitz-Slyozov model in the case of incoming characteristics. The condition reads $\lim_{x\to 0} (a(x)u(t)-b(x))f(t,x) = \alpha u(t)^2$ where $f$ is the volume distribution function, solution of the Lifshitz-Slyozov equation, $a$ and $b$ the aggregation and fragmentation rates, $u$ the concentration of free particles and $\alpha$ a nucleation constant emerging from the microscopic model. It is the main novelty of this work and it answers to a question that has been conjectured or suggested by both mathematicians and physicists. We emphasize that this boundary value depends on a particular scaling (as opposed to a modeling choice) and is the result of a separation of time scale and an averaging of fast (fluctuating) variables.Comment: 42 pages, 3 figures, video on supplementary materials at http://yvinec.perso.math.cnrs.fr/video.htm

Ornstein-Zernike equation and Percus-Yevick theory for molecular crystals

We derive the Ornstein-Zernike equation for molecular crystals of axially symmetric particles and apply the Percus-Yevick approximation to this system. The one-particle orientational distribution function has a nontrivial dependence on the orientation and is needed as an input. Despite some differences, the Ornstein-Zernike equation for molecular crystals has a similar structure as for liquids. We solve both equations for hard ellipsoids on a sc lattice. Compared to molecular liquids, the tensorial orientational correlators exhibit less structure. However, depending on the lengths a and b of the rotation axis and the perpendicular axes of the ellipsoids, different behavior is found. For oblate and prolate ellipsoids with b >= 0.35 (units of the lattice constant), damped oscillations in distinct directions of direct space occur for some correlators. They manifest themselves in some correlators in reciprocal space as a maximum at the Brillouin zone edge, accompanied by maxima at the zone center for other correlators. The oscillations indicate alternating orientational fluctuations, while the maxima at the zone center originate from nematic-like orientational fluctuations. For a <= 2.5 and b <= 0.35, the oscillations are weaker. For a >= 3.0 and b <= 0.35, no oscillations occur any longer. For many of the correlators in reciprocal space, an increase of a at fixed b leads to a divergence at the zone center q = 0, consistent with nematic-like long range fluctuations, and for some oblate and prolate systems with b ~< 1.0 a simultaneous tendency to divergence of few other correlators at the zone edge is observed. Comparison with correlators from MC simulations shows satisfactory agreement. We also obtain a phase boundary for order-disorder transitions.Comment: 20 pages, 13 figures, submitted to Phys. Rev.

Dynamic behavior of stochastic gene expression models in the presence of bursting

This paper considers the behavior of discrete and continuous mathematical models for gene expression in the presence of transcriptional/translational bursting. We treat this problem in generality with respect to the distribution of the burst size as well as the frequency of bursting, and our results are applicable to both inducible and repressible expression patterns in prokaryotes and eukaryotes. We have given numerous examples of the applicability of our results, especially in the experimentally observed situation that burst size is geometrically or exponentially distributed.Comment: 22 page

Adiabatic reduction of a model of stochastic gene expression with jump Markov process

This paper considers adiabatic reduction in a model of stochastic gene expression with bursting transcription considered as a jump Markov process. In this model, the process of gene expression with auto-regulation is described by fast/slow dynamics. The production of mRNA is assumed to follow a compound Poisson process occurring at a rate depending on protein levels (the phenomena called bursting in molecular biology) and the production of protein is a linear function of mRNA numbers. When the dynamics of mRNA is assumed to be a fast process (due to faster mRNA degradation than that of protein) we prove that, with appropriate scalings in the burst rate, jump size or translational rate, the bursting phenomena can be transmitted to the slow variable. We show that, depending on the scaling, the reduced equation is either a stochastic differential equation with a jump Poisson process or a deterministic ordinary differential equation. These results are significant because adiabatic reduction techniques seem to have not been rigorously justified for a stochastic differential system containing a jump Markov process. We expect that the results can be generalized to adiabatic methods in more general stochastic hybrid systems.Comment: 17 page

Computing largest circles separating two sets of segments

A circle $C$ separates two planar sets if it encloses one of the sets and its open interior disk does not meet the other set. A separating circle is a largest one if it cannot be locally increased while still separating the two given sets. An Theta(n log n) optimal algorithm is proposed to find all largest circles separating two given sets of line segments when line segments are allowed to meet only at their endpoints. In the general case, when line segments may intersect $\Omega(n^2)$ times, our algorithm can be adapted to work in O(n alpha(n) log n) time and O(n \alpha(n)) space, where alpha(n) represents the extremely slowly growing inverse of the Ackermann function.Comment: 14 pages, 3 figures, abstract presented at 8th Canadian Conference on Computational Geometry, 199

Microscopic theory of glassy dynamics and glass transition for molecular crystals

We derive a microscopic equation of motion for the dynamical orientational correlators of molecular crystals. Our approach is based upon mode coupling theory. Compared to liquids we find four main differences: (i) the memory kernel contains Umklapp processes, (ii) besides the static two-molecule orientational correlators one also needs the static one-molecule orientational density as an input, where the latter is nontrivial, (iii) the static orientational current density correlator does contribute an anisotropic, inertia-independent part to the memory kernel, (iv) if the molecules are assumed to be fixed on a rigid lattice, the tensorial orientational correlators and the memory kernel have vanishing l,l'=0 components. The resulting mode coupling equations are solved for hard ellipsoids of revolution on a rigid sc-lattice. Using the static orientational correlators from Percus-Yevick theory we find an ideal glass transition generated due to precursors of orientational order which depend on X and p, the aspect ratio and packing fraction of the ellipsoids. The glass formation of oblate ellipsoids is enhanced compared to that for prolate ones. For oblate ellipsoids with X <~ 0.7 and prolate ellipsoids with X >~ 4, the critical diagonal nonergodicity parameters in reciprocal space exhibit more or less sharp maxima at the zone center with very small values elsewhere, while for prolate ellipsoids with 2 <~ X <~ 2.5 we have maxima at the zone edge. The off-diagonal nonergodicity parameters are not restricted to positive values and show similar behavior. For 0.7 <~ X <~ 2, no glass transition is found. In the glass phase, the nonergodicity parameters show a pronounced q-dependence.Comment: 17 pages, 12 figures, accepted at Phys. Rev. E. v4 is almost identical to the final paper version. It includes, compared to former versions v2/v3, no new physical content, but only some corrected formulas in the appendices and corrected typos in text. In comparison to version v1, in v2-v4 some new results have been included and text has been change

Quantitative comparisons of forward problems in MEEG.

This document gives comparisons between several methods that solve the forward problem in MEEG by comparing their precision on a 3-layer spherical model. These methods are based on finite elements which either use surfacic meshes with triangles, volumic meshes with tetrahedra, or implicit elements deduced from levelsets

Evaluating Signs of Determinants Using Single-Precision Arithmetic

We propose a method of evaluating signs of 2×2 and 3×3 determinants with b-bit integer entries using only b and (b + 1)-bit arithmetic, respectively. This algorithm has numerous applications in geometric computation and provides a general and practical approach to robustness. The algorithm has been implemented and compared with other exact computation methods