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.

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.

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.

Visualizing the Residue Interaction Landscape of Proteins by Temporal Network Embedding

Characterizing the structural dynamics of proteins with heterogeneous conformational landscapes is crucial to understanding complex biomolecular processes. To this end, dimensionality reduction algorithms are used to produce low-dimensional embeddings of the high-dimensional conformational phase space. However, identifying a compact and informative set of input features for the embedding remains an ongoing challenge. Here, we propose to harness the power of Residue Interaction Networks (RINs) and their centrality measures, established tools to provide a graph theoretical view on molecular structure. Specifically, we combine the closeness centrality, which captures global features of the protein conformation at residue-wise resolution, with EncoderMap, a hybrid neural-network autoencoder/multidimensional-scaling like dimensionality reduction algorithm. We find that the resulting low-dimensional embedding is a meaningful visualization of the residue interaction landscape that resolves structural details of the protein behavior while retaining global interpretability. This feature-based graph embedding of temporal protein graphs makes it possible to apply the general descriptive power of RIN formalisms to the analysis of protein simulations of complex processes such as protein folding and multidomain interactions requiring no protein-specific input. We demonstrate this on simulations of the fast folding protein Trp-Cage and the multidomain signaling protein FAT10. Due to its generality and modularity, the presented approach can easily be transferred to other protein systems.

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.

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.

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 pK_a 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.

Lipid-mediated activation of plasma membrane-localized deubiquitylating enzymes modulate endosomal trafficking

The abundance of plasma membrane-resident receptors and transporters has to be tightly regulated by ubiquitin-mediated endosomal degradation for the proper coordination of environmental stimuli and intracellular signaling. Arabidopsis OVARIAN TUMOR PROTEASE (OTU) 11 and OTU12 are plasma membrane-localized deubiquitylating enzymes (DUBs) that bind to phospholipids through a polybasic motif in the OTU domain. Here we show that the DUB activity of OTU11 and OTU12 towards K63-linked ubiquitin is stimulated by binding to lipid membranes containing anionic lipids. In addition, we show that the DUB activity of OTU11 against K6- and K11-linkages is also stimulated by anionic lipids, and that OTU11 and OTU12 can modulate the endosomal degradation of a model cargo and the auxin efflux transporter PIN2-GFP in vivo. Our results suggest that the catalytic activity of OTU11 and OTU12 is tightly connected to their ability to bind membranes and that OTU11 and OTU12 are involved in the fine-tuning of plasma membrane proteins in Arabidopsis.

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 (D₂O) 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 (H₂O) 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 H₂O or D₂O. 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 D₂O. 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 H₂O and D₂O have significant, distinct effects on mineralization mechanisms, and that care has to be taken when experimental data from D₂O studies are used, e.g., for the development of H₂O-based computer models.

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.

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.

Three Reasons Why Aspartic Acid and Glutamic Acid Sequences Have a Surprisingly Different Influence on Mineralization

Understanding the role of polymers rich in aspartic acid (Asp) and glutamic acid (Glu) is the key to gaining precise control over mineralization processes. Despite their chemical similarity, experiments revealed a surprisingly different influence of Asp and Glu sequences. We conducted molecular dynamics simulations of Asp and Glu peptides in the presence of calcium and chloride ions to elucidate the underlying phenomena. In line with experimental differences, in our simulations, we indeed find strong differences in the way the peptides interact with ions in solution. The investigated Asp pentapeptide tends to pull a lot of ions into its vicinity, and many structures with clusters of calcium and chloride ions on the surface of the peptide can be observed. Under the same conditions, comparatively fewer ions can be found in proximity of the investigated Glu pentapeptide, and the structures are characterized by single calcium ions bound to multiple carboxylate groups. Based on our simulation data, we identified three reasons contributing to these differences, leading to a new level of understanding additive-ion interactions.

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.

Coarse grained simulation of the aggregation and structure control of polyethylene nanocrystals

