Learning Protein-Ligand Unbinding Pathways via Single-Parameter Community Detection.

Understanding the dynamics of biomolecular complexes, e.g., of protein-ligand (un)binding, requires the comprehension of paths such systems take between metastable states.

Nonadiabatic dynamics of molecules interacting with metal surfaces: A quantum-classical approach based on Langevin dynamics and the hierarchical equations of motion.

A novel mixed quantum-classical approach to simulating nonadiabatic dynamics of molecules at metal surfaces is presented. The method combines the numerically exact hierarchical equations of motion approach for the quantum electronic degrees of freedom with Langevin dynamics for the classical degrees of freedom, namely, low-frequency vibrational modes within the molecule. The approach extends previous mixed quantum-classical methods based on Langevin equations to models containing strong electron-electron or quantum electronic-vibrational interactions, while maintaining a nonperturbative and non-Markovian treatment of the molecule-metal coupling.

Basal delamination during mouse gastrulation primes pluripotent cells for differentiation.

The blueprint of the mammalian body plan is laid out during gastrulation, when a trilaminar embryo is formed. This process entails a burst of proliferation, the ingression of embryonic epiblast cells at the primitive streak, and their priming toward primitive streak fates. How these different events are coordinated remains unknown.

Investigation of Rare Protein Conformational Transitions via Dissipation-Corrected Targeted Molecular Dynamics.

To sample rare events, dissipation-corrected targeted molecular dynamics (dcTMD) applies a constant velocity constraint along a one-dimensional reaction coordinate , which drives an atomistic system from an initial state into a target state. Employing a cumulant approximation of Jarzynski's identity, the free energy Δ() is calculated from the mean external work and dissipated work of the process. By calculating the friction coefficient Γ() from the dissipated work, in a second step, the equilibrium dynamics of the process can be studied by propagating a Langevin equation.

The Onset of Molecule-Spanning Dynamics in Heat Shock Protein Hsp90.

Protein dynamics have been investigated on a wide range of time scales. Nano- and picosecond dynamics have been assigned to local fluctuations, while slower dynamics have been attributed to larger conformational changes. However, it is largely unknown how fast (local) fluctuations can lead to slow global (allosteric) changes.

Predicting Protein-Ligand Binding and Unbinding Kinetics with Biased MD Simulations and Coarse-Graining of Dynamics: Current State and Challenges.

The prediction of drug-target binding and unbinding kinetics that occur on time scales between milliseconds and several hours is a prime challenge for biased molecular dynamics simulation approaches. This Perspective gives a concise summary of the theory and the current state-of-the-art of such predictions via biased simulations, of insights into the molecular mechanisms defining binding and unbinding kinetics as well as of the extraordinary challenges predictions of ligand kinetics pose in comparison to binding free energy predictions.

Path separation of dissipation-corrected targeted molecular dynamics simulations of protein-ligand unbinding.

Protein-ligand (un)binding simulations are a recent focus of biased molecular dynamics simulations. Such binding and unbinding can occur via different pathways in and out of a binding site. Here, we present a theoretical framework on how to compute kinetics along separate paths and on how to combine the path-specific rates into global binding and unbinding rates for comparison with experimental results.

Anisotropic Friction in a Ligand-Protein Complex.

The effect of an externally applied directional force on molecular friction is so far poorly understood. Here, we study the force-driven dissociation of the ligand-protein complex biotin-streptavidin and identify anisotropic friction as a not yet described type of molecular friction. Using AFM-based stereographic single molecule force spectroscopy and targeted molecular dynamics simulations, we find that the rupture force and friction for biotin-streptavidin vary with the pulling angle.

Energy Transport and Its Function in Heptahelical Transmembrane Proteins.

Photoproteins such as bacteriorhodopsin (bR) and rhodopsin (Rho) need to effectively dissipate photoinduced excess energy to prevent themselves from damage. Another well-studied seven transmembrane (TM) helices protein is the β adrenergic receptor (βAR), a G protein-coupled receptor for which energy dissipation paths have been linked with allosteric communication. To study the vibrational energy transport in the active and inactive states of these proteins, a master equation approach [, , 045103] is employed, which uses scaling rules that allow us to calculate energy transport rates solely based on the protein structure.

