An automated Bayesian pipeline for rapid analysis of single-molecule binding data

Probing actin-activated ATP turnover kinetics of human cardiac myosin II by single molecule fluorescence

Mechanistic insights into myosin II energy transduction in striated muscle in health and disease would benefit from functional studies of a wide range of point-mutants. This approach is, however, hampered by the slow turnaround of myosin II expression that usually relies on adenoviruses for gene transfer. A recently developed virus-free method is more time effective but would yield too small amounts of myosin for standard biochemical analyses. However, if the fluorescent adenosine triphosphate (ATP) and single molecule (sm) total internal reflection fluorescence microscopy previously used to analyze basal ATP turnover by myosin alone, can be expanded to actin-activated ATP turnover, it would appreciably reduce the required amount of myosin. To that end, we here describe zero-length cross-linking of human cardiac myosin II motor fragments (sub-fragment 1 long [S1L]) to surface-immobilized actin filaments in a configuration with maintained actin-activated ATP turnover. After optimizing the analysis of sm fluorescence events, we show that the amount of myosin produced from C2C12 cells in one 60 mm cell culture plate is sufficient to obtain both the basal myosin ATP turnover rate and the maximum actin-activated rate constant (kcat). Our analysis of many single binding events of fluorescent ATP to many S1L motor fragments revealed processes reflecting basal and actin-activated ATPase, but also a third exponential process consistent with non-specific ATP-binding outside the active site.

Quantifying the minimum localization uncertainty of image scanning localization microscopy

Modulation enhanced single-molecule localization microscopy (meSMLM), where emitters are sparsely activated with sequentially applied patterned illumination, increases the localization precision over single-molecule localization microscopy (SMLM). The precision improvement of modulation enhanced SMLM is derived from retrieving the position of an emitter relative to individual illumination patterns, which adds to existing point spread function information from SMLM. Here, we introduce SpinFlux: modulation enhanced localization for spinning disk confocal microscopy. SpinFlux uses a spinning disk with pinholes in its illumination and emission paths, to sequentially illuminate regions in the sample during each measurement. The resulting intensity-modulated emission signal is analyzed for each individual pattern to localize emitters with improved precision. We derive a statistical image formation model for SpinFlux and we quantify the theoretical minimum localization uncertainty in terms of the Cramér-Rao lower bound. Using the theoretical minimum uncertainty, we compare SpinFlux to localization on Fourier reweighted image scanning microscopy reconstructions. We find that localization on image scanning microscopy reconstructions with Fourier reweighting ideally results in a global precision improvement of 2.1 over SMLM. When SpinFlux is used for sequential illumination with three patterns around the emitter position, the localization precision improvement over SMLM is twofold when patterns are focused around the emitter position. If four donut-shaped illumination patterns are used for SpinFlux, the maximum local precision improvement over SMLM is increased to 3.5. Localization of image scanning microscopy reconstructions thus has the largest potential for global improvements of the localization precision, where SpinFlux is the method of choice for local refinements.

Structural role for DNA Ligase IV in promoting the fidelity of non-homologous end joining

Nonhomologous end joining (NHEJ), the primary pathway of vertebrate DNA double-strand-break (DSB) repair, directly re-ligates broken DNA ends. Damaged DSB ends that cannot be immediately re-ligated are modified by NHEJ processing enzymes, including error-prone polymerases and nucleases, to enable ligation. However, DSB ends that are initially compatible for re-ligation are typically joined without end processing. As both ligation and end processing occur in the short-range (SR) synaptic complex that closely aligns DNA ends, it remains unclear how ligation of compatible ends is prioritized over end processing. In this study, we identify structural interactions of the NHEJ-specific DNA Ligase IV (Lig4) within the SR complex that prioritize ligation and promote NHEJ fidelity. Mutational analysis demonstrates that Lig4 must bind DNA ends to form the SR complex. Furthermore, single-molecule experiments show that a single Lig4 binds both DNA ends at the instant of SR synapsis. Thus, Lig4 is poised to ligate compatible ends upon initial formation of the SR complex before error-prone processing. Our results provide a molecular basis for the fidelity of NHEJ.

