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.

Titin kinase ubiquitination aligns autophagy receptors with mechanical signals in the sarcomere.

Striated muscle undergoes remodelling in response to mechanical and physiological stress, but little is known about the integration of such varied signals in the myofibril. The interaction of the elastic kinase region from sarcomeric titin (A168-M1) with the autophagy receptors Nbr1/p62 and MuRF E3 ubiquitin ligases is well suited to link mechanosensing with the trophic response of the myofibril. To investigate the mechanisms of signal cross-talk at this titin node, we elucidated its 3D structure, analysed its response to stretch using steered molecular dynamics simulations and explored its functional relation to MuRF1 and Nbr1/p62 using cellular assays. We found that MuRF1-mediated ubiquitination of titin kinase promotes its scaffolding of Nbr1/p62 and that the process can be dynamically down-regulated by the mechanical unfolding of a linker sequence joining titin kinase with the MuRF1 receptor site in titin. We propose that titin ubiquitination is sensitive to the mechanical state of the sarcomere, the regulation of sarcomere targeting by Nbr1/p62 being a functional outcome. We conclude that MuRF1/Titin Kinase/Nbr1/p62 constitutes a distinct assembly that predictably promotes sarcomere breakdown in inactive muscle.

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.

On the Binding Mechanisms of Calcium Ions to Polycarboxylates: Effects of Molecular Weight, Side Chain, and Backbone Chemistry.

We experimentally determined the characteristics and Langmuir parameters of the binding of calcium ions to different polycarboxylates. By using potentiometric titrations and isothermal titration calorimetry, the effects of side chain chemistry, pH value, and chain length were systematically investigated using the linear polymers poly(aspartic acid), poly(glutamic acid), and poly(acrylic acid). We demonstrate that for polymers with high polymerization degrees, the binding process is governed by higher-order effects, such as the change of apparent pKa of carboxyl groups, and contributions arising from the whole polymer chain while the chemistry of the monomer unit constituting the polymer plays a subordinate role. In addition, primary binding sites need to be present in the polymer, thus rendering the abundance and sequential arrangement of protonated and deprotonated groups important. The detection of higher-order effects contradicts the assumptions posed by the Langmuir model of noninteracting binding sites and puts a question mark on whether ion binding to polycarboxylates can be described using solely a Langmuir binding model. No single uniform mechanism fits all investigated systems, and the whole polymer chain, including terminal groups, needs to be considered for the interpretation of binding data. Therefore, one needs to be careful when explaining ion binding to polymers solely based on studies on monomers or oligomers.

Solvent-mediated isotope effects strongly influence the early stages of calcium carbonate formation: exploring D2O vs. H2O in a combined computational and experimental approach.

In experimental studies, heavy water (D2O) is employed, e.g., so as to shift the spectroscopic solvent background, but any potential effects of this solvent exchange on reaction pathways are often neglected. While the important role of light water (H2O) during the early stages of calcium carbonate formation has been realized, studies into the actual effects of aqueous solvent exchanges are scarce. Here, we present a combined computational and experimental approach to start to fill this gap. We extended a suitable force field for molecular dynamics (MD) simulations. Experimentally, we utilised advanced titration assays and time-resolved attenuated total reflection Fourier transform infrared (ATR-FTIR) spectroscopy. We find distinct effects in various mixtures of the two aqueous solvents, and in pure H2O or D2O. Disagreements between the computational results and experimental data regarding the stabilities of ion associates might be due to the unexplored role of HDO, or an unprobed complex phase behaviour of the solvent mixtures in the simulations. Altogether, however, our data suggest that calcium carbonate formation might proceed "more classically" in D2O. Also, there are indications for the formation of new structures in amorphous and crystalline calcium carbonates. There is huge potential towards further improving the understanding of mineralization mechanisms by studying solvent-mediated isotope effects, also beyond calcium carbonate. Last, it must be appreciated that H2O and D2O have significant, distinct effects on mineralization mechanisms, and that care has to be taken when experimental data from D2O studies are used, e.g., for the development of H2O-based computer models.

Interactions Determining the Structural Integrity of the Trimer of Plant Light Harvesting Complex in Lipid Membranes.

