Histone dynamics mediate DNA unwrapping and sliding in nucleosomes
Nucleosomes are elementary building blocks of chromatin in eukaryotes. They tightly wrap ∼147 DNA base pairs around an octamer of histone proteins. How nucleosome structural dynamics affect genome functioning is not completely clear. Here we report all-atom molecular dynamics simulations of nucleosome core particles at a timescale of 15 micro- seconds. At this timescale, functional modes of nucleosome dynamics such as spontaneous nucleosomal DNA breathing, unwrapping, twisting, and sliding were observed. We identified atomistic mechanisms of these processes by analyzing the accompanying structural rear- rangements of the histone octamer and histone-DNA contacts. Octamer dynamics and plasticity were found to enable DNA unwrapping and sliding. Through multi-scale modeling, we showed that nucleosomal DNA dynamics contribute to significant conformational variability of the chromatin fiber at the supranucleosomal level. Our study further supports mechanistic coupling between fine details of histone dynamics and chromatin functioning, provides a framework for understanding the effects of various chromatin modifications.
To gain success in the evolutionary “arms race”, venomous animals such as scorpions produce diverse neurotoxins selected to hit targets in the nervous system of prey. Scorpion α-toxins affect insect and/or mammalian voltage-gated sodium channels (Nav’s) and thereby modify the excitability of muscle and nerve cells. Although more than a hundred α-toxins are known and a number of them have been studied into detail, the molecular mechanism of their interaction with Nav’s is still poorly understood. Here, we employ extensive molecular dynamics simulations and spatial mapping of hydrophobic/hydrophilic properties distributed over the molecular surface of α-toxins. It is revealed that in spite of the small size and relatively rigid structure, these toxins possess modular organization from structural, functional and evolutionary perspectives. The more conserved and rigid “core module” is supplemented with the “specificity module” (SM) that is comparatively flexible and variable, and determines the taxon (mammal vs. insect) specificity of α-toxin activity. We further show that SMs in mammal toxins are more flexible and hydrophilic than in insect toxins. Concomitant sequence-based analysis of Nav’s extracellular loops suggests that α-toxins recognize the channels using both modules. We propose that the core module binds to the voltage-sensing domain of repeat IV, whereas the more versatile SM interacts with the pore domain in repeat I of Nav’s. These findings corroborate and expand the hypothesis on different functional epitopes of toxins that has been reported previously. In effect, we propose that the modular structure in toxins evolved to match the domain architecture of Nav’s.
Non-B DNA structures have a great potential to form and influence various genomic processes including transcription. One of the mechanisms of transcription regulation is nucleo- some positioning. Even though only B-DNA can be wrapped around a nucleosome, non-B DNA structures can compete with a nucleosome for a genomic location. Here we used perman- ganate/S1 nuclease footprinting data on non-B DNA structures, such as Z-DNA, H-DNA, G- quadruplexes and stress-induced duplex destabilization (SIDD) sites, together with MNase-seq data on nucleosome positioning in the mouse genome. We found three types of patterns of nucleosome positioning around non-B DNA structures: a structure is surrounded by nucleo- somes from both sides, from one side, or nucleosome free region. Machine learning models based on random forest and XGBoost algorithms were constructed to recognize DNA regions of 1kB length containing a particular pattern of nucleosome positioning for four types of DNA structures (Z-DNA, H-DNA, G-quadruplexes and SIDD sites) based on statistics of di- and tri- nucleotides. The best performance (94% of accuracy) was reached for G-quadruplexes while for other types of structures the accuracy was under 70%. We conclude that 1kB regions con- taining G-quadruplexes have distinct compositional properties, and this fact points to preferen- tial locations of such pattern in the genome and requires further investigation. For other DNA structures a region composition is not a sufficient predictive factor and one should take into account other physical and structural DNA properties to improve nucleosome-DNA-structure pattern recognition.
Non-B DNA structures have a great potential to form and influence various genomic processes including transcription. One of the mechanisms of transcription regulation is nucleosome positioning. Even though only B-DNA can be wrapped around a nucleosome, non-B DNA structures can compete with a nucleosome for a genomic location. Here we used permanganate/S1 nuclease footprinting data on non-B DNA structures, such as Z-DNA, H-DNA, G-quadruplexes and stress-induced duplex destabilization (SIDD) sites, together with MNase-seq data on nucleosome positioning in the mouse genome. We found three types of patterns of nucleosome positioning around non-B DNA structures: a structure is surrounded by nucleosomes from both sides, from one side, or nucleosome free region. Machine learning models based on random forest and XGBoost algorithms were constructed to recognize DNA regions of 1kB length containing a particular pattern of nucleosome positioning for four types of DNA structures (Z-DNA, H-DNA, G-quadruplexes and SIDD sites) based on statistics of di- and tri-nucleotides. The best performance (94% of accuracy) was reached for Gquadruplexes while for other types of structures the accuracy was under 70%. We conclude that 1kB regions containing Gquadruplexes have distinct compositional properties, and this fact points to preferential locations of such pattern in the genome and requires further investigation. Gene ontology analysis revealed that the genes intersecting with the discovered patterns are enriched in channel and transmembrane activity, transcription factor and receptor binding. The direction for further research is to study the distribution of the discovered patterns in different tissues to identify well-positioned and dynamic nucleosomes and reveal genes, regulated via DNA structures and nucleosome positioning.
Mutations in transmembrane (TM) domains of receptor tyrosine kinases are shown to cause a number of inherited diseases and cancer development. Here, we use a combined molecular modeling approach to understand molecular mechanism of effect of G370R and A391E mutations on dimerization of TM domains of human fibroblast growth factor receptor 3 (FGFR3). According to results of Monte-Carlo conformational search in the implicit membrane and further molecular dynamics simulations, TM dimer of this receptor is able to form a number of various conformations, which differ significantly by the free energy of association in a full-atom model bilayer. The aforementioned mutations affect dimerization efficiency of TM segments and lead to repopulation of conformational ensemble for the dimer. Particularly, both mutations do not change the dimerization free energy of the predominant (putative “non-active”) symmetric conformation of TM dimer, while affect dimerization efficiency of its asymmetric (“intermediate”) and alternative symmetric (putative “active”) models. Results of our simulations provide novel atomistic prospective of the role of G370 and A391E mutations in dimerization of TM domains of FGFR3 and their consecutive contributions to the activation pathway of the receptor.
Ribosomal protein S2 is an essential component of translation machinery, and its viable mutated variants conferring distinct phenotypes serve as a valuable tool in studying the role of S2 in translation regulation. One of a few available rpsB mutants, rpsB1, shows thermosensitivity and ensures enhanced expression of leaderless mRNAs. In this study, we identified the nature of the rpsB1 mutation. Sequencing of the rpsB1 allele revealed a G-to-A transition in the part of the rpsB gene which encodes a coiled-coil domain of S2. The resulting E132K substitution resides in a highly conserved site, TKKE, a so-called N-terminal capping box, at the beginning of the second alpha helix. The protruding coiled-coil domain of S2 is known to provide binding with 16S rRNA in the head of the 30S subunit and, in addition, to interact with a key mRNA binding protein, S1. Molecular dynamics simulations revealed a detrimental impact of the E132K mutation on the coiled-coil structure and thereby on the interactions between S2 and 16S rRNA, providing a clue for the thermosensitivity of the rpsB1 mutant. Using a strain producing a leaderless lacZ transcript from the chromosomal lac promoter, we demonstrated that not only the rpsB1 mutation generating S2/S1-deficient ribosomes but also the rpsA::IS10 mutation leading to partial deficiency in S1 alone increased translation efficiency of the leaderless mRNA by about 10-fold. Moderate overexpression of S1 relieved all these effects and, moreover, suppressed the thermosensitive phenotype of rpsB1, indicating the role of S1 as an extragenic suppressor of the E132K mutation.
A model for organizing cargo transportation between two node stations connected by a railway line which contains a certain number of intermediate stations is considered. The movement of cargo is in one direction. Such a situation may occur, for example, if one of the node stations is located in a region which produce raw material for manufacturing industry located in another region, and there is another node station. The organization of freight traﬃc is performed by means of a number of technologies. These technologies determine the rules for taking on cargo at the initial node station, the rules of interaction between neighboring stations, as well as the rule of distribution of cargo to the ﬁnal node stations. The process of cargo transportation is followed by the set rule of control. For such a model, one must determine possible modes of cargo transportation and describe their properties. This model is described by a ﬁnite-dimensional system of diﬀerential equations with nonlocal linear restrictions. The class of the solution satisfying nonlocal linear restrictions is extremely narrow. It results in the need for the “correct” extension of solutions of a system of diﬀerential equations to a class of quasi-solutions having the distinctive feature of gaps in a countable number of points. It was possible numerically using the Runge–Kutta method of the fourth order to build these quasi-solutions and determine their rate of growth. Let us note that in the technical plan the main complexity consisted in obtaining quasi-solutions satisfying the nonlocal linear restrictions. Furthermore, we investigated the dependence of quasi-solutions and, in particular, sizes of gaps (jumps) of solutions on a number of parameters of the model characterizing a rule of control, technologies for transportation of cargo and intensity of giving of cargo on a node station.
Event logs collected by modern information and technical systems usually contain enough data for automated process models discovery. A variety of algorithms was developed for process models discovery, conformance checking, log to model alignment, comparison of process models, etc., nevertheless a quick analysis of ad-hoc selected parts of a journal still have not get a full-fledged implementation. This paper describes an ROLAP-based method of multidimensional event logs storage for process mining. The result of the analysis of the journal is visualized as directed graph representing the union of all possible event sequences, ranked by their occurrence probability. Our implementation allows the analyst to discover process models for sublogs defined by ad-hoc selection of criteria and value of occurrence probability
Existing approaches suggest that IT strategy should be a reflection of business strategy. However, actually organisations do not often follow business strategy even if it is formally declared. In these conditions, IT strategy can be viewed not as a plan, but as an organisational shared view on the role of information systems. This approach generally reflects only a top-down perspective of IT strategy. So, it can be supplemented by a strategic behaviour pattern (i.e., more or less standard response to a changes that is formed as result of previous experience) to implement bottom-up approach. Two components that can help to establish effective reaction regarding new initiatives in IT are proposed here: model of IT-related decision making, and efficiency measurement metric to estimate maturity of business processes and appropriate IT. Usage of proposed tools is demonstrated in practical cases.