Physical modeling of nucleosome clustering in euchromatin resulting from interactions between epigenetic reader proteins.

Euchromatin is an accessible phase of genetic material containing genes that encode proteins with increased expression levels. The structure of euchromatin in vitro has been described as a 30-nm fiber formed from ordered nucleosome arrays. However, recent advances in microscopy have revealed an in vivo euchromatin architecture that is much more disordered, characterized by variable-length linker DNA and sporadic nucleosome clusters. In this work, we develop a theoretical model to elucidate factors contributing to the disordered in vivo architecture of euchromatin. We begin by developing a 1D model of nucleosome positioning that captures the interactions between bound epigenetic reader proteins to predict the distribution of DNA linker lengths between adjacent nucleosomes. We then use the predicted linker lengths to construct 3D chromatin configurations consistent with the physical properties of DNA within the nucleosome array, and we evaluate the distribution of nucleosome cluster sizes in those configurations. Our model reproduces experimental cluster-size distributions, which are dramatically influenced by the local pattern of epigenetic marks and the concentration of reader proteins. Based on our model, we attribute the disordered arrangement of euchromatin to the heterogeneous binding of reader proteins and subsequent short-range interactions between bound reader proteins on adjacent nucleosomes. By replicating experimental results with our physics-based model, we propose a mechanism for euchromatin organization in the nucleus that impacts gene regulation and the maintenance of epigenetic marks.

Closing the loop between microstructure and charge transport in conjugated polymers by combining microscopy and simulation.

A grand challenge in materials science is to identify the impact of molecular composition and structure across a range of length scales on macroscopic properties. We demonstrate a unified experimental-theoretical framework that coordinates experimental measurements of mesoscale structure with molecular-level physical modeling to bridge multiple scales of physical behavior. Here we apply this framework to understand charge transport in a semiconducting polymer. Spatially-resolved nanodiffraction in a transmission electron microscope is combined with a self-consistent framework of the polymer chain statistics to yield a detailed picture of the polymer microstructure ranging from the molecular to device relevant scale. Using these data as inputs for charge transport calculations, the combined multiscale approach highlights the underrepresented role of defects in existing transport models. Short-range transport is shown to be more chaotic than is often pictured, with the drift velocity accounting for a small portion of overall charge motion. Local transport is sensitive to the alignment and geometry of polymer chains. At longer length scales, large domains and gradual grain boundaries funnel charges preferentially to certain regions, creating inhomogeneous charge distributions. While alignment generally improves mobility, these funneling effects negatively impact mobility. The microstructure is modified in silico to explore possible design rules, showing chain stiffness and alignment to be beneficial while local homogeneity has no positive effect. This combined approach creates a flexible and extensible pipeline for analyzing multiscale functional properties and a general strategy for extending the accesible length scales of experimental and theoretical probes by harnessing their combined strengths.

Diffusion and distal linkages govern interchromosomal dynamics during meiotic prophase.

SignificanceEssential for sexual reproduction, meiosis is a specialized cell division required for the production of haploid gametes. Critical to this process are the pairing, recombination, and segregation of homologous chromosomes (homologs). While pairing and recombination are linked, it is not known how many linkages are sufficient to hold homologs in proximity. Here, we reveal that random diffusion and the placement of a small number of linkages are sufficient to establish the apparent "pairing" of homologs. We also show that colocalization between any two loci is more dynamic than anticipated. Our study provides observations of live interchromosomal dynamics during meiosis and illustrates the power of combining single-cell measurements with theoretical polymer modeling.

Local changes in protein filament properties drive large-scale membrane transformations involved in endosome tethering and fusion.

Large-scale cellular transformations are triggered by subtle physical and structural changes to individual biomacromolecular and membrane components. A prototypical example of such an event is the orchestrated fusion of membranes within an endosome that enables transport of cargo and processing of biochemical moieties. In this work, we demonstrate how protein filaments on the endosomal membrane surface can leverage a rigid-to-flexible transformation to elicit a large-scale change in membrane flexibility to enable membrane fusion. We develop a polymer field-theoretic model that captures molecular alignment arising from nematic interactions with varying surface density and fraction of flexible filaments, which are biologically controlled within the endosomal membrane. We then predict the collective elasticity of the filament brush in response to changes in the filament alignment, predicting a greater than 20-fold increase of the effective membrane elasticity over the bare membrane elasticity that is triggered by filament alignment. These results show that the endosome can modulate the filament properties to orchestrate membrane fluidization that facilitates vesicle fusion, providing an example of how active processes that modulate local molecular properties can result in large-scale transformations that are essential to cellular survival.