Nonequilibrium Modeling of the Elementary Step in PDZ3 Allosteric Communication.

While allostery is of paramount importance for protein signaling and regulation, the underlying dynamical process of allosteric communication is not well understood. The PDZ3 domain represents a prime example of an allosteric single-domain protein, as it features a well-established long-range coupling between the C-terminal α-helix and ligand binding. In an intriguing experiment, Hamm and co-workers employed photoswitching of the α-helix to initiate a conformational change of PDZ3 that propagates from the C-terminus to the bound ligand within 200 ns.

Ligand Unbinding Pathway and Mechanism Analysis Assisted by Machine Learning and Graph Methods.

We present two methods to reveal protein-ligand unbinding mechanisms in biased unbinding simulations by clustering trajectories into ensembles representing unbinding paths. The first approach is based on a contact principal component analysis for reducing the dimensionality of the input data, followed by identification of unbinding paths and training a machine learning model for trajectory clustering. The second approach clusters trajectories according to their pairwise mean Euclidean distance employing the neighbor-net algorithm, which takes into account input data bias in the distances set and is superior to dendrogram construction.

Cooperative Protein Allosteric Transition Mediated by a Fluctuating Transmission Network.

Allosteric communication between distant protein sites represents a key mechanism of biomolecular regulation and signal transduction. Compared to other processes such as protein folding, however, the dynamical evolution of allosteric transitions is still not well understood. As an example of allosteric coupling between distant protein regions, we consider the global open-closed motion of the two domains of T4 lysozyme, which is triggered by local motion in the hinge region.

Molecular Origin of Driving-Dependent Friction in Fluids.

The friction coefficient of fluids may become a function of the velocity at increased external driving. This non-Newtonian behavior is of general theoretical interest and of great practical importance, for example, for the design of lubricants. Although the effect has been observed in large-scale atomistic simulations of bulk liquids, its theoretical formulation and microscopic origin are not well understood.

Temporal control of the integrated stress response by a stochastic molecular switch.

Stress granules (SGs) are formed in the cytosol as an acute response to environmental cues and activation of the integrated stress response (ISR), a central signaling pathway controlling protein synthesis. Using chronic virus infection as stress model, we previously uncovered a unique temporal control of the ISR resulting in recurrent phases of SG assembly and disassembly. Here, we elucidate the molecular network generating this fluctuating stress response by integrating quantitative experiments with mathematical modeling and find that the ISR operates as a stochastic switch.

Predicting Ion Channel Conductance via Dissipation-Corrected Targeted Molecular Dynamics and Langevin Equation Simulations.

Ion channels are important proteins for physiological information transfer and functional control. To predict the microscopic origins of their voltage-conductance characteristics, here we applied dissipation-corrected targeted molecular dynamics in combination with Langevin equation simulations to potassium diffusion through the gramicidin A channel as a test system. Performing a nonequilibrium principal component analysis on backbone dihedral angles, we find coupled protein-ion dynamics to occur during ion transfer.

finDr: A web server for D-peptide ligand identification.

In the rapidly expanding field of peptide therapeutics, the short half-life of peptides represents a considerable limitation for drug action. D-peptides, consisting entirely of the dextrorotatory enantiomers of naturally occurring levorotatory amino acids (AAs), do not suffer from these shortcomings as they are intrinsically resistant to proteolytic degradation, resulting in a favourable pharmacokinetic profile. To experimentally identify D-peptide binders to interesting therapeutic targets, so-called mirror-image phage display is typically performed, whereby the target is synthesized in D-form and L-peptide binders are screened as in conventional phage display.

Current approaches to fate mapping and lineage tracing using image data.

Visualizing, tracking and reconstructing cell lineages in developing embryos has been an ongoing effort for well over a century. Recent advances in light microscopy, labelling strategies and computational methods to analyse complex image datasets have enabled detailed investigations into the fates of cells. Combined with powerful new advances in genomics and single-cell transcriptomics, the field of developmental biology is able to describe the formation of the embryo like never before.

