Dissecting DNA-histone interactions in the nucleosome by molecular dynamics simulations of DNA unwrapping.

The nucleosome complex of DNA wrapped around a histone protein octamer organizes the genome of eukaryotes and regulates the access of protein factors to the DNA. We performed molecular dynamics simulations of the nucleosome in explicit water to study the dynamics of its histone-DNA interactions. A high-resolution histone-DNA interaction map was derived that revealed a five-nucleotide periodicity, in which the two DNA strands of the double helix made alternating contacts. On the 100-ns timescale, the histone tails mostly maintained their initial positions relative to the DNA, and the spontaneous unwrapping of DNA was limited to 1-2 basepairs. In steered molecular dynamics simulations, external forces were applied to the linker DNA to investigate the unwrapping pathway of the nucleosomal DNA. In comparison with a nucleosome without the unstructured N-terminal histone tails, the following findings were obtained: 1), Two main barriers during unwrapping were identified at DNA position ±70 and ±45 basepairs relative to the central DNA basepair at the dyad axis. 2), DNA interactions of the histone H3 N-terminus and the histone H2A C-terminus opposed the initiation of unwrapping. 3), The N-terminal tails of H2A, H2B, and H4 counteracted the unwrapping process at later stages and were essential determinants of nucleosome dynamics. Our detailed analysis of DNA-histone interactions revealed molecular mechanisms for modulating access to nucleosomal DNA via conformational rearrangements of its structure.

Force spectroscopy of chromatin fibers: extracting energetics and structural information from Monte Carlo simulations.

The folding of the nucleosome chain into a chromatin fiber is a central factor for controlling the DNA access of protein factors involved in transcription, DNA replication and repair. Force spectroscopy experiments with chromatin fibers are ideally suited to dissect the interactions that drive this process, and to probe the underlying fiber conformation. However, the interpretation of the experimental data is fraught with difficulties due to the complex interplay of the nucleosome geometry and the different energy terms involved. Here, we apply a Monte Carlo simulation approach to derive virtual chromatin fiber force spectroscopy curves. In the simulations, the effect of the nucleosome geometry, repeat length, nucleosome-nucleosome interaction potential, and the unwrapping of the DNA from the histone protein core on the shape of the force-extension curves was investigated. These simulations provide a framework for the evaluation of experimental data sets. We demonstrate how the relative contributions of DNA bending and twisting, nucleosome unstacking and unwrapping the nucleosomal DNA from the histone octamer can be dissected for a given fiber geometry.

Parallel high-performance grid computing: capabilities and opportunities of a novel demanding service and business class allowing highest resource efficiency.

Especially in the life-science and the health-care sectors the huge IT requirements are imminent due to the large and complex systems to be analysed and simulated. Grid infrastructures play here a rapidly increasing role for research, diagnostics, and treatment, since they provide the necessary large-scale resources efficiently. Whereas grids were first used for huge number crunching of trivially parallelizable problems, increasingly parallel high-performance computing is required. Here, we show for the prime example of molecular dynamic simulations how the presence of large grid clusters including very fast network interconnects within grid infrastructures allows now parallel high-performance grid computing efficiently and thus combines the benefits of dedicated super-computing centres and grid infrastructures. The demands for this service class are the highest since the user group has very heterogeneous requirements: i) two to many thousands of CPUs, ii) different memory architectures, iii) huge storage capabilities, and iv) fast communication via network interconnects, are all needed in different combinations and must be considered in a highly dedicated manner to reach highest performance efficiency. Beyond, advanced and dedicated i) interaction with users, ii) the management of jobs, iii) accounting, and iv) billing, not only combines classic with parallel high-performance grid usage, but more importantly is also able to increase the efficiency of IT resource providers. Consequently, the mere "yes-we-can" becomes a huge opportunity like e.g. the life-science and health-care sectors as well as grid infrastructures by reaching higher level of resource efficiency.

Exploring the conformational space of chromatin fibers and their stability by numerical dynamic phase diagrams.

