Transposition of HOPPLA in siRNA-deficient plants suggests a limited effect of the environment on retrotransposon mobility in Brachypodium distachyon.

Long terminal repeat retrotransposons (LTR-RTs) are powerful mutagens regarded as a major source of genetic novelty and important drivers of evolution. Yet, the uncontrolled and potentially selfish proliferation of LTR-RTs can lead to deleterious mutations and genome instability, with large fitness costs for their host. While population genomics data suggest that an ongoing LTR-RT mobility is common in many species, the understanding of their dual role in evolution is limited. Here, we harness the genetic diversity of 320 sequenced natural accessions of the Mediterranean grass Brachypodium distachyon to characterize how genetic and environmental factors influence plant LTR-RT dynamics in the wild. When combining a coverage-based approach to estimate global LTR-RT copy number variations with mobilome-sequencing of nine accessions exposed to eight different stresses, we find little evidence for a major role of environmental factors in LTR-RT accumulations in B. distachyon natural accessions. Instead, we show that loss of RNA polymerase IV (Pol IV), which mediates RNA-directed DNA methylation in plants, results in high transcriptional and transpositional activities of RLC_BdisC024 (HOPPLA) LTR-RT family elements, and that these effects are not stress-specific. This work supports findings indicating an ongoing mobility in B. distachyon and reveals that host RNA-directed DNA methylation rather than environmental factors controls their mobility in this wild grass model.

Dynamical Reweighting for Biased Rare Event Simulations.

Dynamical reweighting techniques aim to recover the correct molecular dynamics from a simulation at a modified potential energy surface. They are important for unbiasing enhanced sampling simulations of molecular rare events. Here, we review the theoretical frameworks of dynamical reweighting for modified potentials. Based on an overview of kinetic models with increasing level of detail, we discuss techniques to reweight two-state dynamics, multistate dynamics, and path integrals. We explore the natural link to transition path sampling and how the effect of nonequilibrium forces can be reweighted. We end by providing an outlook on how dynamical reweighting integrates with techniques for optimizing collective variables and with modern potential energy surfaces.

Thermal isomerization rates in retinal analogues using Ab-Initio molecular dynamics.

For a detailed understanding of chemical processes in nature and industry, we need accurate models of chemical reactions in complex environments. While Eyring transition state theory is commonly used for modeling chemical reactions, it is most accurate for small molecules in the gas phase. A wide range of alternative rate theories exist that can better capture reactions involving complex molecules and environmental effects. However, they require that the chemical reaction is sampled by molecular dynamics simulations. This is a formidable challenge since the accessible simulation timescales are many orders of magnitude smaller than typical timescales of chemical reactions. To overcome these limitations, rare event methods involving enhanced molecular dynamics sampling are employed. In this work, thermal isomerization of retinal is studied using tight-binding density functional theory. Results from transition state theory are compared to those obtained from enhanced sampling. Rates obtained from dynamical reweighting using infrequent metadynamics simulations were in close agreement with those from transition state theory. Meanwhile, rates obtained from application of Kramers' rate equation to a sampled free energy profile along a torsional dihedral reaction coordinate were found to be up to three orders of magnitude higher. This discrepancy raises concerns about applying rate methods to one-dimensional reaction coordinates in chemical reactions.

How Chromophore Labels Shape the Structure and Dynamics of a Peptide Hydrogel.

Biocompatible and functionalizable hydrogels have a wide range of (potential) medicinal applications. The hydrogelation process, particularly for systems with very low polymer weight percentages (<1 wt %), remains poorly understood, making it challenging to predict the self-assembly of a given molecular building block into a hydrogel. This severely hinders the rational design of self-assembled hydrogels. In this study, we demonstrate the impact of an N-terminal group on the self-assembly and rheology of the peptide hydrogel hFF03 (hydrogelating, fibril forming peptide 03) using molecular dynamics simulations, oscillatory shear rheology, and circular dichroism spectroscopy. We find that the chromophore and even its specific regioisomers have a significant influence on the microscopic structure and dynamics of the self-assembled fibril, and on the macroscopic mechanical properties. This is because the chromophore influences the possible salt bridges, which form and stabilize the fibril formation. Furthermore, we find that the solvation shell fibrils by itself cannot explain the viscoelasticity of hFF03 hydrogels. Our atomistic model of the hFF03 fibril formation enables a more rational design of these hydrogels. In particular, altering the N-terminal chromophore emerges as a design strategy to tune the mechanic properties of these self-assembled peptide hydrogels.

