Role of diffusion and reaction of the constituents in spreading of histone modification marks.

Cells switch genes ON or OFF by altering the state of chromatin via histone modifications at specific regulatory locations along the chromatin polymer. These gene regulation processes are carried out by a network of reactions in which the histone marks spread to neighboring regions with the help of enzymes. In the literature, this spreading has been studied as a purely kinetic, non-diffusive process considering the interactions between neighboring nucleosomes. In this work, we go beyond this framework and study the spreading of modifications using a reaction-diffusion (RD) model accounting for the diffusion of the constituents. We quantitatively segregate the modification profiles generated from kinetic and RD models. The diffusion and degradation of enzymes set a natural length scale for limiting the domain size of modification spreading, and the resulting enzyme limitation is inherent in our model. We also demonstrate the emergence of confined modification domains without the explicit requirement of a nucleation site. We explore polymer compaction effects on spreading and show that single-cell domains may differ from averaged profiles. We find that the modification profiles from our model are comparable with existing H3K9me3 data of S. pombe.

Nonequilibrium switching of segmental states can influence compaction of chromatin.

Knowledge about the dynamic nature of chromatin organization is essential to understand the regulation of processes like DNA transcription and repair. The existing models of chromatin assume that protein organization and chemical states along chromatin are static and the 3D organization is purely a result of protein-mediated intra-chromatin interactions. Here we present a new hypothesis that certain nonequilibrium processes, such as switching of chemical and physical states due to nucleosome assembly/disassembly or gene repression/activation, can also simultaneously influence chromatin configurations. To understand the implications of this inherent nonequilibrium switching, we present a block copolymer model of chromatin, with switching of its segmental states between two states, mimicking active/repressed or protein unbound/bound states. We show that competition between switching timescale Tt, polymer relaxation timescale τp, and segmental relaxation timescale τs can lead to non-trivial changes in chromatin organization, leading to changes in local compaction and contact probabilities. As a function of the switching timescale, the radius of gyration of chromatin shows a non-monotonic behavior with a prominent minimum when Tt ≈ τp and a maximum when Tt ≈ τs. We find that polymers with a small segment length exhibit a more compact structure than those with larger segment lengths. We also find that the switching can lead to higher contact probability and better mixing of far-away segments. Our study also shows that the nature of the distribution of chromatin clusters varies widely as we change the switching rate.

Intermolecular interactions underlie protein/peptide phase separation irrespective of sequence and structure at crowded milieu.

Liquid-liquid phase separation (LLPS) has emerged as a crucial biological phenomenon underlying the sequestration of macromolecules (such as proteins and nucleic acids) into membraneless organelles in cells. Unstructured and intrinsically disordered domains are known to facilitate multivalent interactions driving protein LLPS. We hypothesized that LLPS could be an intrinsic property of proteins/polypeptides but with distinct phase regimes irrespective of their sequence and structure. To examine this, we studied many (a total of 23) proteins/polypeptides with different structures and sequences for LLPS study in the presence and absence of molecular crowder, polyethylene glycol (PEG-8000). We showed that all proteins and even highly charged polypeptides (under study) can undergo liquid condensate formation, however with different phase regimes and intermolecular interactions. We further demonstrated that electrostatic, hydrophobic, and H-bonding or a combination of such intermolecular interactions plays a crucial role in individual protein/peptide LLPS.

Predicting scale-dependent chromatin polymer properties from systematic coarse-graining.

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.

SUMO1 hinders α-Synuclein fibrillation by inducing structural compaction.