Polyethylene (PE) telechelics with carboxylate functional groups at both ends have been shown to assemble into hexagonal nanocrystal platelets with a height defined by their chain length in basic CsOH-solution. In this coarse grained (CG) simulation study we show how properties of the functional groups alter the aggregation and crystallization behavior of those telechelics. Systematic variation of the parameters of the CG model showed that important factors which control nanoparticle stability and structure are the PE chain length and the hydrophilicity and the steric demand of the head groups. To characterize the aggregation process we analyzed the number and size of the obtained aggregates as well as intramolecular order and intermolecular alignment of the polymer chains. By comparison of CG and atomistic simulation data, it could be shown that atomistic simulations representing different chemical systems can be emulated with specific, different CG parameter sets. Thus, the results from the (generic) CG simulation models can be used to explain the effect of different head groups and different counterions on the aggregation of PE telechelics and the order of the obtained nanocrystals.

Deciphering molecular details of the RAC-ribosome interaction by EPR spectroscopy

The eukaryotic ribosome-associated complex (RAC) plays a significant role in de novo protein folding. Its unique interaction with the ribosome, comprising contacts to both ribosomal subunits, suggests a RAC-mediated coordination between translation elongation and co-translational protein folding. Here, we apply electron paramagnetic resonance (EPR) spectroscopy combined with site-directed spin labeling (SDSL) to gain deeper insights into a RAC-ribosome contact affecting translational accuracy. We identified a local contact point of RAC to the ribosome. The data provide the first experimental evidence for the existence of a four-helix bundle as well as a long α-helix in full-length RAC, in solution as well as on the ribosome. Additionally, we complemented the structural picture of the region mediating this functionally important contact on the 40S ribosomal subunit. In sum, this study constitutes the first application of SDSL-EPR spectroscopy to elucidate the molecular details of the interaction between the 3.3 MDa translation machinery and a chaperone complex.

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.

Coarse-Grained Simulation of the Adsorption of Water on Au(111) Surfaces Using a Modified Stillinger-Weber Potential

For reproducing the behavior of water molecules adsorbed on gold surfaces in terms of density of both bulk and interfacial water and in terms of structuring of water on top of gold atoms, the implementation of a multibody potential is necessary, thus the Stillinger-Weber potential was tested. The goal is using a single nonbonded potential for coarse-grained models, without the usage of explicit charges. In order to modify the angular part of the Stillinger-Weber potential from a single cosine to a piecewise function accounting for multiple equilibrium angles, employed for Au-Au-Au and Au-Au-water triplets, it is necessary to create a version of the simulation package LAMMPS that supports the assignment of multiple favored angles. This novel approach is able to reproduce the data obtained using quantum mechanical calculations and density profiles of both bulk and adsorbed water molecules obtained using classical polarizable force fields.

Machine Learning Force Fields and Coarse-Grained Variables in Molecular Dynamics: Application to Materials and Biological Systems

Machine learning encompasses tools and algorithms that are now becoming popular in almost all scientific and technological fields. This is true for molecular dynamics as well, where machine learning offers promises of extracting valuable information from the enormous amounts of data generated by simulation of complex systems. We provide here a review of our current understanding of goals, benefits, and limitations of machine learning techniques for computational studies on atomistic systems, focusing on the construction of empirical force fields from ab initio databases and the determination of reaction coordinates for free energy computation and enhanced sampling.

Getting to the bottom of bycatch: a GIS-based toolbox to assess the risk of marine mammal bycatch

Marine mammal bycatch poses a particular challenge in developing countries, where data to document bycatch and its effects are often lacking. Using the Bycatch Risk Assessment (ByRA) toolkit, based on InVEST open-source models, we chose 4 field sites in Southeast Asia with varying amounts of data on marine mammals and fishing occurrence: Trat province in the eastern Gulf of Thailand, the Sibu-Tinggi Islands and Kuching Bay, Malaysia, and Kien Giang Biosphere Reserve in southwestern Vietnam. These field sites have similar species of coastal marine mammals, small-scale and commercial fisheries, and support for research from universities and/or management. In Thailand and Kuching, results showed changing patterns of fishing and Irrawaddy dolphin Orcaella brevirostris habitat use across seasons, showing how bycatch risk could change throughout the year. Risk maps for dugongs Dugong dugon in peninsular Malaysia highlighted patterns of bycatch risk concentrated around a mainland fishing pier, and revealed high risk in a northern subregion. In Vietnam, first maps of bycatch risk for the Irrawaddy dolphin showed the highest risk driven by intensive use of gillnets and trawling gear. ByRA pinpointed areas of spatial and seasonal bycatch exposure, and estimated the consequence of bycatch on local species, providing managers with critical information on where to focus bycatch mitigation and meet new global standards for US Marine Mammal Protection Act and other international regulation (e.g. Official Journal of the European Union 2019; Regulation 2019/1241) compliance. The toolbox, a transferable open-source tool, can be used to guide fisheries management, marine mammal conservation, spatial planning, and further research.

