Enhanced diffusion of tracer particles in nonreciprocal mixtures.

We study the diffusivity of a tagged particle in a binary mixture of Brownian particles with nonreciprocal interactions. Numerical simulations reveal that, for a broad class of interaction potentials, nonreciprocity can significantly increase the long-time diffusion coefficient of tracer particles and that this diffusion enhancement is associated with a breakdown of the Einstein relation. These observations are quantified and confirmed via two different and complementary analytical approaches: (i) a linearized stochastic density field theory, which is particularly accurate in the limit of soft interactions, and (ii) a reduced two-body description, which is exact at leading order in the density of particles. The latter reveals that diffusion enhancement can be attributed to the formation of transiently propelled dimers of particles, whose cohesion and speed are controlled by the nonreciprocal interactions.

Isotropic active colloids: explicit vs. implicit descriptions of propulsion mechanisms.

Modeling the couplings between active particles often neglects the possible many-body effects that control the propulsion mechanism. Accounting for such effects requires the explicit modeling of the molecular details at the origin of activity. Here, we take advantage of a recent two-dimensional model of isotropic active particles whose propulsion originates from the interactions between solute particles in the bath. The colloid catalyzes a chemical reaction in its vicinity, which results in a local phase separation of solute particles, and the density fluctuations of solute particles cause the enhanced diffusion of the colloid. In this paper, we investigate an assembly of such active particles, using (i) an explicit model, where the microscopic dynamics of the solute particles is accounted for; and (ii) an implicit model, whose parameters are inferred from the explicit model at infinite dilution. In the explicit solute model, the long-time diffusion coefficient of the active colloids strongly decreases with density, an effect which is not captured by the derived implicit model. This suggests that classical models, which usually decouple pair interactions from activity, fail to describe collective dynamics in active colloidal systems driven by solute-solute interactions.

Confined Electrolytes Show Bulk Dynamics Modulated by Hydrodynamic Couplings with the Walls.

We use numerical simulations at the mesoscopic scale, namely, multiparticle collision dynamics (MPCD), to investigate the properties of electrolyte solutions in a charged slit pore. The solution is described within the primitive model of electrolytes, where ions are charged hard spheres embedded in a dielectric medium. Hydrodynamic couplings between ions and with the charged walls are precisely accounted for by the MPCD algorithm. We show that the dynamic properties of ions in this situation strongly differ from the behavior at infinite dilution (ideal case), contrary to what is usually assumed in the usual Poisson-Nernst-Planck description of this kind of systems. As a consequence of confinement, the diffusion coefficients of ions unexpectedly increase with the average ionic density in the systems. This is due to a decrease of the proportion of ions that are slowed down by the wall. Moreover, nonequilibrium simulations are used to estimate the electrical conductivity of these confined electrolytes. We show that the simulation results can be accounted for quantitatively by combining bulk descriptions of the electrical conductivity of electrolytes with a simple description of the hydrodynamics of ions in a slit pore.

Conditions for the propulsion of a colloid surrounded by a mesoscale phase separation.

We study a two-dimensional model of an active isotropic colloid whose propulsion is linked to the interactions between solute particles of the bath. The colloid catalyzes a chemical reaction in its vicinity, that yields a local phase separation of solute particles. The density fluctuations of solute particles result in the enhanced diffusion of the colloid. Using numerical simulations, we thoroughly investigate the conditions under which activity occurs, and we establish a state diagram for the activity of the colloid as a function of the parameters of the model. We use the generated data to unravel a key observable that controls the existence and the intensity of activity: The filling fraction of the reaction area. Remarkably, we finally show that propulsion also occurs in three-dimensional geometries, which confirms the interest of this mechanism for experimental applications.

Electroosmotic Flow Induced Lift Forces on Polymer Chains in Nanochannels.

A major objective of research in nanofluidics is to achieve better selectivity in manipulating the fluxes of nano-objects and in particular of biopolymers. Numerical simulations allow one to better understand the physical mechanisms at play in such situations. We performed hybrid mesoscale simulations to investigate the properties of polymers under flows in slit pores at the nanoscale. We use multiparticle collision dynamics, an algorithm that includes hydrodynamics and thermal fluctuations, to investigate the properties of fully flexible and stiff polymers under several types of flow, showing that Poiseuille flows and electroosmotic flows can lead to quantitatively and qualitatively different behaviors of the chain. In particular, a counterintuitive phenomenon occurs in the presence of an electroosmotic flow: When the monomers are attracted by the solid surfaces through van der Waals forces, shear-induced forces lead to a stronger repulsion of the polymers from these surfaces. Such focusing of the chain in the middle of the channel increases its flowing velocity, a phenomenon that may be exploited to separate different types of polymers.

