Source author record

Jean-Charles Walter

Jean-Charles Walter appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

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

15 published item(s)

preprint2022arXiv

Modelling the effect of ribosome mobility on the rate of protein synthesis

Translation is one of the main steps in the synthesis of proteins. It consists of ribosomes that translate sequences of nucleotides encoded on mRNA into polypeptide sequences of amino acids. Ribosomes bound to mRNA move unidirectionally, while unbound ribosomes diffuse in the cytoplasm. It has been hypothesized that finite diffusion of ribosomes plays an important role in ribosome recycling and that mRNA circularization enhances the efficiency of translation. In order to estimate the effect of cytoplasmic diffusion on the rate of translation, we consider a Totally Asymmetric Simple Exclusion Process (TASEP) coupled to a finite diffusive reservoir, which we call the Ribosome Transport model with Diffusion (RTD). In this model, we derive an analytical expression for the rate of protein synthesis as a function of the diffusion constant of ribosomes, which is corroborated with results from continuous-time Monte Carlo simulations. Using a wide range of biological relevant parameters, we conclude that diffusion in biological cells is fast enough so that it does not play a role in controlling the rate of translation initiation.

preprint2022arXiv

Physical modeling of ribosomes along messenger RNA: estimating kinetic parameters from ribosome profiling experiments using a ballistic model

Gene expression consists in the synthesis of proteins from the information encoded on DNA. One of the two main steps of gene expression is the translation of messenger RNA (mRNA) into polypeptide sequences of amino acids. Here, by taking into account mRNA degradation, we model the motion of ribosomes along mRNA with a ballistic model where particles advance along a filament without excluded volume interactions. Unidirectional models of transport have previously been used to fit the average density of ribosomes obtained by the experimental ribo-sequencing (Ribo-seq) technique. In this case an inverse fit gives access to the kinetic rates: the position-dependent speeds and the entry rate of ribosomes onto mRNA. The degradation rate is not, however, accounted for and experimental data from different experiments are needed to have enough parameters for the fit. Here, we propose an entirely novel experimental setup and theoretical framework consisting in splitting the mRNAs into categories depending on the number of ribosomes from one to four. We solve analytically the ballistic model for a fixed number of ribosomes per mRNA, study the different regimes of degradation, and propose a criteria for the quality of the inverse fit. The proposed method provides a high sensitivity to the mRNA degradation rate. The additional equations coming from using the monosome (single ribosome) and polysome (arbitrary number) ribo-seq profiles enable us to determine all the kinetic rates in terms of the experimentally accessible mRNA degradation rate.

preprint2020arXiv

Modeling supercoiled DNA interacting with an anchored cluster of proteins: towards a quantitative estimation of chromosomal DNA supercoiling

We investigate the measurement of DNA supercoiling density ($σ$) along chromosomes using interaction frequencies between DNA and DNA-anchored clusters of proteins. Specifically, we show how the physics of DNA supercoiling leads, in bacteria, to the quantitative modeling of binding properties of ParB proteins around their centromere-like site, {\it parS}. Using this framework, we provide an upper bound for $σ$ in the {\it Escherichia coli} chromosome, consistent with plasmid values, and offer a proof of concept for a high accuracy measurement. To reach these conclusions, we revisit the problem of the formation of ParB clusters. We predict, in particular, that they result from a non-equilibrium, stationary balance between an influx of produced proteins and an outflux of excess proteins, i.e., they behave like liquid-like protein condensates with unconventional ``leaky'' boundaries.

preprint2016arXiv

Torque-Induced Rotational Dynamics in Polymers: Torsional Blobs and Thinning

By using the blob theory and computer simulations, we investigate the properties of a linear polymer performing a stationary rotational motion around a long impenetrable rod. In particular, in the simulations the rotation is induced by a torque applied to the end of the polymer that is tethered to the rod. Three different regimes are found, in close analogy with the case of polymers pulled by a constant force at one end. For low torques the polymer rotates maintaining its equilibrium conformation. At intermediate torques the polymer assumes a trumpet shape, being composed by blobs of increasing size. At even larger torques the polymer is partially wrapped around the rod. We derive several scaling relations between various quantities as angular velocity, elongation and torque. The analytical predictions match the simulation data well. Interestingly, we find a "thinning" regime where the torque has a very weak (logarithmic) dependence on the angular velocity. We discuss the origin of this behavior, which has no counterpart in polymers pulled by an applied force.

preprint2014arXiv

An introduction to Monte Carlo methods

