Combining molecular dynamics simulations and scoring method to computationally model ubiquitylated linker histones in chromatosomes.

The chromatin in eukaryotic cells plays a fundamental role in all processes during a cell's life cycle. This nucleoprotein is normally tightly packed but needs to be unpacked for expression and division. The linker histones are critical for such packaging processes and while most experimental and simulation works recognize their crucial importance, the focus is nearly always set on the nucleosome as the basic chromatin building block. Linker histones can undergo several modifications, but only few studies on their ubiquitylation have been conducted. Mono-ubiquitylated linker histones (HUb), while poorly understood, are expected to influence DNA compaction. The size of ubiquitin and the globular domain of the linker histone are comparable and one would expect an increased disorder upon ubiquitylation of the linker histone. However, the formation of higher order chromatin is not hindered and ubiquitylation of the linker histone may even promote gene expression. Structural data on chromatosomes is rare and HUb has never been modeled in a chromatosome so far. Descriptions of the chromatin complex with HUb would greatly benefit from computational structural data. In this study we generate molecular dynamics simulation data for six differently linked HUb variants with the help of a sampling scheme tailored to drive the exploration of phase space. We identify conformational sub-states of the six HUb variants using the sketch-map algorithm for dimensionality reduction and iterative HDBSCAN for clustering on the excessively sampled, shallow free energy landscapes. We present a highly efficient geometric scoring method to identify sub-states of HUb that fit into the nucleosome. We predict HUb conformations inside a nucleosome using on-dyad and off-dyad chromatosome structures as reference and show that unbiased simulations of HUb produce significantly more fitting than non-fitting HUb conformations. A tetranucleosome array is used to show that ubiquitylation can even occur in chromatin without too much steric clashes.

Deep Learning to Analyze Sliding Drops.

State-of-the-art contact angle measurements usually involve image analysis of sessile drops. The drops are symmetric and images can be taken at high resolution. The analysis of videos of drops sliding down a tilted plate is hampered due to the low resolution of the cutout area where the drop is visible. The challenge is to analyze all video images automatically, while the drops are not symmetric anymore and contact angles change while sliding down the tilted plate. To increase the accuracy of contact angles, we present a 4-segment super-resolution optimized-fitting (4S-SROF) method. We developed a deep learning-based super-resolution model with an upscale ratio of 3; i.e., the trained model is able to enlarge drop images 9 times accurately (PSNR = 36.39). In addition, a systematic experiment using synthetic images was conducted to determine the best parameters for polynomial fitting of contact angles. Our method improved the accuracy by 21% for contact angles lower than 90° and by 33% for contact angles higher than 90°.

Determining glass transition in all-atom acrylic polymeric melt simulations using machine learning.

The functionality of many polymeric materials depends on their glass transition temperatures (Tg). In computer simulations, Tg is often calculated from the gradual change in macroscopic properties. Precise determination of this change depends on the fitting protocols. We previously proposed a robust data-driven approach to determine Tg from the molecular dynamics simulation data of a coarse-grained semiflexible polymer model. In contrast to the global macroscopic properties, our method relies on high-resolution microscopic details. Here, we demonstrate the generality of our approach by using various dimensionality reduction and clustering methods and apply it to an atomistic model of acrylic polymers. Our study reveals the explicit contribution of the side chain and backbone residues in influencing the determination of the glass transition temperature.

Fast conformational clustering of extensive molecular dynamics simulation data.

We present an unsupervised data processing workflow that is specifically designed to obtain a fast conformational clustering of long molecular dynamics simulation trajectories. In this approach, we combine two dimensionality reduction algorithms (cc_analysis and encodermap) with a density-based spatial clustering algorithm (hierarchical density-based spatial clustering of applications with noise). The proposed scheme benefits from the strengths of the three algorithms while avoiding most of the drawbacks of the individual methods. Here, the cc_analysis algorithm is applied for the first time to molecular simulation data. The encodermap algorithm complements cc_analysis by providing an efficient way to process and assign large amounts of data to clusters. The main goal of the procedure is to maximize the number of assigned frames of a given trajectory while keeping a clear conformational identity of the clusters that are found. In practice, we achieve this by using an iterative clustering approach and a tunable root-mean-square-deviation-based criterion in the final cluster assignment. This allows us to find clusters of different densities and different degrees of structural identity. With the help of four protein systems, we illustrate the capability and performance of this clustering workflow: wild-type and thermostable mutant of the Trp-cage protein (TC5b and TC10b), NTL9, and Protein B. Each of these test systems poses their individual challenges to the scheme, which, in total, give a nice overview of the advantages and potential difficulties that can arise when using the proposed method.

