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. In this Review, we discuss some of the different strategies and applications to lineage tracing in live-imaging data and outline software methodologies that can be applied to various cell-tracking challenges.

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. To demonstrate the approach, nonequilibrium transport observables are calculated for a molecular nanojunction containing strong interactions.

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. This observation holds true for friction extracted from Kramers' rate expression and by dissipation-corrected targeted molecular dynamics simulations based on Jarzynski's identity. We rule out ligand solvation and protein-internal friction as sources of the angle-dependent friction. Instead, we observe a heterogeneity in free energy barriers along an experimentally uncontrolled orientation parameter, which increases the rupture force variance and therefore the overall friction. We anticipate that anisotropic friction needs to be accounted for in a complete understanding of friction in biomolecular dynamics and anisotropic mechanical environments.

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. Using dissipation-corrected targeted molecular dynamics in combination with temperature-boosted Langevin equation simulations [S. Wolf et al., Nat. Commun. 11, 2918 (2020)] applied to a two-dimensional model and the trypsin-benzamidine complex as test systems, we assess the robustness of the procedure and discuss the aspects of its practical applicability to predict multisecond kinetics of complex biomolecular systems.

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. Sensitive, specific and quantitative serological assays are urgently needed for a better understanding of humoral immune response against the virus as a basis for developing public health strategies to control viral spread. The procedure described here has been used for clinical studies and provides a general framework for the application of quantitative high-throughput microscopy to rapidly develop serological assays for emerging virus infections.

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. The subsequent response of the protein, covering four decades of time, ranging from ∼1 ns to ∼µs, can be rationalized by a remodeling of its rugged free-energy landscape, with very subtle shifts in the populations of a small number of structurally well-defined states. It is proposed that structurally and dynamically driven allostery, often discussed as limiting scenarios of allosteric communication, actually go hand-in-hand, allowing the protein to adapt its free-energy landscape to incoming signals.

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. Finally, we describe a more complex case where the reaction coordinate relevant for path identification is a single intraligand hydrogen bond, highlighting the challenges involved in unbinding path reaction coordinate detection.

An objective comparison of cell-tracking algorithms.

We present a combined report on the results of three editions of the Cell Tracking Challenge, an ongoing initiative aimed at promoting the development and objective evaluation of cell segmentation and tracking algorithms. With 21 participating algorithms and a data repository consisting of 13 data sets from various microscopy modalities, the challenge displays today's state-of-the-art methodology in the field. We analyzed the challenge results using performance measures for segmentation and tracking that rank all participating methods. We also analyzed the performance of all of the algorithms in terms of biological measures and practical usability. Although some methods scored high in all technical aspects, none obtained fully correct solutions. We found that methods that either take prior information into account using learning strategies or analyze cells in a global spatiotemporal video context performed better than other methods under the segmentation and tracking scenarios included in the challenge.

Photocontrolling Protein-Peptide Interactions: From Minimal Perturbation to Complete Unbinding.

An azobenzene-derived photoswitch has been covalently cross-linked to two sites of the S-peptide in the RNase S complex in a manner that the α-helical content of the S-peptide reduces upon cis-to-trans isomerization of the photoswitch. Three complementary experimental techniques have been employed, isothermal titration calorimetry, circular dichroism spectroscopy and intrinsic tyrosine fluorescence quenching, to determine the binding affinity of the S-peptide to the S-protein in the two states of the photoswitch. Five mutants with the photoswitch attached to different sites of the S-peptide have been explored, with the goal to maximize the change in binding affinity upon photoswitching, and to identify the mechanisms that determine the binding affinity. With regard to the first goal, one mutant has been identified, which binds with reasonable affinity in the one state of the photoswitch, while specific binding is completely switched off in the other state. With regard to the second goal, accompanying molecular dynamics simulations combined with a quantitative structure activity relationship revealed that the α-helicity of the S-peptide in the binding pocket correlates surprisingly well with measured dissociation constants. Moreover, the simulations show that both configurations of all S-peptides exhibit quite well-defined structures, even in apparently disordered states.

Estimation of Protein-Ligand Unbinding Kinetics Using Non-Equilibrium Targeted Molecular Dynamics Simulations.