The structural basis for the stability of the trimeric form of the light harvesting complex (LHCII), a pigmented protein from green plants pivotal for photosynthesis, remains elusive till date. The protein embedded in a dipalmitoylphosphatidylcholine (DPPC) lipid membrane is investigated using all-atom molecular dynamics simulations to find out the interactions responsible for the structural integrity of the trimer and its relation to antenna function. Central association of chlorophyll a (CLA) molecules near the LHCII chains is attributed to a conserved coordination between the Mg of CLA and the oxygen of a specific residue of the first helix of a chain. The residue forms a salt-bridge with the fourth helix of the same chain of the trimer, not of the monomer. In an earlier experiment, three residues (WYR) at each chain of the trimer have been found indispensable for the trimerization and referred to as trimerization motif. We find that the residues of the trimerization motif are connected to the lipids or pigments by a chain of interactions rather than a direct contact. Synergistic effects of sequentially located hydrogen bonds and salt-bridges within monomers of the trimer keep the trimer conformation stable in association with the pigments or the lipids. These interactions are exclusively present in the pigmented trimer and not present in the monomer or in the unpigmented trimer. Thus, our results provide a molecular basis for the inherent stability of the LHCII trimer in a lipid membrane and explain many pre-existing experimental data.

Snapshots of a modified nucleotide moving through the confines of a DNA polymerase.

DNA polymerases have evolved to process the four canonical nucleotides accurately. Nevertheless, these enzymes are also known to process modified nucleotides, which is the key to numerous core biotechnology applications. Processing of modified nucleotides includes incorporation of the modified nucleotide and postincorporation elongation to proceed with the synthesis of the nascent DNA strand. The structural basis for postincorporation elongation is currently unknown. We addressed this issue and successfully crystallized KlenTaq DNA polymerase in six closed ternary complexes containing the enzyme, the modified DNA substrate, and the incoming nucleotide. Each structure shows a high-resolution snapshot of the elongation of a modified primer, where the modification "moves" from the 3'-primer terminus upstream to the sixth nucleotide in the primer strand. Combining these data with quantum mechanics/molecular mechanics calculations and biochemical studies elucidates how the enzyme and the modified substrate mutually modulate their conformations without compromising the enzyme's activity significantly. The study highlights the plasticity of the system as origin of the broad substrate properties of DNA polymerases and facilitates the design of improved systems.

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.

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.

EncoderMap(II): Visualizing Important Molecular Motions with Improved Generation of Protein Conformations.

Dimensionality reduction can be used to project high-dimensional molecular data into a simplified, low-dimensional map. One feature of our recently introduced dimensionality reduction technique EncoderMap, which relies on the combination of an autoencoder with multidimensional scaling, is its ability to do the reverse. It is able to generate conformations for any selected points in the low-dimensional map. This transfers the simplified, low-dimensional map back into the high-dimensional conformational space. Although the output is again high-dimensional, certain aspects of the simplification are preserved. The generated conformations only mirror the most dominant conformational differences that determine the positions of conformational states in the low-dimensional map. This allows depicting such differences and-in consequence-visualizing molecular motions and gives a unique perspective on high-dimensional conformational data. In our previous work, protein conformations described in backbone dihedral angle space were used as the input for EncoderMap, and conformations were also generated in this space. For large proteins, however, the generation of conformations is inaccurate with this approach due to the local character of backbone dihedral angles. Here, we present an improved variant of EncoderMap which is able to generate large protein conformations that are accurate in short-range and long-range orders. This is achieved by differentiable reconstruction of Cartesian coordinates from the generated dihedrals, which allows adding a contribution to the cost function that monitors the accuracy of all pairwise distances between the Cα-atoms of the generated conformations. The improved capabilities to generate conformations of large, even multidomain, proteins are demonstrated for two examples: diubiquitin and a part of the Ssa1 Hsp70 yeast chaperone. We show that the improved variant of EncoderMap can nicely visualize motions of protein domains relative to each other but is also able to highlight important conformational changes within the individual domains.

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.

Single molecule translocation in smectics illustrates the challenge for time-mapping in simulations on multiple scales.

Understanding the connections between the characteristic dynamical time scales associated with a coarse-grained (CG) and a detailed representation is central to the applicability of the coarse-graining methods to understand molecular processes. The process of coarse graining leads to an accelerated dynamics, owing to the smoothening of the underlying free-energy landscapes. Often a single time-mapping factor is used to relate the time scales associated with the two representations. We critically examine this idea using a model system ideally suited for this purpose. Single molecular transport properties are studied via molecular dynamics simulations of the CG and atomistic representations of a liquid crystalline, azobenzene containing mesogen, simulated in the smectic and the isotropic phases. The out-of-plane dynamics in the smectic phase occurs via molecular hops from one smectic layer to the next. Hopping can occur via two mechanisms, with and without significant reorientation. The out-of-plane transport can be understood as a superposition of two (one associated with each mode of transport) independent continuous time random walks for which a single time-mapping factor would be rather inadequate. A comparison of the free-energy surfaces, relevant to the out-of-plane transport, qualitatively supports the above observations. Thus, this work underlines the need for building CG models that exhibit both structural and dynamical consistency to the underlying atomistic model.

