Simulating chromatin is crucial for predicting genome organization and dynamics. Although coarse-grained bead-spring polymer models are commonly used to describe chromatin, the relevant bead ...dimensions, elastic properties, and the nature of inter-bead potentials are unknown. Using nucleosome-resolution contact probability (Micro-C) data, we systematically coarse-grain chromatin and predict quantities essential for polymer representation of chromatin. We compute size distributions of chromatin beads for different coarse-graining scales, quantify fluctuations and distributions of bond lengths between neighboring regions, and derive effective spring constant values. Unlike the prevalent notion, our findings argue that coarse-grained chromatin beads must be considered as soft particles that can overlap, and we derive an effective inter-bead soft potential and quantify an overlap parameter. We also compute angle distributions giving insights into intrinsic folding and local bendability of chromatin. While the nucleosome-linker DNA bond angle naturally emerges from our work, we show two populations of local structural states. The bead sizes, bond lengths, and bond angles show different mean behavior at Topologically Associating Domain (TAD) boundaries and TAD interiors. We integrate our findings into a coarse-grained polymer model and provide quantitative estimates of all model parameters, which can serve as a foundational basis for all future coarse-grained chromatin simulations.
Collapsed conformations of chromatin have been long suspected of being mediated by interactions with multivalent binding proteins, which can bring together distant sections of the chromatin fiber. In ...this study, we use Langevin dynamics simulation of a coarse grained chromatin polymer to show that the role of binding proteins can be more nuanced than previously suspected. In particular, for chromatin polymer in confinement, entropic forces can drive reswelling of collapsed chromatin with increasing binder concentrations, and this reswelling transition happens at physiologically relevant binder concentrations. Both the extent of collapse, and also of reswelling depends on the strength of confinement. We also study the kinetics of collapse and reswelling and show that both processes occur in similar timescales. We characterise this reswelling of chromatin in biologically relevant regimes and discuss the non-trivial role of multivalent binding proteins in mediating the spatial organisation of the genome.
Multivalent binding proteins can drive reswelling of collapsed chromatin beyond a critical binder concentration. This reswelling is driven by volume exclusion and entropic forces, and has implications for spatial organisation of chromatin inside the nucleus.
An important question in the context of the three-dimensional organization of chromosomes is the mechanism of formation of large loops between distant basepairs. Recent experiments suggest that the ...formation of loops might be mediated by loop extrusion factor proteins such as cohesin. Experiments on cohesin have shown that cohesins walk diffusively on the DNA and that nucleosomes act as obstacles to the diffusion, lowering the permeability and hence reducing the effective diffusion constant. An estimation of the times required to form the loops of typical sizes seen in Hi-C experiments using these low-effective-diffusion constants leads to times that are unphysically large. The puzzle then is the following: how does a cohesin molecule diffusing on the DNA backbone achieve speeds necessary to form the large loops seen in experiments? We propose a simple answer to this puzzle and show that although at low densities, nucleosomes act as barriers to cohesin diffusion, beyond a certain concentration they can reduce loop formation times because of a subtle interplay between the nucleosome size and the mean linker length. This effect is further enhanced on considering stochastic binding kinetics of nucleosomes on the DNA backbone and leads to predictions of lower loop formation times than might be expected from a naive obstacle picture of nucleosomes.
Dynein motors exhibit catch bonding, where the unbinding rate of the motors from microtubule filaments decreases with increasing opposing load. The implications of this catch bond on the transport ...properties of dynein-driven cargo are yet to be fully understood. In this context, optical trapping assays constitute an important means of accurately measuring the forces generated by molecular motor proteins. We investigate, using theory and stochastic simulations, the transport properties of cargo transported by catch bonded dynein molecular motors - both singly and in teams - in a harmonic potential, which mimics the variable force experienced by cargo in an optical trap. We estimate the biologically relevant measures of first passage time - the time during which the cargo remains bound to the microtubule and detachment force - the force at which the cargo unbinds from the microtubule, using both two-dimensional and one-dimensional force balance frameworks. Our results suggest that even for cargo transported by a single motor, catch bonding may play a role depending on the force scale which marks the onset of the catch bond. By comparing with experimental measurements on single dynein-driven transport, we estimate realistic bounds of this catch bond force scale. Generically, catch bonding results in increased persistent motion, and can also generate non-monotonic behaviour of first passage times. For cargo transported by multiple motors, emergent collective effects due to catch bonding can result in non-trivial re-entrant phenomena wherein average first passage times and detachment forces exhibit non-monotonic behaviour as a function of the stall force and the motor velocity.
A theoretical study of cargo transport by dynein motors in an optical trap shows non-trivial transport behaviour due to the dynein catch bond. We compare our model with experimental data to estimate the force scale at which the catch bond sets in.
We propose a simple model for chromatin organization based on the interaction of the chromatin fibers with lamin proteins along the nuclear membrane. Lamin proteins are known to be a major factor ...that influences chromatin organization and hence gene expression in the cells. We provide a quantitative understanding of lamin-associated chromatin organization in a crowded macromolecular environment by systematically varying the heteropolymer segment distribution and the strength of the lamin-chromatin attractive interaction. Our minimal polymer model reproduces the formation of lamin-associated-domains and provides an in silico tool for quantifying domain length distributions for different distributions of heteropolymer segments. We show that a Gaussian distribution of heteropolymer segments, coupled with strong lamin-chromatin interactions, can qualitatively reproduce observed length distributions of lamin-associated-domains. Further, lamin-mediated interaction can enhance the formation of chromosome territories as well as the organization of chromatin into tightly packed heterochromatin and the loosely packed gene-rich euchromatin regions.
Generation of mechanical oscillations is ubiquitous to a wide variety of intracellular processes, ranging from activity of muscle fibers to oscillations of the mitotic spindle. The activity of motors ...plays a vital role in maintaining the integrity of the mitotic spindle structure and generating spontaneous oscillations. Although the structural features and properties of the individual motors are well characterized, their implications on the functional behavior of motor-filament complexes are more involved. We show that force-induced allosteric deformations in dynein, which result in catchbonding behavior, provide a generic mechanism to generate spontaneous oscillations in motor-cytoskeletal filament complexes. The resultant phase diagram of such motor-filament systems—characterized by force-induced allosteric deformations—exhibits bistability and sustained limit-cycle oscillations in biologically relevant regimes, such as for catchbonded dynein. The results reported here elucidate the central role of this mechanism in fashioning a distinctive stability behavior and oscillations in motor-filament complexes such as mitotic spindles.
Abstract
In an attempt to understand the role of dysregulated circadian rhythm in glioma, our recent findings highlighted the existence of a feed-forward loop between tumour metabolite lactate, ...pro-inflammatory cytokine IL-1β and circadian CLOCK. To further elucidate the implication of this complex interplay, we developed a mathematical model that quantitatively describes this lactate dehydrogenase A (LDHA)-IL-1β-CLOCK/BMAL1 circuit and predicts potential therapeutic targets. The model was calibrated on quantitative western blotting data in two glioma cell lines in response to either lactate inhibition or IL-1β stimulation. The calibrated model described the experimental data well and most of the parameters were identifiable, thus the model was predictive. Sensitivity analysis identified IL-1β and LDHA as potential intervention points. Mathematical models described here can be useful to understand the complex interrelationship between metabolism, inflammation and circadian rhythm, and in designing effective therapeutic strategies. Our findings underscore the importance of including the circadian clock when developing pharmacological approaches that target aberrant tumour metabolism and inflammation.
Insight box
The complex interplay of metabolism-inflammation-circadian rhythm in tumours is not well understood. Our recent findings provided evidence of a feed-forward loop between tumour metabolite lactate, pro-inflammatory cytokine IL-1β and circadian CLOCK/BMAL1 in glioma. To elucidate the implication of this complex interplay, we developed a mathematical model that quantitatively describes this LDHA-IL-1β-CLOCK/BMAL1 circuit and integrates experimental data to predict potential therapeutic targets. The study employed a multi-start optimization strategy and profile likelihood estimations for parameter estimation and assessing identifiability. The simulations are in reasonable agreement with the experimental data. Sensitivity analysis found LDHA and IL-1β as potential therapeutic points. Mathematical models described here can provide insights to understand the complex interrelationship between metabolism, inflammation and circadian rhythm, and in identifying effective therapeutic targets.
Syncytial cells contain multiple nuclei and have local distribution and function of cellular components despite being synthesized in a common cytoplasm. The syncytial Drosophila blastoderm embryo ...shows reduced spread of organelle and plasma membrane-associated proteins between adjacent nucleo-cytoplasmic domains. Anchoring to the cytoarchitecture within a nucleo-cytoplasmic domain is likely to decrease the spread of molecules; however, its role in restricting this spread has not been assessed. In order to analyze the cellular mechanisms that regulate the rate of spread of plasma membrane-associated molecules in the syncytial Drosophila embryos, we express a pleckstrin homology (PH) domain in a localized manner at the anterior of the embryo by tagging it with the bicoid mRNA localization signal. Anteriorly expressed PH domain forms an exponential gradient in the anteroposterior axis with a longer length scale compared with Bicoid. Using a combination of experiments and theoretical modeling, we find that the characteristic distribution and length scale emerge due to plasma membrane sequestration and restriction within an energid. Loss of plasma membrane remodeling to form pseudocleavage furrows shows an enhanced spread of PH domain but not Bicoid. Modeling analysis suggests that the enhanced spread of the PH domain occurs due to the increased spread of the cytoplasmic population of the PH domain in pseudocleavage furrow mutants. Our analysis of cytoarchitecture interaction in regulating plasma membrane protein distribution and constraining its spread has implications on the mechanisms of spread of various molecules, such as morphogens in syncytial cells.
We present a theory for polyelectrolyte gels that allow the effective charge of the polymer backbone to self-regulate. Using a variational approach, we obtain an expression for the free energy of ...gels that accounts for the gel elasticity, free energy of mixing, counterion adsorption, local dielectric constant, electrostatic interaction among polymer segments, electrolyte ion correlations, and self-consistent charge regularization on the polymer strands. This free energy is then minimized to predict the behavior of the system as characterized by the gel volume fraction as a function of external variables such as temperature and salt concentration. We present results for the volume transition of polyelectrolyte gels in salt-free solvents, solvents with monovalent salts, and solvents with divalent salts. The results of our theoretical analysis capture the essential features of existing experimental results and also provide predictions for further experimentation. Our analysis highlights the importance of the self-regularization of the effective charge for the volume transition of gels in particular, and for charged polymer systems in general. Our analysis also enables us to identify the dominant free energy contributions for charged polymer networks and provides a framework for further investigation of specific experimental systems.
We report a two-layer microfluidic device to study the combined effect of confinement and chemical gradient on the motility of wild-type E. coli . We track individual E. coli in 50 μm and 10 μm wide ...microchannels, with a channel height of 2 μm, to generate quasi-2D conditions. We find that contrary to expectations, bacterial trajectories are superdiffusive even in the absence of a chemical (glucose) gradient. The superdiffusive behaviour becomes more pronounced upon introducing a chemical gradient or strengthening the lateral confinement. Run length distributions for weak lateral confinement in the absence of chemical gradients follow an exponential distribution. Both confinement and chemoattraction induce deviations from this behaviour, with the run length distributions approaching a power-law form under these conditions. Both confinement and chemoattraction suppress large-angle tumbles as well. Our results suggest that wild-type E. coli modulates both its runs and tumbles in a similar manner under physical confinement and chemical gradient. Our findings have implications for understanding how bacteria modulate their motility behaviour in natural habitats.