We here report on nonequilibrium targeted molecular dynamics simulations as a tool for the estimation of protein-ligand unbinding kinetics. Correlating simulations with experimental data from SPR kinetics measurements and X-ray crystallography on two small molecule compound libraries bound to the N-terminal domain of the chaperone Hsp90, we show that the mean nonequilibrium work computed in an ensemble of trajectories of enforced ligand unbinding is a promising predictor for ligand unbinding rates. We furthermore investigate the molecular basis determining unbinding rates within the compound libraries. We propose ligand conformational changes and protein-ligand nonbonded interactions to impact on unbinding rates. Ligands may remain longer at the protein if they exhibit strong electrostatic and/or van der Waals interactions with the target. In the case of ligands with a rigid chemical scaffold that exhibit longer residence times, transient electrostatic interactions with the protein appear to facilitate unbinding. Our results imply that understanding the unbinding pathway and the protein-ligand interactions along this path is crucial for the prediction of small molecule ligands with defined unbinding kinetics.

Principal component analysis of nonequilibrium molecular dynamics simulations.

Principal component analysis (PCA) represents a standard approach to identify collective variables {xi} = x, which can be used to construct the free energy landscape ΔG(x) of a molecular system. While PCA is routinely applied to equilibrium molecular dynamics (MD) simulations, it is less obvious as to how to extend the approach to nonequilibrium simulation techniques. This includes, e.g., the definition of the statistical averages employed in PCA as well as the relation between the equilibrium free energy landscape ΔG(x) and the energy landscapes ΔG(x) obtained from nonequilibrium MD. As an example for a nonequilibrium method, "targeted MD" is considered which employs a moving distance constraint to enforce rare transitions along some biasing coordinate s. The introduced bias can be described by a weighting function P(s), which provides a direct relation between equilibrium and nonequilibrium data, and thus establishes a well-defined way to perform PCA on nonequilibrium data. While the resulting distribution P(x) and energy ΔG∝lnP will not reflect the equilibrium state of the system, the nonequilibrium energy landscape ΔG(x) may directly reveal the molecular reaction mechanism. Applied to targeted MD simulations of the unfolding of decaalanine, for example, a PCA performed on backbone dihedral angles is shown to discriminate several unfolding pathways. Although the formulation is in principle exact, its practical use depends critically on the choice of the biasing coordinate s, which should account for a naturally occurring motion between two well-defined end-states of the system.

2D-IR Spectroscopy of an AHA Labeled Photoswitchable PDZ2 Domain.

We explore the capability of the non-natural amino acid azidohomoalanine (AHA) as an IR label to sense relatively small structural changes in proteins with the help of 2D IR difference spectroscopy. To that end, we AHA-labeled an allosteric protein (the PDZ2 domain from human tyrosine-phosphatase 1E) and furthermore covalently linked it to an azobenzene-derived photoswitch as to mimic its conformational transition upon ligand binding. To determine the strengths and limitations of the AHA label, in total six mutants have been investigated with the label at sites with varying properties. Only one mutant revealed a measurable 2D IR difference signal. In contrast to the commonly observed frequency shifts that report on the degree of solvation, in this case we observe an intensity change. To understand this spectral response, we performed classical MD simulations, evaluating local contacts of the AHA labels to water molecules and protein side chains and calculating the vibrational frequency on the basis of an electrostatic model. Although these simulations revealed in part significant and complex changes of the number of intraprotein and water contacts upon trans-cis photoisomerization, they could not provide a clear explanation of why this one label would stick out. Subsequent quantum-chemistry calculations suggest that the response is the result of an electronic interaction involving charge transfer of the azido group with sulfonate groups from the photoswitch. To the best of our knowledge, such an effect has not been described before.

Radioablation of liver malignancies with interstitial high-dose-rate brachytherapy : Complications and risk factors.

BACKGROUND: To evaluate complications and identify risk factors for adverse events in patients undergoing high-dose-rate interstitial brachytherapy (iBT). MATERIAL AND METHODS: Data from 192 patients treated in 343 CT- or MRI-guided interventions from 2006-2009 at our institution were analyzed. In 41 %, the largest tumor treated was ≥ 5 cm, 6 % of the patients had tumors ≥ 10 cm. Prior to iBT, 60 % of the patients had chemotherapy, 22 % liver resection, 19 % thermoablation or transarterial chemoembolization (TACE). Safety was the primary endpoint; survival data were obtained as the secondary endpoints. During follow-up, MRI or CT imaging was performed and clinical and laboratory parameters were obtained. RESULTS: The rate of major complications was below 5 %. Five major bleedings (1.5 %) occurred. The frequency of severe bleeding was significantly higher in patients with advanced liver cirrhosis. One patient developed signs of a nonclassic radiation-induced liver disease. In 3 patients, symptomatic gastrointestinal (GI) ulcers were detected. A dose exposure to the GI wall above 14 Gy/ml was a reliable threshold to predict ulcer formation. A combination of C-reactive protein ≥ 165 mg/l and/or leukocyte count ≥ 12.7 Gpt/l on the second day after the intervention predicted infection (sensitivity 90.0 %; specificity 92.8 %.) Two patients (0.6 %) died within 30 days. Median overall survival after the first liver treatment was 20.1 months for all patients and the local recurrence-free surviving proportion was 89 % after 12 months. CONCLUSIONS: Image-guided iBT yields a low rate of major complications and is effective.