Small Ubiquitin-like Modifier 1 (SUMO1) is an essential protein for many cellular functions, including regulation, signaling, etc., achieved by a process known as SUMOylation, which involves covalent attachment of SUMO1 to target proteins. SUMO1 also regulates the function of several proteins via non-covalent interactions involving the hydrophobic patch in the target protein identified as SUMO Binding or Interacting Motif (SBM/SIM). Here, we demonstrate a crucial functional potential of SUMO1 mediated by its non-covalent interactions with α-Synuclein, a protein responsible for many neurodegenerative diseases called α-Synucleinopathies. SUMO1 hinders the fibrillation of α-Synuclein, an intrinsically disordered protein (IDP) that undergoes a transition to ß-structures during the fibrillation process. Using a plethora of biophysical techniques, we show that SUMO1 transiently binds to the N-terminus region of α-Synuclein non-covalently and causes structural compaction, which hinders the self-association process and thereby delays the fibrillation process. On the one hand, this study demonstrates an essential functional role of SUMO1 protein concerning neurodegeneration; it also illustrates the commonly stated mechanism that IDPs carry out multiple functions by structural adaptation to suit specific target proteins, on the other. Residue-level details about the SUMO1-α-Synuclein interaction obtained here also serve as a reliable approach for investigating the detailed mechanisms of IDP functions.

Disruption of polyhomeotic polymerization decreases nucleosome occupancy and alters genome accessibility.

Chromatin attains its three-dimensional (3D) conformation by establishing contacts between different noncontiguous regions. Sterile Alpha Motif (SAM)-mediated polymerization of the polyhomeotic (PH) protein regulates subnuclear clustering of Polycomb Repressive Complex 1 (PRC1) and chromatin topology. The mutations that perturb the ability of the PH to polymerize, disrupt long-range chromatin contacts, alter Hox gene expression, and lead to developmental defects. To understand the underlying mechanism, we combined the experiments and theory to investigate the effect of this SAM domain mutation on nucleosome occupancy and accessibility on a genome wide scale. Our data show that disruption of PH polymerization because of SAM domain mutation decreases nucleosome occupancy and alters accessibility. Polymer simulations investigating the interplay between distant chromatin contacts and nucleosome occupancy, both of which are regulated by PH polymerization, suggest that nucleosome density increases when contacts between different regions of chromatin are established. Taken together, it appears that SAM domain-mediated PH polymerization biomechanically regulates the organization of chromatin at multiple scales from nucleosomes to chromosomes and we suggest that higher order organization can have a top-down causation effect on nucleosome occupancy.

Nucleosome sliding can influence the spreading of histone modifications.

Nucleosomes are the fundamental building blocks of chromatin that not only help in the folding of chromatin, but also in carrying epigenetic information. It is known that nucleosome sliding is responsible for dynamically organizing chromatin structure and the resulting gene regulation. Since sliding can move two neighboring nucleosomes physically close or away, can it play a role in the spreading of histone modifications? We investigate this by simulating a stochastic model that couples nucleosome dynamics with the kinetics of histone modifications. We show that the sliding of nucleosomes can affect the modification pattern as well as the time it takes to modify a given region of chromatin. Exploring different nucleosome densities and modification kinetic parameters, we show that nucleosome sliding can be important for creating histone modification domains. Our model predicts that nucleosome density coupled with sliding dynamics can create an asymmetric histone modification profile around regulatory regions. We also compute the probability distribution of modified nucleosomes and relaxation kinetics of modifications. Our predictions are comparable with known experimental results.

Direct Demonstration of Seed Size-Dependent α-Synuclein Amyloid Amplification.

The size of amyloid seeds is known to modulate their autocatalytic amplification and cellular toxicity. However, the seed size-dependent secondary nucleation mechanism, toxicity, and disease-associated biological processes mediated by α-synuclein (α-Syn) fibrils are largely unknown. Using the cellular model and in vitro reconstitution, we showed that the size of α-Syn fibril seeds dictates not only their cellular internalization and associated cell death but also the distinct mechanisms of fibril amplification pathways involved in the pathological conformational change of α-Syn. Specifically, small fibril seeds showed elongation possibly through monomer addition at the fibril termini, whereas longer fibrils template the fibril amplification by surface-mediated nucleation as demonstrated by super-resolution microscopy. The distinct mechanism of fibril amplification and cellular uptake along with toxicity suggest that breakage of fibrils into seeds of different sizes determines the underlying pathological outcome of synucleinopathies.