Effect of local active fluctuations on structure and dynamics of flexible biopolymers.

Active fluctuations play a significant role in the structure and dynamics of biopolymers (e.g. chromatin and cytoskeletal proteins) that are instrumental in the functioning of living cells. For a large range of experimentally accessible length and time scales, these polymers can be represented as flexible chains that are subjected to spatially and temporally varying fluctuating forces. In this work, we introduce a mathematical framework that correlates the spatial and temporal patterns of the fluctuations to different observables that describe the dynamics and conformations of the polymer. We demonstrate the power of this approach by analyzing the case of a point fluctuation on the polymer with an exponential decay of correlation in time with a finite time constant. Specifically, we identify the length and time scale over which the behavior of the polymer exhibits a significant departure from the behavior of a Rouse chain and the range of impact of the fluctuation along the chain. Furthermore, we show that the conformation of the polymer retains the memory of the active fluctuation from earlier times. Altogether, this work sets the basis for understanding and interpreting the role of spatio-temporal patterns of fluctuations in the dynamics, conformation, and functionality of biopolymers in living cells.

Leveraging polymer modeling to reconstruct chromatin connectivity from live images.

Chromosomal dynamics plays a central role in a number of critical biological processes, such as transcriptional regulation, genetic recombination, and DNA replication. However, visualization of chromatin is generally limited to live imaging of a few fluorescently labeled chromosomal loci or high-resolution reconstruction of multiple loci from a single time frame. To aid in mapping the underlying chromosomal structure based on parsimonious experimental measurements, we present an exact analytical expression for the evolution of the polymer configuration based on a flexible-polymer model, and we propose an algorithm that tracks the polymer configuration from live images of chromatin marked with several fluorescent marks. Our theory identifies the resolution of microscopy needed to achieve high-accuracy tracking for a given spacing of markers, establishing the statistical confidence in the assignment of genome identity to the visualized marks. We then leverage experimental data of locus-tracking measurements to demonstrate the validity of our modeling approach and to establish a basis for the design of experiments with a desired resolution. Altogether, this work provides a computational approach founded on polymer physics that vastly improves the interpretation of in vivo measurements of biopolymer dynamics.

Physical modeling of the heritability and maintenance of epigenetic modifications.

We develop a predictive theoretical model of the physical mechanisms that govern the heritability and maintenance of epigenetic modifications. This model focuses on a particular modification, methylation of lysine-9 of histone H3 (H3K9), which is one of the most representative and critical epigenetic marks that affects chromatin organization and gene expression. Our model combines the effect of segregation and compaction on chromosomal organization with the effect of the interaction between proteins that compact the chromatin (heterochromatin protein 1) and the methyltransferases that affect methyl spreading. Our chromatin model demonstrates that a block of H3K9 methylations in the epigenetic sequence determines the compaction state at any particular location in the chromatin. Using our predictive model for chromatin compaction, we develop a methylation model to address the reestablishment of the methylation sequence following DNA replication. Our model reliably maintains methylation over generations, thereby establishing the robustness of the epigenetic code.

Coarse-grained modeling reveals the impact of supercoiling and loop length in DNA looping kinetics.

Measurements of protein-mediated DNA looping reveal that in vivo conditions favor the formation of loops shorter than those that occur in vitro, yet the precise physical mechanisms underlying this shift remain unclear. To understand the extent to which in vivo supercoiling may explain these shifts, we develop a theoretical model based on coarse-grained molecular simulation and analytical transition state theory, enabling us to map out looping energetics and kinetics as a function of two key biophysical parameters: superhelical density and loop length. We show that loops on the scale of a persistence length respond to supercoiling over a much wider range of superhelical densities and to a larger extent than longer loops. This effect arises from a tendency for loops to be centered on the plectonemic end region, which bends progressively more tightly with superhelical density. This trend reveals a mechanism by which supercoiling favors shorter loop lengths. In addition, our model predicts a complex kinetic response to supercoiling for a given loop length, governed by a competition between an enhanced rate of looping due to torsional buckling and a reduction in looping rate due to chain straightening as the plectoneme tightens at higher superhelical densities. Together, these effects lead to a flattening of the kinetic response to supercoiling within the physiological range for all but the shortest loops. Using experimental estimates for in vivo superhelical densities, we discuss our model's ability to explain available looping data, highlighting both the importance of supercoiling as a regulatory force in genetics and the additional complexities of looping phenomena in vivo.