The three-dimensional structure of chromatin affects DNA accessibility and is therefore a key regulator of gene expression. However, the path of the DNA between consecutive nucleosomes, and the resulting chromatin fiber organization remain controversial. The conformational space available for the folding of the nucleosome chain has been analytically described by phase diagrams with a two-angle model, which describes the chain trajectory by a DNA entry-exit angle at the nucleosome and a torsion angle between consecutive nucleosomes. Here, a novel type of numerical phase diagrams is introduced that relates the geometric phase space to the energy associated with a given chromatin conformation. The resulting phase diagrams revealed differences in the energy landscape that reflect the probability of a given conformation to form in thermal equilibrium. Furthermore, we investigated the effects of entropy and additional degrees of freedom in the dynamic phase diagrams by performing Monte Carlo simulations of the initial chain trajectories. Using our approach, we were able to demonstrate that conformations that initially were geometrically impossible could evolve into energetically favorable states in thermal equilibrium due to DNA bending and torsion. In addition, dynamic phase diagrams were applied to identify chromatin fibers that reflect certain experimentally determined features.

DNA-DNA interactions in tight supercoils are described by a small effective charge density.

DNA-DNA interactions are important for genome compaction and transcription regulation. In studies of such complex processes, DNA is often modeled as a homogeneously charged cylinder and its electrostatic interactions are calculated within the framework of the Poisson-Boltzmann equation. Commonly, a charge adaptation factor is used to address limitations of this theoretical approach. Despite considerable theoretical and experimental efforts, a rigorous quantitative assessment of this parameter is lacking. Here, we comprehensively characterized DNA-DNA interactions in the presence of monovalent ions by analyzing the supercoiling behavior of single DNA molecules held under constant tension. Both a theoretical model and coarse-grained simulations of this process revealed a surprisingly small effective DNA charge of 40% of the nominal charge density, which was additionally supported by all-atom molecular dynamics simulations.

The effect of internucleosomal interaction on folding of the chromatin fiber.

The folding of the nucleosome chain into a chromatin fiber modulates DNA accessibility and is therefore an important factor for the control of gene expression. The fiber conformation depends crucially on the interaction between individual nucleosomes. However, this parameter has not been accurately determined experimentally, and it is affected by posttranslational histone modifications and binding of chromosomal proteins. Here, the effect of different internucleosomal interaction strengths on the fiber conformation was investigated by Monte Carlo computer simulations. The fiber geometry was modeled to fit that of chicken erythrocyte chromatin, which has been examined in numerous experimental studies. In the Monte Carlo simulation, the nucleosome shape was described as an oblate spherocylinder, and a replica exchange protocol was developed to reach thermal equilibrium for a broad range of internucleosomal interaction energies. The simulations revealed the large impact of the nucleosome geometry and the nucleosome repeat length on the compaction of the chromatin fiber. At high internucleosomal interaction energies, a lateral self-association of distant fiber parts and an interdigitation of nucleosomes were apparent. These results identify key factors for the control of the compaction and higher order folding of the chromatin fiber.

Nucleosome geometry and internucleosomal interactions control the chromatin fiber conformation.

Based on model structures with atomic resolution, a coarse-grained model for the nucleosome geometry was implemented. The dependence of the chromatin fiber conformation on the spatial orientation of nucleosomes and the path and length of the linker DNA was systematically explored by Monte Carlo simulations. Two fiber types were analyzed in detail that represent nucleosome chains without and with linker histones, respectively: two-start helices with crossed-linker DNA (CL conformation) and interdigitated one-start helices (ID conformation) with different nucleosome tilt angles. The CL conformation was derived from a tetranucleosome crystal structure that was extended into a fiber. At thermal equilibrium, the fiber shape persisted but relaxed into a structure with a somewhat lower linear mass density of 3.1 +/- 0.1 nucleosomes/11 nm fiber. Stable ID fibers required local nucleosome tilt angles between 40 degrees and 60 degrees. For these configurations, much higher mass densities of up to 7.9 +/- 0.2 nucleosomes/11 nm fiber were obtained. A model is proposed, in which the transition between a CL and ID fiber is mediated by relatively small changes of the local nucleosome geometry. These were found to be in very good agreement with changes induced by linker histone H1 binding as predicted from the high resolution model structures.