Guanidine-II aptamer conformations and ligand binding modes through the lens of molecular simulation.

Regulation of gene expression via riboswitches is a widespread mechanism in bacteria. Here, we investigate ligand binding of a member of the guanidine sensing riboswitch family, the guanidine-II riboswitch (Gd-II). It consists of two stem-loops forming a dimer upon ligand binding. Using extensive molecular dynamics simulations we have identified conformational states corresponding to ligand-bound and unbound states in a monomeric stem-loop of Gd-II and studied the selectivity of this binding. To characterize these states and ligand-dependent conformational changes we applied a combination of dimensionality reduction, clustering, and feature selection methods. In absence of a ligand, the shape of the binding pocket alternates between the conformation observed in presence of guanidinium and a collapsed conformation, which is associated with a deformation of the dimerization interface. Furthermore, the structural features responsible for the ability to discriminate against closely related analogs of guanidine are resolved. Based on these insights, we propose a mechanism that couples ligand binding to aptamer dimerization in the Gd-II system, demonstrating the value of computational methods in the field of nucleic acids research.

Towards a molecular basis of ubiquitin signaling: A dual-scale simulation study of ubiquitin dimers.

Covalent modification of proteins by ubiquitin or ubiquitin chains is one of the most prevalent post-translational modifications in eukaryotes. Different types of ubiquitin chains are assumed to selectively signal respectively modified proteins for different fates. In support of this hypothesis, structural studies have shown that the eight possible ubiquitin dimers adopt different conformations. However, at least in some cases, these structures cannot sufficiently explain the molecular basis of the selective signaling mechanisms. This indicates that the available structures represent only a few distinct conformations within the entire conformational space adopted by a ubiquitin dimer. Here, molecular simulations on different levels of resolution can complement the structural information. We have combined exhaustive coarse grained and atomistic simulations of all eight possible ubiquitin dimers with a suitable dimensionality reduction technique and a new method to characterize protein-protein interfaces and the conformational landscape of protein conjugates. We found that ubiquitin dimers exhibit characteristic linkage type-dependent properties in solution, such as interface stability and the character of contacts between the subunits, which can be directly correlated with experimentally observed linkage-specific properties.

Back-mapping based sampling: Coarse grained free energy landscapes as a guideline for atomistic exploration.

One ongoing topic of research in MD simulations is how to enable sampling to chemically and biologically relevant time scales. We address this question by introducing a back-mapping based sampling (BMBS) that combines multiple aspects of different sampling techniques. BMBS uses coarse grained (CG) free energy surfaces (FESs) and dimensionality reduction to initiate new atomistic simulations. These new simulations are started from atomistic conformations that were back-mapped from CG points all over the FES in order to sample the entire accessible phase space as fast as possible. In the context of BMBS, we address relevant back-mapping related questions like where to start the back-mapping from and how to judge the atomistic ensemble that results from the BMBS. The latter is done with the use of the earth mover's distance, which allows us to quantitatively compare distributions of CG and atomistic ensembles. By using this metric, we can also show that the BMBS is able to correct inaccuracies of the CG model. In this paper, BMBS is applied to a just recently introduced neural network (NN) based approach for a radical coarse graining to predict free energy surfaces for oligopeptides. The BMBS scheme back-maps these FESs to the atomistic scale, justifying and complementing the proposed NN based CG approach. The efficiency benefit of the algorithm scales with the length of the oligomer. Already for the heptamers, the algorithm is about one order of magnitude faster in sampling compared to a standard MD simulation.

Estimating sliding drop width via side-view features using recurrent neural networks.

High speed side-view videos of sliding drops enable researchers to investigate drop dynamics and surface properties. However, understanding the physics of sliding requires knowledge of the drop width. A front-view perspective of the drop is necessary. In particular, the drop's width is a crucial parameter owing to its association with the friction force. Incorporating extra cameras or mirrors to monitor changes in the width of drops from a front-view perspective is cumbersome and limits the viewing area. This limitation impedes a comprehensive analysis of sliding drops, especially when they interact with surface defects. Our study explores the use of various regression and multivariate sequence analysis (MSA) models to estimate the drop width at a solid surface solely from side-view videos. This approach eliminates the need to incorporate additional equipment into the experimental setup. In addition, it ensures an unlimited viewing area of sliding drops. The Long Short Term Memory (LSTM) model with a 20 sliding window size has the best performance with the lowest root mean square error (RMSE) of 67 µm. Within the spectrum of drop widths in our dataset, ranging from 1.6 to 4.4 mm, this RMSE indicates that we can predict the width of sliding drops with an error of 2.4%. Furthermore, the applied LSTM model provides a drop width across the whole sliding length of 5 cm, previously unattainable.