Active and thermal fluctuations in multi-scale polymer structure and dynamics.

The presence of athermal noise or biological fluctuations control and maintain crucial life-processes. In this work, we present an exact analytical treatment of the dynamic behavior of a flexible polymer chain that is subjected to both thermal and active forces. Our model for active forces incorporates temporal correlation associated with the characteristic time scale and processivity of enzymatic function (driven by ATP hydrolysis), leading to an active-force time scale that competes with relaxation processes within the polymer chain. We analyze the structure and dynamics of an active-Brownian polymer using our exact results for the dynamic structure factor and the looping time for the chain ends. The spectrum of relaxation times within a polymer chain implies two different behaviors at small and large length scales. Small length-scale relaxation is faster than the active-force time scale, and the dynamic and structural behavior at these scales are oblivious to active forces and, are thus governed by the true thermal temperature. Large length-scale behavior is governed by relaxation times that are much longer than the active-force time scale, resulting in an effective active-Brownian temperature that dramatically alters structural and dynamic behavior. These complex multi-scale effects imply a time-dependent temperature that governs living and non-equilibrium systems, serving as a unifying concept for interpreting and predicting their physical behavior.

Semiflexible polymer solutions. II. Fluctuations and Frank elastic constants.

We study the collective elastic behavior of semiflexible polymer solutions in a nematic liquid-crystalline state using polymer field theory. Our polymer field-theoretic model of semiflexible polymer solutions is extended to include second-order fluctuation corrections to the free energy, permitting the evaluation of the Frank elastic constants based on orientational order fluctuations in the nematic state. Our exact treatment of wormlike chain statistics permits the evaluation of behavior from the nematic state, thus accurately capturing the impact of single-chain behavior on collective elastic response. Results for the Frank elastic constants are presented as a function of aligning field strength and chain length, and we explore the impact of conformation fluctuations and hairpin defects on the twist, splay, and bend moduli. Our results indicate that the twist elastic constant Ktwist is smaller than both bend and splay constants (Kbend and Ksplay, respectively) for the entire range of polymer rigidity. Splay and bend elastic constants exhibit regimes of dominance over the range of chain stiffness, where Ksplay > Kbend for flexible polymers (large-N limit) while the opposite is true for rigid polymers. Theoretical analysis also suggests the splay modulus tracks exactly to that of the end-to-end distance in the transverse direction for semiflexible polymers at intermediate to large-N. These results provide insight into the role of conformation fluctuations and hairpin defects on the collective response of polymer solutions.

Impact of chromosomal organization on epigenetic drift and domain stability revealed by physics-based simulations.

We examine the relationship between the size of domains of epigenetic marks and the stability of those domains using our theoretical model that captures the physical mechanisms governing the maintenance of epigenetic modifications. We focus our study on histone H3 lysine-9 trimethylation, one of the most common and consequential epigenetic marks with roles in chromatin compaction and gene repression. Our model combines the effects of methyl spreading by methyltransferases and chromatin segregation into heterochromatin and euchromatin because of preferential heterochromatin protein 1 (HP1) binding. Our model indicates that, although large methylated domains are passed successfully from one chromatin generation to the next, small alterations to the methylation sequence are not maintained during chromatin replication. Using our predictive model, we investigate the size required for an epigenetic domain to persist over chromatin generations while surrounded by a much larger domain of opposite methylation and compaction state. We find that there is a critical size threshold in the hundreds-of-nucleosomes scale above which an epigenetic domain will be reliably maintained over generations. The precise size of the threshold differs for heterochromatic and euchromatic domains. Our results are consistent with natural alterations to the epigenetic sequence occurring during embryonic development and due to age-related epigenetic drift.

A fluorogenic array for temporally unlimited single-molecule tracking.

