Exploring the Free Energy Landscape of Nucleosomes

Nucleosomes as allosteric scaffolds for genetic regulation

Nucleosomes are stable yet highly dynamic complexes exhibiting diverse types of motions, such as sliding, DNA unwrapping, and disassembly, encoding a landscape with a large number of metastable states. In this review, describing recent studies on these nucleosome structure changes, we propose that the nucleosome can be viewed as an ideal allosteric scaffold: regulated by effector molecules such as transcription factors and chromatin remodelers, the nucleosome controls the downstream gene activity. Binding of transcription factors to the nucleosome can enhance DNA unwrapping or slide the DNA, altering either the binding or the unbinding of other transcription factors to nearby sites. ATP-dependent chromatin remodelers induce a series of DNA deformations, which allosterically propagate throughout the nucleosome to induce DNA sliding or histone exchange.

Ensembles of Breathing Nucleosomes: A Computational Study

About three-fourths of the human DNA molecules are wrapped into nucleosomes, protein spools with DNA. Nucleosomes are highly dynamic, transiently exposing their DNA through spontaneous unspooling. Recent experiments allowed to observe the DNA of an ensemble of such breathing nucleosomes through x-ray diffraction with contrast matching between the solvent and the protein core. In this study, we calculate such an ensemble through a Monte Carlo simulation of a coarse-grained nucleosome model with sequence-dependent DNA mechanics. Our analysis gives detailed insights into the sequence dependence of nucleosome breathing observed in the experiment and allows us to determine the adsorption energy of the DNA bound to the protein core as a function of the ionic strength. Moreover, we predict the breathing behavior of other potentially interesting sequences and compare the findings to earlier related experiments.

"Direct" and "Indirect" Effects of Histone Modifications: Modulation of Sterical Bulk as a Novel Source of Functionality

The chromatin-regulatory principles of histone post-translational modifications (PTMs) are discussed with a focus on the potential alterations in chromatin functional state due to steric and mechanical constraints imposed by bulky histone modifications such as ubiquitin and SUMO. In the classical view, PTMs operate as recruitment platforms for histone "readers," and as determinants of chromatin array compaction. Alterations of histone charges by "small" chemical modifications (e.g., acetylation, phosphorylation) could regulate nucleosome spontaneous dynamics without globally affecting nucleosome structure. These fluctuations in nucleosome wrapping can be exploited by chromatin-processing machinery. In contrast, ubiquitin and SUMO are comparable in size to histones, and it seems logical that these PTMs could conflict with canonical nucleosome organization. An experimentally testable hypothesis that by adding sterical bulk these PTMs can robustly alter nucleosome primary structure is proposed. The model presented here stresses the diversity of mechanisms by which histone PTMs regulate chromatin dynamics, primary structure and, hence, functionality.

Mean first passage times reconstruct the slowest relaxations in potential energy landscapes of nanoclusters

Relaxation modes are the collective modes in which all probability deviations from equilibrium states decay with the same relaxation rates. In contrast, a first passage time is the required time for arriving for the first time from one state to another. In this paper, we discuss how and why the slowest relaxation rates of relaxation modes are reconstructed from the first passage times. As an illustrative model, we use a continuous-time Markov state model of vacancy diffusion in KCl nanoclusters. Using this model, we reveal that all characteristics of the relaxations in KCl nanoclusters come from the fact that they are hybrids of two kinetically different regions of the fast surface and slow bulk diffusions. The origin of the different diffusivities turns out to come from the heterogeneity of the activation energies on the potential energy landscapes. We also develop a stationary population method to compute the mean first passage times as mean times required for pair annihilations of particle-hole pairs, which enables us to obtain the symmetric results of relaxation rates under the exchange of the sinks and the sources. With this symmetric method, we finally show why the slowest relaxation times can be reconstructed from the mean first passage times.

A free-energy landscape for the glucagon-like peptide 1 receptor GLP1R

G-protein-coupled receptors (GPCRs) are among the most important receptors in human physiology and pathology. They serve as master regulators of numerous key processes and are involved in as well as cause debilitating diseases. Consequently, GPCRs are among the most attractive targets for drug design and pharmaceutical interventions (>30% of drugs on the market). The glucagon-like peptide 1 (GLP-1) hormone receptor GLP1R is closely involved in insulin secretion by pancreatic β-cells and constitutes a major druggable target for the development of anti-diabetes and obesity agents. GLP1R structure was recently solved, with ligands, allosteric modulators and as part of a complex with its cognate G protein. However, the translation of this structural data into structure/function understanding remains limited. The current study functionally characterizes GLP1R with special emphasis on ligand and cellular partner binding interactions and presents a free-energy landscape as well as a functional model of the activation cycle of GLP1R. Our results should facilitate a deeper understanding of the molecular mechanism underlying GLP1R activation, forming a basis for improved development of targeted therapeutics for diabetes and related disorders.