Data-Driven Identification and Analysis of the Glass Transition in Polymer Melts.

Understanding the nature of glass transition, as well as the precise estimation of the glass transition temperature for polymeric materials, remains open questions in both experimental and theoretical polymer sciences. We propose a data-driven approach, which utilizes the high-resolution details accessible through the molecular dynamics simulation and considers the structural information on individual chains. It clearly identifies the glass transition temperature of polymer melts of weakly semiflexible chains. By combining principal component analysis and clustering, we identify the glass transition temperature in the asymptotic limit even from relatively short time trajectories, which just reach into the Rouse-like monomer displacement regime. We demonstrate that fluctuations captured by the principal component analysis reflect the change in a chain's behavior: from conformational rearrangement above to small fluctuations below the glass transition temperature. Our approach is straightforward to apply and should be applicable to other polymeric glass-forming liquids.

Generating a conformational landscape of ubiquitin chains at atomistic resolution by back-mapping based sampling.

Ubiquitin chains are flexible multidomain proteins that have important biological functions in cellular signalling. Computational studies with all-atom molecular dynamics simulations of the conformational spaces of polyubiquitins can be challenging due to the system size and a multitude of long-lived meta-stable states. Coarse graining is an efficient approach to overcome this problem-at the cost of losing high-resolution details. Recently, we proposed the back-mapping based sampling (BMBS) approach that reintroduces atomistic information into a given coarse grained (CG) sampling based on a two-dimensional (2D) projection of the conformational landscape, produces an atomistic ensemble and allows to systematically compare the ensembles at the two levels of resolution. Here, we apply BMBS to K48-linked tri-ubiquitin, showing its applicability to larger systems than those it was originally introduced on and demonstrating that the algorithm scales very well with system size. In an extension of the original BMBS we test three different seeding strategies, i.e. different approaches from where in the CG landscape atomistic trajectories are initiated. Furthermore, we apply a recently introduced conformational clustering algorithm to the back-mapped atomistic ensemble. Thus, we obtain insight into the structural composition of the 2D landscape and illustrate that the dimensionality reduction algorithm separates different conformational characteristics very well into different regions of the map. This cluster analysis allows us to show how atomistic trajectories sample conformational states, move through the projection space and in sum converge to an atomistic conformational landscape that slightly differs from the original CG map, indicating a correction of flaws in the CG template.

Efficient Sampling and Characterization of Free Energy Landscapes of Ion-Peptide Systems.

Proteins that influence nucleation, growth, or polymorph selection during biomineralization processes are often rich in glutamic- or aspartic acid. Here, the interactions between carboxylate side chains and ions lead to an interplay of peptide conformations and ion structuring in solution. Molecular dynamics simulations are an ideal tool to mechanistically investigate these processes. Unfortunately, the formation of strong ion-peptide contacts and ion bridges drastically impedes structural reorganization of ionic bonds and conformational transitions of the polymers. Thus, to obtain a complete thermodynamical picture of such systems, enhanced sampling techniques become necessary as well as the methods to characterize the conformational states of these partially disordered polymer-ion systems. Here, we propose a new set of Hamiltonian replica exchange (HRE) parameters for efficient simulations of peptide-ion systems, with an aspartic acid trimer in the presence of Ca2+ and Cl- ions as a test system. We introduce dimensionality reduction and clustering strategies to characterize the states of such a multicomponent system and to analyze the outcome of the proposed HRE with different reweighting methods.

Using Dimensionality Reduction to Systematically Expand Conformational Sampling of Intrinsically Disordered Peptides.

One of the approaches to improve our ability to characterize biologically important processes and to map out an underlying free energy landscape is to direct MD simulations to explore molecular conformational phase space faster. Intrinsically disordered systems with shallow free energy landscapes of a huge number of metastable minima pose a particular challenge in this regard. Both characterization of the often ill-defined conformational states as well as the assessment of the degree of convergence of phase space exploration are problematic. We have used a multidimensional scaling-like embedding (sketch-map) to describe the energetically accessible regions of phase space for a peptide fragment of the intrinsically disordered protein α-synuclein. Using sketch-map coordinates from a short initial simulation, we guided additional MD simulations to efficiently expand sampling of the conformational space. The sketch-map projections are very well suited to detect rare but possibly functionally relevant events, metastable intermediates, and transition states in the vast amount of data.