Single-molecule tracking of Nanog and Oct4 in living mouse embryonic stem cells uncovers a feedback mechanism of pluripotency maintenance

Nanog and Oct4 are core transcription factors that form part of a gene regulatory network to regulate hundreds of target genes for pluripotency maintenance in mouse embryonic stem cells (ESCs). To understand their function in the pluripotency maintenance, we visualised and quantified the dynamics of single molecules of Nanog and Oct4 in a mouse ESCs during pluripotency loss. Interestingly, Nanog interacted longer with its target loci upon reduced expression or at the onset of differentiation, suggesting a feedback mechanism to maintain the pluripotent state. The expression level and interaction time of Nanog and Oct4 correlate with their fluctuation and interaction frequency, respectively, which in turn depend on the ESC differentiation status. The DNA viscoelasticity near the Oct4 target locus remained flexible during differentiation, supporting its role either in chromatin opening or a preferred binding to uncondensed chromatin regions. Based on these results, we propose a new negative feedback mechanism for pluripotency maintenance via the DNA condensation state-dependent interplay of Nanog and Oct4.

ExTrack characterizes transition kinetics and diffusion in noisy single-particle tracks

Single-particle tracking microscopy is a powerful technique to investigate how proteins dynamically interact with their environment in live cells. However, the analysis of tracks is confounded by noisy molecule localization, short tracks, and rapid transitions between different motion states, notably between immobile and diffusive states. Here, we propose a probabilistic method termed ExTrack that uses the full spatio-temporal information of tracks to extract global model parameters, to calculate state probabilities at every time point, to reveal distributions of state durations, and to refine the positions of bound molecules. ExTrack works for a wide range of diffusion coefficients and transition rates, even if experimental data deviate from model assumptions. We demonstrate its capacity by applying it to slowly diffusing and rapidly transitioning bacterial envelope proteins. ExTrack greatly increases the regime of computationally analyzable noisy single-particle tracks. The ExTrack package is available in ImageJ and Python.

ExTrack characterizes transition kinetics and diffusion in noisy single-particle tracks

Single-particle tracking microscopy is a powerful technique to investigate how proteins dynamically interact with their environment in live cells. However, the analysis of tracks is confounded by noisy molecule localization, short tracks, and rapid transitions between different motion states, notably between immobile and diffusive states. Here, we propose a probabilistic method termed ExTrack that uses the full spatio-temporal information of tracks to extract global model parameters, to calculate state probabilities at every time point, to reveal distributions of state durations, and to refine the positions of bound molecules. ExTrack works for a wide range of diffusion coefficients and transition rates, even if experimental data deviate from model assumptions. We demonstrate its capacity by applying it to slowly diffusing and rapidly transitioning bacterial envelope proteins. ExTrack greatly increases the regime of computationally analyzable noisy single-particle tracks. The ExTrack package is available in ImageJ and Python.

SOLEIL: single-objective lens inclined light sheet localization microscopy

High-NA light sheet illumination can improve the resolution of single-molecule localization microscopy (SMLM) by reducing the background fluorescence. These approaches currently require custom-made sample holders or additional specialized objectives, which makes the sample mounting or the optical system complex and therefore reduces the usability of these approaches. Here, we developed a single-objective lens-inclined light sheet microscope (SOLEIL) that is capable of 2D and 3D SMLM in thick samples. SOLEIL combines oblique illumination with point spread function PSF engineering to enable dSTORM imaging in a wide variety of samples. SOLEIL is compatible with standard sample holders and off-the-shelve optics and standard high NA objectives. To accomplish optimal optical sectioning we show that there is an ideal oblique angle and sheet thickness. Furthermore, to show what optical sectioning delivers for SMLM we benchmark SOLEIL against widefield and HILO microscopy with several biological samples. SOLEIL delivers in 15 μm thick Caco2-BBE cells a 374% higher intensity to background ratio and a 54% improvement in the estimated CRLB compared to widefield illumination, and a 184% higher intensity to background ratio and a 20% improvement in the estimated CRLB compared to HILO illumination.