Prebound State Discovered in the Unbinding Pathway of Fluorinated Variants of the Trypsin-BPTI Complex Using Random Acceleration Molecular Dynamics Simulations.

The serine protease trypsin forms a tightly bound inhibitor complex with the bovine pancreatic trypsin inhibitor (BPTI). The complex is stabilized by the P1 residue Lys15, which interacts with negatively charged amino acids at the bottom of the S1 pocket. Truncating the P1 residue of wildtype BPTI to α-aminobutyric acid (Abu) leaves a complex with moderate inhibitor strength, which is held in place by additional hydrogen bonds at the protein-protein interface. Fluorination of the Abu residue partially restores the inhibitor strength. The mechanism with which fluorination can restore the inhibitor strength is unknown, and accurate computational investigation requires knowledge of the binding and unbinding pathways. The preferred unbinding pathway is likely to be complex, as encounter states have been described before, and unrestrained umbrella sampling simulations of these complexes suggest additional energetic minima. Here, we use random acceleration molecular dynamics to find a new metastable state in the unbinding pathway of Abu-BPTI variants and wildtype BPTI from trypsin, which we call the prebound state. The prebound state and the fully bound state differ by a substantial shift in the position, a slight shift in the orientation of the BPTI variants, and changes in the interaction pattern. Particularly important is the breaking of three hydrogen bonds around Arg17. Fluorination of the P1 residue lowers the energy barrier of the transition between the fully bound state and prebound state and also lowers the energy minimum of the prebound state. While the effect of fluorination is in general difficult to quantify, here, it is in part caused by favorable stabilization of a hydrogen bond between Gln194 and Cys14. The interaction pattern of the prebound state offers insights into the inhibitory mechanism of BPTI and might add valuable information for the design of serine protease inhibitors.

Impact of glycan nature on structure and viscoelastic properties of glycopeptide hydrogels.

Mucus is a complex biological hydrogel that acts as a barrier for almost everything entering or exiting the body. It is therefore of emerging interest for biomedical and pharmaceutical applications. Besides water, the most abundant components are the large and densely glycosylated mucins, glycoproteins of up to 20 MDa and carbohydrate content of up to 80 wt%. Here, we designed and explored a library of glycosylated peptides to deconstruct the complexity of mucus. Using the well-characterized hFF03 coiled-coil system as a hydrogel-forming peptide scaffold, we systematically probed the contribution of single glycans to the secondary structure as well as the formation and viscoelastic properties of the resulting hydrogels. We show that glycan-decoration does not affect α-helix and coiled-coil formation while it alters gel stiffness. By using oscillatory macrorheology, dynamic light scattering microrheology, and fluorescence lifetime-based nanorheology, we characterized the glycopeptide materials over several length scales. Molecular simulations revealed that the glycosylated linker may extend into the solvent, but more frequently interacts with the peptide, thereby likely modifying the stability of the self-assembled fibers. This systematic study highlights the interplay between glycan structure and hydrogel properties and may guide the development of synthetic mucus mimetics.

Implementation of Girsanov Reweighting in OpenMM and Deeptime.

Classical molecular dynamics (MD) simulations provide invaluable insights into complex molecular systems but face limitations in capturing phenomena occurring on time scales beyond their reach. To bridge this gap, various enhanced sampling techniques have been developed, which are complemented by reweighting techniques to recover the unbiased dynamics. Girsanov reweighting is a reweighting technique that reweights simulation paths, generated by a stochastic MD integrator, without evoking an effective model of the dynamics. Instead, it calculates the relative path probability density at the time resolution of the MD integrator. Efficient implementation of Girsanov reweighting requires that the reweighting factors are calculated on-the-fly during the simulations and thus needs to be implemented within the MD integrator. Here, we present a comprehensive guide for implementing Girsanov reweighting into MD simulations. We demonstrate the implementation in the MD simulation package OpenMM by extending the library openmmtools. Additionally, we implemented a reweighted Markov state model estimator within the time series analysis package Deeptime.