Monte Carlo simulations are methods for simulating statistical systems. The aim is to generate a representative ensemble of configurations to access thermodynamical quantities without the need to solve the system analytically or to perform an exact enumeration. The main principles of Monte Carlo simulations are ergodicity and detailed balance. The Ising model is a lattice spin system with nearest neighbor interactions that is appropriate to illustrate different examples of Monte Carlo simulations. It displays a second order phase transition between a disordered (high temperature) and ordered (low temperature) phases, leading to different strategies of simulations. The Metropolis algorithm and the Glauber dynamics are efficient at high temperature. Close to the critical temperature, where the spins display long range correlations, cluster algorithms are more efficient. We introduce the rejection free (or continuous time) algorithm and describe in details an interesting alternative representation of the Ising model using graphs instead of spins with the Worm algorithm. We conclude with an important discussion of the dynamical effects such as thermalization and correlation time.

preprint2014arXiv

Hyperscaling above the upper critical dimension

Above the upper critical dimension, the breakdown of hyperscaling is associated with dangerous irrelevant variables in the renormalization group formalism at least for systems with periodic boundary conditions. While these have been extensively studied, there have been only a few analyses of finite-size scaling with free boundary conditions. The conventional expectation there is that, in contrast to periodic geometries, finite-size scaling is Gaussian, governed by a correlation length commensurate with the lattice extent. Here, detailed numerical studies of the five-dimensional Ising model indicate that this expectation is unsupported, both at the infinite-volume critical point and at the pseudocritical point where the finite-size susceptibility peaks. Instead the evidence indicates that finite-size scaling at the pseudocritical point is similar to that in the periodic case. An analytic explanation is offered which allows hyperscaling to be extended beyond the upper critical dimension.

preprint2014arXiv

Motor proteins traffic regulation by supply-demand balance of resources

In cells and in vitro assays the number of motor proteins involved in biological transport processes is far from being unlimited. The cytoskeletal binding sites are in contact with the same finite reservoir of motors (either the cytosol or the flow chamber) and hence compete for recruiting the available motors, potentially depleting the reservoir and affecting cytoskeletal transport. In this work we provide a theoretical framework to study, analytically and numerically, how motor density profiles and crowding along cytoskeletal filaments depend on the competition of motors for their binding sites. We propose two models in which finite processive motor proteins actively advance along cytoskeletal filaments and are continuously exchanged with the motor pool. We first look at homogeneous reservoirs and then examine the effects of free motor diffusion in the surrounding medium. We consider as a reference situation recent in vitro experimental setups of kinesin-8 motors binding and moving along microtubule filaments in a flow chamber. We investigate how the crowding of linear motor proteins moving on a filament can be regulated by the balance between supply (concentration of motor proteins in the flow chamber) and demand (total number of polymerised tubulin heterodimers). We present analytical results for the density profiles of bound motors, the reservoir depletion, and propose novel phase diagrams that present the formation of jams of motor proteins on the filament as a function of two tuneable experimental parameters: the motor protein concentration and the concentration of tubulins polymerized into cytoskeletal filaments. Extensive numerical simulations corroborate the analytical results for parameters in the experimental range and also address the effects of diffusion of motor proteins in the reservoir.

preprint2014arXiv

Rotational dynamics of entangled polymers

Some recent results on the rotational dynamics of polymers are reviewed and extended. We focus here on the relaxation of a polymer, either flexible or semiflexible, initially wrapped around a rigid rod. We also study the steady polymer rotation generated by a constant torque on the rod. The interplay of frictional and entropic forces leads to a complex dynamical behavior characterized by non-trivial universal exponents. The results are based on extensive simulations of polymers undergoing Rouse dynamics and on an analytical approach using force balance and scaling arguments. The analytical results are in general in good agreement with the simulations, showing how a simplified approach can correctly capture the complex dynamical behavior of rotating polymers.

preprint2014arXiv

Unwinding Dynamics of a Helically Wrapped Polymer

We study the rotational dynamics of a flexible polymer initially wrapped around a rigid rod and unwinding from it. This dynamics is of interest in several problems in biology and constitutes a fundamental instance of polymer relaxation from a state of minimal entropy. We investigate the dynamics of several quantities such as the total and local winding angles and metric quantities. The results of simulations performed in two and three dimensions, with and without self-avoidance, are explained by a theory based on scaling arguments and on a balance between frictional and entropic forces. The early stage of the dynamics is particularly rich, being characterized by three coexisting phases.

preprint2013arXiv

Unwinding relaxation dynamics of polymers

The relaxation dynamics of a polymer wound around a fixed obstacle constitutes a fundamental instance of polymer with twist and torque and it is of relevance also for DNA denaturation dynamics. We investigate it by simulations and Langevin equation analysis. The latter predicts a relaxation time scaling as a power of the polymer length times a logarithmic correction related to the equilibrium fluctuations of the winding angle. The numerical data support this result and show that at short times the winding angle decreases as a power-law. This is also in agreement with the Langevin equation provided a winding-dependent friction is used, suggesting that such reduced description of the system captures the basic features of the problem.