We describe three optical tags, ArrayG, ArrayD and ArrayG/N, for intracellular tracking of single molecules over milliseconds to hours. ArrayG is a fluorogenic tag composed of a green fluorescent protein-nanobody array and monomeric wild-type green fluorescent protein binders that are initially dim but brighten ~26-fold on binding with the array. By balancing the rates of binder production, photobleaching and stochastic binder exchange, we achieve temporally unlimited tracking of single molecules. High-speed tracking of ArrayG-tagged kinesins and integrins for thousands of frames reveals novel dynamical features. Tracking of single histones at 0.5 Hz for >1 hour with the import competent ArrayG/N tag shows that chromosomal loci behave as Rouse polymers with visco-elastic memory and exhibit a non-Gaussian displacement distribution. ArrayD, based on a dihydrofolate reductase nanobody array and dihydrofolate reductase-fluorophore binder, enables dual-color imaging. The arrays combine brightness, fluorogenicity, fluorescence replenishment and extended fluorophore choice, opening new avenues for tracking single molecules in living cells.

Dynamic light scattering microrheology for soft and living materials.

We present a method for using dynamic light scattering in the single-scattering limit to measure the viscoelastic moduli of soft materials. This microrheology technique only requires a small sample volume of 12 µL to measure up to six decades in time of rheological behavior. We demonstrate the use of dynamic light scattering microrheology (DLSµR) on a variety of soft materials, including dilute polymer solutions, covalently-crosslinked polymer gels, and active, biological fluids. In this work, we detail the procedure for applying the technique to new materials and discuss the critical considerations for implementing the technique, including a custom analysis script for analyzing data output. We focus on the advantages of applying DLSµR to biologically relevant materials: breast cancer cells encapsulated in a collagen gel and cystic fibrosis sputum. DLSµR is an easy, efficient, and economical rheological technique that can guide the design of new polymeric materials and facilitate the understanding of the underlying physics governing behavior of naturally derived materials.

Bottom-up modeling of chromatin segregation due to epigenetic modifications.

We use a chromosome-scale simulation to show that the preferential binding of heterochromatin protein 1 (HP1) to regions high in histone methylation (specifically H3K9me3) results in phase segregation and reproduces features of the observed Hi-C contact map. Specifically, we perform Monte Carlo simulations with one computational bead per nucleosome and an H3K9me3 pattern based on published ChIP-seq signals. We implement a binding model in which HP1 preferentially binds to trimethylated histone tails and then oligomerizes to bridge together nucleosomes. We observe a phase reminiscent of heterochromatin-dense and high in H3K9me3-and another reminiscent of euchromatin-less dense and lacking H3K9me3. This segregation results in a plaid contact probability map that matches the general shape and position of published Hi-C data. Analysis suggests that a roughly 20-kb segment of H3K9me3 enrichment is required to drive segregation into the heterochromatic phase.

Chromatin Compaction Leads to a Preference for Peripheral Heterochromatin.

A layer of dense heterochromatin is found at the periphery of the nucleus. Because this peripheral heterochromatin functions as a repressive phase, mechanisms that relocate genes to the periphery play an important role in regulating transcription. Using Monte Carlo simulations, we show that an interaction that attracts euchromatin and heterochromatin equally to the nuclear envelope will still preferentially locate heterochromatin to the nuclear periphery. This observation considerably broadens the class of possible interactions that result in peripheral positioning to include boundary interactions that either weakly attract all chromatin or strongly bind to a randomly chosen 0.05% of nucleosomes. The key distinguishing feature of heterochromatin is its high chromatin density with respect to euchromatin. In our model, this densification is caused by heterochromatin protein 1's preferential binding to histone H3 tails with a methylated lysine at the ninth residue, a hallmark of heterochromatin. We find that a global rearrangement of chromatin to place heterochromatin at the nuclear periphery can be accomplished by attaching a small subset of loci, even if these loci are uncorrelated with heterochromatin. Hence, factors that densify chromatin determine which genomic regions condense to form peripheral heterochromatin.

Chromosome Structural Mechanics Dictates the Local Spreading of Epigenetic Marks.