α-Synuclein Aggregation Intermediates form Fibril Polymorphs with Distinct Prion-like Properties.

α-Synuclein (α-Syn) amyloids in synucleinopathies are suggested to be structurally and functionally diverse, reminiscent of prion-like strains. The mechanism of how the aggregation of the same precursor protein results in the formation of fibril polymorphs remains elusive. Here, we demonstrate the structure-function relationship of two polymorphs, pre-matured fibrils (PMFs) and helix-matured fibrils (HMFs), based on α-Syn aggregation intermediates. These polymorphs display the structural differences as demonstrated by solid-state NMR and mass spectrometry studies and also possess different cellular activities such as seeding, internalization, and cell-to-cell transfer of aggregates. HMFs, with a compact core structure, exhibit low seeding potency but readily internalize and transfer from one cell to another. The less structured PMFs lack transcellular transfer ability but induce abundant α-Syn pathology and trigger the formation of aggresomes in cells. Overall, the study highlights that the conformational heterogeneity in the aggregation pathway may lead to fibril polymorphs with distinct prion-like behavior.

Heterogeneous interactions and polymer entropy decide organization and dynamics of chromatin domains.

Chromatin is known to be organized into multiple domains of varying sizes and compaction. While these domains are often imagined as static structures, they are highly dynamic and show cell-to-cell variability. Since processes such as gene regulation and DNA replication occur in the context of these domains, it is important to understand their organization, fluctuation, and dynamics. To simulate chromatin domains, one requires knowledge of interaction strengths among chromatin segments. Here, we derive interaction-strength parameters from experimentally known contact maps and use them to predict chromatin organization and dynamics. Taking two domains on the human chromosome as examples, we investigate its three-dimensional organization, size/shape fluctuations, and dynamics of different segments within a domain, accounting for hydrodynamic effects. Considering different cell types, we quantify changes in interaction strengths and chromatin shape fluctuations in different epigenetic states. Perturbing the interaction strengths systematically, we further investigate how epigenetic-like changes can alter the spatio-temporal nature of the domains. Our results show that heterogeneous weak interactions are crucial in determining the organization of the domains. Computing effective stiffness and relaxation times, we investigate how perturbations in interactions affect the solid- and liquid-like nature of chromatin domains. Quantifying dynamics of chromatin segments within a domain, we show how the competition between polymer entropy and interaction energy influence the timescales of loop formation and maintenance of stable loops.

Role of non-specific interactions in the phase-separation and maturation of macromolecules.

Phase separation of biomolecules could be mediated by both specific and non-specific interactions. How the interplay between non-specific and specific interactions along with polymer entropy influences phase separation is an open question. We address this question by simulating self-associating molecules as polymer chains with a short core stretch that forms the specifically interacting functional interface and longer non-core regions that participate in non-specific/promiscuous interactions. Our results show that the interplay of specific (strength, ϵsp) and non-specific interactions (strength, ϵns) could result in phase separation of polymers and its transition to solid-like aggregates (mature state). In the absence of ϵns, the polymer chains do not dwell long enough in the vicinity of each other to undergo phase separation and transition into a mature state. On the other hand, in the limit of strong ϵns, the assemblies cannot transition into the mature state and form a non-specific assembly, suggesting an optimal range of interactions favoring mature multimers. In the scenario where only a fraction (Nfrac) of the non-core regions participate in attractive interactions, we find that slight modifications to either ϵns or Nfrac can result in dramatically altered self-assembled states. Using a combination of heterogeneous and homogeneous mix of polymers, we establish how this interplay between interaction energies dictates the propensity of biomolecules to find the correct binding partner at dilute concentrations in crowded environments.

High fidelity epigenetic inheritance: Information theoretic model predicts threshold filling of histone modifications post replication.