Folding profiles of antimicrobial scorpion venom-derived peptides on hydrophobic surfaces: a molecular dynamics study

Most helical antimicrobial peptides (AMPs) are usually unfolded in aqueous solution; however they acquire their secondary structure in the presence of a hydrophobic environment such as lipid membranes. Being the biological membranes the main target of many AMPs it is necessary to understand their way of action. Pandinin 2 (Pin2) is an alpha-helical AMP isolated from the venom of the African scorpion Pandinus imperator which shows high antimicrobial activity against Gram-positive bacteria and it is less active against Gram-negative bacteria, nevertheless, it has strong hemolytic activity. Its chemically synthesized Pin2GVG analog has low hemolytic activity while keeping its antimicrobial activity. With the aim of exploring the partition and subsequent folding of these peptides, in this work we report the results of extensive molecular dynamics simulations of Pin2 and Pin2GVG peptides in the presence of 2 hydrophobic environments such as dodecyl-phosphocholine (DPC) micelle and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocoline (POPC) membrane. Our results indicate that Pin2 folds in DPC with a 79% of alpha-helical content, which is in agreement with the experimental results, while in POPC it has 62.5% of alpha-helical content. On the other hand, Pin2GVG presents a higher percentage of alpha-helical structure in POPC and a smaller content in DPC when compared with Pin2. These results can help to better choose the starting structures in future molecular dynamics simulations of AMPs, because these peptides can adopt slightly different conformations depending on the hydrophobic environment.Communicated by Ramaswamy H. Sarma.

Theory of active chromatin remodeling.. [DOI: 10.1101/687145] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

Nucleosome positioning controls the accessible regions of chromatin and plays essential roles in DNA-templated processes. ATP driven remodeling enzymes are known to be crucial for its establishment in vivo, but their non-equilibrium nature has hindered the development of a unified theoretical framework for nucleosome positioning. Using a perturbation theory, we show that the effect of these enzymes can be well approximated by effective equilibrium models with rescaled temperatures and interactions. Numerical simulations support the accuracy of the theory in predicting both kinetic and steady-state quantities, including the effective temperature and the radial distribution function, in biologically relevant regimes. The energy landscape view emerging from our study provides an intuitive understanding for the impact of remodeling enzymes in either reinforcing or overwriting intrinsic signals for nucleosome positioning, and may help improve the accuracy of computational models for its prediction in silico.

1CPN: A coarse-grained multi-scale model of chromatin

A central question in epigenetics is how histone modifications influence the 3D structure of eukaryotic genomes and, ultimately, how this 3D structure is manifested in gene expression. The wide range of length scales that influence the 3D genome structure presents important challenges; epigenetic modifications to histones occur on scales of angstroms, yet the resulting effects of these modifications on genome structure can span micrometers. There is a scarcity of computational tools capable of providing a mechanistic picture of how molecular information from individual histones is propagated up to large regions of the genome. In this work, a new molecular model of chromatin is presented that provides such a picture. This new model, referred to as 1CPN, is structured around a rigorous multiscale approach, whereby free energies from an established and extensively validated model of the nucleosome are mapped onto a reduced coarse-grained topology. As such, 1CPN incorporates detailed physics from the nucleosome, such as histone modifications and DNA sequence, while maintaining the computational efficiency that is required to permit kilobase-scale simulations of genomic DNA. The 1CPN model reproduces the free energies and dynamics of both single nucleosomes and short chromatin fibers, and it is shown to be compatible with recently developed models of the linker histone. It is applied here to examine the effects of the linker DNA on the free energies of chromatin assembly and to demonstrate that these free energies are strongly dependent on the linker DNA length, pitch, and even DNA sequence. The 1CPN model is implemented in the LAMMPS simulation package and is distributed freely for public use.

Critical role of histone tail entropy in nucleosome unwinding

The nucleosome is the fundamental packaging unit for the genome. It must remain tightly wound to ensure genome stability while simultaneously being flexible enough to keep the DNA molecule accessible for genome function. The set of physicochemical interactions responsible for the delicate balance between these naturally opposed processes have not been determined due to challenges in resolving partially unwound nucleosome configurations at atomic resolution. Using a near atomistic protein-DNA model and advanced sampling techniques, we calculate the free energy cost of nucleosome DNA unwinding. Our simulations identify a large energetic barrier that decouples the outer and the inner DNA unwinding into two separate processes, occurring on different time scales. This dynamical decoupling allows the exposure of outer DNA at a modest cost to ensure accessibility while keeping the inner DNA and the histone core intact to maintain stability. We also reveal that this energetic barrier arises from a delayed loss of contacts between disordered histone tails and the DNA and is, surprisingly, largely offset by an entropic contribution from these tails. Implications of this enthalpy entropy compensation for the regulation of nucleosome stability and genome function are discussed.