Tipping the Scale from Disorder to Alpha-helix: Folding of Amphiphilic Peptides in the Presence of Macroscopic and Molecular Interfaces.

Secondary amphiphilicity is inherent to the secondary structural elements of proteins. By forming energetically favorable contacts with each other these amphiphilic building blocks give rise to the formation of a tertiary structure. Small proteins and peptides, on the other hand, are usually too short to form multiple structural elements and cannot stabilize them internally. Therefore, these molecules are often found to be structurally ambiguous up to the point of a large degree of intrinsic disorder in solution. Consequently, their conformational preference is particularly susceptible to environmental conditions such as pH, salts, or presence of interfaces. In this study we use molecular dynamics simulations to analyze the conformational behavior of two synthetic peptides, LKKLLKLLKKLLKL (LK) and EAALAEALAEALAE (EALA), with built-in secondary amphiphilicity upon forming an alpha-helix. We use these model peptides to systematically study their aggregation and the influence of macroscopic and molecular interfaces on their conformational preferences. We show that the peptides are neither random coils in bulk water nor fully formed alpha helices, but adopt multiple conformations and secondary structure elements with short lifetimes. These provide a basis for conformation-selection and population-shift upon environmental changes. Differences in these peptides' response to macroscopic and molecular interfaces (presented by an aggregation partner) can be linked to their inherent alpha-helical tendencies in bulk water. We find that the peptides' aggregation behavior is also strongly affected by presence or absence of an interface, and rather subtly depends on their surface charge and hydrophobicity.

Derivation of coarse-grained simulation models of chlorophyll molecules in lipid bilayers for applications in light harvesting systems.

The correct interplay of interactions between protein, pigment and lipid molecules is highly relevant for our understanding of the association behavior of the light harvesting complex (LHCII) of green plants. To cover the relevant time and length scales in this multicomponent system, a multi-scale simulation ansatz is employed that subsequently uses a classical all atomistic (AA) model to derive a suitable coarse grained (CG) model which can be backmapped into the AA resolution, aiming for a seamless conversion between two scales. Such an approach requires a faithful description of not only the protein and lipid components, but also the interaction functions for the indispensable pigment molecules, chlorophyll b and chlorophyll a (referred to as chl b/chl a). In this paper we develop a CG model for chl b and chl a in a dipalmitoylphosphatidyl choline (DPPC) bilayer system. The structural properties and the distribution behavior of chl within the lipid bilayer in the CG simulations are consistent with those of AA reference simulations. The non-bonded potentials are parameterized such that they fit to the thermodynamics based MARTINI force-field for the lipid bilayer and the protein. The CG simulation shows chl aggregation in the lipid bilayer which is supported by fluorescence quenching experiments. It is shown that the derived chl model is well suited for CG simulations of stable, structurally consistent, trimeric LHCII and can in the future be used to study their large scale aggregation behavior.

Relative resolution: A hybrid formalism for fluid mixtures.

We show here that molecular resolution is inherently hybrid in terms of relative separation. While nearest neighbors are characterized by a fine-grained (geometrically detailed) model, other neighbors are characterized by a coarse-grained (isotropically simplified) model. We notably present an analytical expression for relating the two models via energy conservation. This hybrid framework is correspondingly capable of retrieving the structural and thermal behavior of various multi-component and multi-phase fluids across state space.

Molecular dynamics simulations of peptides at the air-water interface: influencing factors on peptide-templated mineralization.