Diffusion of a tracer in a dense mixture of soft particles connected to different thermostats.

We study the dynamics of a tracer in a dense mixture of particles connected to different thermostats. Starting from the overdamped Langevin equations that describe the evolution of the system, we derive the expression of the self-diffusion coefficient of a tagged particle in the suspension, in the limit of soft interactions between the particles. Our derivation, which relies on the linearization of the Dean-Kawasaki equations obeyed by the density fields and on a path-integral representation of the dynamics of the tracer, extends previous derivations that held for tracers in contact with a single bath. Our analytical result is confronted to results from Brownian dynamics simulations. The agreement with numerical simulations is very good even for high densities. We show how the diffusivity of tracers can be affected by the activity of a dense environment of soft particles that may represent polymer coils-a result that could be of relevance in the interpretation of measurements of diffusivity in biological media. Finally, our analytical result is general and can be applied to the diffusion of tracers coupled to different types of fluctuating environments, provided that their evolution equations are linear and that the coupling between the tracer and the bath is weak.

Spontaneous propulsion of an isotropic colloid in a phase-separating environment.

The motion of active colloids is generally achieved through their anisotropy, as exemplified by Janus colloids. Recently, there was a growing interest in the propulsion of isotropic colloids, which requires some local symmetry breaking. Although several mechanisms for such propulsion were proposed, little is known about the role played by the interactions within the environment of the colloid, which can have a dramatic effect on its propulsion. Here, we propose a minimal model of an isotropic colloid in a bath of solute particles that interact with each other. These interactions lead to a spontaneous phase transition close to the colloid, to directed motion of the colloid over very long timescales and to significantly enhanced diffusion, in spite of the crowding induced by solute particles. We determine the range of parameters where this effect is observable in the model, and we propose an effective Langevin equation that accounts for it and allows one to determine the different contributions at stake in self-propulsion and enhanced diffusion.

Coarse-Grained Models of Aqueous Solutions of Polyelectrolytes: Significance of Explicit Charges.

The structure of polyelectrolytes is highly sensitive to small changes in interactions between their monomers. In particular, interactions mediated by counterions play a significant role and are affected by both specific molecular effects and generic concentration effects. The ability of coarse-grained models to reproduce the structural properties of an atomic model is thus a challenging task. Our present study compares the ability of different kinds of coarse-grained models: (i) to reproduce the structure of an atomistic model of a polyelectrolyte (the sodium polyacrylate) and (ii) to reproduce the variations of this structure with the number of monomers and with the concentration of different species. We show that adequate scalings of the gyration radius of the polymer Rg with the number of monomers N and with the box size Lbox are only obtained, first, if the monomer charges and the counterions are explicitly described and, second, if an attractive Lennard-Jones contribution is added to the interaction between distant monomers. Also, we show that implicit ion models are relevant only to the high electrostatic screening regime.

Rouse model with transient intramolecular contacts on a timescale of seconds recapitulates folding and fluctuation of yeast chromosomes.

DNA folding and dynamics along with major nuclear functions are determined by chromosome structural properties, which remain, thus far, elusive in vivo. Here, we combine polymer modeling and single particle tracking experiments to determine the physico-chemical parameters of chromatin in vitro and in living yeast. We find that the motion of reconstituted chromatin fibers can be recapitulated by the Rouse model using mechanical parameters of nucleosome arrays deduced from structural simulations. Conversely, we report that the Rouse model shows some inconsistencies to analyze the motion and structural properties inferred from yeast chromosomes determined with chromosome conformation capture techniques (specifically, Hi-C). We hence introduce the Rouse model with Transient Internal Contacts (RouseTIC), in which random association and dissociation occurs along the chromosome contour. The parametrization of this model by fitting motion and Hi-C data allows us to measure the kinetic parameters of the contact formation reaction. Chromosome contacts appear to be transient; associated to a lifetime of seconds and characterized by an attractive energy of -0.3 to -0.5 kBT. We suggest attributing this energy to the occurrence of histone tail-DNA contacts and notice that its amplitude sets chromosomes in 'theta' conditions, in which they are poised for compartmentalization and phase separation.