The Oligomerization Landscape of Histones

In eukaryotes, DNA is packaged within nucleosomes. The DNA of each nucleosome is typically centered around an octameric histone protein core: one central tetramer plus two separate dimers. Studying the assembly mechanisms of histones is essential for understanding the dynamics of entire nucleosomes and higher-order DNA packaging. Here, we investigate canonical histone assembly and that of the centromere-specific histone variant, centromere protein A (CENP-A), using molecular dynamics simulations. We quantitatively characterize their thermodynamical and dynamical features, showing that two H3/H4 dimers form a structurally floppy, weakly bound complex, the latter exhibiting large instability around the central interface manifested via a swiveling motion of two halves. This finding is consistent with the recently observed DNA handedness flipping of the tetrasome. In contrast, the variant CENP-A encodes distinctive stability to its tetramer with a rigid but twisted interface compared to the crystal structure, implying diverse structural possibilities of the histone variant. Interestingly, the observed tetramer dynamics alter significantly and appear to reach a new balance when H2A/H2B dimers are present. Furthermore, we found that the preferred structure for the (CENP-A/H4)₂ tetramer is incongruent with the octameric structure, explaining many of the unusual dynamical behaviors of the CENP-A nucleosome. In all, these data reveal key mechanistic insights and structural details for the assembly of canonical and variant histone tetramers and octamers, providing theoretical quantifications and physical interpretations for longstanding and recent experimental observations. Based on these findings, we propose different chaperone-assisted binding and nucleosome assembly mechanisms for the canonical and CENP-A histone oligomers.

The Free Energy Landscape of Internucleosome Interactions and Its Relation to Chromatin Fiber Structure

The supramolecular chromatin fiber is governed by molecular scale energetics and interactions. Such energetics originate from the fiber's building block, the nucleosome core particle (NCP). In recent years, the chromatin fiber has been examined through perturbative methods in attempts to extract the energetics of nucleosome association in the fiber. This body of work has led to different results from experiments and simulations concerning the nucleosome-nucleosome energetics. Here, we expand on previous experiments and use coarse-grained simulations to evaluate the energetics inherent to nucleosomes across a variety of parameters in configurational and environmental space. Through this effort, we are able to uncover molecular processes that are critical to understanding the 30 nm chromatin fiber structure. In particular, we describe the NCP-NCP interactions by relying on an anisotropic energetic landscape, rather than a single potential energy value. The attractions in that landscape arise predominantly from the highly anisotropic interactions provided by the NCP histone N-terminal domain (NTD) tails. Our results are found to be in good agreement with recent nucleosome interaction experiments that suggest a maximum interaction energy of 2.69k _B T. Furthermore, we examine the influence of crucial epigenetic modifications, such as acetylation of the H4 tail, and how they modify the underlying landscape. Our results for acetylated NCP interactions are also in agreement with experiment. We additionally find an induced chirality in NCP-NCP interactions upon acetylation that reduces interactions which would correspond to a left-handed superhelical chromatin fiber.

Disparity in anomalous diffusion of proteins searching for their target DNA sites in a crowded medium is controlled by the size, shape and mobility of macromolecular crowders

Using extensive computer simulations, we analyzed the role of physical properties of molecular crowding agents such as size, shape and mobility in the target search dynamics of DNA binding proteins. Our main result is that the sub-diffusive dynamics of a protein inside a crowded medium strongly depends on the crowder properties and also on the protein's mode of diffusion. For instance, while scanning the DNA one-dimensionally, the protein dynamics does not vary with the change in crowder properties. Conversely, the diffusion exponent varies non-monotonically during 3D diffusion and is maximally affected when the crowders match the protein physically. The investigation shows that the effect stems from the ruggedness of the associated potential energy landscape, which is regulated by the protein-crowder and DNA-crowder interactions. Our findings have broad significance in understanding the target search dynamics of proteins on DNA in crowded cellular milieu and selecting appropriate crowding agents when designing in vitro experiments.

Improving Coarse-Grained Protein Force Fields with Small-Angle X-ray Scattering Data