Precision in iterative modulation enhanced single-molecule localization microscopy

Modulation enhanced single-molecule localization microscopy (meSMLM) methods improve the localization precision by using patterned illumination to encode additional position information. Iterative meSMLM (imeSMLM) methods iteratively generate prior information on emitter positions, used to locally improve the localization precision during subsequent iterations. The Cramér-Rao lower bound (CRLB) cannot incorporate prior information to bound the best achievable localization precision, because it requires estimators to be unbiased. By treating estimands as random variables with a known prior distribution, the Van Trees inequality (VTI) can be used to bound the best possible localization precision of imeSMLM methods. An imeSMLM method is considered, where the positions of in-plane standing wave illumination patterns are controlled over the course of multiple iterations. Using the VTI, we analytically approximate a lower bound on the maximum localization precision of imeSMLM methods that make use of standing wave illumination patterns. In addition, we evaluate the maximally achievable localization precision for different illumination pattern placement strategies using Monte Carlo simulations. We show that in the absence of background and under perfect modulation, the information content of signal photons increases exponentially as a function of the iteration count. However, the information increase is no longer exponential as a function of the iteration count under non-zero background, imperfect modulation or limited mechanical resolution of the illumination positioning system. As a result, imeSMLM with two iterations reaches at most a 5-fold improvement over SMLM at 8 expected background photons per pixel and 95% modulation contrast. Moreover, the information increase from imeSMLM is balanced by a reduced signal photon rate. Therefore, SMLM outperforms imeSMLM when considering an equal measurement time and illumination power per iteration. Finally, the VTI is an excellent tool for the assessment of the performance of illumination control and is therefore the method of choice for optimal design and control of imeSMLM methods.

Principles and pitfalls of high-throughput analysis of microRNA-binding thermodynamics and kinetics by RNA Bind-n-Seq

RNA Bind-n-Seq (RBNS) is a cost-effective, high-throughput method capable of identifying the sequence preferences of RNA-binding proteins and of qualitatively defining relative dissociation constants. Although RBNS is often described as an unbiased method, several factors may influence the outcome of the analysis. Here, we discuss these biases and present an analytical strategy to estimate absolute binding affinities from RBNS data, extend RBNS to kinetic studies, and develop a framework to compute relative association and dissociation rate constants. As proof of principle, we measured the equilibrium binding properties of mammalian Argonaute2 (AGO2) guided by eight microRNAs (miRNAs) and kinetic parameters for let-7a. The miRNA-binding site repertoires, dissociation constants, and kinetic parameters calculated from RBNS data using our methods correlate well with values measured by traditional ensemble and single-molecule approaches. Our data provide additional quantitative measurements for Argonaute-bound miRNA binding that should facilitate development of quantitative targeting rules for individual miRNAs.

High-throughput biochemical profiling reveals functional adaptation of a bacterial Argonaute

Argonautes are nucleic acid-guided proteins that perform numerous cellular functions across all domains of life. Little is known about how distinct evolutionary pressures have shaped each Argonaute's biophysical properties. We applied high-throughput biochemistry to characterize how Thermus thermophilus Argonaute (TtAgo), a DNA-guided DNA endonuclease, finds, binds, and cleaves its targets. We found that TtAgo uses biophysical adaptations similar to those of eukaryotic Argonautes for rapid association but requires more extensive complementarity to achieve high-affinity target binding. Using these data, we constructed models for TtAgo association rates and equilibrium binding affinities that estimate the nucleic acid- and protein-mediated components of the target interaction energies. Finally, we showed that TtAgo cleavage rates vary widely based on the DNA guide, suggesting that only a subset of guides cleaves targets on physiologically relevant timescales.

Bayesian machine learning analysis of single-molecule fluorescence colocalization images