Polymer coil-globule phase transition is a universal folding principle of Drosophila epigenetic domains.

BACKGROUND: Localized functional domains within chromosomes, known as topologically associating domains (TADs), have been recently highlighted. In Drosophila, TADs are biochemically defined by epigenetic marks, this suggesting that the 3D arrangement may be the "missing link" between epigenetics and gene activity. Recent observations (Boettiger et al. in Nature 529(7586):418-422, 2016) provide access to structural features of these domains with unprecedented resolution thanks to super-resolution experiments. In particular, they give access to the distribution of the radii of gyration for domains of different linear length and associated with different transcriptional activity states: active, inactive or repressed. Intriguingly, the observed scaling laws lack consistent interpretation in polymer physics. RESULTS: We develop a new methodology conceived to extract the best information from such super-resolution data by exploiting the whole distribution of gyration radii, and to place these experimental results on a theoretical framework. We show that the experimental data are compatible with the finite-size behavior of a self-attracting polymer. The same generic polymer model leads to quantitative differences between active, inactive and repressed domains. Active domains behave as pure polymer coils, while inactive and repressed domains both lie at the coil-globule crossover. For the first time, the "color-specificity" of both the persistence length and the mean interaction energy are estimated, leading to important differences between epigenetic states. CONCLUSION: These results point toward a crucial role of criticality to enhance the system responsivity, resulting in both energy transitions and structural rearrangements. We get strong indications that epigenetically induced changes in nucleosome-nucleosome interaction can cause chromatin to shift between different activity states.

Computation of the Hydrodynamic Radius of Charged Nanoparticles from Nonequilibrium Molecular Dynamics.

We have used nonequilibrium molecular dynamics to simulate the flow of water molecules around a charged nanoparticle described at the atomic scale. These nonequilibrium simulations allowed us to compute the friction coefficient of the nanoparticle and then to deduce its hydrodynamic radius. We have compared two different strategies to thermostat the simulation box, since the low symmetry of the flow field renders the control of temperature non trivial. We show that both lead to an adequate control of the temperature of the system. To deduce the hydrodynamic radius of the nanoparticle we have employed a partial thermostat, which exploits the cylindrical symmetry of the flow field. Thereby, only a part of the simulation box far from the nanoparticle is thermostated. We have taken into account the finite concentration of the nanoparticle when calculating the friction force acting on it. We have focused on the case of polyoxometalate ions, which are inorganic charged nanoparticles. It appears that, for a given structure of the nanoparticle at the atomic level, the hydrodynamic radius significantly increases with the nanoparticles charge, a phenomenon that had not been quantified so far using molecular dynamics. The presence of an added salt only slightly modifies the hydrodynamic radius.

Comparison of different coupling schemes between counterions and charged nanoparticles in multiparticle collision dynamics.

We applied the multiparticle collision dynamics (MPC) simulation technique to highly asymmetric electrolytes in solution, i.e., charged nanoparticles and their counterions in a solvent. These systems belong to a domain of solute size which ranges between the electrolyte and the colloidal domains, where most analytical theories are expected to fail, and efficient simulation techniques are still missing. MPC is a mesoscopic simulation method which mimics hydrodynamics properties of a fluid, includes thermal fluctuations, and can be coupled to a molecular dynamics of solutes. We took advantage of the size asymmetry between nanoparticles and counterions to treat the coupling between solutes and the solvent bath within the MPC method. Counterions were coupled to the solvent bath during the collision step and nanoparticles either through a direct interaction force or with stochastic rotation rules which mimic stick boundary conditions. Moreover, we adapted the simulation procedure to address the issue of the strong electrostatic interactions between solutes of opposite charges. We show that the short-ranged repulsion between counterions and nanoparticles can be modeled by stochastic reflection rules. This simulation scheme is very efficient from a computational point of view. We have also computed the transport coefficients for various densities. The diffusion of counterions was found in one case to increase slightly with the volume fraction of nanoparticles. The deviation of the electric conductivity from the ideal behavior (solutes at infinite dilution without any direct interactions) is found to be strong.