We present a theoretical model that demonstrates the integral role chromosome organization and structural mechanics play in the spreading of histone modifications involved in epigenetic regulation. Our model shows that heterogeneous nucleosome positioning, and the resulting position-dependent mechanical properties, must be included to reproduce several qualitative features of experimental data of histone methylation spreading around an artificially induced "nucleation site." We show that our model recreates both the extent of spreading and the presence of a subdominant peak upstream of the transcription start site. Our model indicates that the spreading of epigenetic modifications is sensitive to heterogeneity in chromatin organization and the resulting variability in the chromatin's mechanical properties, suggesting that nucleosome spacing can directly control the conferral of epigenetic marks by modifying the structural mechanics of the chromosome. It further illustrates how the physical organization of the DNA polymer may play a significant role in re-establishing the epigenetic code upon cell division.

Geometrical Heterogeneity Dominates Thermal Fluctuations in Facilitating Chromatin Contacts.

Within a living cell, the myriad of proteins that bind DNA introduce heterogeneously spaced kinks into an otherwise semiflexible DNA double helix. To investigate the effects of heterogeneous nucleosome binding on chromatin organization, we extend the wormlike chain model to include statistically spaced, rigid kinks. On timescales where nucleosome positions are fixed, we find that the probability of chromatin loop formation can vary by up to six orders of magnitude between two sets of nucleosome positions drawn from the same distribution. On longer timescales, we show that continuous rerandomization due to nucleosome turnover results in chromatin tracing out an effective WLC with a dramatically smaller Kuhn length than bare DNA. Together, these observations demonstrate that nucleosome spacing acts as the primary source of the structural heterogeneity that dominates local and global chromatin organization.

Polymer physics across scales: Modeling the multiscale behavior of functional soft materials and biological systems.

Polymeric materials are ubiquitous in our daily lives, and they play a significant role in many technological applications. The general predictive framework for the behavior of soft polymeric materials can be divided into two vastly different approaches. Highly coarse-grained models capture polymers as flexible random walks, resulting in general predictions of physical behavior but lack chemical specificity. Detailed atomistic models contain molecular detail but are frequently computationally intractable for exhaustive materials discovery. In this perspective, we discuss theoretical models that successfully bridge these disparate approaches. We identify intermediate-scale physical models that are amenable to theoretical analyses while containing sufficient granular detail to capture a range of molecular-level processes. We then provide several problems in materials engineering and biological physics where multiscale physics is essential in their behavior.

Polymer Semiflexibility Induces Nonuniversal Phase Transitions in Diblock Copolymers.

The order-disorder phase transition and the associated phase diagrams of semiflexible diblock copolymers are investigated using the wormlike chain model, incorporating concentration fluctuations. The free energy up to quartic order in concentration fluctuations is developed with chain-rigidity-dependent coefficients, evaluated using our exact results for the wormlike chain model, and a one-loop renormalization treatment is used to account for fluctuation effects. The chain length N and the monomer aspect ratio α directly control the strength of immiscibility (defined by the Flory-Huggins parameter χ) at the order-disorder transition and the resulting microstructures at different chemical compositions f_{A}. When monomers are infinitely thin (i.e., large aspect ratio α), the finite chain length N lowers the χN at the phase transition. However, fluctuation effects become important when chains have a finite radius, and a decrease in the chain length N elevates the χN at the phase transition. Phase diagrams of diblock copolymers over a wide range of N and α are calculated based on our fluctuation theory. We find that both finite N and α enhance the stability of the lamellar phase above the order-disorder transition. Our results demonstrate that polymer semiflexibility plays a dramatic role in the phase behavior, even for large chain lengths (e.g., N≈100).

Active DNA Olympic Hydrogels Driven by Topoisomerase Activity.

Biological systems are equipped with a diverse repertoire of proteins that regulate DNA topology with precision that is beyond the reach of conventional polymer chemistry. Here, we harness the unique properties of topoisomerases to synthesize Olympic hydrogels formed by topologically interlinked DNA rings. Using dynamic light scattering microrheology to probe the viscoelasticity of DNA topological networks, we show that topoisomerase II enables the facile preparation of active, adenosine triphosphate-driven Olympic hydrogels that can be switched between liquid and solid states on demand. Our results provide a versatile system for engineering switchable topological materials that may be broadly leveraged to model the impact of topological constraints and active dynamics in the physics of chromosomes and other polymeric materials.