BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Christopher P. Calderon

Christopher P. Calderon contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
Quantitative Methodscond-mat.mes-hallcond-mat.softcond-mat.stat-mechApplicationsBiological PhysicsSubcellular Processes

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
8works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

8 published item(s)

preprint2016arXiv

Motion Blur Filtering: A Statistical Approach for Extracting Confinement Forces and Diffusivity from a Single Blurred Trajectory

Single Particle Tracking (SPT) can aid in understanding complex spatio-temporal processes. However, quantifying diffusivity and forces from individual live cell trajectories is complicated by inter- & intra-trajectory kinetic heterogeneity, thermal fluctuations, and statistical temporal dependence inherent to the underlying molecule's time correlated confined dynamics experienced in the cell. Experimental artifacts such as localization uncertainty and motion blur also obscure the data. We introduce a new maximum likelihood estimation (MLE) technique that decouples the above noise sources and systematically treats temporal correlation via a likelihood function (permitting more reliable extraction of effective forces from position vs. time data). Our estimator is demonstrated to be consistent over a wide range of exposure times, diffusion coefficients, and confinement "radii". The algorithm and corresponding software can reliably extract motion parameters independent of exposure time in trajectories exhibiting confined and/or non-stationary dynamics and will aid in directly comparing trajectories obtained from different imaging modalities.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Inferring Latent States and Refining Force Estimates via Hierarchical Dirichlet Process Modeling in Single Particle Tracking Experiments

Optical microscopy provides rich spatio-temporal information characterizing in vivo molecular motion. However, effective forces and other parameters used to summarize molecular motion change over time in live cells due to latent state changes, e.g., changes induced by dynamic micro-environments, photobleaching, and other heterogeneity inherent in biological processes. This study focuses on techniques for analyzing Single Particle Tracking (SPT) data experiencing abrupt state changes. We demonstrate the approach on GFP tagged chromatids experiencing metaphase in yeast cells and probe the effective forces resulting from dynamic interactions that reflect the sum of a number of physical phenomena. State changes are induced by factors such as microtubule dynamics exerting force through the centromere, thermal polymer fluctuations, etc. Simulations are used to demonstrate the relevance of the approach in more general SPT data analyses. Refined force estimates are obtained by adopting and modifying a nonparametric Bayesian modeling technique, the Hierarchical Dirichlet Process Switching Linear Dynamical System (HDP-SLDS), for SPT applications. The HDP-SLDS method shows promise in systematically identifying dynamical regime changes induced by unobserved state changes when the number of underlying states is unknown in advance (a common problem in SPT applications). We expand on the relevance of the HDP-SLDS approach, review the relevant background of Hierarchical Dirichlet Processes, show how to map discrete time HDP-SLDS models to classic SPT models, and discuss limitations of the approach. In addition, we demonstrate new computational techniques for tuning hyperparameters and for checking the statistical consistency of model assumptions directly against individual experimental trajectories; the techniques circumvent the need for "ground-truth" and subjective information.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Correcting for Bias of Molecular Confinement Parameters Induced by Small Time Series Sample Sizes in Single-Molecule Trajectories Containing Measurement Noise

Several single-molecule studies aim to reliably extract parameters characterizing molecular confinement or transient kinetic trapping from experimental observations. Pioneering works from single particle tracking in membrane diffusion studies [Kusumi et al., Biophysical J., 1993] appealed to Mean Square Displacement tools for extracting diffusivity and other parameters quantifying the degree of confinement. More recently, the practical utility of systematically treating multiple noise sources (including noise induced by random photon counts) through likelihood techniques have been more broadly realized in the SPT community. However, bias induced by finite time series sample sizes has not received great attention. Mitigating parameter bias induced by finite sampling is important to any scientific endeavor aiming for high accuracy, but correcting for bias is also often an important step in the construction of optimal parameter estimates. In this article, it is demonstrated how a popular model of confinement can be corrected for finite sample bias in situations where the underlying data exhibits Brownian diffusion and observations are measured with non-negligible experimental noise (e.g., noise induced by finite photon counts). The work of Tang and Chen [J. Econometrics, 2009] is extended to correct for bias in the estimated corral radius (a parameter commonly used to quantify confinement in SPT studies) in the presence of measurement noise. It is shown that the approach presented is capable of reliably extracting the corral radius using only hundreds of discretely sampled observations in situations where other methods (including MSD and Bayesian techniques) would encounter serious difficulties. The ability to accurately statistically characterize transient confinement suggests new techniques for quantifying confined and/or hop diffusion in complex environments.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Robust Hypothesis Tests for Detecting Statistical Evidence of 2D and 3D Interactions in Single-Molecule Measurements

A variety of experimental techniques have improved the 2D and 3D spatial resolution that can be extracted from \emph{in vivo} single-molecule measurements. This enables researchers to quantitatively infer the magnitude and directionality of forces experienced by biomolecules in their native cellular environments. Situations where such forces are biologically relevant range from mitosis to directed transport of protein cargo along cytoskeletal structures. Models commonly applied to quantify single-molecule dynamics assume that effective forces and velocity in the $x,y$ (or $x,y,z$) directions are statistically independent, but this assumption is physically unrealistic in many situations. We present a hypothesis testing approach capable of determining if there is evidence of statistical dependence between positional coordinates in experimentally measured trajectories; if the hypothesis of independence between spatial coordinates is rejected, then a new model accounting for 2D (3D) interactions should be considered to more faithfully represent the underlying experimental kinetics. The technique is robust in the sense that 2D (3D) interactions can be detected via statistical hypothesis testing even if there is substantial inconsistency between the physical particle's actual noise sources and the simplified model's assumed noise structure. For example, 2D (3D) interactions can be reliably detected even if the researcher assumes normal diffusion, but the experimental data experiences "anomalous diffusion" and/or is subjected to a measurement noise characterized by a distribution differing from that assumed by the fitted model. The approach is demonstrated on control simulations and on experimental data (IFT88 directed transport in the primary cilium).

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Uncertainty Quantification of Discrete Association Problems in Image Sequence-based Tracking