Shape manipulation of ion irradiated Ag nanoparticles embedded in lithium niobate.

Spherical silver nanoparticles were prepared by means of ion beam synthesis in lithium niobate. The embedded nanoparticles were then irradiated with energetic (84)Kr and (197)Au ions, resulting in different electronic energy losses between 8.1 and 27.5 keV nm(-1) in the top layer of the samples. Due to the high electronic energy losses of the irradiating ions, molten ion tracks are formed inside the lithium niobate in which the elongated Ag nanoparticles are formed. This process is strongly dependent on the initial particle size and leads to a broad aspect ratio distribution. Extinction spectra of the samples feature the extinction maximum with shoulders on either side. While the maximum is caused by numerous remaining spherical nanoparticles, the shoulders can be attributed to elongated particles. The latter could be verified by COMSOL simulations. The extinction spectra are thus a superposition of the spectra of all individual particles.

7TM Domain Structure of Adhesion GPCRs.

Schematic presentation of the overall adhesion G Protein-Coupled Receptor (aGPCR) structure and functional domains, covering an extracellular N-terminal fragment (NTF), a membrane-spanning C-terminal fragment (CTF) and a GPCR proteolysis site (GPS). (Left side) aGPCR model constructed based on the seven-transmembrane (7TM) structure (blue) of secretin family glucagon receptor (GCGR) (PDB, 4L6R) [11] and the GPCR autoproteolysis inducing (GAIN) domain (magenta) structure of latrophilin 1 (PDB, 4DLQ) [9]. The ß-13 strand residues are depicted in green. (Right side) The experimentally validated full-length secretin family GCGR structure combining structural and experimental information from the GCGR 7TM crystal structure (PDB, 4L6R) (blue), the GCGR extracellular domain (ECD) structure (PDB, 4ERS) (magenta) and the ECD structure of glucagon-like peptide-1 (GLP-1)-bound glucagon-like peptide-1 receptor (GLP-1R) (PDB, 3IOL) (green), complemented by site-directed mutagenesis, electron microscopy (EM), hydrogen-deuterium exchange (HDX) and cross-linking studies [11-13]) Despite the recent breakthroughs in the elucidation of the three-dimensional structures of the seven transmembrane (7TM) domain of the G protein-coupled receptor (GPCR) superfamily, a corresponding structure of a member of the adhesion GPCR (aGPCR) family has not yet been solved. In this chapter, we give an overview of the current knowledge of the 7TM domain of aGPCRs by comparative structure-based sequence similarity analyses between aGPCRs and GPCRs with known crystal structure. Of the GPCR superfamily, only the secretin family shares some sequence similarity with aGPCRs. This chapter will therefore emphasize on the comparison of these two GPCR families. Two 7TM domain structures of secretin family GPCRs are known that provide insight into the structure-function relationships of conserved sequence motifs that play important roles and are also present in most aGPCRs. This suggests that the 7TM domains of aGPCRs and secretin family GPCRs share a similar structural fold and that the conserved residues in both families may be involved in similar intermolecular interaction networks and facilitate similar conformational changes. Comparison of the residues that line the large peptide hormone binding pocket in the 7TM domain of secretin family GPCRs with corresponding residues in aGPCRs indicates that in the latter, the corresponding pocket in the 7TM domain is relatively hydrophobic and may be even larger. Improved knowledge on these conserved sequence motifs will help to understand the interactions of the aGPCR 7TM domain with ligands and gain insight into the activation mechanism of aGPCRs.

The role of protein-bound water molecules in microbial rhodopsins.