Machine Learning Driven Analysis of Large Scale Simulations Reveals Conformational Characteristics of Ubiquitin Chains

Understanding the conformational characteristics of protein complexes in solution is crucial for a deeper insight in their biological function. Molecular dynamics simulations performed on high performance computing plants and with modern simulation techniques can be used to obtain large data sets that contain conformational and thermodynamic information about biomolecular systems. While this can in principle give a detailed picture of protein-protein interactions in solution and therefore complement experimental data, it also raises the challenge of processing exceedingly large high-dimensional data sets with several million samples. Here we present a novel method for the characterization of protein-protein interactions, which combines a neural network based dimensionality reduction technique to obtain a two-dimensional representation of the conformational space with a density based clustering algorithm for state detection and a metric which assesses the (dis)similarity between different conformational spaces. This method is highly scalable and therefore makes the analysis of massive data sets computationally tractable. We demonstrate the power of this approach to large scale data analysis by characterizing highly dynamic conformational phase spaces of differently linked ubiquitin (Ub) oligomers from coarse-grained simulations. We are able to extract a protein-protein interaction model for two unlinked Ub proteins which is then used to determine how the Ub-Ub interaction pattern is altered in Ub oligomers by the introduction of a covalent linkage. We find that the Ub chain conformational ensemble depends highly on the linkage type and for some cases also on the Ub chain length. By this, we obtain insight into the conformational characteristics of different Ub chains and how this may contribute to linkage type and chain length specific recognition.

Lipoprotein apheresis is an optimal therapeutic option to reduce increased Lp(a) levels

BACKGROUND

Lipoprotein(a) (Lp(a)) is a genetic risk factor for cardiovascular disease (CVD) and is associated with the induction and sustaining of atherosclerotic cardiovascular diseases (ASCVD). Since 2008 Lp(a) along with progressive CVD has been approved as an indication for regular lipoprotein apheresis (LA) in Germany. The German Lipoprotein Apheresis Registry (GLAR) has been initiated to provide statistical evidence for the assessment of extracorporeal procedures to treat dyslipidemia for both LDL-cholesterol (LDL-C) and Lp(a). The GLAR now allows prospective investigations over a 5-year period about annual incidence rates of cardiovascular events. Here Lp(a) patients (LDL-C < 100 mg/dl; Lp(a) > 60 mg/dl or >120 nmol/l) showed the same reduction of major coronary (83%) and non-coronary events (63%) as had been formerly shown in the Pro(a)LiFe study. However, Lp(a) is not only an apolipoprotein(a) (apo(a)) and LDL-C containing particle, which is covalently bound to a LDL-C core by a disulphide bridge. The composition of this particle, inter alia containing oxidized phospholipids, gives pro-atherosclerotic, pro-inflammatory, and pro-thrombotic properties, inducing atherosclerotic processes mainly in the arterial wall. However, recent investigations have shown that a reduction of inflammatory settings without LDL-C or Lp(a) reduction may reduce ASCVD events. Lipoprotein apheresis (LA) could not only reduce LDL-C and Lp(a) in parallel, but also different inflammatory and coagulation parameters. In summary lipoprotein apheresis is not only anti-atherosclerotic, but also anti-inflammatory and anti-thrombotic and therefore an ideal treatment option with respect to the shown reduction of major adverse coronary events (MACE) and major adverse non-coronary events (MANCE) by reducing Lp(a) levels.

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.

Simulating and analysing configurational landscapes of protein-protein contact formation