During cell devision, maintaining the epigenetic information encoded in histone modification patterns is crucial for survival and identity of cells. The faithful inheritance of the histone marks from the parental to the daughter strands is a puzzle, given that each strand gets only half of the parental nucleosomes. Mapping DNA replication and reconstruction of modifications to equivalent problems in communication of information, we ask how well enzymes can recover the parental modifications, if they were ideal computing machines. Studying a parameter regime where realistic enzymes can function, our analysis predicts that enzymes may implement a critical threshold filling algorithm which fills unmodified regions of length at most k. This algorithm, motivated from communication theory, is derived from the maximum à posteriori probability (MAP) decoding which identifies the most probable modification sequence based on available observations. Simulations using our method produce modification patterns similar to what has been observed in recent experiments. We also show that our results can be naturally extended to explain inheritance of spatially distinct antagonistic modifications.

Continuous variable responses and signal gating form kinetic bases for pulsatile insulin signaling and emergence of resistance.

Understanding kinetic control of biological processes is as important as identifying components that constitute pathways. Insulin signaling is central for almost all metazoans, and its perturbations are associated with various developmental disorders, metabolic diseases, and aging. While temporal phosphorylation changes and kinetic constants have provided some insights, constant or variable parameters that establish and maintain signal topology are poorly understood. Here, we report kinetic parameters that encode insulin concentration and nutrient-dependent flow of information using iterative experimental and mathematical simulation-based approaches. Our results illustrate how dynamics of distinct phosphorylation events collectively contribute to selective kinetic gating of signals and maximum connectivity of the signaling cascade under normo-insulinemic but not hyper-insulinemic states. In addition to identifying parameters that provide predictive value for maintaining the balance between metabolic and growth-factor arms, we posit a kinetic basis for the emergence of insulin resistance. Given that pulsatile insulin secretion during a fasted state precedes a fed response, our findings reveal rewiring of insulin signaling akin to memory and anticipation, which was hitherto unknown. Striking disparate temporal behavior of key phosphorylation events that destroy the topology under hyper-insulinemic states underscores the importance of unraveling regulatory components that act as bandwidth filters. In conclusion, besides providing fundamental insights, our study will help in identifying therapeutic strategies that conserve coupling between metabolic and growth-factor arms, which is lost in diseases and conditions of hyper-insulinemia.

Orc4 spatiotemporally stabilizes centromeric chromatin.

The establishment of centromeric chromatin and its propagation by the centromere-specific histone CENPA is mediated by epigenetic mechanisms in most eukaryotes. DNA replication origins, origin binding proteins, and replication timing of centromere DNA are important determinants of centromere function. The epigenetically regulated regional centromeres in the budding yeast Candida albicans have unique DNA sequences that replicate earliest in every chromosome and are clustered throughout the cell cycle. In this study, the genome-wide occupancy of the replication initiation protein Orc4 reveals its abundance at all centromeres in C. albicans Orc4 is associated with four different DNA sequence motifs, one of which coincides with tRNA genes (tDNA) that replicate early and cluster together in space. Hi-C combined with genome-wide replication timing analyses identify that early replicating Orc4-bound regions interact with themselves stronger than with late replicating Orc4-bound regions. We simulate a polymer model of chromosomes of C. albicans and propose that the early replicating and highly enriched Orc4-bound sites preferentially localize around the clustered kinetochores. We also observe that Orc4 is constitutively localized to centromeres, and both Orc4 and the helicase Mcm2 are essential for cell viability and CENPA stability in C. albicans Finally, we show that new molecules of CENPA are recruited to centromeres during late anaphase/telophase, which coincides with the stage at which the CENPA-specific chaperone Scm3 localizes to the kinetochore. We propose that the spatiotemporal localization of Orc4 within the nucleus, in collaboration with Mcm2 and Scm3, maintains centromeric chromatin stability and CENPA recruitment in C. albicans.

The Accidental Ally: Nucleosome Barriers Can Accelerate Cohesin-Mediated Loop Formation in Chromatin.

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.