Self-diffusion and activity coefficients of ions in charged disordered media.

Self-diffusion and single ion activity coefficients of ions of size symmetric electrolytes were studied in the presence of a collection of charged obstacles (called matrix) within a "soft" version of the primitive model of electrolyte solutions. The matrix subsystem possesses a net charge, depending on the concentration and charge of obstacles. The brownian dynamics method was used to calculate the self-diffusion coefficients of mobile species. The replica Ornstein-Zernike theory for the partly quenched systems was used to calculate the individual activity coefficients of mobile ionic species. The results reflect the competition between attractive (obstacle-counterion, co-ion-counterion), and repulsive (obstacle-co-ion) interactions in these model systems. For the simplest possible system of symmetric monovalent ions the latter effect wins: Co-ions are excluded from the area around obstacles, and this slows down their diffusion compared to that of counterions. Introduction of divalent charges into the system changes this result when the concentration of obstacles is low. We compare these results to those obtained for the corresponding fully annealed systems, i.e., where all the species are mobile. In most cases the self-diffusion and activity coefficients of counterions and co-ions in the presence of charged obstacles follow the trends of the fully annealed solution, which are dictated by the composition of the mixture. In few situations, however, the presence of charged obstacles modifies these trends. Our study allows us to clearly identify the effects due to obstacles, and to separate them from those arising from the composition of the solution. In the case of charge and size symmetric systems, the results for the individual activity coefficients fully support the hypothesis of the "electrostatic excluded volume". Thermodynamic and dynamic results are consistent in explaining the behavior of the systems studied.

Coordinate linkage of HIV evolution reveals regions of immunological vulnerability.

Cellular immune control of HIV is mediated, in part, by induction of single amino acid mutations that reduce viral fitness, but compensatory mutations limit this effect. Here, we sought to determine if higher order constraints on viral evolution exist, because some coordinately linked combinations of mutations may hurt viability. Immune targeting of multiple sites in such a multidimensionally conserved region might render the virus particularly vulnerable, because viable escape pathways would be greatly restricted. We analyzed available HIV sequences using a method from physics to reveal distinct groups of amino acids whose mutations are collectively coordinated ("HIV sectors"). From the standpoint of mutations at individual sites, one such group in Gag is as conserved as other collectively coevolving groups of sites in Gag. However, it exhibits higher order conservation indicating constraints on the viability of viral strains with multiple mutations. Mapping amino acids from this group onto protein structures shows that combined mutations likely destabilize multiprotein structural interactions critical for viral function. Persons who durably control HIV without medications preferentially target the sector in Gag predicted to be most vulnerable. By sequencing circulating viruses from these individuals, we find that individual mutations occur with similar frequency in this sector as in other targeted Gag sectors. However, multiple mutations within this sector are very rare, indicating previously unrecognized multidimensional constraints on HIV evolution. Targeting such regions with higher order evolutionary constraints provides a novel approach to immunogen design for a vaccine against HIV and other rapidly mutating viruses.

Effective interaction between charged nanoparticles and DNA.

We investigate the effective interaction mediated by salt ions between charged nanoparticles (NPs) and DNA. DNA is modeled as an infinite cylinder with a constant surface charge in an implicit solvent. Monte Carlo simulations are used to compute the free energy of the system described in the framework of the primitive model of electrolytes, which accounts for excluded volumes of salt ions. A mean-field Poisson-Boltzmann theory also allows us to compute the free energy and provides us with explicit formulae for its main characteristics (position and depth of the minimum). We intend here to identify the physical parameters that have a major impact on the NP-DNA interaction, in an attempt to evaluate physico-chemical properties which could play a role in genotoxicity or, which could be exploited for therapeutic use. Thus, we investigate the influence on the effective interaction of: the shape of the nanoparticle, the magnitude of the nanoparticle charge and its distribution, the value of the pH of the solution, the magnitude of Van der Waals interactions depending on the nature of the constitutive material of the NP (metal vs. dielectric). We show that for positively charged concave NPs the effective interaction is repulsive at short distance, so that it presents a minimum at distance from the DNA. This short-range repulsion is specific to indented particles and is a robust property that holds for a large range of materials and charge densities.

