Ensembles of Breathing Nucleosomes: A Computational Study.
Ontology highlight
ABSTRACT:
SUBMITTER: van Deelen K
PROVIDER: S-EPMC7202939 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Genomic DNA is packaged in chromatin, a dynamic fiber variable in size and compaction. In chromatin, repeating nucleosome units wrap 145-147 DNA basepairs around histone proteins. Genetic and epigenetic regulation of genes relies on structural transitions in chromatin which are driven by intra- and inter-nucleosome dynamics and modulated by chemical modifications of the unstructured terminal tails of histones. Here we demonstrate how the interplay between histone H3 and H2A tails control ample nucleosome breathing motions. We monitored large openings of two genomic nucleosomes, and only moderate breathing of an engineered nucleosome in atomistic molecular simulations amounting to 24 μs. Transitions between open and closed nucleosome conformations were mediated by the displacement and changes in compaction of the two histone tails. These motions involved changes in the DNA interaction profiles of clusters of epigenetic regulatory aminoacids in the tails. Removing the histone tails resulted in a large increase of the amplitude of nucleosome breathing but did not change the sequence dependent pattern of the motions. Histone tail modulated nucleosome breathing is a key mechanism of chromatin dynamics with important implications for epigenetic regulation.
| S-EPMC8174689 | biostudies-literature
Project description:Prediction problems in biomedical sciences are generally quite difficult, partially due to incomplete knowledge of how the phenomenon of interest is influenced by the variables and measurements used for prediction, as well as a lack of consensus regarding the ideal predictor(s) for specific problems. In these situations, a powerful approach to improving prediction performance is to construct ensembles that combine the outputs of many individual base predictors, which have been successful for many biomedical prediction tasks. Moreover, selecting a parsimonious ensemble can be of even greater value for biomedical sciences, where it is not only important to learn an accurate predictor, but also to interpret what novel knowledge it can provide about the target problem. Ensemble selection is a promising approach for this task because of its ability to select a collectively predictive subset, often a relatively small one, of all input base predictors. One of the most well-known algorithms for ensemble selection, CES (Caruana et al.'s Ensemble Selection), generally performs well in practice, but faces several challenges due to the difficulty of choosing the right values of its various parameters. Since the choices made for these parameters are usually ad-hoc, good performance of CES is difficult to guarantee for a variety of problems or datasets. To address these challenges with CES and other such algorithms, we propose a novel heterogeneous ensemble selection approach based on the paradigm of reinforcement learning (RL), which offers a more systematic and mathematically sound methodology for exploring the many possible combinations of base predictors that can be selected into an ensemble. We develop three RL-based strategies for constructing ensembles and analyze their results on two unbalanced computational genomics problems, namely the prediction of protein function and splice sites in eukaryotic genomes. We show that the resultant ensembles are indeed substantially more parsimonious as compared to the full set of base predictors, yet still offer almost the same classification power, especially for larger datasets. The RL ensembles also yield a better combination of parsimony and predictive performance as compared to CES.
| S-EPMC5147733 | biostudies-literature
Project description:Although poorly positioned nucleosomes are ubiquitous in the eukaryotic genome, they are difficult to identify with existing nucleosome identification methods. Recently available enhanced high-throughput chromatin conformation capture techniques such as Micro-C, DNase Hi-C, and Hi-CO characterize nucleosome-level chromatin proximity, probing the positions of mono-nucleosomes and the spacing between nucleosome pairs at the same time, enabling nucleosome profiling in poorly positioned regions. Here we develop a novel computational approach, NucleoMap, to identify nucleosome positioning from ultra-high resolution chromatin contact maps. By integrating nucleosome read density, contact distances, and binding preferences, NucleoMap precisely locates nucleosomes in both prokaryotic and eukaryotic genomes and outperforms existing nucleosome identification methods in both precision and recall. We rigorously characterize genome-wide association in eukaryotes between the spatial organization of mono-nucleosomes and their corresponding histone modifications, protein binding activities, and higher-order chromatin functions. We also find evidence of two tetra-nucleosome folding structures in human embryonic stem cells and analyze their association with multiple structural and functional regions. Based on the identified nucleosomes, nucleosome contact maps are constructed, reflecting the inter-nucleosome distances and preserving the contact distance profiles in original contact maps.
| S-EPMC9321407 | biostudies-literature
Project description:The human nose serves vital physiological functions, including warming, filtration, humidification, and olfaction. These functions are based on transport phenomena that depend on nasal airflow patterns and turbulence. Accurate prediction of these airflow properties requires careful selection of computational fluid dynamics models and rigorous validation. The validation studies in the past have been limited by poor representations of the complex nasal geometry, lack of detailed airflow comparisons, and restricted ranges of flow rate. The objective of this study is to validate various numerical methods based on an anatomically accurate nasal model against published experimentally measured data under breathing flow rates from 180 to 1100ml/s. The numerical results of velocity profiles and turbulence intensities were obtained using the laminar model, four widely used Reynolds-averaged Navier-Stokes (RANS) turbulence models (i.e., k-?, standard k-?, Shear Stress Transport k-?, and Reynolds Stress Model), large eddy simulation (LES) model, and direct numerical simulation (DNS). It was found that, despite certain irregularity in the flow field, the laminar model achieved good agreement with experimental results under restful breathing condition (180ml/s) and performed better than the RANS models. As the breathing flow rate increased, the RANS models achieved more accurate predictions but still performed worse than LES and DNS. As expected, LES and DNS can provide accurate predictions of the nasal airflow under all flow conditions but have an approximately 100-fold higher computational cost. Among all the RANS models tested, the standard k-? model agrees most closely with the experimental values in terms of velocity profile and turbulence intensity.
| S-EPMC5694356 | biostudies-literature
Project description:The characterization of intrinsically disordered protein (IDP) ensembles is complicated both by inherent heterogeneity and by the fact that many common experimental techniques function poorly when applied to IDPs. For this reason, the development of alternative structural tools for probing IDP ensembles has attracted considerable attention. Here we describe our recent work in developing experimental and computational tools for characterizing IDP ensembles using Amide I (backbone carbonyl stretch) vibrational spectroscopy. In this approach, the infrared (IR) absorption frequencies of isotope-labeled amide bonds probe their local electrostatic environments and structures. Empirical frequency maps allow us to use this spectroscopic data as a direct experimental test of atomistic structural models. We apply these methods to a family of short elastin-like peptides (ELPs), fragments of the elastin protein based around the Pro-Gly turn motif characteristic of the elastomeric segments of the full protein. Using a maximum entropy analysis that applies constraints from experimental spectra to weighting predicted spectra from molecular dynamics (MD) ensembles, we find that peptides with Ala or Val side chains preceding the Pro-Gly turn unit exhibit a stronger tendency toward extended structures than do Gly-Pro-Gly motifs, suggesting an important role for steric interactions in tuning the molecular properties of elastin.
| S-EPMC5312701 | biostudies-literature
Project description:The stability, activity, and solubility of a protein sequence are determined by a delicate balance of molecular interactions in a variety of conformational states. Even so, most computational protein design methods model sequences in the context of a single native conformation. Simulations that model the native state as an ensemble have been mostly neglected due to the lack of sufficiently powerful optimization algorithms for multistate design. Here, we have applied our multistate design algorithm to study the potential utility of various forms of input structural data for design. To facilitate a more thorough analysis, we developed new methods for the design and high-throughput stability determination of combinatorial mutation libraries based on protein design calculations. The application of these methods to the core design of a small model system produced many variants with improved thermodynamic stability and showed that multistate design methods can be readily applied to large structural ensembles. We found that exhaustive screening of our designed libraries helped to clarify several sources of simulation error that would have otherwise been difficult to ascertain. Interestingly, the lack of correlation between our simulated and experimentally measured stability values shows clearly that a design procedure need not reproduce experimental data exactly to achieve success. This surprising result suggests potentially fruitful directions for the improvement of computational protein design technology.
| S-EPMC2993350 | biostudies-literature
Project description:BackgroundThere is on-going controversy regarding the potential for increased respiratory effort to generate patient self-inflicted lung injury (P-SILI) in spontaneously breathing patients with COVID-19 acute hypoxaemic respiratory failure. However, direct clinical evidence linking increased inspiratory effort to lung injury is scarce. We adapted a computational simulator of cardiopulmonary pathophysiology to quantify the mechanical forces that could lead to P-SILI at different levels of respiratory effort. In accordance with recent data, the simulator parameters were manually adjusted to generate a population of 10 patients that recapitulate clinical features exhibited by certain COVID-19 patients, i.e., severe hypoxaemia combined with relatively well-preserved lung mechanics, being treated with supplemental oxygen.ResultsSimulations were conducted at tidal volumes (VT) and respiratory rates (RR) of 7 ml/kg and 14 breaths/min (representing normal respiratory effort) and at VT/RR of 7/20, 7/30, 10/14, 10/20 and 10/30 ml/kg / breaths/min. While oxygenation improved with higher respiratory efforts, significant increases in multiple indicators of the potential for lung injury were observed at all higher VT/RR combinations tested. Pleural pressure swing increased from 12.0 ± 0.3 cmH2O at baseline to 33.8 ± 0.4 cmH2O at VT/RR of 7 ml/kg/30 breaths/min and to 46.2 ± 0.5 cmH2O at 10 ml/kg/30 breaths/min. Transpulmonary pressure swing increased from 4.7 ± 0.1 cmH2O at baseline to 17.9 ± 0.3 cmH2O at VT/RR of 7 ml/kg/30 breaths/min and to 24.2 ± 0.3 cmH2O at 10 ml/kg/30 breaths/min. Total lung strain increased from 0.29 ± 0.006 at baseline to 0.65 ± 0.016 at 10 ml/kg/30 breaths/min. Mechanical power increased from 1.6 ± 0.1 J/min at baseline to 12.9 ± 0.2 J/min at VT/RR of 7 ml/kg/30 breaths/min, and to 24.9 ± 0.3 J/min at 10 ml/kg/30 breaths/min. Driving pressure increased from 7.7 ± 0.2 cmH2O at baseline to 19.6 ± 0.2 cmH2O at VT/RR of 7 ml/kg/30 breaths/min, and to 26.9 ± 0.3 cmH2O at 10 ml/kg/30 breaths/min.ConclusionsOur results suggest that the forces generated by increased inspiratory effort commonly seen in COVID-19 acute hypoxaemic respiratory failure are comparable with those that have been associated with ventilator-induced lung injury during mechanical ventilation. Respiratory efforts in these patients should be carefully monitored and controlled to minimise the risk of lung injury.
| S-EPMC8276227 | biostudies-literature
Project description:Formaldehyde has long been recognized as a hazardous environmental agent highly reactive with DNA. Recently, it has been realized that due to the activity of histone demethylation enzymes within the cell nucleus, formaldehyde is produced endogenously, in direct vicinity of genomic DNA. Should it lead to extensive DNA damage? We address this question with the aid of a computational mapping method, analogous to X-ray and nuclear magnetic resonance techniques for observing weakly specific interactions of small organic compounds with a macromolecule in order to establish important functional sites. We concentrate on the leading reaction of formaldehyde with free bases: hydroxymethylation of cytosine amino groups. Our results show that in B-DNA, cytosine amino groups are totally inaccessible for the formaldehyde attack. Then, we explore the effect of recently discovered transient flipping of Watson-Crick (WC) pairs into Hoogsteen (HG) pairs (HG breathing). Our results show that the HG base pair formation dramatically affects the accessibility for formaldehyde of cytosine amino nitrogens within WC base pairs adjacent to HG base pairs. The extensive literature on DNA interaction with formaldehyde is analyzed in light of the new findings. The obtained data emphasize the significance of DNA HG breathing.
| S-EPMC3439909 | biostudies-literature
Project description:Accommodating backbone flexibility continues to be the most difficult challenge in computational docking of protein-protein complexes. Towards that end, we simulate four distinct biophysical models of protein binding in RosettaDock, a multiscale Monte-Carlo-based algorithm that uses a quasi-kinetic search process to emulate the diffusional encounter of two proteins and to identify low-energy complexes. The four binding models are as follows: (1) key-lock (KL) model, using rigid-backbone docking; (2) conformer selection (CS) model, using a novel ensemble docking algorithm; (3) induced fit (IF) model, using energy-gradient-based backbone minimization; and (4) combined conformer selection/induced fit (CS/IF) model. Backbone flexibility was limited to the smaller partner of the complex, structural ensembles were generated using Rosetta refinement methods, and docking consisted of local perturbations around the complexed conformation using unbound component crystal structures for a set of 21 target complexes. The lowest-energy structure contained >30% of the native residue-residue contacts for 9, 13, 13, and 14 targets for KL, CS, IF, and CS/IF docking, respectively. When applied to 15 targets using nuclear magnetic resonance ensembles of the smaller protein, the lowest-energy structure recovered at least 30% native residue contacts in 3, 8, 4, and 8 targets for KL, CS, IF, and CS/IF docking, respectively. CS/IF docking of the nuclear magnetic resonance ensemble performed equally well or better than KL docking with the unbound crystal structure in 10 of 15 cases. The marked success of CS and CS/IF docking shows that ensemble docking can be a versatile and effective method for accommodating conformational plasticity in docking and serves as a demonstration for the CS theory--that binding-competent conformers exist in the unbound ensemble and can be selected based on their favorable binding energies.
| S-EPMC2573042 | biostudies-literature
Project description:Micrococcal nuclease (MNase) is widely used to map nucleosomes. However, nucleosomes are highly dynamic and susceptible to experimental conditions, resulting in extreme variability across nucleosome maps, which complicates the generation of accurate nucleosome organization data. We mapped nucleosomes from different individuals using improved MNase-seq. The improvements included setting different digestion levels (low, medium, high) and naked DNA correction to remove the noise caused by experimental manipulation and comparing maps to obtain the accurate position and occupancy of strong nucleosomes (SNs) in the whole genome. In addition, the characteristics of SNs were further excavated. SNs were enriched in Alu elements and near the centromere of Chr12. SNs contain some specific sequences, and the GC content of SNs is different from that of dynamic nucleosomes. The findings suggest that nucleosome location in the genome and the DNA sequence may affect nucleosome stability.
| S-EPMC9249913 | biostudies-literature