Multi-wavelength single-molecule fluorescence colocalization (CoSMoS) methods allow elucidation of complex biochemical reaction mechanisms. However, analysis of CoSMoS data is intrinsically challenging because of low image signal-to-noise ratios, non-specific surface binding of the fluorescent molecules, and analysis methods that require subjective inputs to achieve accurate results. Here, we use Bayesian probabilistic programming to implement Tapqir, an unsupervised machine learning method that incorporates a holistic, physics-based causal model of CoSMoS data. This method accounts for uncertainties in image analysis due to photon and camera noise, optical non-uniformities, non-specific binding, and spot detection. Rather than merely producing a binary 'spot/no spot' classification of unspecified reliability, Tapqir objectively assigns spot classification probabilities that allow accurate downstream analysis of molecular dynamics, thermodynamics, and kinetics. We both quantitatively validate Tapqir performance against simulated CoSMoS image data with known properties and also demonstrate that it implements fully objective, automated analysis of experiment-derived data sets with a wide range of signal, noise, and non-specific binding characteristics.

Drift correction in localization microscopy using entropy minimization

Localization microscopy offers resolutions down to a single nanometer but currently requires additional dedicated hardware or fiducial markers to reduce resolution loss from the drift of the sample. Drift estimation without fiducial markers is typically implemented using redundant cross correlation (RCC). We show that RCC has sub-optimal precision and bias, which leaves room for improvement. Here, we minimize a bound on the entropy of the obtained localizations to efficiently compute a precise drift estimate. Within practical compute-time constraints, simulations show a 5x improvement in drift estimation precision over the widely used RCC algorithm. The algorithm operates directly on fluorophore localizations and is tested on simulated and experimental datasets in 2D and 3D. An open source implementation is provided, implemented in Python and C++, and can utilize a GPU if available.

KERA: analysis tool for multi-process, multi-state single-molecule data

Molecular machines within cells dynamically assemble, disassemble and reorganize. Molecular interactions between their components can be observed at the single-molecule level and quantified using colocalization single-molecule spectroscopy, in which individual labeled molecules are seen transiently associating with a surface-tethered partner, or other total internal reflection fluorescence microscopy approaches in which the interactions elicit changes in fluorescence in the labeled surface-tethered partner. When multiple interacting partners can form ternary, quaternary and higher order complexes, the types of spatial and temporal organization of these complexes can be deduced from the order of appearance and reorganization of the components. Time evolution of complex architectures can be followed by changes in the fluorescence behavior in multiple channels. Here, we describe the kinetic event resolving algorithm (KERA), a software tool for organizing and sorting the discretized fluorescent trajectories from a range of single-molecule experiments. KERA organizes the data in groups by transition patterns, and displays exhaustive dwell time data for each interaction sequence. Enumerating and quantifying sequences of molecular interactions provides important information regarding the underlying mechanism of the assembly, dynamics and architecture of the macromolecular complexes. We demonstrate KERA's utility by analyzing conformational dynamics of two DNA binding proteins: replication protein A and xeroderma pigmentosum complementation group D helicase.

Bayesian Inference: The Comprehensive Approach to Analyzing Single-Molecule Experiments

Biophysics experiments performed at single-molecule resolution provide exceptional insight into the structural details and dynamic behavior of biological systems. However, extracting this information from the corresponding experimental data unequivocally requires applying a biophysical model. In this review, we discuss how to use probability theory to apply these models to single-molecule data. Many current single-molecule data analysis methods apply parts of probability theory, sometimes unknowingly, and thus miss out on the full set of benefits provided by this self-consistent framework. The full application of probability theory involves a process called Bayesian inference that fully accounts for the uncertainties inherent to single-molecule experiments. Additionally, using Bayesian inference provides a scientifically rigorous method of incorporating information from multiple experiments into a single analysis and finding the best biophysical model for an experiment without the risk of overfitting the data. These benefits make the Bayesian approach ideal for analyzing any type of single-molecule experiment.

Thermus thermophilus Argonaute Functions in the Completion of DNA Replication