Applications, ranging from tracking molecular motion within cells to analyzing complex animal foraging behavior, require algorithms for associating a collection of spot-like particles in one image with particles contained in another image. These associations are often made via network flow algorithms. However, it is often the case that many candidate association solutions (the output of network flow algorithms) have nearly optimal scores; in this case, the optimal assignment solution is of dubious quality. Algorithms for reliably computing the uncertainty of candidate association solutions are under-developed in situations where many particles are tracked over multiple frames of data. This is due in part to the fact that exact uncertainty quantification (UQ) in large association problems is computationally intractable because the exact computation exhibits exponential dependence on the number of particles tracked. We introduce a technique that can accurately and efficiently quantify association ambiguity (i.e., UQ for discrete association problems) without requiring the evaluation of the cost of each feasible association solution. Our method can readily be wrapped around existing tracking algorithms and can efficiently handle a variety of 2D association problems. The applications presented are focused on tracking molecules in live cells. Our method is validated via both simulations and experiments. The experimental applications aim to accurately form tracks and quantify diffusion of quantum dot labeled proteins from \emph{in vivo} measurements; here association problems involving cost matrices possessing hundreds to thousands of rows/columns are encountered. For such large-scale problems, we discuss how our approach can efficiently and accurately quantify inherent uncertainty in candidate data associations.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A Random Force is a Force, of Course, of Coarse: Decomposing Complex Enzyme Kinetics with Surrogate Models

The temporal autocorrelation (AC) function associated with monitoring order parameters characterizing conformational fluctuations of an enzyme is analyzed using a collection of surrogate models. The surrogates considered are phenomenological stochastic differential equation (SDE) models. It is demonstrated how an ensemble of such surrogate models, each surrogate being calibrated from a single trajectory, indirectly contains information about unresolved conformational degrees of freedom. This ensemble can be used to construct complex temporal ACs associated with a "non-Markovian" process. The ensemble of surrogates approach allows researchers to consider models more flexible than a mixture of exponentials to describe relaxation times and at the same time gain physical information about the system. The relevance of this type of analysis to matching single-molecule experiments to computer simulations and how more complex stochastic processes can emerge from a mixture of simpler processes is also discussed. The ideas are illustrated on a toy SDE model and on molecular dynamics simulations of the enzyme dihydrofolate reductase.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

On the Use of Local Diffusion Models for Path Ensemble Averaging in Potential of Mean Force Computations

We use a constant velocity steered molecular dynamics (SMD) simulation of the stretching of deca-alanine in vacuum to demonstrate a technique that can be used to create surrogate stochastic processes using the time series that come out of SMD simulations. The surrogate processes are constructed by first estimating a sequence of local parametric models along a SMD trajectory and then a single global model is constructed by piecing the local models together through smoothing splines (estimation is made computationally feasible by likelihood function approximations). The calibrated surrogate models are then "bootstrapped" in order to simulate the large number of work paths typically needed to construct a potential of mean force (PMF) by appealing to Jarzynski's work theorem. When this procedure is repeated for a small number of SMD paths, it is shown that the global models appear to come from a single family of closely related diffusion processes. Possible techniques for exploiting this observation are also briefly discussed. The findings of this paper have potential relevance to computationally expensive computer simulations and experimental works involving optical tweezers where it difficult to collect a large number of samples, but possible to sample accurately and frequently in time.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Fitting Effective Diffusion Models to Data Associated with a "Glassy Potential": Estimation, Classical Inference Procedures and Some Heuristics

A variety of researchers have successfully obtained the parameters of low dimensional diffusion models using the data that comes out of atomistic simulations. This naturally raises a variety of questions about efficient estimation, goodness-of-fit tests, and confidence interval estimation. The first part of this article uses maximum likelihood estimation to obtain the parameters of a diffusion model from a scalar time series. I address numerical issues associated with attempting to realize asymptotic statistics results with moderate sample sizes in the presence of exact and approximated transition densities. Approximate transition densities are used because the analytic solution of a transition density associated with a parametric diffusion model is often unknown.I am primarily interested in how well the deterministic transition density expansions of Ait-Sahalia capture the curvature of the transition density in (idealized) situations that occur when one carries out simulations in the presence of a "glassy" interaction potential. Accurate approximation of the curvature of the transition density is desirable because it can be used to quantify the goodness-of-fit of the model and to calculate asymptotic confidence intervals of the estimated parameters. The second part of this paper contributes a heuristic estimation technique for approximating a nonlinear diffusion model. A "global" nonlinear model is obtained by taking a batch of time series and applying simple local models to portions of the data. I demonstrate the technique on a diffusion model with a known transition density and on data generated by the Stochastic Simulation Algorithm.

0 citations0 reviewsNo code linkedOpen access