Fluorinated Protein-Ligand Complexes: A Computational Perspective.

Fluorine is an element renowned for its unique properties. Its powerful capability to modulate molecular properties makes it an attractive substituent for protein binding ligands; however, the rational design of fluorination can be challenging with effects on interactions and binding energies being difficult to predict. In this Perspective, we highlight how computational methods help us to understand the role of fluorine in protein-ligand binding with a focus on molecular simulation. We underline the importance of an accurate force field, present fluoride channels as a showcase for biomolecular interactions with fluorine, and discuss fluorine specific interactions like the ability to form hydrogen bonds and interactions with aryl groups. We put special emphasis on the disruption of water networks and entropic effects.

Nutritional support in hospitalised patients with diabetes and risk for malnutrition: a secondary analysis of an investigator-initiated, Swiss, randomised controlled multicentre trial.

OBJECTIVES: The main objective of this study was to investigate the effects of nutritional support on mortality in hospitalised patients with diabetes and nutritional risk participating in the Effect of early nutritional support on Frailty, Functional Outcomes, and Recovery of malnourished medical inpatients Trial (EFFORT) trial. DESIGN: Secondary analysis of a Swiss-wide multicentre, randomised controlled trial. PARTICIPANTS: Patients with diabetes and risk for malnutrition. INTERVENTIONS: Individualised nutritional support versus usual care. PRIMARY OUTCOME MEASURE: 30-day all-cause mortality. RESULTS: Of the 2028 patients included in the original trial, 445 patients were diagnosed with diabetes and included in this analysis. In terms of efficacy of nutritional therapy, there was a 25% lower risk for mortality in patients with diabetes receiving nutritional support compared with controls (7% vs 10%, adjusted HR 0.75 (95% CI 0.39 to 1.43)), a finding that was not statistically significant but similar to the overall trial effects with no evidence of interaction (p=0.92). Regarding safety of nutritional therapy, there was no increase in diabetes-specific complications associated with nutritional support, particularly there was no increase in risk for hyperglycaemia (adjusted OR 0.97, 95% CI 0.56 to 1.67 p=0.90). CONCLUSION: Patients with diabetes and malnutrition in the hospital setting have a particularly high risk for adverse outcomes and mortality. Individualised nutritional support reduced mortality in this secondary analysis of a randomized trial, but this effect was not significant calling for further large-scale trials in this vhighly ulnerable patient population. TRIAL REGISTRATION NUMBER: NCT02517476.

Accelerated generation of gene-engineered monoclonal CHO cell lines using FluidFM nanoinjection and CRISPR/Cas9.

Chinese hamster ovary (CHO) cells are the commonly used mammalian host system to manufacture recombinant proteins including monoclonal antibodies. However unfavorable non-human glycoprofile displayed on CHO-produced monoclonal antibodies have negative impacts on product quality, pharmacokinetics, and therapeutic efficiency. Glycoengineering such as genetic elimination of genes involved in glycosylation pathway in CHO cells is a viable solution but constrained due to longer timeline and laborious workflow. Here, in this proof-of-concept (PoC) study, we present a novel approach coined CellEDIT to engineer CHO cells by intranuclear delivery of the CRISPR components to single cells using the FluidFM technology. Co-injection of CRISPR system targeting BAX, DHFR, and FUT8 directly into the nucleus of single cells, enabled us to generate triple knockout CHO-K1 cell lines within a short time frame. The proposed technique assures the origin of monoclonality without the requirement of limiting dilution, cell sorting or positive selection. Furthermore, the approach is compatible to develop both single and multiple knockout clones (FUT8, BAX, and DHFR) in CHO cells. Further analyses on single and multiple knockout clones confirmed the targeted genetic disruption and altered protein expression. The knockout CHO-K1 clones showed the persistence of gene editing during the subsequent passages, compatible with serum free chemically defined media and showed equivalent transgene expression like parental clone.