Biomineralization is the intricate, biomedically highly relevant process by which living organisms deposit minerals on biological matrices to stiffen tissues and build skeletal structures and shells. Rapaport and coworkers ( J. Am. Chem. Soc. 2000 , 122 , 12523 ; Adv. Funct. Mater. 2008 , 18 , 2889 ; Acta Biomater. 2012 , 8 , 2466 ) have designed a class of self-assembling amphiphilic peptides that are capable of forming hydrogels and attracting ions from the environment, generating structures akin to the extracellular matrix and promoting bone regeneration. The air-water interface serves both in experiment and in simulations as a model hydrophobic surface to mimic the cell's organic-aqueous interface and to investigate the organization of the peptide matrix into ordered ß-pleated monolayers and the subsequent onset of biomineral formation. To obtain insight into the underlying molecular mechanism, we have used molecular dynamics simulations to study the effect of peptide sequence on aggregate stability and ion-peptide interactions. We find-in excellent agreement with experimental observations-that the nature of the peptide termini (proline vs phenylalanine) affect the aggregate order, while the nature of the acidic side chains (aspartic vs glutamic acid) affect the aggregate's stability in the presence of ions. These simulations provide valuable microscopic insight into the way ions and peptide templates mutually affect each other during the early stages of biomineralization preceding nucleation.

Sticky water surfaces: helix-coil transitions suppressed in a cell-penetrating peptide at the air-water interface.

GALA is a 30 amino acid synthetic peptide consisting of a Glu-Ala-Leu-Ala repeat and is known to undergo a reversible structural transition from a disordered to an α-helical structure when changing the pH from basic to acidic values. In its helical state GALA can insert into and disintegrate lipid membranes. This effect has generated much interest in GALA as a candidate for pH triggered, targeted drug delivery. GALA also serves as a well-defined model system to understand cell penetration mechanisms and protein folding triggered by external stimuli. Structural transitions of GALA in solution have been studied extensively. However, cell penetration is an interfacial effect and potential biomedical applications of GALA would involve a variety of surfaces, e.g., nanoparticles, lipid membranes, tubing, and liquid-gas interfaces. Despite the apparent importance of interfaces in the functioning of GALA, the effect of surfaces on the reversible folding of GALA has not yet been studied. Here, we use sum frequency generation vibrational spectroscopy (SFG) to probe the structural response of GALA at the air-water interface and IR spectroscopy to follow GALA folding in bulk solution. We combine the SFG data with molecular dynamics simulations to obtain a molecular-level picture of the interaction of GALA with the air-water interface. Surprisingly, while the fully reversible structural transition was observed in solution, at the water-air interface, a large fraction of the GALA population remained helical at high pH. This "stickiness" of the air-water interface can be explained by the stabilizing interactions of hydrophobic leucine and alanine side chains with the water surface.

Cryo-EM of soft-landed ß-galactosidase: Gas-phase and native structures are remarkably similar.

Native mass spectrometry (MS) has become widely accepted in structural biology, providing information on stoichiometry, interactions, homogeneity, and shape of protein complexes. Yet, the fundamental assumption that proteins inside the mass spectrometer retain a structure faithful to native proteins in solution remains a matter of intense debate. Here, we reveal the gas-phase structure of ß-galactosidase using single-particle cryo-electron microscopy (cryo-EM) down to 2.6-Å resolution, enabled by soft landing of mass-selected protein complexes onto cold transmission electron microscopy (TEM) grids followed by in situ ice coating. We find that large parts of the secondary and tertiary structure are retained from the solution. Dehydration-driven subunit reorientation leads to consistent compaction in the gas phase. By providing a direct link between high-resolution imaging and the capability to handle and select protein complexes that behave problematically in conventional sample preparation, the approach has the potential to expand the scope of both native mass spectrometry and cryo-EM.

A transferable coarse-grained model for diphenylalanine: how to represent an environment driven conformational transition.

One of the major challenges in the development of coarse grained (CG) simulation models that aim at biomolecular structure formation processes is the correct representation of an environment-driven conformational change, for example, a folding/unfolding event upon interaction with an interface or upon aggregation. In the present study, we investigate this transferability challenge for a CG model using the example of diphenylalanine. This dipeptide displays a transition from a trans-like to a cis-like conformation upon aggregation as well as upon transfer from bulk water to the cyclohexane/water interface. Here, we show that one can construct a single CG model that can reproduce both the bulk and interface conformational behavior and the segregation between hydrophobic/hydrophilic medium. While the general strategy to obtain nonbonded interactions in the present CG model is to reproduce solvation free energies of small molecules representing the CG beads in the respective solvents, the success of the model strongly depends on nontrivial decisions one has to make to capture the delicate balance between the bonded and nonbonded interactions. In particular, we found that the peptide's conformational behavior is qualitatively affected by the cyclohexane/water interaction potential, an interaction that does not directly involve the peptide at all but merely influences the properties of the hydrophobic/hydrophilic interface. Furthermore, we show that a small modification to improve the structural/conformational properties of the CG model could dramatically alter the thermodynamic properties.