Small-angle X-ray scattering (SAXS) experiments provide valuable structural data for biomolecules in solution. We develop a highly efficient maximum entropy approach to fit SAXS data by introducing minimal biases to a coarse-grained protein force field, the associative memory, water mediated, structure, and energy model (AWSEM). We demonstrate that the resulting force field, AWSEM-SAXS, succeeds in reproducing scattering profiles and models protein structures with shapes that are in much better agreement with experimental results. Quantitative metrics further reveal a modest, but consistent, improvement in the accuracy of modeled structures when SAXS data are incorporated into the force field. Additionally, when applied to a multiconformational protein, we find that AWSEM-SAXS is able to recover the population of different protein conformations from SAXS data alone. We, therefore, conclude that the maximum entropy approach is effective in fine-tuning the force field to better characterize both protein structure and conformational fluctuation.

Molecular Mechanism of Spontaneous Nucleosome Unraveling

Meters of DNA wrap around histone proteins to form nucleosomes and fit inside the micron-diameter nucleus. For the genetic information encoded in the DNA to become available for transcription, replication, and repair, the DNA-histone assembly must be disrupted. Experiment has indicated that the outer stretches of nucleosomal DNA "breathe" by spontaneously detaching from and reattaching to the histone core. Here, we report direct observation of spontaneous DNA breathing in atomistic molecular dynamics simulations, detailing a microscopic mechanism of the DNA breathing process. According to our simulations, the outer stretches of nucleosomal DNA detach in discrete steps involving 5 or 10 base pairs, with the detachment process being orchestrated by the motion of several conserved histone residues. The inner stretches of nucleosomal DNA are found to be more stably associated with the histone core by more abundant nonspecific DNA-protein contacts, providing a microscopic interpretation of nucleosome unraveling experiments. The CG content of nucleosomal DNA is found to anticorrelate with the extent of unwrapping, supporting the possibility that AT-rich segments may signal the start of transcription by forming less stable nucleosomes.

AWSEM-IDP: A Coarse-Grained Force Field for Intrinsically Disordered Proteins

The associative memory, water-mediated, structure and energy model (AWSEM) has been successfully used to study protein folding, binding, and aggregation problems. In this work, we introduce AWSEM-IDP, a new AWSEM branch for simulating intrinsically disordered proteins (IDPs), where the weights of the potentials determining secondary structure formation have been finely tuned, and a novel potential is introduced that helps to precisely control both the average extent of protein chain collapse and the chain's fluctuations in size. AWSEM-IDP can efficiently sample large conformational spaces, while retaining sufficient molecular accuracy to realistically model proteins. We applied this new model to two IDPs, demonstrating that AWSEM-IDP can reasonably well reproduce higher-resolution reference data, thus providing the foundation for a transferable IDP force field. Finally, we used thermodynamic perturbation theory to show that, in general, the conformational ensembles of IDPs are highly sensitive to fine-tuning of force field parameters.

Extensile motor activity drives coherent motions in a model of interphase chromatin

The 3D spatiotemporal organization of the human genome inside the cell nucleus remains a major open question in cellular biology. In the time between two cell divisions, chromatin-the functional form of DNA in cells-fills the nucleus in its uncondensed polymeric form. Recent in vivo imaging experiments reveal that the chromatin moves coherently, having displacements with long-ranged correlations on the scale of micrometers and lasting for seconds. To elucidate the mechanism(s) behind these motions, we develop a coarse-grained active polymer model where chromatin is represented as a confined flexible chain acted upon by molecular motors that drive fluid flows by exerting dipolar forces on the system. Numerical simulations of this model account for steric and hydrodynamic interactions as well as internal chain mechanics. These demonstrate that coherent motions emerge in systems involving extensile dipoles and are accompanied by large-scale chain reconfigurations and nematic ordering. Comparisons with experiments show good qualitative agreement and support the hypothesis that self-organizing long-ranged hydrodynamic couplings between chromatin-associated active motor proteins are responsible for the observed coherent dynamics.

Chromatin remodelers couple inchworm motion with twist-defect formation to slide nucleosomal DNA

ATP-dependent chromatin remodelers are molecular machines that control genome organization by repositioning, ejecting, or editing nucleosomes, activities that confer them essential regulatory roles on gene expression and DNA replication. Here, we investigate the molecular mechanism of active nucleosome sliding by means of molecular dynamics simulations of the Snf2 remodeler translocase in complex with a nucleosome. During its inchworm motion driven by ATP consumption, the translocase overwrites the original nucleosome energy landscape via steric and electrostatic interactions to induce sliding of nucleosomal DNA unidirectionally. The sliding is initiated at the remodeler binding location via the generation of a pair of twist defects, which then spontaneously propagate to complete sliding throughout the entire nucleosome. We also reveal how remodeler mutations and DNA sequence control active nucleosome repositioning, explaining several past experimental observations. These results offer a detailed mechanistic picture of remodeling important for the complete understanding of these key biological processes.

