Understanding the interactions between viruses and surfaces or interfaces is important, as they provide the principles underpinning the cleaning and disinfection of contaminated surfaces. Yet, the physics of such interactions is currently poorly unde
rstood. For instance, there are longstanding experimental observations suggesting that the presence of air-water interfaces can generically inactivate and kill viruses, yet the mechanism underlying this phenomenon remains unknown. Here we use theory and simulations to show that electrostatics provides one such mechanism, and that this is very general. Thus, we predict that the free energy of an RNA virus should increase by several thousands of $k_BT$ as the virion breaches an air-water interface. We also show that the fate of a virus approaching a generic liquid-liquid interface depends strongly on the detailed balance between interfacial and electrostatic forces, which can be tuned, for instance, by choosing different media to contact a virus-laden respiratory droplet. We propose that these results can be used to design effective strategies for surface disinfection. Intriguingly, tunability requires electrostatic and interfacial forces to scale similarly with viral size, which naturally occurs when charges are arranged in a double-shell distribution as in RNA viruses like influenza and all coronaviruses.
We discuss a polymer model for the 3D organization of human chromosomes. A chromosome is represented by a string of beads, with each bead being colored according to 1D bioinformatic data (e.g., chromatin state, histone modification, GC content). Indi
vidual spheres (representing bi- and multi-valent transcription factors) can bind reversibly and selectively to beads with the appropriate color. During molecular dynamics simulations, the factors bind, and the string spontaneously folds into loops, rosettes, and topologically-associating domains (TADs). This organization occurs in the absence of any specified interactions between distant DNA segments, or between transcription factors. A comparison with Hi-C data shows that simulations predict the location of most boundaries between TADs correctly. The model is fitting-free in the sense that it does not use Hi-C data as an input; consequently, one of its strengths is that it can -- in principle -- be used to predict the 3D organization of any region of interest, or whole chromosome, in a given organism, or cell line, in the absence of existing Hi-C data. We discuss how this simple model might be refined to include more transcription factors and binding sites, and to correctly predict contacts between convergent CTCF binding sites.
We present a generic coarse-grained model to describe molecular motors acting on polymer substrates, mimicking, for example, RNA polymerase on DNA or kinesin on microtubules. The polymer is modeled as a connected chain of beads; motors are represente
d as freely diffusing beads which, upon encountering the substrate, bind to it through a short-ranged attractive potential. When bound, motors and polymer beads experience an equal and opposite active force, directed tangential to the polymer; this leads to motion of the motors along the polymer contour. The inclusion of explicit motors differentiates our model from other recent active polymer models. We study, by means of Langevin dynamics simulations, the effect of the motor activity on both the conformational and dynamical properties of the substrate. We find that activity leads, in addition to the expected enhancement of polymer diffusion, to an effective reduction of its persistence length. We discover that this effective softening is a consequence of the emergence of double-folded branches, or hairpins, and that it can be tuned by changing the number of motors or the force they generate. Finally, we investigate the effect of the motors on the probability of knot formation. Counter-intuitively our simulations reveal that, even though at equilibrium a more flexible substrate would show an increased knotting probability, motor activity leads to a marked decrease in the occurrence of knotted conformations with respect to equilibrium.
Chromatin loop extrusion is a popular model for the formation of CTCF loops and topological domains. Recent HiC data have revealed a strong bias in favour of a particular arrangement of the CTCF binding motifs that stabilize loops, and extrusion is t
he only model to date which can explain this. However, the model requires a motor to generate the loops, and although cohesin is a strong candidate for the extruding factor, a suitable motor protein (or a motor activity in cohesin itself) has yet to be found. Here we explore a new hypothesis: that there is no motor, and thermal motion within the nucleus drives extrusion. Using theoretical modelling and computer simulations we ask whether such diffusive extrusion could feasibly generate loops. Our simulations uncover an interesting ratchet effect (where an osmotic pressure promotes loop growth), and suggest, by comparison to recent in vitro and in vivo measurements, that diffusive extrusion can in principle generate loops of the size observed in the data. Extra View on : C. A. Brackley, J. Johnson, D. Michieletto, A. N. Morozov, M. Nicodemi, P. R. Cook, and D. Marenduzzo Non-equilibrium chromosome looping via molecular slip-links, Physical Review Letters 119, 138101 (2017)
There are many proteins or protein complexes which have multiple DNA binding domains. This allows them to bind to multiple points on a DNA molecule (or chromatin fibre) at the same time. There are also many proteins which have been found to be able t
o compact DNA in vitro, and many others have been observed in foci or puncta when fluorescently labelled and imaged in vivo. In this work we study, using coarse-grained Langevin dynamics simulations, the compaction of polymers by simple model proteins and a phenomenon known as the bridging-induced attraction. The latter is a mechanism observed in previous simulations [Brackley et al., Proc. Natl. Acad. Sci. USA 110 (2013)], where proteins modelled as spheres form clusters via their multivalent interactions with a polymer, even in the absence of any explicit protein-protein attractive interactions. Here we extend this concept to consider more detailed model proteins, represented as simple patchy particles interacting with a semi-flexible bead-and-spring polymer. We find that both the compacting ability and the effect of the bridging-induced attraction depend on the valence of the model proteins. These effects also depend on the shape of the protein, which determines its ability to form bridges.
We study the effect of transcription on the kinetics of DNA supercoiling in 3D by means of Brownian dynamics simulations of a single nucleotide resolution coarse-grained model for double stranded DNA. By accounting for the action of a transcribing RN
A polymerase (RNAP), we characterise the geometry and non equilibrium dynamics of the twin supercoiling domains forming on each side of the RNAP. Textbook pictures depict such domains as symmetric, with plectonemes (writhed DNA) appearing close to the RNAP. On the contrary, we find that the twist generated by transcription results in asymmetric domains, with plectonemes formed far from the RNAP. We show that this translates into an action-at-a-distance on DNA-binding proteins: for instance, positive supercoils downstream of an elongating RNAP destabilise nucleosomes long before the transcriptional machinery reaches the histone octamer. To understand these observations we use our framework to quantitatively analyse the relaxation dynamics of supercoiled DNA. We find a striking separation of timescales between twist diffusion, which is a simple and fast process, and writhe relaxation, which is slow and entails multiple steps.
We study the compression and extension dynamics of a DNA-like polymer interacting with non-DNA binding and DNA-binding proteins, by means of computer simulations. The geometry we consider is inspired by recent experiments probing the compressional el
asticity of the bacterial nucleoid (DNA plus associated proteins), where DNA is confined into a cylindrical container and subjected to the action of a piston - a spherical bead to which an external force is applied. We quantify the effect of steric interactions (excluded volume) on the force-extension curves as the polymer is compressed. We find that non-DNA-binding proteins, even at low densities, exert an osmotic force which can be a lot larger than the entropic force exerted by the compressed DNA. The trends we observe are qualitatively robust with respect to changes in protein size, and are similar for neutral and charged proteins (and DNA). We also quantify the dynamics of DNA expansion following removal of the piston: while the expansion is well fitted by power laws, the apparent exponent depends on protein concentration, and protein-DNA interaction in a significant way. We further highlight an interesting kinetic process which we observe during the expansion of DNA interacting with DNA-binding proteins when the interaction strength is intermediate: the proteins bind while the DNA is packaged by the compression force, but they pop-off one-by-one as the force is removed, leading to a slow unzipping kinetics. Finally, we quantify the importance of supercoiling, which is an important feature of bacterial DNA in vivo.
We propose a model for the formation of chromatin loops based on the diffusive sliding of a DNA-bound factor which can dimerise to form a molecular slip-link. Our slip-links mimic the behaviour of cohesin-like molecules, which, along with the CTCF pr
otein, stabilize loops which organize the genome. By combining 3D Brownian dynamics simulations and 1D exactly solvable non-equilibrium models, we show that diffusive sliding is sufficient to account for the strong bias in favour of convergent CTCF-mediated chromosome loops observed experimentally. Importantly, our model does not require any underlying, and energetically costly, motor activity of cohesin. We also find that the diffusive motion of multiple slip-links along chromatin may be rectified by an intriguing ratchet effect that arises if slip-links bind to the chromatin at a preferred loading site. This emergent collective behaviour is driven by a 1D osmotic pressure which is set up near the loading point, and favours the extrusion of loops which are much larger than the ones formed by single slip-links.
Fluorescence microscopy reveals that the contents of many (membrane-free) nuclear bodies exchange rapidly with the soluble pool whilst the underlying structure persists; such observations await a satisfactory biophysical explanation. To shed light on
this, we perform large-scale Brownian dynamics simulations of a chromatin fiber interacting with an ensemble of (multivalent) DNA-binding proteins; these proteins switch between two states -- active (binding) and inactive (non-binding). This system provides a model for any DNA-binding protein that can be modified post-translationally to change its affinity for DNA (e.g., like the phosphorylation of a transcription factor). Due to this out-of-equilibrium process, proteins spontaneously assemble into clusters of self-limiting size, as individual proteins in a cluster exchange with the soluble pool with kinetics like those seen in photo-bleaching experiments. This behavior contrasts sharply with that exhibited by equilibrium, or non-switching, proteins that exist only in the binding state; when these bind to DNA non-specifically, they form clusters that grow indefinitely in size. Our results point to post-translational modification of chromatin-bridging proteins as a generic mechanism driving the self-assembly of highly dynamic, non-equilibrium, protein clusters with the properties of nuclear bodies. Such active modification also reshapes intra-chromatin contacts to give networks resembling those seen in topologically-associating domains, as switching markedly favors local (short-range) contacts over distant ones.
We propose a stochastic model for gene transcription coupled to DNA supercoiling, where we incorporate the experimental observation that polymerases create supercoiling as they unwind the DNA helix, and that these enzymes bind more favourably to regi
ons where the genome is unwound. Within this model, we show that when the transcriptionally induced flux of supercoiling increases, there is a sharp crossover from a regime where torsional stresses relax quickly and gene transcription is random, to one where gene expression is highly correlated and tightly regulated by supercoiling. In the latter regime, the model displays transcriptional bursts, waves of supercoiling, and up-regulation of divergent or bidirectional genes. It also predicts that topological enzymes which relax twist and writhe should provide a pathway to down-regulate transcription. This article has been published in Physical Review Letters, May 2016.