Protein-bound internal water molecules are essential features of the structure and function of microbial rhodopsins. Besides structural stabilization, they act as proton conductors and even proton storage sites. Currently, the most understood model system exhibiting such features is bacteriorhodopsin (bR). During the last 20 years, the importance of water molecules for proton transport has been revealed through this protein. It has been shown that water molecules are as essential as amino acids for proton transport and biological function. In this review, we present an overview of the historical development of this research on bR. We furthermore summarize the recently discovered protein-bound water features associated with proton transport. Specifically, we discuss a pentameric water/amino acid arrangement close to the protonated Schiff base as central proton-binding site, a protonated water cluster as proton storage site at the proton-release site, and a transient linear water chain at the proton uptake site. We highlight how protein conformational changes reposition or reorient internal water molecules, thereby guiding proton transport. Last, we compare the water positions in bR with those in other microbial rhodopsins to elucidate how protein-bound water molecules guide the function of microbial rhodopsins. This article is part of a Special Issue entitled: Retinal Proteins - You can teach an old dog new tricks.

A delocalized proton-binding site within a membrane protein.

The role of protein-bound water molecules in protein function and catalysis is an emerging topic. Here, we studied the solvation of an excess proton by protein-bound water molecules and the contribution of the surrounding amino acid residues at the proton release site of the membrane protein bacteriorhodopsin. It hosts an excess proton within a protein-bound water cluster, which is hydrogen bonded to several surrounding amino acids. Indicative of delocalization is a broad continuum absorbance experimentally observed by time-resolved Fourier transform infrared spectroscopy. In combination with site-directed mutagenesis, the involvement of several amino acids (especially Glu-194 and Glu-204) in the delocalization was elaborated. Details regarding the contributions of the glutamates and water molecules to the delocalization mode in biomolecular simulations are controversial. We carried out quantum mechanics/molecular mechanics (QM/MM) self-consistent charge density functional tight-binding simulations for all amino acids that have been experimentally shown to be involved in solvation of the excess proton, and systematically investigated the influence of the quantum box size. We compared calculated theoretical infrared spectra with experimental ones as a measure for the correct description of excess proton delocalization. A continuum absorbance can only be observed for small quantum boxes containing few amino acids and/or water molecules. Larger quantum boxes, including all experimentally shown involved amino acids, resulted in narrow absorbance bands, indicating protonation of a single binding site in contradiction to experimental results. We conclude that small quantum boxes seem to reproduce representative extreme cases of proton delocalization modes: proton delocalization only on water molecules or only between Glu-194 and Glu-204. Extending the experimental spectral region to lower wave numbers, a water-delocalized proton reproduces the observed continuum absorbance better than a glutamate-shared delocalized proton. However, a full agreement between QM simulations and experimental results on the delocalized excess proton will require a larger quantum box as well as more sophisticated QM/MM methods.

Utilizing dynamic annealing during ion implantation: synthesis of silver nanoparticles in crystalline lithium niobate.

Silver nanoparticles (NPs) embedded in lithium niobate were fabricated via ion beam synthesis and are suitable for various plasmonic applications, e.g. enhancement of optical nonlinear effects. After room temperature silver implantation, annealing in the temperature range of 400-600 °C was performed in order to recrystallize the damaged lithium niobate surface layer. The shape of the silver NPs, their optical properties as well as the structural properties of their surrounding matrix have been analyzed for various annealing steps. TEM investigations show that annealing at 400 °C does not lead to recrystallization of the damaged lithium niobate. A recrystallization occurs upon increasing the annealing temperature to 500 or 600 °C, but simultaneously a second phase consisting of lithium triniobate forms. This is additionally supported by XRD measurements. By utilizing dynamic annealing, i.e. implanting silver at elevated temperatures of 400 °C, it is shown that the LiNbO3 matrix stays single crystalline during ion implantation and no LiNb3O8 is formed. This is additionally verified by comparing the positions of the surface plasmon resonances with calculations based on Mie's scattering theory.

Infrared spectral marker bands characterizing a transient water wire inside a hydrophobic membrane protein.

Proton conduction along protein-bound "water wires" is an essential feature in membrane proteins. Here, we analyze in detail a transient water wire, which conducts protons via a hydrophobic barrier within a membrane protein to create a proton gradient. It is formed only for a millisecond out of three water molecules distributed at inactive positions in a polar environment in the ground state. The movement into a hydrophobic environment causes characteristic shifts of the water bands reflecting their different chemical properties. These band shifts are identified by time-resolved Fourier Transform Infrared difference spectroscopy and analyzed by biomolecular Quantum Mechanical/Molecular Mechanical simulations. A non-hydrogen bonded ("dangling") O-H stretching vibration band and a broad continuum absorbance caused by a combined vibration along the water wire are identified as characteristic marker bands of such water wires in a hydrophobic environment. The results provide a basic understanding of water wires in hydrophobic environments.