Nucleosomes are the protein-DNA complexes underlying Eukaryotic genome organization, and serve as regulators of gene expression by occluding DNA to other proteins. This regulation requires the precise positioning of nucleosomes along DNA. Chromatin remodelers are the molecular machines that consume ATP to slide nucleosome at their correct locations, but the mechanisms of remodeling are still unclear. Based on the static structural information of a remodeler bound on nucleosome, we performed molecular dynamics computer simulations revealing the details of how remodelers slide nucleosomal DNA: the inchworm-like motion of remodelers create small DNA deformations called twist defects, which then spontaneously propagate throughout the nucleosome to induce sliding. These simulations explain several past experimental findings and are important for our understanding of genome organization.

Exploring the free-energy landscape of GPCR activation

G-protein-coupled receptors (GPCRs) are a large group of membrane-bound receptor proteins that are involved in a plethora of diverse processes (e.g., vision, hormone response). In mammals, and particularly in humans, GPCRs are involved in many signal transduction pathways and, as such, are heavily studied for their immense pharmaceutical potential. Indeed, a large fraction of drugs target various GPCRs, and drug-development is often aimed at GPCRs. Therefore, understanding the activation of GPCRs is a challenge of major importance both from fundamental and practical considerations. And yet, despite the remarkable progress in structural understanding, we still do not have a translation of the structural information to an energy-based picture. Here we use coarse-grained (CG) modeling to chart the free-energy landscape of the activation process of the β-2 adrenergic receptor (β₂AR) as a representative GPCR. The landscape provides the needed tool for analyzing the processes that lead to activation of the receptor upon binding of the ligand (adrenaline) while limiting constitutive activation. Our results pave the way to better understand the biological mechanisms of action of the β₂AR and GPCRs, from a physical chemistry point of view rather than simply by observing the receptor's behavior physiologically.

Interactions of HP1 Bound to H3K9me3 Dinucleosome by Molecular Simulations and Biochemical Assays

Heterochromatin protein 1 (HP1), associated with heterochromatin formation, recognizes an epigenetically repressive marker, trimethylated lysine 9 in histone H3 (H3K9me3), and generally contributes to long-term silencing. How HP1 induces heterochromatin is not fully understood. Recent experiments suggested that not one, but two nucleosomes provide a platform for this recognition. Integrating previous and new biochemical assays with computational modeling, we provide near-atomic structural models for HP1 binding to the dinucleosomes. We found that the dimeric HP1α tends to bind two H3K9me3s that are in adjacent nucleosomes, thus bridging two nucleosomes. We identified, to our knowledge, a novel DNA binding motif in the hinge region that is specific to HP1α and is essential for recognizing the H3K9me3 sites of two nucleosomes. An HP1 isoform, HP1γ, does not easily bridge two nucleosomes in extended conformations because of the absence of the above binding motif and its shorter hinge region. We propose a molecular mechanism for chromatin structural changes caused by HP1.

DNA sliding in nucleosomes via twist defect propagation revealed by molecular simulations

While nucleosomes are highly stable structures as fundamental units of chromatin, they also slide along the DNA, either spontaneously or by active remodelers. Here, we investigate the microscopic mechanisms of nucleosome sliding by multiscale molecular simulations, characterizing how the screw-like motion of DNA proceeds via the formation and propagation of twist defects. Firstly, coarse-grained molecular simulations reveal that the sliding dynamics is highly dependent on DNA sequence. Depending on the sequence and the nucleosome super-helical location, we find two distinct types of twist defects: a locally under-twisted DNA region, previously observed in crystal structures, and a locally over-twisted DNA, an unprecedented feature. The stability of the over-twist defect was confirmed via all-atom simulations. Analysis of our trajectories via Markov state modeling highlights how the sequence-dependence of the sliding dynamics is due to the different twist defect energy costs, and in particular how nucleosome regions where defects cannot easily form introduce the kinetic bottlenecks slowing down repositioning. Twist defects can also mediate sliding of nucleosomes made with strong positioning sequences, albeit at a much lower diffusion coefficient, due to a high-energy intermediate state. Finally, we discuss how chromatin remodelers may exploit these spontaneous fluctuations to induce unidirectional sliding of nucleosomes.

Slowest kinetic modes revealed by metabasin renormalization