preprint2012arXiv

Fractional Brownian motion and the critical dynamics of zipping polymers

We consider two complementary polymer strands of length $L$ attached by a common end monomer. The two strands bind through complementary monomers and at low temperatures form a double stranded conformation (zipping), while at high temperature they dissociate (unzipping). This is a simple model of DNA (or RNA) hairpin formation. Here we investigate the dynamics of the strands at the equilibrium critical temperature $T=T_c$ using Monte Carlo Rouse dynamics. We find that the dynamics is anomalous, with a characteristic time scaling as $τ\sim L^{2.26(2)}$, exceeding the Rouse time $\sim L^{2.18}$. We investigate the probability distribution function, the velocity autocorrelation function, the survival probability and boundary behaviour of the underlying stochastic process. These quantities scale as expected from a fractional Brownian motion with a Hurst exponent $H=0.44(1)$. We discuss similarities and differences with unbiased polymer translocation.

preprint2012arXiv

Probing Hybridization parameters from microarray experiments: nearest neighbor model and beyond

In this article it is shown how optimized and dedicated microarray experiments can be used to study the thermodynamics of DNA hybridization for a large number of different conformations in a highly parallel fashion. In particular, free energy penalties for mismatches are obtained in two independent ways and are shown to be correlated with values from melting experiments in solution reported in the literature. The additivity principle, which is at the basis of the nearest-neighbor model, and according to which the penalty for two isolated mismatches is equal to the sum of the independent penalties, is thoroughly tested. Additivity is shown to break down for a mismatch distance below 5 nt. The behavior of mismatches in the vicinity of the helix edges, and the behavior of tandem mismatches are also investigated. Finally, some thermodynamic outlying sequences are observed and highlighted. These sequences contain combinations of GA mismatches. The analysis of the microarray data reported in this article provides new insights on the DNA hybridization parameters and can help to increase the accuracy of hybridization-based technologies.

preprint2011arXiv

Equilibrium winding angle of a polymer around a bar

The winding angle probability distribution of a planar self-avoiding walk has been known exactly since a long time: it has a gaussian shape with a variance growing as $<θ^2>\sim \ln L$. For the three-dimensional case of a walk winding around a bar, the same scaling is suggested, based on a first-order epsilon-expansion. We tested this three-dimensional case by means of Monte Carlo simulations up to length $L\approx25\,000$ and using exact enumeration data for sizes $L\le20$. We find that the variance of the winding angle scales as $<θ^2>\sim (\ln L)^{2α}$, with $α=0.75(1)$. The ratio $γ= <θ^4>/<θ^2>^2=3.74(5)$ is incompatible with the gaussian value $γ=3$, but consistent with the observation that the tail of the probability distribution function $p(θ)$ is found to decrease slower than a gaussian function. These findings are at odds with the existing first-order $ε$-expansion results.

preprint2011arXiv

Nonequilibrium effects in DNA microarrays: a multiplatform study

It has recently been shown that in some DNA microarrays the time needed to reach thermal equilibrium may largely exceed the typical experimental time, which is about 15h in standard protocols (Hooyberghs et al. Phys. Rev. E 81, 012901 (2010)). In this paper we discuss how this breakdown of thermodynamic equilibrium could be detected in microarray experiments without resorting to real time hybridization data, which are difficult to implement in standard experimental conditions. The method is based on the analysis of the distribution of fluorescence intensities I from different spots for probes carrying base mismatches. In thermal equilibrium and at sufficiently low concentrations, log I is expected to be linearly related to the hybridization free energy $ΔG$ with a slope equal to $1/RT_{exp}$, where $T_{exp}$ is the experimental temperature and R is the gas constant. The breakdown of equilibrium results in the deviation from this law. A model for hybridization kinetics explaining the observed experimental behavior is discussed, the so-called 3-state model. It predicts that deviations from equilibrium yield a proportionality of $\log I$ to $ΔG/RT_{eff}$. Here, $T_{eff}$ is an effective temperature, higher than the experimental one. This behavior is indeed observed in some experiments on Agilent arrays. We analyze experimental data from two other microarray platforms and discuss, on the basis of the results, the attainment of equilibrium in these cases. Interestingly, the same 3-state model predicts a (dynamical) saturation of the signal at values below the expected one at equilibrium.

preprint2011arXiv

Relaxation at finite temperature in Fully-Frustrated Ising Models

We consider by means of Monte Carlo simulations the relaxation in the paramagnetic phase of the anti-ferromagnetic Ising model on a triangular lattice and of a fully-frustrated Ising model on a square lattice. In contradistinction to previous studies of the second model, we show that spin-spin correlation functions do not decay with a stretched-exponential law at low temperature but that both models display an exponential decay with logarithmic corrections that are interpreted as the signature of topological defects.