Interacting proteins can form aggregates and protein-protein interfaces with multiple patterns and different stabilities. Using molecular simulation one would like to understand the formation of these aggregates and which of the observed states are relevant for protein function and recognition. To characterize the complex configurational ensemble of protein aggregates, one needs a quantitative measure for the similarity of structures. We present well-suited descriptors that capture the essential features of non-covalent protein contact formation and domain motion. This set of collective variables is used with a nonlinear multi-dimensional scaling-based dimensionality reduction technique to obtain a low-dimensional representation of the configurational landscape of two ubiquitin proteins from coarse-grained simulations. We show that this two-dimensional representation is a powerful basis to identify meaningful states in the ensemble of aggregated structures and to calculate distributions and free energy landscapes for different sets of simulations. By using a measure to quantitatively compare free energy landscapes we can show how the introduction of a covalent bond between two ubiquitin proteins at different positions alters the configurational states of these dimers.

Heart rate, respiratory rate, apnoeas and peripheral arterial oxygen saturation in healthy term neonates during quiet sleep

AIM

This study compiled percentiles for cardiorespiratory parameters in healthy term neonates during quiet sleep.

METHODS

We enrolled 215 healthy term neonates born at Hannover Medical School, Germany, between October 2011 and March 2013. They were prospectively observed on the maternity ward at a median age of two days using six-hour recordings of pulse oximeter plethysmography, oxygen saturation, thoracic breathing movements and electrocardiogram during sleep in a supine position. We examined their heart rate, respiratory rate and oxygen saturation during quiet sleep, plus bradycardias, apnoeas lasting at least four-seconds and desaturations below 85%.

RESULTS

The 3rd, 50th and 97th percentiles were calculated as follows: heart rate 87, 112 and 133 beats per minute, respiratory rate 32, 44 and 57 per minute and oxygen saturation 94, 98 and 100%. Desaturations, apnoeas and bradycardias below 80 beats per minute were common and recorded in 54%, 98% and 30% of participants. In contrast, only 7% experienced bradycardias of less than two-thirds of the baseline heart rate and 5% experienced apnoeas exceeding 15 seconds.

CONCLUSION

Our results will facilitate the evidence-based valuation of cardiorespiratory parameters in term neonates and help validate the significance of cardiorespiratory events in preterm infants at discharge.

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.

Post-translational modification of proteins by covalent attachment of ubiquitin is a key cellular process, regulating for example the fate and recycling of proteins. We present a new method to combine multiscale simulation with advanced analysis methods to characterize the states of ubiquitin-ubiquitin conjugates. We found that the linkage position affects the conformational space of ubiquitin dimers, determining the number and stability of relevant states, the character of subunit contacts and the nature of the surface exposed to possible binding partners.

Author Correction: The structure of the ubiquitin-like modifier FAT10 reveals an alternative targeting mechanism for proteasomal degradation

The original version of the Supplementary Information associated with this Article inadvertently omitted Supplementary Table 3. The HTML version of the Article has been updated to include a corrected version of the Supplementary Information.

The structure of the ubiquitin-like modifier FAT10 reveals an alternative targeting mechanism for proteasomal degradation

FAT10 is a ubiquitin-like modifier that directly targets proteins for proteasomal degradation. Here, we report the high-resolution structures of the two individual ubiquitin-like domains (UBD) of FAT10 that are joined by a flexible linker. While the UBDs of FAT10 show the typical ubiquitin-fold, their surfaces are entirely different from each other and from ubiquitin explaining their unique binding specificities. Deletion of the linker abrogates FAT10-conjugation while its mutation blocks auto-FAT10ylation of the FAT10-conjugating enzyme USE1 but not bulk conjugate formation. FAT10- but not ubiquitin-mediated degradation is independent of the segregase VCP/p97 in the presence but not the absence of FAT10’s unstructured N-terminal heptapeptide. Stabilization of the FAT10 UBDs strongly decelerates degradation suggesting that the intrinsic instability of FAT10 together with its disordered N-terminus enables the rapid, joint degradation of FAT10 and its substrates without the need for FAT10 de-conjugation and partial substrate unfolding.

The ubiquitin-like modifier FAT10 is composed of two ubiquitin-like domains (UBDs). Here the authors present the FAT10 UBD structures and show that the unstructured FAT10 N-terminal heptapeptide together with the poor stability of FAT10 facilitate the rapid proteasomal targeting of FAT10 along with its substrates.

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.