Understanding the slowest relaxations of complex systems, such as relaxation of glass-forming materials, diffusion in nanoclusters, and folding of biomolecules, is important for physics, chemistry, and biology. For a kinetic system, the relaxation modes are determined by diagonalizing its transition rate matrix. However, for realistic systems of interest, numerical diagonalization, as well as extracting physical understanding from the diagonalization results, is difficult due to the high dimensionality. Here, we develop an alternative and generally applicable method of extracting the long-time scale relaxation dynamics by combining the metabasin analysis of Okushima et al. [Phys. Rev. E 80, 036112 (2009)PLEEE81539-375510.1103/PhysRevE.80.036112] and a Jacobi method. We test the method on an illustrative model of a four-funnel model, for which we obtain a renormalized kinematic equation of much lower dimension sufficient for determining slow relaxation modes precisely. The method is successfully applied to the vacancy transport problem in ionic nanoparticles [Niiyama et al., Chem. Phys. Lett. 654, 52 (2016)CHPLBC0009-261410.1016/j.cplett.2016.04.088], allowing a clear physical interpretation that the final relaxation consists of two successive, characteristic processes.

Modulation of nucleosomal DNA accessibility via charge-altering post-translational modifications in histone core

BACKGROUND

Controlled modulation of nucleosomal DNA accessibility via post-translational modifications (PTM) is a critical component to many cellular functions. Charge-altering PTMs in the globular histone core-including acetylation, phosphorylation, crotonylation, propionylation, butyrylation, formylation, and citrullination-can alter the strong electrostatic interactions between the oppositely charged nucleosomal DNA and the histone proteins and thus modulate accessibility of the nucleosomal DNA, affecting processes that depend on access to the genetic information, such as transcription. However, direct experimental investigation of the effects of these PTMs is very difficult. Theoretical models can rationalize existing observations, suggest working hypotheses for future experiments, and provide a unifying framework for connecting PTMs with the observed effects.

RESULTS

A physics-based framework is proposed that predicts the effect of charge-altering PTMs in the histone core, quantitatively for several types of lysine charge-neutralizing PTMs including acetylation, and qualitatively for all phosphorylations, on the nucleosome stability and subsequent changes in DNA accessibility, making a connection to resulting biological phenotypes. The framework takes into account multiple partially assembled states of the nucleosome at the atomic resolution. The framework is validated against experimentally known nucleosome stability changes due to the acetylation of specific lysines, and their effect on transcription. The predicted effect of charge-altering PTMs on DNA accessibility can vary dramatically, from virtually none to a strong, region-dependent increase in accessibility of the nucleosomal DNA; in some cases, e.g., H4K44, H2AK75, and H2BK57, the effect is significantly stronger than that of the extensively studied acetylation sites such H3K56, H3K115 or H3K122. Proximity to the DNA is suggestive of the strength of the PTM effect, but there are many exceptions. For the vast majority of charge-altering PTMs, the predicted increase in the DNA accessibility should be large enough to result in a measurable modulation of transcription. However, a few possible PTMs, such as acetylation of H4K77, counterintuitively decrease the DNA accessibility, suggestive of the repressed chromatin. A structural explanation for the phenomenon is provided. For the majority of charge-altering PTMs, the effect on DNA accessibility is simply additive (noncooperative), but there are exceptions, e.g., simultaneous acetylation of H4K79 and H3K122, where the combined effect is amplified. The amplification is a direct consequence of the nucleosome-DNA complex having more than two structural states. The effect of individual PTMs is classified based on changes in the accessibility of various regions throughout the nucleosomal DNA. The PTM's resulting imprint on the DNA accessibility, "PTMprint," is used to predict effects of many yet unexplored PTMs. For example, acetylation of H4K44 yields a PTMprint similar to the PTMprint of H3K56, and thus acetylation of H4K44 is predicted to lead to a wide range of strong biological effects.

CONCLUSION

Charge-altering post-translational modifications in the relatively unexplored globular histone core may provide a precision mechanism for controlling accessibility to the nucleosomal DNA.

Reconstruction of Atomistic Structures from Coarse-Grained Models for Protein-DNA Complexes

While coarse-grained (CG) simulations have widely been used to accelerate structure sampling of large biomolecular complexes, they are unavoidably less accurate and thus the reconstruction of all-atom (AA) structures and the subsequent refinement is desirable. In this study we developed an efficient method to reconstruct AA structures from sampled CG protein-DNA complex models, which attempts to model the protein-DNA interface accurately. First we developed a method to reconstruct atomic details of DNA structures from a three-site per nucleotide CG model, which uses a DNA fragment library. Next, for the protein-DNA interface, we referred to the side chain orientations in the known structure of the target interface when available. The other parts are modeled by existing tools. We confirmed the accuracy of the protocol in various aspects including the structure deviation in the self-reproduction, the base pair reproducibility, atomic contacts at the protein-DNA interface, and feasibility of the posterior AA simulations.

