PMC:5861423
Annnotations
{"target":"https://pubannotation.org/docs/sourcedb/PMC/sourceid/5861423","sourcedb":"PMC","sourceid":"5861423","source_url":"https://www.ncbi.nlm.nih.gov/pmc/5861423","text":"Single-particle trajectories reveal two-state diffusion-kinetics of hOGG1 proteins on DNA\n\nAbstract\nAbstract\nWe reanalyze trajectories of hOGG1 repair proteins diffusing on DNA. A previous analysis of these trajectories with the popular mean-squared-displacement approach revealed only simple diffusion. Here, a new optimal estimator of diffusion coefficients reveals two-state kinetics of the protein. A simple, solvable model, in which the protein randomly switches between a loosely bound, highly mobile state and a tightly bound, less mobile state is the simplest possible dynamic model consistent with the data. It yields accurate estimates of hOGG1’s (i) diffusivity in each state, uncorrupted by experimental errors arising from shot noise, motion blur and thermal fluctuations of the DNA; (ii) rates of switching between states and (iii) rate of detachment from the DNA. The protein spends roughly equal time in each state. It detaches only from the loosely bound state, with a rate that depends on pH and the salt concentration in solution, while its rates for switching between states are insensitive to both. The diffusivity in the loosely bound state depends primarily on pH and is three to ten times higher than in the tightly bound state. We propose and discuss some new experiments that take full advantage of the new tools of analysis presented here.\n\nINTRODUCTION\nMolecular fluorescent labels and super-resolution microscopy allow us to track single biomolecules in cells (1). There, diffusion is ubiquitous, as many cellular processes rely on diffusion for transport (2). A precise understanding of such processes requires a precise determination of diffusion constants. Less than that may miss process-specific details by lumping them into one, simple diffusive process, we show below.\nRecent examples of experimental measurements of diffusion in biological systems include supercoils on DNA (3), proteins on biopolymers such as DNA (4–9) or microtubules (10), on surfaces (11), in natural (12–14) and artificial (15) lipid membranes, in films (16) and inside cells (17,18), all recorded with time-lapse photography.\nDiffusion of proteins on DNA captured the attention of biophysicists nearly half a century ago when it was observed that the Lac repressor was able to locate its cognate site in vitro even faster than theoretically predicted for 3D diffusion in bulk (19). That such rates can even be approached under conditions, where many protein molecules bind DNA nonspecifically is all the more impressive. Such rapid target-binding activity is explained by fast 1D sliding along the DNA in many transcription factor and DNA repair protein systems. But high-speed transport is not a sufficient explanation since these proteins must probe the DNA to identify their targets. For some proteins, target-containing DNA is only subtly different from non-specific DNA, making the search all the more difficult. This is particularly the case for human oxoguanine DNA glycosylase 1 (hOGG1) repair protein (20–22), the subject of the present study. hOGG1 scans DNA to identify oxidized guanine bases and thus must diffuse fast along nonspecific DNA and also recognize and bind stably to oxoguanine bases, which may require local conformational change.\nThe apparently dichotomous requirements of rapid transport on nonspecific DNA and effective probing/stable target binding through DNA shape-sensing and/or base-specific interactions (which almost certainly increases free energy barriers for translocation) is known as the ‘speed-stability paradox’. Recently, it was shown that while the Lac repressor recognizes its cognate sites with low probability in a single pass, the redundancy of 1D diffusional search leads to a high probability of site recognition in the course of a single binding event (23).\nOne proposed solution to the speed-stability paradox is the ability of the searching protein to adopt multiple states: a fast-sliding 1D ‘search’ state that is minimally affected by the presence of targets interconverting with a slower-sliding ‘recognition’ state that recognizes targets with significant probability and has sacrificed speed-of-search in favor of the increased DNA interaction necessary for recognition (24).\nIt is a major technical challenge to determine the dynamics of search and recognition states quantitatively, as they are predicted to interconvert rapidly and stochastically to search DNA effectively. Consequently, the two states are not resolved, but averaged over, in standard ensemble biochemical assays.\nSingle-molecule-tracking assays promise to resolve individual biomolecule dynamics on DNA. Trajectories are mostly short, however, and position data are contaminated by multiple sources of experimental error. In particular: localization errors due to diffraction in microscope optics and limited numbers of recorded photons; motion blur due to particle movement during the camera exposure time; and thermal fluctuations of the DNA on which the particle diffuses. We need to account for these errors or risk severely biased results (25). The analysis is further complicated by the fact that we cannot simply average over multiple trajectories to reduce statistical error since we are interested in resolving individual molecular dynamics (26), and we can neither see nor model the moving substrate: It is invisible, and no solvable theory exists for its motion. Thus we face a sixtuple of challenges that reinforce each other: (i) resolve individual particle dynamics (ii) from mostly short trajectories (iii) recorded with considerable localization errors and motion blur (iv) on a moving substrate (v) that is invisible and (vi) its motion uncharted by theory.\nHere, we show how optimal estimators that treat noise sources in the single-molecule tracking data rigorously enable identification and quantitative characterization of two distinct states of searching of hOGG1 molecules along non-specific DNA. We analyze the diffusion of individual hOGG1 proteins that were tracked in vitro while diffusing on λ DNA stretched in a shear flow over a cover slip (Figure 1A–E).\nFigure 1. Experimental measurements of diffusion coefficients of hOGG1 proteins on flow-stretched λ DNA. (A) Workflow for estimation of diffusion coefficients from experimental tracking data of diffusing particles on DNA. (B) Experimental setup (not to scale): Proteins (green) on flow-stretched DNA (black) attached to cover slip at one end and fluctuating in a shear flow. (C) Image of several fluorescent hOGG1 molecules bound to and diffusing on a DNA molecule at higher density than during data acquisition for illustration (27). Scale-bar = 1 μm. (D) Transverse coordinates y of the trajectories of three hOGG1 proteins diffusing on DNA and recorded with time-lapse Δt = 11 ms. The mean residence time of proteins on DNA ranges from 50 to 500 ms, depending on solution conditions. (E) Longitudinal coordinates x of the same three proteins. (F) Periodogram of the transverse coordinate of a protein diffusing on DNA. The transverse fluctuations fit a Lorentzian plus a constant. The corner frequency (3 dB frequency) fc of the Lorentzian is 6.7 Hz. (G) Periodogram of longitudinal displacements of the same protein. The periodogram of displacements Δx is used, because diffusion is an unbounded process, making the periodogram of x a bad statistics. The expected value of , the power spectrum PΔx, is the sum of a diffusion term, a white-noise term, and a single Lorentzian term describing longitudinal DNA fluctuations. The corner frequency of these fluctuations is fc, x = 2fc, in agreement with our assumption that DNA fluctuations in the y- and z-directions contribute equally to DNA fluctuations in the x-direction. Shown values for data in F, G are block averages, each over 20 periodogram values, and the grey areas mark the 68% confidence interval (CI) for the block averages. We characterize the motion of the invisible substrate using data for the motion of the protein. That motion is recorded in two dimensions, while the substrate, the DNA, is one-dimensional and fairly stretched. Using this, we overcome the absence of a specific theory with a generic yet precise phenomenology that all realistic theories must share.\nUsing this phenomenology, we have constructed optimal estimators of diffusion coefficients of particles diffusing on a fluctuating substrate, e.g. DNA (25). Applying these estimators here, we find diffusion constants that, together with the distribution of residence times of hOGG1 on DNA, point unambiguously to a two-state kinetics of hOGG1 on DNA.\nWe propose an analytically solvable model for the kinetics of hOGG1: a protein binds to DNA in a loosely bound state and switches stochastically between this state and another more tightly bound state, until it detaches again from the DNA, from the loosely bound state. We provide accurate estimates of the kinetic parameters of our model, showing that the loosely bound state has much higher diffusivity than the tightly bound state. We hypothesize that the loosely bound and highly mobile state allows the protein to efficiently cover the length of the DNA, i.e. acting as a 1D travel and ‘search’ state, while the tightly bound and less mobile state allows it to probe for and recognize oxidative damage to DNA bases.\nThe experiment’s pH and salt concentration have been shown to affect the diffusivity and binding times of hOGG1 on DNA (27,28). Our estimates of the kinetic parameters of our model show that pH and salt concentration do not significantly influence the proteins’ transition rates between the two states and thus do not influence the time it spends in each state within the range of conditions tested here. Instead, pH and salt concentration affect hOGG1’s diffusion coefficient in the loosely bound state and its rate of detachment from the DNA from this state.\nSome of the data analyzed here were previously analyzed using a standard method, a straight-line fit to an ensemble of mean squared displacements (MSDs) of trajectories (27). This previous analysis assumed that substrate motion did not contribute to observed longitudinal displacements and that all such displacements were described by a single diffusion coefficient. The present analysis thus demonstrates the usefulness of optimized statistical methods, such as the CVE and MLE (25). They make the analysis of individual trajectories possible and that is crucial to fully exploit the single-molecule resolution offered by single-molecule experiments.\n\nMATERIALS AND METHODS\n\nExperimental data\nThe experimental setup and data acquisition are described in detail in (27). It produced five different data sets consisting of time-series of (x, y)-coordinates of thousands of fluorescently marked hOGG1 proteins that diffused on single flow-stretched λ-phage DNA molecules (48.5 kb = 16 μm long) at a range of different pH-values, 6.6–7.8, and salt concentrations, 0.01–0.1 M. Each data-set consisted of hundreds to thousands of time-series recorded at the experimental conditions listed in Table 1. Three of the data-sets (corresponding to DNA molecules 2–5) were previously analyzed in (27), while two datasets (DNA molecules 1 and 6) have not been analyzed previously.\nTable 1. Goodness-of-fit of the one-state and the two-state models. Column 1: Identity of DNA molecule(s). Column 2: pH-value of buffer. Column 3: Salt concentration of buffer. Column 4: Time lapse and shutter time of movie; the shutter time was set equal to the time lapse to maximize the number of recorded photons (45). Column 5: Number of time-series. Column 6: P-value (and number of degrees of freedom (dof) fitted) for the one-state model’s exponential distribution to the distribution of measured residence times and its predicted constant average diffusion coefficient as function of residence time to measured diffusion coefficients. Column 7: Akaike weight wAIC for the one-state model. Column 8: P-value (and number of degrees of freedom (dof) fitted) for the two-state model’s double-exponential distribution to measured residence times and its predicted average diffusion coefficient as function of residence time to experimental data. Column 9: Akaike weight wAIC for the two-state model; a more complex model that allows detachment from the tightly bound state shows no improvement of the fits, and is thus rejected by the Akaike information criterion (Supplementary Table S3). A low P-value means that the model can be rejected, while a high P-value means that the data supports the model. P ≥ 0.05 means the model is supported by any meaningful significance level. Akaike weights measure the relative probability of the models given the experimental data (47)\n1 2 3 4 5 6 7 8 9\nData set One-state model Two-state model\nDNA no. pH [NaCl] Δt n P (dof) w AIC P (dof) w AIC\n1 6.6 100 mM 25 ms 228 0.11 (11) 0.002 0.64 (8) 1.00\n2 \u0026 3 7.0 10 mM 28.5 ms 254 0.004 (12) 0.10 0.49 (9) 0.90\n4 7.0 75 mM 7 ms 55 3 × 10−6 (10) 2 × 10−8 0.15 (7) 1.00\n5 7.5 10 mM 11 ms 246 6 × 10−10 (13) 2 × 10−12 0.24 (10) 1.00\n6 7.8 50 mM 15 ms 92 8 × 10−8 (12) 6 × 10−12 0.61 (9) 1.00\nAll − − − 875 \u003c10−16 (58) \u003c10−16 0.38 (43) 1.00 In each experiment, a DNA molecule (two in one experiment) was end-biotinylated and fused to a coverslip (Figure 1B). The DNA was stretched to ∼75% of its contour length by a shear flow (the flow speed at the DNA was 100 μm s−1). If the free end of the DNA, where it tends to curl up, is excluded, the remaining DNA is stretched 90% (27). hOGG1 proteins fluorescently marked with Cy3B diffused on the DNA molecule and were filmed using total internal reflection microscopy (TIRF) by an EMCCD camera with pixel width corresponding to 250 nm. The proteins were tracked in the resulting movie until they detached from the DNA.\nTrajectories of protein positions were estimated by fitting 2D Gaussians plus constant backgrounds to the recorded point-spread functions as described in (27). Trajectories were previously analyzed by fitting a straight line to the average over the MSDs of all trajectories in each data-set (27).\nHere, we estimated diffusion coefficients, the variance of localization errors, and parameters characterizing the DNA’s motion statistically from these same trajectories using the optimal MLE and CVE methods developed in (25), as described below and in Supplementary Section S3.\nIn Supplementary Section S3F we compare the CVE and MLE to MSD-based methods and to recent Bayesian methods for inferring multi-state diffusion.\n\nA priori knowledge used to separate the protein’s diffusion along DNA from DNA motion\nFor each tracked protein, we know its trajectory measured in lab coordinates, apart from localization errors and motion blur. It is its trajectory on the DNA, however, that relates to its diffusion coefficient on the DNA. The difference between the two trajectories is given by the motion of the invisible DNA in the lab. This may seem an impasse because the DNA’s motion was not observed and we have no complete theory for its motion.\nWe do, however, have knowledge of general mathematical properties that any linear model of the DNA’s motion must satisfy. This leads to a phenomenological model for the motion of local segments of the DNA, and that is sufficient, since any visiting protein visits only a local segment of the DNA. This model dictates a protocol for how to analyze trajectory data in order to separate the motion of a protein along DNA from the motion of the DNA in the lab and thus obtain accurate estimates of protein diffusion coefficients. We sketch this phenomenological model in the present subsection and give the protocol in the next subsection. Details are given in (25).\n\nWorm-like Chain model in free-draining flow\nThe ultimate way to do describe the DNA’s motion, is to solve a realistic dynamical model for it. This is impossible, however: The simplest realistic model would treat the DNA as a massless, semi-flexible, unstretchable fiber (the Kratky-Porod model, a.k.a. the Worm-Like Chain model) in a free-draining shear flow over the coverslip to which one end of the DNA is attached. The result is a non-linear partial differential equation for the over-damped thermal Brownian motion of the fiber in the flow. The boundary condition on that motion imposed by the coverslip adds to the intractability of this model by exact analytical means.\n\nSpectral theory invoked\nThe non-linear nature of this simplest realistic model means that general results from mathematics’ spectral theory do not apply directly to it. They may apply indirectly, however, e.g., if our theory is well approximated by a linearized version of it. We expect this to be the case since the DNA is stretched taut and its fluctuations consequently are relatively small.\nSo, despite the model’s analytical intractability, qualitative mathematical results exist that may well apply to it. These results are of such a nature that experimental data will reveal whether they apply or not. Thus we need not argue their case theoretically. Right or wrong, our explanation does not affect the correctness of our ensuing data analysis, because that analysis relies solely on the phenomenological spectral theory described here and dictated by experimental data.\nIn any linear theory, the fiber’s motion can be written as a linear superposition of spatial eigen-modes. The amplitudes multiplying each spatial eigen-mode are time-dependent; each has its own thermal dynamics equal to that of a massively overdamped harmonic oscillator at finite temperature, and each mode’s amplitude has its own characteristic relaxation time.\nIt is well known from spectral theory that the spectrum of eigen-values, here the relaxation rates of eigen-modes, is bounded from below and discrete when the system is compact in the mathematical sense of that word. DNA of finite extent is compact. This means that there is a slowest relaxation time. Since DNA is also one-dimensional, we expect the discrete spectral values to be well separated; i.e. we expect the next-slowest relaxation time to be many times faster than the slowest one. Thus, with the limited time-resolution of time-lapse recordings, we may only resolve the motion of the slowest or a few of the slowest spatial eigen-modes even if we could see and track the motion of the DNA. All higher modes show up in our measurements as uncorrelated (white) noise, and simply add to the localization error that we observe.\nThus the task of describing the DNA’s invisible motion is simplified by realizing the manner in which this motion affects our data. Spectral theory facilitates this realization and formulates its result.\nFor illustration of spectral theory, we have analyzed a linear mean-field model for the DNA’s motion that is analytically tractable and discuss two other linear, exactly solvable models for DNA as well (Supplementary Section S1). Their solutions show that the DNA’s motion is composed of a sum of independent, orthogonal modes, each having its own characteristic relaxation time and its own characteristic mean amplitude. They also confirm that contributions from higher modes quickly become negligible: our mean-field model for DNA stretched by a shear flow shows that both the relaxation time (Supplementary Figure S1I and Supplementary Eq. (S31)) and the mean amplitude (Supplementary Figure S1J and Supplementary Eq. (S5)) of the second mode are six times lower than those of the lowest mode; our two other models, DNA stretched by a plug flow and by pulling at its ends, have second modes that have relaxation times and amplitudes that are 5.3 and 4 times lower than the respective values for their respective lowest modes (Supplementary Figures S1A, B, E and F).\nReturning to our experiment, the body of a priori knowledge we have invoked here has reduced our task to a point where we can deduce the DNA’s motion from the information we have at hand, when supplemented with yet another piece of a priori knowledge: two different length-scales occur in the experiment. The large ratio between them invites a highly simplifying approximation.\n\nSeparation of length scales\nA hOGG1 protein diffuses only for a while (its residence time) on a flow-stretched DNA molecule (Figure 1B–E). Thus, its diffusion length while in residence is small (of the order of 1 μm or smaller) compared to the length of the DNA (16 μm). The DNA is stretched out in the x-direction and its dominant fluctuations are large-wave-length according to spectral theory. Consequently, the protein’s diffusion on the DNA changes the protein’s y-coordinate insignificantly compared to changes in the protein’s y-coordinate due to the DNA’s motion in the lab. Thus, the y-coordinate of a protein’s trajectory is, effectively, the trajectory of a fixed, physical point on the DNA, e.g., the mean position, , of the protein, where s is distance along the DNA measured from its tethered end (e.g. in number of base-pairs). The y-coordinate’s periodogram (Figure 1F) bears this out, as discussed in the Results section below.\n\nTethered and unstretchable, the DNA’s unobservable longitudinal motion follows from its observable transversal motion\nSince the DNA is unstretchable, its longitudinal (x-direction) fluctuations are completely determined by its upstream fluctuations in the transverse y- and z- directions through the relation , where the prime on a variable denotes the derivative with respect to s. Since the DNA is taut, this expression may be approximated by x′ ≈ 1 − (y′)2/2 − (z′)2/2. This shows that longitudinal (x-direction) DNA fluctuations are determined by transverse fluctuations through a quadratic dependence. This implies that the power spectrum of is a Lorentzian with a corner frequency that is twice that of the power spectrum of , fc, x = 2fc. Here we have assumed that the statistics of DNA motion is the same in the two transverse directions, y and z (Supplementary Sections S1 and S2).\nWe cannot observe , which is why we make this assumption, so we cannot verify the assumption directly. It is a reasonable assumption if the presence of the cover slip mainly affects the average value of with little effect on fluctuations in about this shifted average. Our mean-field model for the DNA’s motion indicates that this is the case (Supplementary Section S1C), and data are consistent with this assumption (Figure 1G). Our protocol for data analysis exploits this to increase the precision of our statistical description of longitudinal DNA fluctuations based on the observed transverse motion.\n\nStatistical independence of the three contributions to changes in the protein’s x-coordinate\nFinally, x-displacements of the protein, Δxn, as measured in lab coordinates, are a sum of three statistically independent contributions: one from the protein’s actual diffusion along the DNA, another from localization errors, and a third from motion of the DNA in the x-direction. Localization errors are by nature statistically independent of the two actual, physical displacements. The latter two are independent of each other because the separation of length scales makes the observed diffusion of the protein on the DNA independent of the DNA’s motion and vice versa.\nConsequently, the covariance of the measured displacements of a protein’s x-coordinate is the sum of three contributions, one from each of those three statistically independent terms: the first, proportional to the protein’s diffusion coefficient D, which we want to determine, another parameterized by the variance of localization errors, σ2, which we will have to determine simultaneously with D if it is not known otherwise, and a third, the autocovariance of DNA fluctuations, , which we know, up to a prefactor, from substrate motion in the y-direction. Supplementary Equations (S39)–(S41) express this mathematically.\nSupplementary Equations (S39)–(S40) show that the displacements of a protein’s experimentally measured x-coordinate directly give estimates for D and σ2 as solutions to two coupled linear equations: The necessary input is the covariance of the measured displacements and the covariances , which also are determined from experimental data.\n\nProtocol for data analysis\nIn actual practice, we determine the values of D, σ2, and the parameters of simultaneously, using Maximum Likelihood Estimation (MLE) applied simultaneously to the periodogram of the measured displacements along the x-coordinate, and the periodogram of the y-coordinate of the recorded position (Supplementary Eqs. (S44) and (S45)).\nThis section walks through our protocol, motivating each step, while pointing to relevant sections of the Supplementary Information for details. Points 1–5 below describe our procedure for analyzing individual trajectories. Points 6–9 describe our procedure for investigating the results of this analysis and for comparing it to theoretical models for hOGG1’s dynamics on DNA. We discarded time-series that displayed anomalies or were too short to be tested reliably for anomalies (Supplementary Section S3A).\nMotivation: To ensure that the time-series we analyzed correspond to mobile DNA-bound proteins that experienced constant conditions along the DNA.\nFor each time-series longer than N = 25, we used MLE to estimate the protein’s diffusion coefficient D, the variance σ2 of localization errors, and the parameters characterizing DNA fluctuations locally in 1 μm long patches along the DNA with the phenomenological model of substrate motion developed in (25) and sketched above (Supplementary Section S3C).\nMotivation: For long time-series, MLE is unbiased and hence can reliably estimate D and substrate motion simultaneously. So we do that for the relatively few long time-series we have in order to characterize the motion of the DNA locally, where each long time-series visited. We have enough long time-series to characterize the local motion of the entire stretch of DNA visited by the short time-series we analyse below. With the motion of the DNA thus described, we can correct the CVE for DNA motion below.\nFor time-series shorter than N = 25, we used the CVE corrected for substrate motion, using weighted averages of the parameters that had been estimated using MLE for the several long time-series in the same patch (Supplementary Section S3D).\nMotivation: For time-series shorter than N = 25, the MLE may be significantly biased (Supplementary Figure S5). Thus, for such short time-series, recorded in a patch on the substrate that already had been characterized with the help of several long time-series, we used the CVE corrected for substrate motion.\nFor time-series longer than N = 25 in the same patch, we also estimated diffusion coefficients using the uncorrected CVE.\nMotivation: We compared these estimates to those obtained with MLE; their difference is an experimental estimate of the bias bD of the CVE, which we compare to our theoretical prediction for this bias given by Equation (1) below. This step provides a consistency check of our procedure (Figure 2F and Supplementary Section S3E).\nDiffusion coefficients for diffusion along the DNA were obtained by multiplying all estimates with the tortuosity of the DNA, which was constant, equal to 1.2, across the segment of the flow-stretched DNA that was used in our analysis (Supplementary Figure S1J and Supplementary Sections S2 and S3).\nMotivation: Estimated diffusion coefficients describe the projection on the x-axis of diffusion along the DNA. Since the DNA is stretched to 90% along the x-axis, diffusion coefficients for the proteins’ motion along the DNA contour are 1/(0.9)2 ≈ 1.2 times higher than those we measure for the motion projected along the x-axis.\nNote that although this tortuosity factor is a source of error on estimated values for the diffusion coefficients, it does not affect their ratios. So it does not affect any of the conclusions that we draw from our analysis, such as the statistical support for proposed models. Nor does it affect the relative importance of bias due to substrate fluctuations. Optimally, the tortuosity could be measured directly, e.g. by using optical tweezers to stretch the DNA.\nWe plotted the estimated diffusion coefficients of proteins on DNA against the measured residence time of proteins, as scatter plots (Figure 3A–E) and as binned on the residence-time axis (Figure 3F–J). Histograms of measured residence times were also plotted to investigate the distribution of residence times (Figure 3K–O).\nMotivation: These plots were used, with standard statistical testing as described below, to compare experimental data with theoretical models for the underlying kinetics (Supplementary Section S4).\nMaximum likelihood estimates of parameters in the two-state model presented in Figure 4B–G (Supplementary Section S4C) were obtained by a combined fit to the block-averaged estimates of diffusion coefficients (Figure 3F–J) and to the full distribution of residence times (Figure 3K–O).\nMotivation: The two-state model predicts correlated values for a protein’s residence time and the protein’s residence-time-averaged diffusion coefficient. Thus, combined simultaneous fitting of model predictions to experimental residence times as well as experimental diffusion coefficients (which are time-averages, each an average over one residency) is the logical and statistically optimal way to fit to these data.\nWe also fitted the two-state model directly to diffusion coefficient estimates obtained using the uncorrected CVE (Supplementary Figure S2 and Table S4).\nMotivation: This is another consistency test of our results: it leads to the same two-state model for hOGG1’s kinetics, albeit with biased estimates of the diffusivities (Supplementary Table S4). It confirms that the observed correlation between diffusion coefficients and residence times is not a spurious effect induced by a possible bias in the MLE and the corrected CVE due to finite statistics.\nThe next four sections, two on tests and two on Monte Carlo simulations, are peripheral parts of our protocol.\nMotivation: One can reproduce our results without the next four sections, but cannot reproduce our confidence in our results without them. The confidence rests on our understanding of the steps described above as well as on the reflections over self-consistency, precision, and accuracy that the next four sections describe.\nFigure 2. Estimated parameter values as functions of the protein’s position on DNA. Results for molecule no. 5 (pH 7.5 and [NaCl] = 0.01 M). Averages are weighted means and error bars are weighted s.e.m., except for estimates of diffusion coefficients; their averages are simple mean values. Data close to the DNA ends (open gray symbols) are excluded in the following analysis to avoid bias (see discussion in Supplementary Section S3A). (A) MLE of the corner frequency fc of transverse (y-direction) DNA motion. Measured values vary by a factor two between the tethered and the free end. (The hypothesis that fc is constant is refuted with P = 6 × 10−7.) The faster dynamics near the tethered end is due to larger tension in the DNA there. (The corner frequency of the DNA’s longitudinal motion is twice the corner frequency of its transverse motion.) (B) The diffusivity of longitudinal DNA fluctuations is constant along the DNA, corresponding to increasing amplitude in the downstream direction in consequence of decreasing tension; (P = 0.93). (C) The diffusivity of transverse DNA fluctuations is also constant along the DNA; (P = 0.19). (D) The CVE overestimates σ2, the variance of localization errors, by almost a factor two because it does not account for DNA motion. The MLE of σ2 increases slightly towards the free end. The assumption that it is constant has negligible support: P = 0.003. (E) Estimated diffusion coefficients. The uncorrected CVE overestimates diffusion coefficients significantly, because it does not account for DNA motion. The MLE does not, and shows that diffusion coefficients do not depend on the protein’s position on the DNA (P = 0.10). (F) Experimental estimates of the bias bD of the CVE, calculated as , and theoretical estimates, calculated from Equation (1) using weighted means of estimates of fc and Kx (Materials and Methods and Supplementary Section S3D). The theory generally agrees excellently with experiments. The bias increases near the DNA’s free end, where DNA fluctuations are larger and slower.\nFigure 3. Estimated diffusion coefficients and residence times for proteins on DNA. Each row shows results for the experimental conditions listed in the top of the left panels. (A–E) Estimated diffusion coefficients versus residence times for proteins on DNA. (F–J) Block averages of diffusion coefficient estimates in A–E binned on the time axis. (A, F, K) DNA molecule no. 1 (pH 6.6 and [NaCl] = 0.1 M); (B, G, L) DNA molecules nos. 2 and 3 (pH 7.0 and [NaCl] = 0.01 M); (C, H, M) DNA molecule no. 4 (pH 7.0 and [NaCl] = 0.075 M); (D, I, N) DNA molecule no. 5 (pH 7.5 and [NaCl] = 0.01 M); (E, J, O) DNA molecule no. 6 (pH 7.8 and [NaCl] = 0.05 M). A clear dependence on the residence time is seen in the measured diffusion coefficients, quite contrary to what one finds for a simple diffusion process. (K–O) Distribution of protein residence times on DNA. The distributions are not simple exponentials, so the rate of detachment of protein from DNA is not constant but decreases with the time bound. Dashed blue lines and full black lines in F–O mark combined ML fits of the one- and two-state models, respectively, to data in the second and third columns (P-values are given in legends and in Table 1, and estimated parameter values are given in Table 2).\nFigure 4. Two-state model for hOGG1’s diffusion on DNA. (A) Schematic free energy landscape for the two-state model for diffusion on DNA. Maximum likelihood fits of this two-state model to data shown in Figure 3 give P-values which support the two-state model and Akaike weights which overwhelmingly favor the two-state model over the one-state model (Table 1). (B–G) Maximum likelihood estimates of parameters in the model as a function of pH-value (B, D, F) and salt concentration (C, E, G). (B, C) Diffusion coefficients in the loosely bound state, D1, and in the tightly bound state, D2. *The low value of D1 found for [NaCl] = 0.1 M is explained by the low value of the pH (pH 6.6) for this DNA molecule. (D, E) Rates of transitions, b from the loosely to the tightly bound state and r for returning. These rates for changing between states do not show significant dependence on pH (D) or salt concentration (E). (F, G) Detachment rate d1 from the loosely bound state. **Salt concentration is varied 8-fold between the two measurements at pH 7.0, which explains the difference in detachment rates.\n\nTesting for drift of DNA-attached proteins with surrounding flow\nThe estimators for the diffusion coefficient were derived under the assumption of negligible drift. To test for drift of the proteins along the DNA in response to the drag-force from the surrounding flow, we calculated their mean displacement per time-lapse (Supplementary Section S3A). There was no discernible drift, neither locally on the DNA (Supplementary Figure S4D) nor on average (Supplementary Table S5).\n\nStatistical testing\nTheoretical models of the distribution of residence times were tested against experimental data using Pearson’s chi-square goodness-of-fit test with each bin containing at least five observed counts, and the number of bins n ≥ 7. For other tests we used a chi-squared test for variance to test whether averages of estimated parameters agreed with the expected values. Estimates were divided into m bins, and averages were calculated for each bin. The number of estimates in each bin was on average eleven or more to assure that averages were approximately Gaussian distributed. The one- and two-state models for hOGG1 on DNA predict both a distribution of protein residence times and a mean diffusion coefficient for a given residence time. They were fitted to both types of data in a combined fit. Consequently, the goodness-of-fit of each model was determined from a corresponding total chi-squared value: the sum of a Pearson’s chi-squared value for its agreement with the distribution of residence times, and a variance-based chi-squared value for its agreement with measured diffusion coefficients as function of residence time. The total degrees of freedom is n + m − p − 1 for each model, where p is the number of parameters fitted.\n\nMonte Carlo simulations of diffusion on a crowded DNA strand\nIn (27), experimental conditions were chosen that result in a low average density of proteins on the DNA (fewer than three hOGG1 molecules at a time). However, even if the TIRF setup limits bleaching of proteins before they bind to DNA, bleached or unlabeled proteins could also bind to DNA, so the number of proteins bound to DNA might have been higher than what we observed. Thus a protein may occasionally have prevented another protein from diffusing in a given direction. We therefore Monte Carlo simulated diffusing proteins at concentrations up to twenty times larger than in experiments. At these concentrations, estimated diffusion coefficients differ only negligibly from their values at zero concentration and can account neither for the observed dispersion in diffusion coefficients nor the correlation between residence times and diffusion coefficients that we see in experiments (Supplementary Table S6 and Supplementary Section S4A). It is easy to see why: Protein encounters are rare, and the rest of the time diffusion goes on as at zero concentration. This is in agreement with a recent study of the effects of crowding on single-file diffusion (29). Those authors introduced the order parameter τp = c2DΔt, where c is the concentration of diffusers on the DNA. They concluded that for τp ≪ 1, the effect of crowding is negligible; for our experimental conditions we find τp ∼ 10−5.\n\nMonte Carlo simulations of random walkers in a quenched potential landscape\nDNA is not a perfectly homogeneous medium since AT and CG base-pairs differ physically. To investigate whether this could cause our observed nontrivial distributions of residence times and diffusion coefficients, we Monte Carlo simulated the model proposed in (24), which assumes that proteins perform random walks in a quenched energy-landscape with uncorrelated and Gaussian distributed binding energies (Supplementary Figure S6A and Supplementary Section S4B). The time-lapse of our experimental measurements is much longer than the molecular time-scale set by the mean time between consecutive steps in the random walk of a protein. Consequently, a protein diffusing on DNA effectively averages over the DNA’s quenched energy-landscape during a single time-lapse. Thus, the Monte Carlo simulation shows no correlations between diffusion coefficients and residence times (Supplementary Figure S6C), and its distribution of residence times is a simple exponential (Supplementary Figure S6B), which means that detachment is a simple Poisson process with a constant rate of detachment. In summary, the model proposed in (24) predicts that the simple one-state model should explain observations. Since it does not, we have eliminated the model proposed in (24) as an explanation.\n\nRESULTS\n\nCrucial experimental result for the DNA’s motion\nThe y-coordinate’s periodogram (30) (Figure 1F) fits a Lorentzian plus a constant (25). (This Lorentzian is aliased and low-pass filtered due to limited time-resolution and motion blur, respectively (25), Supplementary Eq. (S44).) This fit means that the y-motion of the tracked point on the DNA is that of a massively over-damped harmonic oscillator at finite temperature (31,32) tracked with a white-noise localization error. The only physical explanation of this observation is that locally on the DNA, at the segment whose y-coordinate was tracked, only a single mode contributes discernibly to the DNA’s recorded motion, in agreement with spectral theory (Materials and Methods).\nThis is an extremely useful result. The DNA has very many degrees of freedom, but now we found phenomenologically that we can describe the DNA’s motion locally with a single degree of freedom which, to boot, displays the simplest possible dynamics.\nThe mode with the slowest relaxation time is the dominant contributor to the motion of the DNA at most positions along the DNA. The other modes contribute, but with corner frequencies (= 3 dB frequencies) that increase and amplitudes that decrease so rapidly with mode-number that already the second mode is indistinguishable from white noise in our data: The Nyquist frequency fNyq = 1/(2Δt), where Δt is the time-lapse between measurements, is the maximal frequency that we can resolve in our time-lapse recorded data. It is 45 Hz in Figure 1F. Our mean field approximation makes the relaxation time of the second mode six times shorter than the relaxation time of the first mode (Supplementary Figure S1I). In Figure 1F, that factor six places the corner frequency of the second mode at 40 Hz, which is so close to the Nyquist frequency that our limited statistics cannot distinguish it from white noise, and hence also cannot distinguish higher modes. This, we believe, is why we see only the Lorentzian of a single mode, the first mode, in Figure 1F. The corner frequency of this Lorentzian is in Figure 1F, which is much smaller than the Nyquist frequency.\nWe consequently use this lowest mode to characterize the DNA’s local transverse fluctuations. We know how the DNA’s longitudinal motion follows from the transverse motion that we now have characterized (Materials and Methods and Supplementary Section S2). Thus, we now can separate the DNA’s longitudinal fluctuations from the observed longitudinal motion of the protein, as described in the Materials and Methods section.\n\nCharacterizing the diffusion of hOGG1 on DNA\nThe key steps in our data analysis are outlined in Figure 1A and in the Materials and Methods section. The experimental setup and data samples are shown in Figure 1B–E. Power spectral analysis of the trajectories reveals that the correlation time τ of transverse DNA fluctuations is longer than the time-lapse between frames Δt (Figure 1F, G and Figure 2A)—compare τ = 1/(2πfc) = 25–50 ms to Δt = 11 ms. Diffusion coefficients (Figure 2E), localization error (Figure 2D), and DNA motion parameters (Figure 2A–C) are Maximum Likelihood estimated from time-series longer than N = 25 (Materials and Methods and Supplementary Section S3C). For time-series shorter than N = 25, estimates of parameters of DNA motion from long time-series are used to calculate and correct the bias of the covariance-based estimator (CVE) to obtain unbiased estimates of diffusion coefficients (Materials and Methods and Supplementary Section S3D).\nComparison between diffusion coefficients estimated with MLE and the uncorrected CVE shows that DNA motion causes a bias of up to 0.4 μm2 s−1 (Figure 2E–F, and Materials and Methods and Supplementary Table S1), depending on the protein’s position on the DNA and the time-lapse according to (1) where parametrizes the amplitude of longitudinal fluctuations of the point on the DNA (25). Equation (1) shows that the bias caused by substrate fluctuations depends linearly on Kx, with a proportionality coefficient that is at most equal to unity (with equality when the relaxation time of longitudinal substrate fluctuations, 1/(4πfc), is much smaller than Δt) and decreases with fcΔt as 1 − 6πfcΔt when fcΔt is small and as (4πfcΔt)−3 when fcΔt is large. (For the trajectory analyzed in Figure 1F, G, the estimated bias is , as compared to an estimated diffusion coefficient of .) We measure diffusion coefficients in the range 0.01–0.5 μm2 s−1 (Figure 3F–J and in Table 2), so the bias could, in the worst case, be many times larger than the actual value of the diffusion coefficient.\nTable 2. Estimated parameter values of the two-state model. Column 1: Identity of DNA molecule(s). Column 2: pH-value of buffer. Column 3: Salt concentration of buffer, [NaCl]. Column 4: Dissociation rate from the loosely bound state, d1. Column 5: Rate for switching from the loosely to the tightly bound state, b. Column 6: Rate for switching from the tightly to the loosely bound state, r. Column 7: Diffusion coefficient in the loosely bound state, D1. Column 8: Diffusion coefficient in the tightly bound state, D2. Column 9: Mean fraction of time proteins spend in the loosely bound state, f1. Column 10: Effective mean dissociation rate from DNA, . Column 11: Mean observed diffusivity . Measured values of for DNA molecules 2, 3, 4, and 5 are smaller than those reported for corresponding samples in (27). The exclusion of trajectories shorter than Nmin = 12 points that show higher average diffusivity, compared to Nmin = 5 in (27), as well as correction for substrate fluctuations, account for this difference.\n1 2 3 4 5 6 7 8 9 10 11\nDNA no. pH [NaCl] d 1 b r D 1 D 2 f 1\n(mM) (Hz) (Hz) (Hz) (μm2 s−1) (μm2 s−1) (Hz) (μm2 s−1)\n1 6.6 100 3 ± 3 1.3 ± 1.6 1.7 ± 0.5 0.02 ± 0.01 0.001 ± 0.004 0.47 ± 0.23 2 ± 2 0.009 ± 0.002\n2 \u0026 3 7.0 10 13 ± 6 0.6 ± 0.5 2.2 ± 0.2 0.23 ± 0.04 0.07 ± 0.02 0.73 ± 0.17 10 ± 5 0.11 ± 0.01\n4 7.0 75 33 ± 10 1.4 ± 0.6 3.3 ± 0.8 0.46 ± 0.11 0.01 ± 0.04 0.65 ± 0.11 23 ± 7 0.09 ± 0.04\n5 7.5 10 5 ± 1 2.1 ± 0.9 2.2 ± 0.4 0.44 ± 0.06 0.07 ± 0.04 0.42 ± 0.09 2 ± 1 0.22 ± 0.05\n6 7.8 50 15 ± 5 1.0 ± 0.4 2.3 ± 0.4 0.44 ± 0.07 0.05 ± 0.03 0.64 ± 0.09 11 ± 2 0.16 ± 0.05 The bias of the uncorrected CVE does not depend on the proteins’ residence time on DNA (Supplementary Figure S2 and (25)), and (unbiased) estimated diffusion coefficients do not depend on the proteins’ position on the DNA (Figure 2E). However, the proteins’ estimated diffusion coefficients show a much higher spread around their mean than what would be expected from pure statistical error, i.e. they are not consistent with the proteins having a single unique diffusion coefficient (Supplementary Table S2). The diffusion coefficients furthermore depend on the proteins’ residence time on DNA (Figure 3A–J). This result invites close inspection of the residence times themselves. They are appealingly simple data, as each residence time is just measured with a clock. Nevertheless, the distributions of the protein residence times are non-trivial: they do not fit the simple exponential distributions that would result from a fixed rate of detachment (Figure 3K–O). hOGG1 proteins on DNA must have more than one state of attachment according to this distribution of residence times.\nThe absence of drift in experiments (Materials and Methods and Supplementary Table S5) implies that the non-trivial distribution of hOGG1 residence times and their correlation with measured diffusion coefficients cannot be explained by intermittent excursions of hOGG1 into bulk, so-called ‘hopping’ (33) (Supplementary Section S4D). Monte Carlo simulations of diffusion on ‘crowded’ DNA indicate that the non-trivial dynamics of hOGG1 cannot be explained by crowding either (Materials and Methods and Supplementary Section S4A). Monte Carlo simulations of diffusion of a one-state protein on DNA with a random potential landscape indicate that such a landscape also cannot explain this non-trivial behavior (Materials and Methods and Supplementary Section S4B).\n\nTwo-state model for hOGG1’s kinetics when bound to DNA\nIn order to explain the non-trivial distribution of residence times and their correlation with measured diffusion coefficients, we propose a minimal two-state model for hOGG1 on DNA: proteins bind to DNA in a loosely bound state (State 1); they switch stochastically to a more tightly bound state (State 2) with rate b and return to the loosely bound state with rate r; proteins detach from the loosely bound state with rate d1; and proteins have diffusion coefficient D1 and D2 in the loosely and tightly bound state, respectively (Figure 4A). We derive analytical expressions for the distribution of protein residence times on DNA and for the mean value of a protein’s diffusion coefficient on DNA as a function of its residence time (Supplementary Section S4C). This allows us to directly fit parameters of the two-state model to experimental data using maximum likelihood estimation (estimated parameter values are listed in Table 2). The model shows excellent agreement with experimental data (Figure 3F–O and Table 1). Proteins diffuse much faster in the loosely bound state than in the tightly bound state, and diffusion coefficients in the loosely bound state depend highly on pH (Figure 4B), while diffusivity does not depend on salt concentration (Figure 4C). Taken with the observation that the diffusivity of the H270A hOGG1 mutant protein is pH-insensitive, this suggests that hOGG1’s His270 residue is important for diffusion in the loosely bound state (27). Proteins spend slightly more time loosely bound than tightly bound (Figure 4D, E and Table 2), and they make only few transitions between states during typical residence times on DNA (compare b and r in Figure 4D, E to d1 in Figure 4F, G).\n\nDISCUSSION\nAt the resolution available experimentally, the thermal motion of a stretched DNA molecule is accurately described locally as a single thermally driven overdamped harmonic oscillator. This simple theory works because higher order modes of the DNA’s motion are too fast and too small in amplitude to be resolved. Using this theory, we characterized the fluctuations of the DNA strands from the longer time-series of proteins diffusing on them. This allowed us to predict the bias caused by DNA fluctuations in covariance-based estimates of diffusion coefficients. With this bias removed, we estimated diffusion coefficients with accuracy and optimal precision for single hOGG1 repair proteins diffusing on λ DNA even from short time-series. Our increased precision revealed a two-state kinetics in hOGG1’s diffusion on DNA, which was missed by the cruder but common method based on the ensemble-averaged mean squared displacement, and allowed us to accurately estimate the proteins’ diffusivity in each state as well as the rates of switching between states.\nA two-state kinetics was earlier proposed to explain the motion of Proliferating Cell Nuclear Antigen (PCNA) on DNA (8). This study, however, relied on (i) PCNA undergoing rotation-coupled sliding in one state and rotation-uncoupled linear diffusion in another state, where the protein topologically entraps the substrate, (ii) elaborate experiments using labels of different sizes, and (iii) inference of the existence of the two states of motion from indirect evidence. Another recent study found evidence of two-state kinetics in TALE proteins’ diffusion along DNA (9). The study relied on subjective visual identification of periods of low diffusivity from the recorded time-series: a method which (i) can only be applied when switching between states is so slow that such periods are long enough to be identified visually in a noisy time-series and (ii) provides no rigorous statistical evidence for the inferred model. Importantly, the above approaches did not allow quantitative estimation of the diffusivity in the two states.\nUsing optimal estimation of individual diffusion coefficients, we have shown the existence of two-state kinetics in the diffusion of hOGG1 on DNA directly from single-molecule tracking data from a single experiment. We have derived analytical expressions for the distribution of the proteins’ residence times on DNA and their expected diffusion coefficient as function of their residence time, enabling rigorous statistical testing and model selection, even when state changes are too short to be detected visually, and we have provided accurate estimates of all of the kinetic parameters of the two-state model. A plausible physical mechanism for this two-state kinetics is that hOGG1 undergoes conformal change between a highly mobile state, where it is loosely bound to the DNA, and a more tightly bound, less mobile state that allows greater engagement with the DNA bases to detect oxidative damage. Our analysis showed that salt concentration in solution solely affected hOGG1’s detachment rate from its loosely bound state and not rates for switching between states for the range of conditions tested. This indicates that electrostatic screening due to free ions does not affect hOGG1 dynamics in the tightly bound state. Furthermore, pH affected hOGG1’s diffusivity in the loosely bound state only. This suggests that the pH-sensitive His270 residue, which was previously shown to play a role in hOGG1’s diffusion along DNA (27), is sensitive to solution pH and/or affects sliding in the loosely bound state only, while it is not involved in facilitating changes in the protein’s state of binding.\nHere, we demonstrated the ability to identify and characterize multiple dynamic states of searching protein molecules from single-molecule tracking data, which will enable refinement of theoretical models explaining how the speed-stability paradox is overcome by hOGG1. The speed-stability paradox is a general physical problem in protein-DNA interaction, and we predict that many proteins may use conformationally and dynamically distinct binding modes to solve the paradox. Resolving dynamic states of DNA-bound proteins from single-molecule diffusion data could be broadly useful as an approach for generating targets for structural studies of each bound state and functional studies characterizing the role of each state in search and/or target binding.\nOur results point to sliding as the dominant contributor to hOGG1’s 1D motion and search for target sites along DNA on the mesoscopic scale, in agreement with other single-molecule studies (7,27), but here shown to consist of two distinct sliding modes. Conversely, ensemble studies indicate much shorter mean sliding lengths, of the order of several to tens of base-pairs; much shorter binding times, on the micro to millisecond timescale (34–38); and a relatively greater contribution of hopping (34,35,37). Much of these differences can be attributed to differences between salt concentrations used in ensemble studies (typically high, comparable to intracellular levels) and single-molecule studies (typically 5- to 20-fold lower). The archetypical sensitivity of nucleic acid proteins’ sequence-nonspecific DNA binding to salt concentration results in a strong shift of the binding equilibrium toward the unbound state and an increased rate of unbinding from the DNA with rising salt concentrations.\nNumerical simulations of the Lac repressor’s target search have shown that if the protein follows a helical path along the DNA, e.g., the phosphate backbone or the major or minor groove, short-range hopping does not contribute significantly to target search in vivo due to geometrical constraints on rebinding (33). This picture agrees with our conclusion that hopping does not contribute to hOGG1’s motion along stretched DNA in vitro based on drift analysis ((27) and here), largely salt concentration independent diffusion constants ((27) and here), and strong rotation coupling (7), suggesting persistant sliding on DNA along a helical path for hOGG1. ‘Molecular clock’ ensemble measurements, however, indicate that hOGG1 can overcome obstacles on the DNA phosphate backbone and transfer between strands without dissociating to bulk (34). Groove-focused sliding in the loosely bound state may explain hOGG1’s rotation-coupled sliding together with an ability to transfer between strands and negotiate backbone-bound obstacles without fully escaping from DNA.\nDynamic NMR spectroscopy measurements have revealed that non-specifically bound human uracil DNA glycosylase (hUNG) undergoes conformational dynamics in the millisecond range (39). Others have suggested that the millisecond timescale of hUNG’s conformational dynamics are tuned to the timescale of DNA breathing (i.e. bubble formation and coalescence) (39,40). This is two orders of magnitude faster than the rates for hOGG1’s switching between the sliding states measured here. It would be interesting to perform spectroscopy measurements to see whether the rates of protein conformational changes of non-specifically bound hOGG1 correspond to the rates of switching between states observed here.\nWe finally conclude that more single-molecule tracking measurements should be done, with hOGG1 on DNA and with other diffusers. New experimental protocols ease large-scale single-molecule data collection (41), and optimized tools of analysis are now available and promise novel insights, as demonstrated here with old, suboptimal data. The increased resolution and much higher throughput will enable rigorous exploration of a larger range of pH and salt concentrations and may help integrate results from single-molecule and ensemble experiments.\nOur results indicate that one may increase the precision of estimates many times by using an experimental setup that applies constant tension along the entire substrate; ∼1 pN would be preferred. One configuration capable of this is a dual optical trap setup (42). Dual optical traps give precise control of the extension of the substrate and can measure DNA fluctuations directly via fluctuations of the two trapped beads holding the DNA. Moreover, modeling of the DNA’s motion is simple in this case; a linear model applies and is fully solvable. This has the very practical consequence that recordings of the motions of the two beads and transverse motion of the DNA measured at various points along the DNA all can be combined to determine the few parameters of the linear model that describes the DNA’s motion, as done in (43) for a similar, less simple setup. The model will thus be over-determined and hence can be falsified, and it enables direct and rigorous analysis of time-series that describe diffusion over length scales that are of the order of the DNA’s extent. Also, the local tortuosity (or, equivalently, the degree of stretching) of the DNA can be calculated from the fitted model (43), so absolute values of diffusion constants can be measured with accuracy and precision. With dual traps one can furthermore investigate rigorously how stretching of the substrate affects kinetic parameters, providing an experimental means to verify numerical predictions for the influence of DNA dynamics on sliding (44) and to link single-molecule measurements made on stretched DNA to ensemble measurements made on DNA in free configuration.\nFlow-stretch assays like the one studied here and in (41) trade the precision of optical tweezers for a simpler experimental setup and the possibility to record on multiple DNA molecules in parallel, providing much more data. Here we recommend to use a strong flow to stretch the DNA strands in order to drive their fluctuations as fast as possible relative to the time-lapse of recordings, while assuring that this does not lead to drift of the attached proteins. Fast fluctuations will ensure that the phenomenological single-mode theory for the DNA’s motion is valid and increases the precision of estimated diffusion coefficients (25). This is all the more important when one wants to choose the time-lapse as small as possible in order optimize the number of recorded positions (45). If the fluctuations are fast enough, they may even be absorbed in the localization error term and do not need to be explicitly accounted for, leading to a simple, and more precise, estimation procedure—this should be checked by plotting the periodograms of the recorded time-series.\nFinally, tracking with MINFLUX promises to increase the acquisition rate of recordings substantially without sacrificing localization precision (46). This would facilitate the separation of DNA fluctuations from diffusive motion and improve the precision of inferred model parameters, and would make it possible to probe faster dynamics of single proteins sliding on DNA.\n\nSupplementary Material\nSupplementary Data Click here for additional data file.\n\nSUPPLEMENTARY DATA\nSupplementary Data are available at NAR online.\n\nACKNOWLEDGEMENTS\nThe experimental data analyzed here were produced in the lab of Xiaoliang Sunney Xie in the Department of Chemistry and Chemical Biology at Harvard University.\n\nFUNDING\nNational Institutes of Health [ST32 GM07598-25 to P.C.B.] through the Harvard University’s Molecular, Cellular and Chemical Biology training program; Burroughs Wellcome Fund through a Career Award at the Scientific Interface (to P.C.B); Human Frontier Science Program Research [GP0054/2009-C to H.F.]. Funding for open access charge: Human Frontier Science Program; Burroughs Wellcome Fund.\nConflict of interest statement. None declared.","divisions":[{"label":"article-title","span":{"begin":0,"end":89}},{"label":"p","span":{"begin":14186,"end":14464}},{"label":"abstract","span":{"begin":100,"end":1366}},{"label":"title","span":{"begin":100,"end":108}},{"label":"p","span":{"begin":109,"end":1366}},{"label":"body","span":{"begin":1368,"end":60367}},{"label":"sec","span":{"begin":1368,"end":10551}},{"label":"title","span":{"begin":1368,"end":1380}},{"label":"p","span":{"begin":1381,"end":1804}},{"label":"p","span":{"begin":1805,"end":2135}},{"label":"p","span":{"begin":2136,"end":3265}},{"label":"p","span":{"begin":3266,"end":3818}},{"label":"p","span":{"begin":3819,"end":4244}},{"label":"p","span":{"begin":4245,"end":4552}},{"label":"p","span":{"begin":4553,"end":5714}},{"label":"p","span":{"begin":5715,"end":6124}},{"label":"figure","span":{"begin":6125,"end":7917}},{"label":"label","span":{"begin":6125,"end":6134}},{"label":"caption","span":{"begin":6136,"end":7917}},{"label":"p","span":{"begin":6136,"end":7917}},{"label":"p","span":{"begin":7918,"end":8265}},{"label":"p","span":{"begin":8266,"end":8616}},{"label":"p","span":{"begin":8617,"end":9337}},{"label":"p","span":{"begin":9338,"end":9898}},{"label":"p","span":{"begin":9899,"end":10551}},{"label":"sec","span":{"begin":10553,"end":39472}},{"label":"title","span":{"begin":10553,"end":10574}},{"label":"sec","span":{"begin":10576,"end":14609}},{"label":"title","span":{"begin":10576,"end":10593}},{"label":"p","span":{"begin":10594,"end":11267}},{"label":"table-wrap","span":{"begin":11268,"end":13264}},{"label":"label","span":{"begin":11268,"end":11276}},{"label":"caption","span":{"begin":11278,"end":12746}},{"label":"title","span":{"begin":11278,"end":12746}},{"label":"table","span":{"begin":12747,"end":13264}},{"label":"tr","span":{"begin":12747,"end":12772}},{"label":"th","span":{"begin":12747,"end":12748}},{"label":"th","span":{"begin":12750,"end":12751}},{"label":"th","span":{"begin":12753,"end":12754}},{"label":"th","span":{"begin":12756,"end":12757}},{"label":"th","span":{"begin":12759,"end":12760}},{"label":"th","span":{"begin":12762,"end":12763}},{"label":"th","span":{"begin":12765,"end":12766}},{"label":"th","span":{"begin":12768,"end":12769}},{"label":"th","span":{"begin":12771,"end":12772}},{"label":"tr","span":{"begin":12773,"end":12815}},{"label":"th","span":{"begin":12773,"end":12781}},{"label":"th","span":{"begin":12783,"end":12798}},{"label":"th","span":{"begin":12800,"end":12815}},{"label":"tr","span":{"begin":12816,"end":12879}},{"label":"th","span":{"begin":12816,"end":12823}},{"label":"th","span":{"begin":12825,"end":12827}},{"label":"th","span":{"begin":12829,"end":12835}},{"label":"th","span":{"begin":12837,"end":12839}},{"label":"th","span":{"begin":12841,"end":12843}},{"label":"th","span":{"begin":12845,"end":12853}},{"label":"th","span":{"begin":12855,"end":12861}},{"label":"th","span":{"begin":12863,"end":12871}},{"label":"th","span":{"begin":12873,"end":12879}},{"label":"tr","span":{"begin":12880,"end":12940}},{"label":"td","span":{"begin":12880,"end":12881}},{"label":"td","span":{"begin":12883,"end":12886}},{"label":"td","span":{"begin":12888,"end":12894}},{"label":"td","span":{"begin":12896,"end":12901}},{"label":"td","span":{"begin":12903,"end":12906}},{"label":"td","span":{"begin":12908,"end":12917}},{"label":"td","span":{"begin":12919,"end":12924}},{"label":"td","span":{"begin":12926,"end":12934}},{"label":"td","span":{"begin":12936,"end":12940}},{"label":"tr","span":{"begin":12941,"end":13006}},{"label":"td","span":{"begin":12941,"end":12946}},{"label":"td","span":{"begin":12948,"end":12951}},{"label":"td","span":{"begin":12953,"end":12958}},{"label":"td","span":{"begin":12960,"end":12967}},{"label":"td","span":{"begin":12969,"end":12972}},{"label":"td","span":{"begin":12974,"end":12984}},{"label":"td","span":{"begin":12986,"end":12990}},{"label":"td","span":{"begin":12992,"end":13000}},{"label":"td","span":{"begin":13002,"end":13006}},{"label":"tr","span":{"begin":13007,"end":13071}},{"label":"td","span":{"begin":13007,"end":13008}},{"label":"td","span":{"begin":13010,"end":13013}},{"label":"td","span":{"begin":13015,"end":13020}},{"label":"td","span":{"begin":13022,"end":13026}},{"label":"td","span":{"begin":13028,"end":13030}},{"label":"td","span":{"begin":13032,"end":13045}},{"label":"td","span":{"begin":13047,"end":13055}},{"label":"td","span":{"begin":13057,"end":13065}},{"label":"td","span":{"begin":13067,"end":13071}},{"label":"tr","span":{"begin":13072,"end":13141}},{"label":"td","span":{"begin":13072,"end":13073}},{"label":"td","span":{"begin":13075,"end":13078}},{"label":"td","span":{"begin":13080,"end":13085}},{"label":"td","span":{"begin":13087,"end":13092}},{"label":"td","span":{"begin":13094,"end":13097}},{"label":"td","span":{"begin":13099,"end":13113}},{"label":"td","span":{"begin":13115,"end":13124}},{"label":"td","span":{"begin":13126,"end":13135}},{"label":"td","span":{"begin":13137,"end":13141}},{"label":"tr","span":{"begin":13142,"end":13208}},{"label":"td","span":{"begin":13142,"end":13143}},{"label":"td","span":{"begin":13145,"end":13148}},{"label":"td","span":{"begin":13150,"end":13155}},{"label":"td","span":{"begin":13157,"end":13162}},{"label":"td","span":{"begin":13164,"end":13166}},{"label":"td","span":{"begin":13168,"end":13181}},{"label":"td","span":{"begin":13183,"end":13192}},{"label":"td","span":{"begin":13194,"end":13202}},{"label":"td","span":{"begin":13204,"end":13208}},{"label":"tr","span":{"begin":13209,"end":13264}},{"label":"td","span":{"begin":13209,"end":13212}},{"label":"td","span":{"begin":13214,"end":13215}},{"label":"td","span":{"begin":13217,"end":13218}},{"label":"td","span":{"begin":13220,"end":13221}},{"label":"td","span":{"begin":13223,"end":13226}},{"label":"td","span":{"begin":13228,"end":13239}},{"label":"td","span":{"begin":13241,"end":13247}},{"label":"td","span":{"begin":13249,"end":13258}},{"label":"td","span":{"begin":13260,"end":13264}},{"label":"p","span":{"begin":13265,"end":13888}},{"label":"p","span":{"begin":13889,"end":14185}},{"label":"p","span":{"begin":14465,"end":14609}},{"label":"sec","span":{"begin":14611,"end":24282}},{"label":"title","span":{"begin":14611,"end":14697}},{"label":"p","span":{"begin":14697,"end":15132}},{"label":"p","span":{"begin":15133,"end":15795}},{"label":"sec","span":{"begin":15797,"end":16473}},{"label":"title","span":{"begin":15797,"end":15840}},{"label":"p","span":{"begin":15841,"end":16473}},{"label":"sec","span":{"begin":16475,"end":20202}},{"label":"title","span":{"begin":16475,"end":16498}},{"label":"p","span":{"begin":16499,"end":16869}},{"label":"p","span":{"begin":16870,"end":17352}},{"label":"p","span":{"begin":17353,"end":17716}},{"label":"p","span":{"begin":17717,"end":18551}},{"label":"p","span":{"begin":18552,"end":18755}},{"label":"p","span":{"begin":18756,"end":19824}},{"label":"p","span":{"begin":19825,"end":20202}},{"label":"sec","span":{"begin":20204,"end":21149}},{"label":"title","span":{"begin":20204,"end":20231}},{"label":"p","span":{"begin":20232,"end":21149}},{"label":"sec","span":{"begin":21151,"end":22650}},{"label":"title","span":{"begin":21151,"end":21268}},{"label":"p","span":{"begin":21269,"end":22042}},{"label":"p","span":{"begin":22043,"end":22650}},{"label":"sec","span":{"begin":22652,"end":24282}},{"label":"title","span":{"begin":22652,"end":22744}},{"label":"p","span":{"begin":22745,"end":23317}},{"label":"p","span":{"begin":23318,"end":23942}},{"label":"p","span":{"begin":23943,"end":24282}},{"label":"sec","span":{"begin":24284,"end":34912}},{"label":"title","span":{"begin":24284,"end":24310}},{"label":"p","span":{"begin":24311,"end":24646}},{"label":"p","span":{"begin":24647,"end":30487}},{"label":"p","span":{"begin":25024,"end":25156}},{"label":"p","span":{"begin":25157,"end":25303}},{"label":"p","span":{"begin":25304,"end":25660}},{"label":"p","span":{"begin":25661,"end":26169}},{"label":"p","span":{"begin":26170,"end":26411}},{"label":"p","span":{"begin":26412,"end":26721}},{"label":"p","span":{"begin":26722,"end":26843}},{"label":"p","span":{"begin":26844,"end":27172}},{"label":"p","span":{"begin":27173,"end":27472}},{"label":"p","span":{"begin":27473,"end":27802}},{"label":"p","span":{"begin":27803,"end":28267}},{"label":"p","span":{"begin":28268,"end":28593}},{"label":"p","span":{"begin":28594,"end":28791}},{"label":"p","span":{"begin":28792,"end":29077}},{"label":"p","span":{"begin":29078,"end":29497}},{"label":"p","span":{"begin":29498,"end":29651}},{"label":"p","span":{"begin":29652,"end":30051}},{"label":"p","span":{"begin":30052,"end":30162}},{"label":"p","span":{"begin":30163,"end":30487}},{"label":"figure","span":{"begin":30488,"end":32546}},{"label":"label","span":{"begin":30488,"end":30497}},{"label":"caption","span":{"begin":30499,"end":32546}},{"label":"p","span":{"begin":30499,"end":32546}},{"label":"figure","span":{"begin":32547,"end":33808}},{"label":"label","span":{"begin":32547,"end":32556}},{"label":"caption","span":{"begin":32558,"end":33808}},{"label":"p","span":{"begin":32558,"end":33808}},{"label":"figure","span":{"begin":33809,"end":34912}},{"label":"label","span":{"begin":33809,"end":33818}},{"label":"caption","span":{"begin":33820,"end":34912}},{"label":"p","span":{"begin":33820,"end":34912}},{"label":"sec","span":{"begin":34914,"end":35392}},{"label":"title","span":{"begin":34914,"end":34978}},{"label":"p","span":{"begin":34979,"end":35392}},{"label":"sec","span":{"begin":35394,"end":36653}},{"label":"title","span":{"begin":35394,"end":35413}},{"label":"p","span":{"begin":35414,"end":36653}},{"label":"sec","span":{"begin":36655,"end":38116}},{"label":"title","span":{"begin":36655,"end":36715}},{"label":"p","span":{"begin":36716,"end":38116}},{"label":"sec","span":{"begin":38118,"end":39472}},{"label":"title","span":{"begin":38118,"end":38193}},{"label":"p","span":{"begin":38194,"end":39472}},{"label":"sec","span":{"begin":39474,"end":49411}},{"label":"title","span":{"begin":39474,"end":39481}},{"label":"sec","span":{"begin":39483,"end":42053}},{"label":"title","span":{"begin":39483,"end":39531}},{"label":"p","span":{"begin":39532,"end":40217}},{"label":"p","span":{"begin":40218,"end":40466}},{"label":"p","span":{"begin":40467,"end":41630}},{"label":"p","span":{"begin":41631,"end":42053}},{"label":"sec","span":{"begin":42055,"end":47642}},{"label":"title","span":{"begin":42055,"end":42099}},{"label":"p","span":{"begin":42100,"end":43025}},{"label":"p","span":{"begin":43026,"end":44111}},{"label":"label","span":{"begin":43316,"end":43319}},{"label":"table-wrap","span":{"begin":44112,"end":45794}},{"label":"label","span":{"begin":44112,"end":44120}},{"label":"caption","span":{"begin":44122,"end":45134}},{"label":"title","span":{"begin":44122,"end":45134}},{"label":"table","span":{"begin":45135,"end":45794}},{"label":"tr","span":{"begin":45135,"end":45168}},{"label":"th","span":{"begin":45135,"end":45136}},{"label":"th","span":{"begin":45138,"end":45139}},{"label":"th","span":{"begin":45141,"end":45142}},{"label":"th","span":{"begin":45144,"end":45145}},{"label":"th","span":{"begin":45147,"end":45148}},{"label":"th","span":{"begin":45150,"end":45151}},{"label":"th","span":{"begin":45153,"end":45154}},{"label":"th","span":{"begin":45156,"end":45157}},{"label":"th","span":{"begin":45159,"end":45160}},{"label":"th","span":{"begin":45162,"end":45164}},{"label":"th","span":{"begin":45166,"end":45168}},{"label":"tr","span":{"begin":45169,"end":45220}},{"label":"th","span":{"begin":45169,"end":45176}},{"label":"th","span":{"begin":45178,"end":45180}},{"label":"th","span":{"begin":45182,"end":45188}},{"label":"th","span":{"begin":45190,"end":45194}},{"label":"th","span":{"begin":45196,"end":45198}},{"label":"th","span":{"begin":45200,"end":45202}},{"label":"th","span":{"begin":45204,"end":45208}},{"label":"th","span":{"begin":45210,"end":45214}},{"label":"th","span":{"begin":45216,"end":45220}},{"label":"tr","span":{"begin":45221,"end":45283}},{"label":"th","span":{"begin":45221,"end":45225}},{"label":"th","span":{"begin":45227,"end":45231}},{"label":"th","span":{"begin":45233,"end":45237}},{"label":"th","span":{"begin":45239,"end":45243}},{"label":"th","span":{"begin":45245,"end":45254}},{"label":"th","span":{"begin":45256,"end":45265}},{"label":"th","span":{"begin":45268,"end":45272}},{"label":"th","span":{"begin":45274,"end":45283}},{"label":"tr","span":{"begin":45284,"end":45387}},{"label":"td","span":{"begin":45284,"end":45285}},{"label":"td","span":{"begin":45287,"end":45290}},{"label":"td","span":{"begin":45292,"end":45295}},{"label":"td","span":{"begin":45297,"end":45302}},{"label":"td","span":{"begin":45304,"end":45313}},{"label":"td","span":{"begin":45315,"end":45324}},{"label":"td","span":{"begin":45326,"end":45337}},{"label":"td","span":{"begin":45339,"end":45352}},{"label":"td","span":{"begin":45354,"end":45365}},{"label":"td","span":{"begin":45367,"end":45372}},{"label":"td","span":{"begin":45374,"end":45387}},{"label":"tr","span":{"begin":45388,"end":45492}},{"label":"td","span":{"begin":45388,"end":45393}},{"label":"td","span":{"begin":45395,"end":45398}},{"label":"td","span":{"begin":45400,"end":45402}},{"label":"td","span":{"begin":45404,"end":45410}},{"label":"td","span":{"begin":45412,"end":45421}},{"label":"td","span":{"begin":45423,"end":45432}},{"label":"td","span":{"begin":45434,"end":45445}},{"label":"td","span":{"begin":45447,"end":45458}},{"label":"td","span":{"begin":45460,"end":45471}},{"label":"td","span":{"begin":45473,"end":45479}},{"label":"td","span":{"begin":45481,"end":45492}},{"label":"tr","span":{"begin":45493,"end":45594}},{"label":"td","span":{"begin":45493,"end":45494}},{"label":"td","span":{"begin":45496,"end":45499}},{"label":"td","span":{"begin":45501,"end":45503}},{"label":"td","span":{"begin":45505,"end":45512}},{"label":"td","span":{"begin":45514,"end":45523}},{"label":"td","span":{"begin":45525,"end":45534}},{"label":"td","span":{"begin":45536,"end":45547}},{"label":"td","span":{"begin":45549,"end":45560}},{"label":"td","span":{"begin":45562,"end":45573}},{"label":"td","span":{"begin":45575,"end":45581}},{"label":"td","span":{"begin":45583,"end":45594}},{"label":"tr","span":{"begin":45595,"end":45693}},{"label":"td","span":{"begin":45595,"end":45596}},{"label":"td","span":{"begin":45598,"end":45601}},{"label":"td","span":{"begin":45603,"end":45605}},{"label":"td","span":{"begin":45607,"end":45612}},{"label":"td","span":{"begin":45614,"end":45623}},{"label":"td","span":{"begin":45625,"end":45634}},{"label":"td","span":{"begin":45636,"end":45647}},{"label":"td","span":{"begin":45649,"end":45660}},{"label":"td","span":{"begin":45662,"end":45673}},{"label":"td","span":{"begin":45675,"end":45680}},{"label":"td","span":{"begin":45682,"end":45693}},{"label":"tr","span":{"begin":45694,"end":45794}},{"label":"td","span":{"begin":45694,"end":45695}},{"label":"td","span":{"begin":45697,"end":45700}},{"label":"td","span":{"begin":45702,"end":45704}},{"label":"td","span":{"begin":45706,"end":45712}},{"label":"td","span":{"begin":45714,"end":45723}},{"label":"td","span":{"begin":45725,"end":45734}},{"label":"td","span":{"begin":45736,"end":45747}},{"label":"td","span":{"begin":45749,"end":45760}},{"label":"td","span":{"begin":45762,"end":45773}},{"label":"td","span":{"begin":45775,"end":45781}},{"label":"td","span":{"begin":45783,"end":45794}},{"label":"p","span":{"begin":45795,"end":46879}},{"label":"p","span":{"begin":46880,"end":47642}},{"label":"sec","span":{"begin":47644,"end":49411}},{"label":"title","span":{"begin":47644,"end":47698}},{"label":"p","span":{"begin":47699,"end":49411}},{"label":"sec","span":{"begin":49413,"end":60286}},{"label":"title","span":{"begin":49413,"end":49423}},{"label":"p","span":{"begin":49424,"end":50481}},{"label":"p","span":{"begin":50482,"end":51516}},{"label":"p","span":{"begin":51517,"end":53121}},{"label":"p","span":{"begin":53122,"end":53879}},{"label":"p","span":{"begin":53880,"end":54884}},{"label":"p","span":{"begin":54885,"end":55947}},{"label":"p","span":{"begin":55948,"end":56645}},{"label":"p","span":{"begin":56646,"end":57192}},{"label":"p","span":{"begin":57193,"end":58842}},{"label":"p","span":{"begin":58843,"end":59914}},{"label":"p","span":{"begin":59915,"end":60286}},{"label":"sec","span":{"begin":60288,"end":60367}},{"label":"title","span":{"begin":60288,"end":60310}},{"label":"label","span":{"begin":60311,"end":60329}},{"label":"caption","span":{"begin":60331,"end":60367}},{"label":"p","span":{"begin":60331,"end":60367}},{"label":"back","span":{"begin":60369,"end":61060}},{"label":"sec","span":{"begin":60369,"end":60435}},{"label":"title","span":{"begin":60369,"end":60387}},{"label":"p","span":{"begin":60388,"end":60435}},{"label":"sec","span":{"begin":60437,"end":60613}},{"label":"title","span":{"begin":60437,"end":60453}},{"label":"p","span":{"begin":60454,"end":60613}},{"label":"sec","span":{"begin":60615,"end":61060}},{"label":"title","span":{"begin":60615,"end":60622}},{"label":"p","span":{"begin":60623,"end":61013}},{"label":"p","span":{"begin":61014,"end":61061}}],"tracks":[{"project":"2_test","denotations":[{"id":"29361033-22289453-77258000","span":{"begin":1490,"end":1491},"obj":"22289453"},{"id":"29361033-22983709-77258001","span":{"begin":1912,"end":1913},"obj":"22983709"},{"id":"29361033-16907618-77258002","span":{"begin":1955,"end":1956},"obj":"16907618"},{"id":"29361033-16672973-77258003","span":{"begin":1975,"end":1977},"obj":"16672973"},{"id":"29361033-19058276-77258004","span":{"begin":1993,"end":1995},"obj":"19058276"},{"id":"29361033-18634880-77258005","span":{"begin":2013,"end":2015},"obj":"18634880"},{"id":"29361033-21980531-77258006","span":{"begin":2033,"end":2035},"obj":"21980531"},{"id":"29361033-18181653-77258007","span":{"begin":2064,"end":2066},"obj":"18181653"},{"id":"29361033-23396281-77258008","span":{"begin":2086,"end":2088},"obj":"23396281"},{"id":"29361033-21961607-77258009","span":{"begin":2089,"end":2091},"obj":"21961607"},{"id":"29361033-4924006-77258010","span":{"begin":2387,"end":2389},"obj":"4924006"},{"id":"29361033-9223305-77258011","span":{"begin":3024,"end":3026},"obj":"9223305"},{"id":"29361033-22723426-77258012","span":{"begin":3814,"end":3816},"obj":"22723426"},{"id":"29361033-15465864-77258013","span":{"begin":4240,"end":4242},"obj":"15465864"},{"id":"29361033-16585517-77258014","span":{"begin":6656,"end":6658},"obj":"16585517"},{"id":"29361033-16585517-77258015","span":{"begin":9458,"end":9460},"obj":"16585517"},{"id":"29361033-27571472-77258016","span":{"begin":9461,"end":9463},"obj":"27571472"},{"id":"29361033-16585517-77258017","span":{"begin":10069,"end":10071},"obj":"16585517"},{"id":"29361033-16585517-77258018","span":{"begin":10666,"end":10668},"obj":"16585517"},{"id":"29361033-16585517-77258019","span":{"begin":11185,"end":11187},"obj":"16585517"},{"id":"29361033-27627329-77258020","span":{"begin":11587,"end":11589},"obj":"27627329"},{"id":"29361033-16585517-77258021","span":{"begin":13599,"end":13601},"obj":"16585517"},{"id":"29361033-16585517-77258022","span":{"begin":14044,"end":14046},"obj":"16585517"},{"id":"29361033-16585517-77258023","span":{"begin":14181,"end":14183},"obj":"16585517"},{"id":"29361033-16585517-77258024","span":{"begin":36720,"end":36722},"obj":"16585517"},{"id":"29361033-15465864-77258025","span":{"begin":38455,"end":38457},"obj":"15465864"},{"id":"29361033-15465864-77258026","span":{"begin":39315,"end":39317},"obj":"15465864"},{"id":"29361033-15465864-77258027","span":{"begin":39450,"end":39452},"obj":"15465864"},{"id":"29361033-16585517-77258028","span":{"begin":44923,"end":44925},"obj":"16585517"},{"id":"29361033-16585517-77258029","span":{"begin":45051,"end":45053},"obj":"16585517"},{"id":"29361033-16585517-77258030","span":{"begin":49167,"end":49169},"obj":"16585517"},{"id":"29361033-19411704-77258031","span":{"begin":50599,"end":50600},"obj":"19411704"},{"id":"29361033-16585517-77258032","span":{"begin":52949,"end":52951},"obj":"16585517"},{"id":"29361033-19898474-77258033","span":{"begin":54070,"end":54071},"obj":"19898474"},{"id":"29361033-16585517-77258034","span":{"begin":54072,"end":54074},"obj":"16585517"},{"id":"29361033-25016526-77258035","span":{"begin":54324,"end":54326},"obj":"25016526"},{"id":"29361033-25016526-77258036","span":{"begin":54379,"end":54381},"obj":"25016526"},{"id":"29361033-22231272-77258037","span":{"begin":54382,"end":54384},"obj":"22231272"},{"id":"29361033-25495964-77258038","span":{"begin":54385,"end":54387},"obj":"25495964"},{"id":"29361033-16585517-77258039","span":{"begin":55347,"end":55349},"obj":"16585517"},{"id":"29361033-16585517-77258040","span":{"begin":55423,"end":55425},"obj":"16585517"},{"id":"29361033-19898474-77258041","span":{"begin":55468,"end":55469},"obj":"19898474"},{"id":"29361033-25016526-77258042","span":{"begin":55723,"end":55725},"obj":"25016526"},{"id":"29361033-19339520-77258043","span":{"begin":56124,"end":56126},"obj":"19339520"},{"id":"29361033-19339520-77258044","span":{"begin":56301,"end":56303},"obj":"19339520"},{"id":"29361033-23506309-77258045","span":{"begin":56304,"end":56306},"obj":"23506309"},{"id":"29361033-9929475-77258046","span":{"begin":57453,"end":57455},"obj":"9929475"},{"id":"29361033-27176327-77258047","span":{"begin":58021,"end":58023},"obj":"27176327"},{"id":"29361033-27176327-77258048","span":{"begin":58395,"end":58397},"obj":"27176327"},{"id":"29361033-26400158-77258049","span":{"begin":58716,"end":58718},"obj":"26400158"},{"id":"29361033-27627329-77258050","span":{"begin":59627,"end":59629},"obj":"27627329"},{"id":"29361033-28008086-77258051","span":{"begin":60060,"end":60062},"obj":"28008086"}],"attributes":[{"subj":"29361033-22289453-77258000","pred":"source","obj":"2_test"},{"subj":"29361033-22983709-77258001","pred":"source","obj":"2_test"},{"subj":"29361033-16907618-77258002","pred":"source","obj":"2_test"},{"subj":"29361033-16672973-77258003","pred":"source","obj":"2_test"},{"subj":"29361033-19058276-77258004","pred":"source","obj":"2_test"},{"subj":"29361033-18634880-77258005","pred":"source","obj":"2_test"},{"subj":"29361033-21980531-77258006","pred":"source","obj":"2_test"},{"subj":"29361033-18181653-77258007","pred":"source","obj":"2_test"},{"subj":"29361033-23396281-77258008","pred":"source","obj":"2_test"},{"subj":"29361033-21961607-77258009","pred":"source","obj":"2_test"},{"subj":"29361033-4924006-77258010","pred":"source","obj":"2_test"},{"subj":"29361033-9223305-77258011","pred":"source","obj":"2_test"},{"subj":"29361033-22723426-77258012","pred":"source","obj":"2_test"},{"subj":"29361033-15465864-77258013","pred":"source","obj":"2_test"},{"subj":"29361033-16585517-77258014","pred":"source","obj":"2_test"},{"subj":"29361033-16585517-77258015","pred":"source","obj":"2_test"},{"subj":"29361033-27571472-77258016","pred":"source","obj":"2_test"},{"subj":"29361033-16585517-77258017","pred":"source","obj":"2_test"},{"subj":"29361033-16585517-77258018","pred":"source","obj":"2_test"},{"subj":"29361033-16585517-77258019","pred":"source","obj":"2_test"},{"subj":"29361033-27627329-77258020","pred":"source","obj":"2_test"},{"subj":"29361033-16585517-77258021","pred":"source","obj":"2_test"},{"subj":"29361033-16585517-77258022","pred":"source","obj":"2_test"},{"subj":"29361033-16585517-77258023","pred":"source","obj":"2_test"},{"subj":"29361033-16585517-77258024","pred":"source","obj":"2_test"},{"subj":"29361033-15465864-77258025","pred":"source","obj":"2_test"},{"subj":"29361033-15465864-77258026","pred":"source","obj":"2_test"},{"subj":"29361033-15465864-77258027","pred":"source","obj":"2_test"},{"subj":"29361033-16585517-77258028","pred":"source","obj":"2_test"},{"subj":"29361033-16585517-77258029","pred":"source","obj":"2_test"},{"subj":"29361033-16585517-77258030","pred":"source","obj":"2_test"},{"subj":"29361033-19411704-77258031","pred":"source","obj":"2_test"},{"subj":"29361033-16585517-77258032","pred":"source","obj":"2_test"},{"subj":"29361033-19898474-77258033","pred":"source","obj":"2_test"},{"subj":"29361033-16585517-77258034","pred":"source","obj":"2_test"},{"subj":"29361033-25016526-77258035","pred":"source","obj":"2_test"},{"subj":"29361033-25016526-77258036","pred":"source","obj":"2_test"},{"subj":"29361033-22231272-77258037","pred":"source","obj":"2_test"},{"subj":"29361033-25495964-77258038","pred":"source","obj":"2_test"},{"subj":"29361033-16585517-77258039","pred":"source","obj":"2_test"},{"subj":"29361033-16585517-77258040","pred":"source","obj":"2_test"},{"subj":"29361033-19898474-77258041","pred":"source","obj":"2_test"},{"subj":"29361033-25016526-77258042","pred":"source","obj":"2_test"},{"subj":"29361033-19339520-77258043","pred":"source","obj":"2_test"},{"subj":"29361033-19339520-77258044","pred":"source","obj":"2_test"},{"subj":"29361033-23506309-77258045","pred":"source","obj":"2_test"},{"subj":"29361033-9929475-77258046","pred":"source","obj":"2_test"},{"subj":"29361033-27176327-77258047","pred":"source","obj":"2_test"},{"subj":"29361033-27176327-77258048","pred":"source","obj":"2_test"},{"subj":"29361033-26400158-77258049","pred":"source","obj":"2_test"},{"subj":"29361033-27627329-77258050","pred":"source","obj":"2_test"},{"subj":"29361033-28008086-77258051","pred":"source","obj":"2_test"}]}],"config":{"attribute types":[{"pred":"source","value type":"selection","values":[{"id":"2_test","color":"#adec93","default":true}]}]}}