Trajectories of microsecond molecular dynamics simulations of nucleosomes and nucleosome core particles.
Ontology highlight
ABSTRACT:
SUBMITTER: Shaytan AK
PROVIDER: S-EPMC4872717 | biostudies-literature | 2016 Jun
REPOSITORIES: biostudies-literature
Publications
Trajectories of microsecond molecular dynamics simulations of nucleosomes and nucleosome core particles.
Shaytan Alexey K AK Armeev Grigoriy A GA Goncearenco Alexander A Zhurkin Victor B VB Landsman David D Panchenko Anna R AR
Data in brief 20160506
We present here raw trajectories of molecular dynamics simulations for nucleosome with linker DNA strands as well as minimalistic nucleosome core particle model. The simulations were done in explicit solvent using CHARMM36 force field. We used this data in the research article Shaytan et al., 2016 [1]. The trajectory files are supplemented by TCL scripts providing advanced visualization capabilities. ...[more]
Similar Datasets
Project description:Molecular dynamics (MD) simulations of membranes are often hindered by the slow lateral diffusion of lipids and the limited time scale of MD. In order to study the dynamics of mixing and characterize the lateral distribution of lipids in converged mixtures, we report microsecond-long all-atom MD simulations performed on the special-purpose machine Anton. Two types of mixed bilayers, POPE:POPG (3:1) and POPC:cholesterol (2:1), as well as a pure POPC bilayer, were each simulated for up to 2 μs. These simulations show that POPE:POPG and POPC:cholesterol are each fully miscible at the simulated conditions, with the final states of the mixed bilayers similar to a random mixture. By simulating three POPE:POPG bilayers at different NaCl concentrations (0, 0.15, and 1 M), we also examined the effect of salt concentration on lipid mixing. While an increase in NaCl concentration is shown to affect the area per lipid, tail order, and lipid lateral diffusion, the final states of mixing remain unaltered, which is explained by the largely uniform increase in Na(+) ions around POPE and POPG. Direct measurement of water permeation reveals that the POPE:POPG bilayer with 1 M NaCl has reduced water permeability compared with those at zero or low salt concentration. Our calculations provide a benchmark to estimate the convergence time scale of all-atom MD simulations of lipid mixing. Additionally, equilibrated structures of POPE:POPG and POPC:cholesterol, which are frequently used to mimic bacterial and mammalian membranes, respectively, can be used as starting points of simulations involving these membranes.
| S-EPMC4196744 | biostudies-literature
Project description:Dissolution of many plant viruses is thought to start with swelling of the capsid caused by calcium removal following infection, but no high-resolution structures of swollen capsids exist. Here we have used microsecond all-atom molecular simulations to describe the dynamics of the capsid of satellite tobacco necrosis virus with and without the 92 structural calcium ions. The capsid expanded 2.5% upon removal of the calcium, in good agreement with experimental estimates. The water permeability of the native capsid was similar to that of a phospholipid membrane, but the permeability increased 10-fold after removing the calcium, predominantly between the 2-fold and 3-fold related subunits. The two calcium binding sites close to the icosahedral 3-fold symmetry axis were pivotal in the expansion and capsid-opening process, while the binding site on the 5-fold axis changed little structurally. These findings suggest that the dissociation of the capsid is initiated at the 3-fold axis.
| S-EPMC3349721 | biostudies-literature
Project description:DNA loop formation on nucleosomes is strongly implicated in chromatin remodeling and occurs spontaneously in nucleosomes subjected to superhelical stress. The nature of such loops depends crucially on the balance between DNA deformation and DNA interaction with the nucleosome core. Currently, no high-resolution structural data on these loops exist. Although uniform rod models have been used to study loop size and shape, these models make assumptions concerning DNA mechanics and DNA-core binding. We present here atomic-scale molecular dynamics simulations for two different loop sizes. The results point to the key role of localized DNA kinking within the loops. Kinks enable the relaxation of DNA bending strain to be coupled with improved DNA-core interactions. Kinks lead to small, irregularly shaped loops that are asymmetrically positioned with respect to the nucleosome core. We also find that loop position can influence the dynamics of the DNA segments at the extremities of the nucleosome.
| S-EPMC4914111 | biostudies-literature
Project description:Extensive all-atom molecular dynamics (?24 ?s total) allowed exploration of configurational space and calculation of lateral diffusion coefficients of the components of a protein-embedded, cholesterol-containing model bilayer. The three model membranes are composed of an ?50/50 (by mole) dipalmitoylphosphatidylcholine (DPPC)/cholesterol bilayer and contained an ?-helical transmembrane protein (HIV-1 gp41 TM). Despite the high concentration of cholesterol, normal Brownian motion was observed and the calculated diffusion coefficients (on the order of 10(-9) cm(2)/s) are consistent with experiments. Diffusion is sensitive to a variety of parameters, and a temperature difference of ?4 K from thermostat artifacts resulted in 2-10-fold differences in diffusion coefficients and significant differences in lipid order, membrane thickness, and unit cell area. Also, the specific peptide sequence likely underlies the consistently observed faster diffusion in one leaflet. Although the simulations here present molecular dynamics (MD) an order of magnitude longer than those from previous studies, the three systems did not approach ergodicity. The distributions of cholesterol and DPPC around the peptides changed on the microsecond time scale, but not significantly enough to thoroughly explore configurational space. These simulations support conclusions of other recent microsecond MD in that even longer time scales are needed for equilibration of model membranes and simulations of more realistic cellular or viral bilayers.
| S-EPMC4254001 | biostudies-literature
Project description:Acetylcholinesterase, with a deep, narrow active-site gorge, attracts enormous interest due to its particularly high catalytic efficiency and its inhibitors used for treatment of Alzheimer's disease. To facilitate the massive pass-through of the substrate and inhibitors, "breathing" motions to modulate the size of the gorge are an important prerequisite. However, the molecular mechanism that governs such motions is not well explored. Here, to systematically investigate intrinsic motions of the enzyme, we performed microsecond molecular dynamics simulations on the monomer and dimer of Torpedo californica acetylcholinesterase (TcAChE) as well as the complex of TcAChE bound with the drug E2020. It has been revealed that protein-ligand interactions and dimerization both keep the gorge in bulk, and opening events of the gorge increase dramatically compared to the monomer. Dynamics of three subdomains, S3, S4 and the Ω-loop, are tightly associated with variations of the gorge size while the dynamics can be changed by ligand binding or protein dimerization. Moreover, high correlations among these subdomains provide a basis for remote residues allosterically modulating the gorge motions. These observations are propitious to expand our understanding of protein structure and function as well as providing clues for performing structure-based drug design.
| S-EPMC5468367 | biostudies-literature
Project description:C5'-Hydrogen atoms are frequently abstracted during DNA oxidation. The oxidized abasic lesion 5'-(2-phosphoryl-1,4-dioxobutane) (DOB) is an electrophilic product of the C5'-radical. DOB is a potent irreversible inhibitor of DNA polymerase β, and forms interstrand cross-links in free DNA. We examined the reactivity of DOB within nucleosomes and nucleosome core particles (NCPs), the monomeric component of chromatin. Depending upon the position at which DOB is generated within a NCP, it is excised from nucleosomal DNA at a rate 275-1500-fold faster than that in free DNA. The half-life of DOB (7.0-16.8 min) in NCPs is shorter than any other abasic lesion. DOB's lifetime in NCPs is also significantly shorter than the estimated lifetime of an abasic site within a cell, suggesting that the observed chemistry would occur intracellularly. Histones also catalyze DOB excision when the lesion is present in the DNA linker region of a nucleosome. Schiff-base formation between DOB and histone proteins is detected in nucleosomes and NCPs, resulting in pyrrolone formation at the lysine residues. The lysines modified by DOB are often post-translationally modified. Consequently, the histone modifications described herein could affect the regulation of gene expression and may provide a chemical basis for the cytotoxicity of the DNA damaging agents that produce this lesion.
| S-EPMC4612368 | biostudies-other
Project description:We present an implementation of generalized Born implicit solvent all-atom classical molecular dynamics (MD) within the AMBER program package that runs entirely on CUDA enabled NVIDIA graphics processing units (GPUs). We discuss the algorithms that are used to exploit the processing power of the GPUs and show the performance that can be achieved in comparison to simulations on conventional CPU clusters. The implementation supports three different precision models in which the contributions to the forces are calculated in single precision floating point arithmetic but accumulated in double precision (SPDP), or everything is computed in single precision (SPSP) or double precision (DPDP). In addition to performance, we have focused on understanding the implications of the different precision models on the outcome of implicit solvent MD simulations. We show results for a range of tests including the accuracy of single point force evaluations and energy conservation as well as structural properties pertainining to protein dynamics. The numerical noise due to rounding errors within the SPSP precision model is sufficiently large to lead to an accumulation of errors which can result in unphysical trajectories for long time scale simulations. We recommend the use of the mixed-precision SPDP model since the numerical results obtained are comparable with those of the full double precision DPDP model and the reference double precision CPU implementation but at significantly reduced computational cost. Our implementation provides performance for GB simulations on a single desktop that is on par with, and in some cases exceeds, that of traditional supercomputers.
| S-EPMC3348677 | biostudies-literature
Project description:Among the multiple effects involved in chromatin condensation and decondensation processes, interactions between nucleosome core particles are suspected to play a crucial role. We analyze them in the absence of linker DNA and added proteins, after the self-assembly of isolated nucleosome core particles under controlled ionic conditions. We describe an original lamellar mesophase forming tubules on the mesoscopic scale. High resolution imaging of cryosections of vitrified samples reveals how nucleosome core particles stack on top of one another into columns which themselves align to form bilayers that repel one another through a solvent layer. We deduce from this structural organization how the particles interact through attractive interactions between top and bottom faces and lateral polar interactions that originate in the heterogeneous charge distribution at the surface of the particle. These interactions, at work under conditions comparable with those found in the living cell, should be of importance in the mechanisms governing chromatin compaction in vivo.
| S-EPMC1301712 | biostudies-other
Project description:The bacterial chaperone trigger factor (TF) is the first chaperone to be encountered by a nascent protein chain as it emerges from the ribosome exit tunnel. Experimental results suggest that TF possesses considerable conformational flexibility, and in an attempt to provide an atomic-level view of this flexibility, we have performed independent 1.5-μs molecular dynamics simulations of TF in explicit solvent using two different simulation force fields (OPLS-AA/L and AMBER ff99SB-ILDN). Both simulations indicate that TF possesses tremendous flexibility, with huge excursions from the crystallographic conformation caused by reorientations of the protein's constituent domains; both simulations also predict the formation of extensive contacts between TF's PPIase domain and the Arm 1 domain that is involved in nascent-chain binding. In the OPLS simulation, however, TF rapidly settles into a very compact conformation that persists for at least 1 μs, whereas in the AMBER simulation, it remains highly dynamic; additional simulations in which the two force fields were swapped suggest that these differences are at least partly attributable to sampling issues. The simulation results provide potential rationalizations of a number of experimental observations regarding TF's conformational behavior and have implications for using simulations to model TF's function on translating ribosomes.
| S-EPMC3736688 | biostudies-literature
Project description:Sphingosine 1-phosphate (S1P) is a lysophospholipid mediator which activates G protein-coupled sphingosine 1-phosphate receptors and thus evokes a variety of cell and tissue responses including lymphocyte trafficking, endothelial development, integrity, and maturation. We performed five all-atom 700 ns molecular dynamics simulations of the sphingosine 1-phosphate receptor 1 (S1P₁) based on recently released crystal structure of that receptor with an antagonist. We found that the initial movements of amino acid residues occurred in the area of highly conserved W269⁶·⁴⁸ in TM6 which is close to the ligand binding location. Those residues located in the central part of the receptor and adjacent to kinks of TM helices comprise of a transmission switch. Side chains movements of those residues were coupled to the movements of water molecules inside the receptor which helped in the gradual opening of intracellular part of the receptor. The most stable parts of the protein were helices TM1 and TM2, while the largest movement was observed for TM7, possibly due to the short intracellular part starting with a helix kink at P⁷·⁵⁰, which might be the first helix to move at the intracellular side. We show for the first time the detailed view of the concerted action of the transmission switch and Trp (W⁶·⁴⁸) rotamer toggle switch leading to redirection of water molecules flow in the central part of the receptor. That event is a prerequisite for subsequent changes in intracellular part of the receptor involving water influx and opening of the receptor structure.
| S-EPMC3789783 | biostudies-literature