Role of Macromolecular Crowding on the Intracellular Diffusion of DNA Binding Proteins

Recent experiments suggest that cellular crowding facilitates the target search dynamics of proteins on DNA, the mechanism of which is not yet known. By using large scale computer simulations, we show that two competing factors, namely the width of the depletion layer that separates the crowder cloud from the DNA molecule and the degree of protein-crowder crosstalk, act in harmony to affect the target search dynamics of proteins. The impacts vary from nonspecific to specific target search regime. During a nonspecific search, dynamics of a protein is only minimally affected, whereas, a significantly different behaviour is observed when the protein starts forming a specific protein-DNA complex. We also find that the severity of impacts largely depends upon physiological crowder concentration and deviation from it leads to attenuation in the binding kinetics. Based on extensive kinetic study and binding energy landscape analysis, we further present a comprehensive molecular description of the search process that allows us to interpret the experimental findings.

Sequence-dependent nucleosome sliding in rotation-coupled and uncoupled modes revealed by molecular simulations

While nucleosome positioning on eukaryotic genome play important roles for genetic regulation, molecular mechanisms of nucleosome positioning and sliding along DNA are not well understood. Here we investigated thermally-activated spontaneous nucleosome sliding mechanisms developing and applying a coarse-grained molecular simulation method that incorporates both long-range electrostatic and short-range hydrogen-bond interactions between histone octamer and DNA. The simulations revealed two distinct sliding modes depending on the nucleosomal DNA sequence. A uniform DNA sequence showed frequent sliding with one base pair step in a rotation-coupled manner, akin to screw-like motions. On the contrary, a strong positioning sequence, the so-called 601 sequence, exhibits rare, abrupt transitions of five and ten base pair steps without rotation. Moreover, we evaluated the importance of hydrogen bond interactions on the sliding mode, finding that strong and weak bonds favor respectively the rotation-coupled and -uncoupled sliding movements.

Electrostatic forces govern the binding mechanism of intrinsically disordered histone chaperones

A unified picture to understand the protein recognition and function must include the native binding complex structure ensembles and the underlying binding mechanisms involved in specific biological processes. However, quantifications of both binding complex structures and dynamical mechanisms are still challenging for IDP. In this study, we have investigated the underlying molecular mechanism of the chaperone Chz1 and histone H2A.Z-H2B association by equilibrium and kinetic stopped-flow fluorescence spectroscopy. The dependence of free energy and kinetic rate constant on electrolyte mean activity coefficient and urea concentration are uncovered. Our results indicate a previous unseen binding kinetic intermediate. An initial conformation selection step of Chz1 is also revealed before the formation of this intermediate state. Based on these observations, a mixed mechanism of three steps including both conformation selection and induced fit is proposed. By combination of the ion- and denaturant-induced experiments, we demonstrate that electrostatic forces play a dominant role in the recognition of bipolar charged intrinsically disordered protein Chz1 to its preferred partner H2A.Z-H2B. Both the intra-chain and inter-chain electrostatic interactions have direct impacts on the native collapsed structure and binding mechanism.

Kilobase Pair Chromatin Fiber Contacts Promoted by Living-System-Like DNA Linker Length Distributions and Nucleosome Depletion

Nucleosome placement, or DNA linker length patterns, are believed to yield specific spatial features in chromatin fibers, but details are unknown. Here we examine by mesoscale modeling how kilobase (kb) range contacts and fiber looping depend on linker lengths ranging from 18 to 45 bp, with values modeled after living systems, including nucleosome free regions (NFRs) and gene encoding segments. We also compare artificial constructs with alternating versus randomly distributed linker lengths in the range of 18-72 bp. We show that nonuniform distributions with NFRs enhance flexibility and encourage kb-range contacts. NFRs between neighboring gene segments diminish short-range contacts between flanking nucleosomes, while enhancing kb-range contacts via hierarchical looping. We also demonstrate that variances in linker lengths enhance such contacts. In particular, moderate sized variations in fiber linker lengths (∼27 bp) encourage long-range contacts in randomly distributed linker length fibers. Our work underscores the importance of linker length patterns, alongside bound proteins, in biological regulation. Contacts formed by kb-range chromatin folding are crucial to gene activity. Because we find that special linker length distributions in living systems promote kb contacts, our work suggests ways to manipulate these patterns for regulation of gene activity.