In many eukaryotes, Argonaute proteins, guided by short RNA sequences, defend cells against transposons and viruses. In the eubacterium Thermus thermophilus, the DNA-guided Argonaute TtAgo defends against transformation by DNA plasmids. Here, we report that TtAgo also participates in DNA replication. In vivo, TtAgo binds 15- to 18-nt DNA guides derived from the chromosomal region where replication terminates and associates with proteins known to act in DNA replication. When gyrase, the sole T. thermophilus type II topoisomerase, is inhibited, TtAgo allows the bacterium to finish replicating its circular genome. In contrast, loss of gyrase and TtAgo activity slows growth and produces long sausage-like filaments in which the individual bacteria are linked by DNA. Finally, wild-type T. thermophilus outcompetes an otherwise isogenic strain lacking TtAgo. We propose that the primary role of TtAgo is to help T. thermophilus disentangle the catenated circular chromosomes generated by DNA replication.

Top-down machine learning approach for high-throughput single-molecule analysis

Single-molecule approaches provide enormous insight into the dynamics of biomolecules, but adequately sampling distributions of states and events often requires extensive sampling. Although emerging experimental techniques can generate such large datasets, existing analysis tools are not suitable to process the large volume of data obtained in high-throughput paradigms. Here, we present a new analysis platform (DISC) that accelerates unsupervised analysis of single-molecule trajectories. By merging model-free statistical learning with the Viterbi algorithm, DISC idealizes single-molecule trajectories up to three orders of magnitude faster with improved accuracy compared to other commonly used algorithms. Further, we demonstrate the utility of DISC algorithm to probe cooperativity between multiple binding events in the cyclic nucleotide binding domains of HCN pacemaker channel. Given the flexible and efficient nature of DISC, we anticipate it will be a powerful tool for unsupervised processing of high-throughput data across a range of single-molecule experiments.

Argonaute bypasses cellular obstacles without hindrance during target search

Argonaute (Ago) proteins are key players in both gene regulation (eukaryotes) and host defense (prokaryotes). Acting on single-stranded nucleic-acid substrates, Ago relies on base pairing between a small nucleic-acid guide and its complementary target sequences for specificity. To efficiently scan nucleic-acid chains for targets, Ago diffuses laterally along the substrate and must bypass secondary structures as well as protein barriers. Using single-molecule FRET in conjunction with kinetic modelling, we reveal that target scanning is mediated through loose protein-nucleic acid interactions, allowing Ago to slide short distances over secondary structures, as well as to bypass protein barriers via intersegmental transfer. Our combined single-molecule experiment and kinetic modelling approach may serve as a platform to dissect search processes and study the effect of sequence on search kinetics for other nucleic acid-guided proteins.

FISIK: Framework for the Inference of In Situ Interaction Kinetics from Single-Molecule Imaging Data

Recent experimental and computational developments have been pushing the limits of live-cell single-molecule imaging, enabling the monitoring of intermolecular interactions in their native environment with high spatiotemporal resolution. However, interactions are captured only for the labeled subset of molecules, which tends to be a small fraction. As a result, it has remained a challenge to calculate molecular interaction kinetics, in particular association rates, from live-cell single-molecule tracking data. To overcome this challenge, we developed a mathematical modeling-based Framework for the Inference of in Situ Interaction Kinetics (FISIK) from single-molecule imaging data with substoichiometric labeling. FISIK consists of (I) devising a mathematical model of molecular movement and interactions, mimicking the biological system and data-acquisition setup, and (II) estimating the unknown model parameters, including molecular association and dissociation rates, by fitting the model to experimental single-molecule data. Due to the stochastic nature of the model and data, we adapted the method of indirect inference for model calibration. We validated FISIK using a series of tests in which we simulated trajectories of diffusing molecules that interact with each other, considering a wide range of model parameters, and including resolution limitations, tracking errors, and mismatches between the model and the biological system it mimics. We found that FISIK has the sensitivity to determine association and dissociation rates, with accuracy and precision depending on the labeled fraction of molecules and the extent of molecule tracking errors. For cases where the labeled fraction is too low (e.g., to afford accurate tracking), combining dynamic but sparse single-molecule imaging data with almost-whole population oligomer distribution data improves FISIK's performance. All in all, FISIK is a promising approach for the derivation of molecular interaction kinetics in their native environment from single-molecule imaging data with substoichiometric labeling.