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.

Specifying conformational heterogeneity of multi-domain proteins at atomic resolution.

The conformational landscape of multi-domain proteins is inherently linked to their specific functions. This also holds for polyubiquitin chains that are assembled by two or more ubiquitin domains connected by a flexible linker thus showing a large interdomain mobility. However, molecular recognition and signal transduction are associated with particular conformational substates that are populated in solution. Here, we apply high-resolution NMR spectroscopy in combination with dual-scale MD simulations to explore the conformational space of K6-, K29-, and K33-linked diubiquitin molecules. The conformational ensembles are evaluated utilizing a paramagnetic cosolute reporting on solvent exposure plus a set of complementary NMR parameters. This approach unravels a conformational heterogeneity of diubiquitins and explains the diversity of structural models that have been determined for K6-, K29-, and K33-linked diubiquitins in free and ligand-bound states so far. We propose a general application of the approach developed here to demystify multi-domain proteins occurring in nature.

Multiscale simulations of protein and membrane systems.

Classical multiscale simulations are perfectly suited to investigate biological soft matter systems. Owing to the bridging between microscopically realistic and lower-resolution models or the integration of a hierarchy of subsystems, one gets access to biologically relevant system sizes and timescales. In recent years, increasingly complex systems and processes have come into focus such as multidomain proteins, phase separation processes in biopolymer solutions, multicomponent biomembranes, or multiprotein complexes up to entire viruses. The review shows factors that have contributed to this progress - from improved models to machine-learning-based analysis and scale-bridging methods.

A Ligand Selection Strategy Identifies Chemical Probes Targeting the Proteases of SARS-CoV-2.

Activity-based probes are valuable tools for chemical biology. However, finding probes that specifically target the active site of an enzyme remains a challenging task. Herein, we present a ligand selection strategy that allows to rapidly tailor electrophilic probes to a target of choice and showcase its application for the two cysteine proteases of SARS-CoV-2 as proof of concept. The resulting probes were specific for the active site labeling of 3CLpro and PLpro with sufficient selectivity in a live cell model as well as in the background of a native human proteome. Exploiting the probes as tools for competitive profiling of a natural product library identified salvianolic acid derivatives as promising 3CLpro inhibitors. We anticipate that our ligand selection strategy will be useful to rapidly develop customized probes and discover inhibitors for a wide range of target proteins also beyond corona virus proteases.

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.