Soluble Oligomeric Nucleants: Simulations of Chain Length, Binding Strength, and Volume Fraction Effects

Recent theories and simulations suggest that molecular additives can bind to the surfaces of nuclei, lower the surface energy, and accelerate nucleation. Experiments have shown that oligomeric and polymeric additives can also modify nucleation rates of proteins, ice, and minerals; however, general design principles for oligomeric or polymeric promoters do not yet exist. Here we investigate oligomeric additives for which each segment of the oligomer can bind to surfaces of nuclei. We use semigrand canonical Monte Carlo simulations in a Potts lattice gas model to study the effects of oligomer chain length, volume fraction, and binding strength. We find that increasing each of those parameters lowers the nucleation barrier. At extremely low oligomer concentrations, the nucleation kinetics can be modeled as though each oligomer is a heterogeneous nucleation site in solution.

Current insights into the German Lipoprotein Apheresis Registry (GLAR) - Almost 5 years on

BACKGROUND

In recent years the Federal Joint Committee (G-BA), a paramount decision-making body of the German health care system required a reassessment of the approval of chronic lipoprotein apheresis therapy for regular reimbursement. Since 2005 an interdisciplinary German apheresis working group has been established by members of both German Societies of Nephrology. In 2009 the working group completed the indication for lipoprotein apheresis with respect to current cardiovascular guidelines and current scientific knowledge for the registry. In 2011 the German Lipoprotein Apheresis Registry (GLAR) was launched and data acquired over nearly 5 years can now be reported.

METHODS AND RESULTS

All data were collected and analyzed during the time period 2012-2015. Over this time interval, 68 German apheresis centers collected retrospective and prospective observational data of 1.283 patients undergoing lipoprotein apheresis (LA) treatment of high LDL-cholesterol (LDL-C) levels and/or high lipoprotein(a) (Lp(a)) levels suffering from progressive cardiovascular disease (CVD). A total of 15,167 documented LA treatments were investigated. All patients treated by LA exhibited a median LDL-C reduction rate of 68.6%, and a median Lp(a) reduction rate of 70.4%. Analogue to the Pro(a)LiFe pattern, patient data were analyzed and compared with respect to the incidence rate of coronary events (MACE) 1 and 2 years before the start of LA treatment (y-2 and y-1) and prospectively one year on LA treatment (y+1). During the first year of LA treatment a MACE reduction of 97% was be observed. In the years considered, LA treatment side effects occurred at a low rate (ca. 5%) and mainly comprised puncture problems.

CONCLUSIONS

For the first time data generated by the GLAR shows that LA lowers the incidence rate of cardiovascular events in patients with high LDL-C and/or high Lp(a) levels, progressive CVD and maximally tolerated lipid lowering medication. In addition LA treatments were found to be safe, exhibiting a low rate of side effects.

A green direct preparation of a magnetic ordered mesoporous carbon catalyst containing Fe–Pd alloys: application to Suzuki–Miyaura reactions in propane-1,2-diol

A mesoporous carbon containing Fe–Pd alloys showed excellent activity for Suzuki–Miyaura couplings using tiny amounts of Pd and afforded Pd-free products after reaction.

Systematic comparison of model polymer nanocomposite mechanics

Polymer nanocomposites render a range of outstanding materials from natural products such as silk, sea shells and bones, to synthesized nanoclay or carbon nanotube reinforced polymer systems. In contrast to the fast expanding interest in this type of material, the fundamental mechanisms of their mixing, phase behavior and reinforcement, especially for higher nanoparticle content as relevant for bio-inorganic composites, are still not fully understood. Although polymer nanocomposites exhibit diverse morphologies, qualitatively their mechanical properties are believed to be governed by a few parameters, namely their internal polymer network topology, nanoparticle volume fraction, particle surface properties and so on. Relating material mechanics to such elementary parameters is the purpose of this work. By taking a coarse-grained molecular modeling approach, we study an range of different polymer nanocomposites. We vary polymer nanoparticle connectivity, surface geometry and volume fraction to systematically study rheological/mechanical properties. Our models cover different materials, and reproduce key characteristics of real nanocomposites, such as phase separation, mechanical reinforcement. The results shed light on establishing elementary structure, property and function relationship of polymer nanocomposites.