Ion-mediated interactions in suspensions of oppositely charged nanoparticles.

The structure of oppositely charged spherical nanoparticles (polyions), dispersed in ionic solutions with continuous solvent (primitive model), is investigated by Monte Carlo (MC) simulations, within explicit and implicit microion representations, over a range of polyion valences and densities, and microion concentrations. Systems with explicit microions are explored by semigrand canonical MC simulations, and allow density-dependent effective polyion pair potentials v(alphabeta) (eff)(r) to be extracted from measured partial pair distribution functions. Implicit microion MC simulations are based on pair potentials of mean force v(alphabeta) ((2))(r) computed by explicit microion simulations of two charged polyions, in the low density limit. In the vicinity of the liquid-gas separation expected for oppositely charged polyions, the implicit microion representation leads to an instability against density fluctuations for polyion valences mid R:Zmid R: significantly below those at which the instability sets in within the exact explicit microion representation. Far from this instability region, the v(alphabeta) ((2))(r) are found to be fairly close to but consistently more repulsive than the effective pair potentials v(alphabeta) (eff)(r). This is corroborated by additional calculations of three-body forces between polyion triplets, which are repulsive when one polyion is of opposite charge to the other two. The explicit microion MC data were exploited to determine the ratio of salt concentrations c and c(o) within the dispersion and the reservoir (Donnan effect). c/c(o) is found to first increase before finally decreasing as a function of the polyion packing fraction.

Nonspecific DNA-protein interaction: why proteins can diffuse along DNA.

Recent single molecule experiments have reported that DNA binding proteins (DNA-BPs) can diffuse along DNA. This suggests that interactions between proteins and DNA play a role during the target search even far from their specific site on DNA. Here we show by means of Monte Carlo simulations and analytical calculations that there is a counterintuitive repulsion between the two oppositely charged macromolecules at a nanometer range. For the concave shape of DNA-BPs, and for realistic protein charge densities, we find that the DNA-protein interaction free energy has a minimum at a finite surface-to-surface separation, in which proteins can easily slide. When a protein encounters its target, the free energy barrier is completely counterbalanced by the H-bond interaction, thus enabling the sequence recognition.

Toward the description of electrostatic interactions between globular proteins: potential of mean force in the primitive model.

Monte Carlo simulations are used to calculate the exact potential of mean force between charged globular proteins in aqueous solution. The aim of the present paper is to study the influence of the ions of the added salt on the effective interaction between these nanoparticles. The charges of the model proteins, either identical or opposite, are either central or distributed on a discrete pattern. Contrarily to Poisson-Boltzmann predictions, attractive, and repulsive direct forces between proteins are not screened similarly. Moreover, it has been shown that the relative orientations of the charge patterns strongly influence salt-mediated interactions. More precisely, for short distances between the proteins, ions enhance the difference of the effective forces between (i) like-charged and oppositely charged proteins, (ii) attractive and repulsive relative orientations of the proteins, which may affect the selectivity of protein/protein recognition. Finally, such results observed with the simplest models are applied to a more elaborate one to demonstrate their generality.

New coarse-graining procedure for the dynamics of charged spherical nanoparticles in solution.

A multiscale strategy based on the Brownian dynamics (BD) simulation method is presented here. It leads to an approximate but realistic reproduction of the dynamics of charged nanoparticles in suspension. This method is particularly suited to systems containing highly dissymmetric electrolytes with added salts, such as micellar suspensions or protein solutions. The coarse-graining procedure leads to a description where only the translational degrees of freedom of the nanoparticles are left, all the degrees of freedom related to the smallest solutes being rigorously averaged out. The authors' contribution aims at quantitatively evaluating the influence of the eliminated forces on the dynamics of the nanoparticles. For this purpose, an effective diffusion coefficient has to be calculated. In practice, this effective diffusion coefficient is taken as an input of a coarse-grained simulation that uses the potential of mean force between nanoparticles. The procedure has been validated by the quantitative comparison between the coarse-grained calculations and BD simulations at the "microscopic" level of description (which explicitly include microions). For a model of aqueous solutions of 10-1 electrolyte with a 1-1 added salt, the agreement is found to be excellent. This new method allows us to compute the diffusion coefficients of nanoparticles with a computation time at least one order of magnitude lower than with explicit BD.