Data-Driven Langevin Modeling of Nonequilibrium Processes.

Given nonstationary data from molecular dynamics simulations, a Markovian Langevin model is constructed that aims to reproduce the time evolution of the underlying process. While at equilibrium the free energy landscape is sampled, nonequilibrium processes can be associated with a biased energy landscape, which accounts for finite sampling effects and external driving. When the data-driven Langevin equation (dLE) approach [.

Hierarchical dynamics in allostery following ATP hydrolysis monitored by single molecule FRET measurements and MD simulations.

We report on a study that combines advanced fluorescence methods with molecular dynamics (MD) simulations to cover timescales from nanoseconds to milliseconds for a large protein. This allows us to delineate how ATP hydrolysis in a protein causes allosteric changes at a distant protein binding site, using the chaperone Hsp90 as test system. The allosteric process occurs hierarchical dynamics involving timescales from nano- to milliseconds and length scales from Ångstroms to several nanometers.

Through bonds or contacts? Mapping protein vibrational energy transfer using non-canonical amino acids.

Vibrational energy transfer (VET) is essential for protein function. It is responsible for efficient energy dissipation in reaction sites, and has been linked to pathways of allosteric communication. While it is understood that VET occurs via backbone as well as via non-covalent contacts, little is known about the competition of these two transport channels, which determines the VET pathways.

Microscopy-based assay for semi-quantitative detection of SARS-CoV-2 specific antibodies in human sera: A semi-quantitative, high throughput, microscopy-based assay expands existing approaches to measure SARS-CoV-2 specific antibody levels in human sera.

Emergence of the novel pathogenic coronavirus SARS-CoV-2 and its rapid pandemic spread presents challenges that demand immediate attention. Here, we describe the development of a semi-quantitative high-content microscopy-based assay for detection of three major classes (IgG, IgA, and IgM) of SARS-CoV-2 specific antibodies in human samples. The possibility to detect antibodies against the entire viral proteome together with a robust semi-automated image analysis workflow resulted in specific, sensitive and unbiased assay that complements the portfolio of SARS-CoV-2 serological assays.

Real-time observation of ligand-induced allosteric transitions in a PDZ domain.

While allostery is of paramount importance for protein regulation, the underlying dynamical process of ligand (un)binding at one site, resulting time evolution of the protein structure, and change of the binding affinity at a remote site are not well understood. Here the ligand-induced conformational transition in a widely studied model system of allostery, the PDZ2 domain, is investigated by transient infrared spectroscopy accompanied by molecular dynamics simulations. To this end, an azobenzene-derived photoswitch is linked to a peptide ligand in a way that its binding affinity to the PDZ2 domain changes upon switching, thus initiating an allosteric transition in the PDZ2 domain protein.

Multisecond ligand dissociation dynamics from atomistic simulations.

Coarse-graining of fully atomistic molecular dynamics simulations is a long-standing goal in order to allow the description of processes occurring on biologically relevant timescales. For example, the prediction of pathways, rates and rate-limiting steps in protein-ligand unbinding is crucial for modern drug discovery. To achieve the enhanced sampling, we perform dissipation-corrected targeted molecular dynamics simulations, which yield free energy and friction profiles of molecular processes under consideration.

Biomechanical comparison of a 3.5 mm anterior locking plate to cannulated screws with anterior tension band wiring in comminuted patellar fractures.

Purpose: To date, surgically treated multifragmentary patellar fractures are still associated with high rates of complications, such as i.e. secondary fracture displacement.

The Mutex Watershed and its Objective: Efficient, Parameter-Free Graph Partitioning.

Image partitioning, or segmentation without semantics, is the task of decomposing an image into distinct segments, or equivalently to detect closed contours. Most prior work either requires seeds, one per segment; or a threshold; or formulates the task as multicut / correlation clustering, an NP-hard problem. Here, we propose an efficient algorithm for graph partitioning, the "Mutex Watershed".