Computational Model for Studying Breakage-Dependent Amyloid Growth.

Amyloid fibrils are typically associated with neurodegenerative diseases. Recent studies have suggested that, similar to prions, many amyloid proteins are infectious in nature and may cause spreading and dissemination of diseases. Typical amyloid infection propagates by recruiting functional proteins into amyloidogenic form and multiplying by breaking the existing fibril. In this study, we model the kinetics of fibril growth through breakage and the subsequent elongation process, similar to the prion infection process. Using kinetic Monte Carlo simulations as well as mathematical counting methods, we show how the measurable quantities like the 50% aggregation time (T50) and the maximum growth rate (Vmax) scale with various parameters in the problem. This study has a direct application where it can be used to understand experiments that amplify the minute amount of amyloid seeds present in biological fluid for early detection of human disease. Using the knowledge from our simulations, we can predict the initial seed concentration, known as the filament kinetics.

Nucleosome positioning and chromatin organization.

In our cells, DNA is folded and packed with the help of many proteins into chromatin whose basic unit is a nucleosome-DNA wrapped around octamer of histone proteins. The chain of nucleosomes is further folded and arranged into many layers and has a dynamic organization. How does the complex chromatin organization emerge from interactions among DNA, histones, and non-histone proteins have been a question of great interest. Here we review recent literature that investigated how nucleosome positioning and nucleosome-mediated interactions drive chromatin organization. Unlike our earlier understanding, chromatin is organized into 3D domains of various sizes having irregularly organized nucleosomes. These domains emerge due to heterogeneous nucleosome positioning and diverse inter-nucleosome interactions that vary in space and time.

α-Synuclein aggregation nucleates through liquid-liquid phase separation.

α-Synuclein (α-Syn) aggregation and amyloid formation is directly linked with Parkinson's disease pathogenesis. However, the early events involved in this process remain unclear. Here, using the in vitro reconstitution and cellular model, we show that liquid-liquid phase separation of α-Syn precedes its aggregation. In particular, in vitro generated α-Syn liquid-like droplets eventually undergo a liquid-to-solid transition and form an amyloid hydrogel that contains oligomers and fibrillar species. Factors known to aggravate α-Syn aggregation, such as low pH, phosphomimetic substitution and familial Parkinson's disease mutations, also promote α-Syn liquid-liquid phase separation and its subsequent maturation. We further demonstrate α-Syn liquid-droplet formation in cells. These cellular α-Syn droplets eventually transform into perinuclear aggresomes, the process regulated by microtubules. This work provides detailed insights into the phase-separation behaviour of natively unstructured α-Syn and its conversion to a disease-associated aggregated state, which is highly relevant in Parkinson's disease pathogenesis.

Computing 3D Chromatin Configurations from Contact Probability Maps by Inverse Brownian Dynamics.

The three-dimensional (3D) organization of chromatin, on the length scale of a few genes, is crucial in determining the functional state-accessibility and amount of gene expression-of the chromatin. Recent advances in chromosome conformation capture experiments provide partial information on the chromatin organization in a cell population, namely the contact count between any segment pairs, but not on the interaction strength that leads to these contact counts. However, given the contact matrix, determining the complete 3D organization of the whole chromatin polymer is an inverse problem. In this work, a novel inverse Brownian dynamics method based on a coarse-grained bead-spring chain model has been proposed to compute the optimal interaction strengths between different segments of chromatin such that the experimentally measured contact count probability constraints are satisfied. Applying this method to the α-globin gene locus in two different cell types, we predict the 3D organizations corresponding to active and repressed states of chromatin at the locus. We show that the average distance between any two segments of the region has a broad distribution and cannot be computed as a simple inverse relation based on the contact probability alone. The results presented for multiple normalization methods suggest that all measurable quantities may crucially depend on the nature of normalization. We argue that by experimentally measuring predicted quantities, one may infer the appropriate form of normalization.