Genomic Energy Landscapes

Energy landscape theory, developed in the context of protein folding, provides, to our knowledge, a new perspective on chromosome architecture. We review what has been learned concerning the topology and structure of both the interphase and mitotic chromosomes from effective energy landscapes constructed using Hi-C data. Energy landscape thinking raises new questions about the nonequilibrium dynamics of the chromosome and gene regulation.

Linking Chromatin Fibers to Gene Folding by Hierarchical Looping

While much is known about DNA structure on the basepair level, this scale represents only a fraction of the structural levels involved in folding the genomic material. With recent advances in experimental and theoretical techniques, a variety of structures have been observed on the fiber, gene, and chromosome levels of genome organization. Here we view chromatin architecture from nucleosomes and fibers to genes and chromosomes, highlighting the rich structural diversity and fiber fluidity emerging from both experimental and theoretical techniques. In this context, we discuss our recently proposed folding mechanism, which we call "hierarchical looping", similar to rope flaking used in mountain climbing, where 10-nm zigzag chromatin fibers are compacted laterally into self-associating loops which then stack and fold in space. We propose that hierarchical looping may act as a bridge between fibers and genes as well as provide a mechanism to relate key features of interphase and metaphase chromosome architecture to genome structural changes. This motif emerged by analysis of ultrastructural internucleosome contact data by electron microscopy-assisted nucleosome interaction capture cross-linking experiments, in combination with mesoscale modeling. We suggest that while the local folding of chromatin can be regulated at the fiber level by adjustment of internal factors such as linker-histone binding affinities, linker DNA lengths, and divalent ion levels, hierarchical looping on the gene level can additionally be controlled by posttranslational modifications and external factors such as polycomb group proteins. From a combination of 3C data and mesoscale modeling, we suggest that hierarchical looping could also play a role in epigenetic gene silencing, as stacked loops may occlude access to transcription start sites. With advances in crystallography, single-molecule in vitro biochemistry, in vivo imaging techniques, and genome-wide contact data experiments, various modeling approaches are allowing for previously unavailable structural interpretation of these data at multiple spatial and temporal scales. An unprecedented level of productivity and opportunity is on the horizon for the chromatin structure field.

Protein dynamics: Conformational footprints

The mechanisms by which proteins evolve new functions can be slow and mysterious. Comprehensive structural analysis of enzyme variants reveal how gradual enrichments of pre-existing populations with the right productive dynamics for new functions can accomplish this aim.

Molecular Mechanism of Facilitated Dissociation of Fis Protein from DNA

Fis protein is a nucleoid-associated protein that plays many roles in transcriptional regulation and DNA site-specific recombination. In contrast to the naïve expectation based on stoichiometry, recent single-molecule studies have shown that the dissociation of Fis protein from DNA is accelerated by increasing the concentration of the Fis protein. Because the detailed molecular mechanism of facilitated dissociation is still not clear, in this study, we employ computational methods to explore the binding landscapes of Fis:DNA complexes with various stoichiometries. When two Fis molecules are present, simulations uncover a ternary complex, where the originally bound Fis protein is partially dissociated from DNA. The simulations support a three-state sequential kinetic model (N ⇄ I → D) for facilitated dissociation, thus explaining the concentration-dependent dissociation.

Single-Molecule Observation Reveals Spontaneous Protein Dynamics in the Nucleosome

Structural dynamics of a protein molecule is often critical to its function. Single-molecule methods provide efficient ways to investigate protein dynamics, although it is very challenging to achieve a millisecond or higher temporal resolution. Here we report spontaneous structural dynamics of the histone protein core in the nucleosome based on a single-molecule method that can reveal submillisecond dynamics by combining maximum likelihood estimation and fluorescence correlation spectroscopy. The nucleosome, comprising ∼147 bp DNA and an octameric histone protein core consisting of H2A, H2B, H3, and H4, is the fundamental packing unit of the eukaryotic genome. The nucleosome imposes a physical barrier that should be overcome during various DNA-templated processes. Structural fluctuation of the nucleosome in the histone core has been hypothesized to be required for nucleosome disassembly but has yet to be directly probed. Our results indicate that at 100 mM NaCl the histone H2A-H2B dimer dissociates from the histone core transiently once every 3.6 ± 0.6 ms and returns to its position within 2.0 ± 0.3 ms. We also found that the motion is facilitated upon H3K56 acetylation and inhibited upon replacing H2A with H2A.Z. These results provide the first direct examples of how a localized post-translational modification or an epigenetic variation affects the kinetic and thermodynamic stabilities of a macromolecular protein complex, which may directly contribute to its functions.