Fitting Force-dependent Kinetic Models to Single-Molecule Dwell Time Records
|
|
- Gavin Henderson
- 5 years ago
- Views:
Transcription
1 Fitting Force-dependent Kinetic Models to Single-Molecule Dwell Time Records Abstract Micromanipulation techniques such as optical tweezers atomic force microscopy allow us to identify characterize force-dependent transitions in molecular machines mechanochemical cycles by studying how dwelltime distributions vary as assisting or hindering loads are applied. We present an new simple maximum likelihood method for estimating the zero-force rate distance to the transition state. As a practical benefit to the experimenter, our method allows for relaxed experimental conditions; one no longer has to perform dwell time measurements repeatedly at very precisely controlled fixed forces, but can use data obtained at widely varying forces as long as the forces are measured or known. Moreover, our method is generalizable to more complicated kinetic pathways including those in which some rates do not vary with force. Introduction The development of micromanipulation techniques that allow measurement control of piconewton-scale forces nanometer-scale displacements with time resolution on the order of milliseconds, such as optical tweezers or atomic force microscopy, has made it possible to probe, under mechanical load, the kinetics of biomolecular processes at the single-molecule level. Such high-resolution measurements in combination with novel analysis methods have enabled determination of the distributions of durations of experimentally observable states within the underlying mechanochemical pathway,,3. Such experiments have proved invaluable in the study of biomolecular machines, including the processive motion of cargo-transporting motor proteins such as kinesin myosin-v, enzyme-substrate interactions involving mechanical deformations, such as initiation by T7 RNA polymerase or DNA cleavage by restriction enzymes. 4,5,6 * " ' 4 * 3 * + %,! & / " ', + " # % " # -. Figure : Illustration of the experimental setup for kinesin attached to a bead, translocating along a microtubule subject to optical of the experimental setup for kinesin attached to a bead, translocating tweezersfigure pulling : on Illustration bead a microtubule microscope slide coverslip, subject tweezers is pulling The along distribution of dwellbound timestoinathe experimentally observable statestoofoptical such systems determined by the underlying kinetic pathway. Force dependence of the rate of a single mechanochemical sub-step in a pathway may backusing on bead. be modeled transition-state theory, in which the rate Gthirteen microtubule s structure consists of a sheet of exp typically or fourteen different parallel k. kb T tubulin protofilaments wrapped to form a cylinder, kinesin remains parallel to individual protofilaments so that it effectively undergoes one-dimensional motion Ray et al., 993. Securing individual microtubules to the surface of a microscope coverslip in a predominantly uniform direction, a biophysicist may record the linear displacement of a bead pulled by a
2 Statistics Case Studies: Single Molecule Dwell Times Math 363 Spring 9 Here k B is Boltzmann s constant T is the absolute temperature. Taking displacement to be our reaction coordinate, the free energy G may be split into the sum of the work done against or assisted by the external force the unloaded free energy, G ± F δ where δ is the distance to the transition state the sign depends on whether the applied force helps or hinders the reaction. We thus have k = k e ±F δ/k BT, where k is the zero-force rate, for the form of the force-dependent rates in the system. Current Method Estimation of rates by fitting force-dependent dwell-time distributions to this model has to date been predominantly done by an ad hoc procedure, wherein data are grouped into force bins, lifetime distributions are fit to each bin, the resulting force-dependent rates are fit to???. This procedure only works well when the underlying kinetics consist of a single rate-limiting step force bins are kept narrow by taking many measurements are taken at each of a series of fixed forces. Problems arise when the data cannot, due to the nature of the experiment, be collected at precisely predetermined forces???. Because the dwell times vary nonlinearly with force, the mean dwell time computed from data taken at a range of forces is not an accurate estimator of the mean dwell time at the force at a force-bin s center. We illustrate this in Figure??, which shows the distribution of mean dwell times estimated using a force bin one piconewton in width containing simulated force-dwell time pairs. To compute each mean dwell time, we drew force values from the uniform distribution from 3.5 pn to 4.5 pn, assigned to each force a rom, exponentially-distributed dwell time drawn from a distribution with density ft τ, δ, F = τ ef δ/k BT exp t τ e F δ/k BT with k B T =4. pn nm, δ=nm τ =3 arbitrary time units. The mean dwell time for the bin was then computed by averaging these rom dwell times. This was repeated ten thous times to generate the distribution shown in the figure. It is clear that the mean dwell time is biased towards higher values than the mean expected at the bin s center. This bias does not diminish if we increase the number of points in a bin even to experimentally unobtainable numbers. Moreover, the distribution of such an estimate with points in the bin is non-normal. Similar bias non-normality not shown were found when the forces of a force-bin were drawn from normal distributions. In both cases, the bias gradually decreased as the width of the uniform or normal distribution was lessened. Such narrow force binning is not always experimentally possible, especially when high-force events are intrinsically rare or difficult to measure. 3 Fitting a Single, Force-Dependent Transition Rather than grouping data into force bins, we propose fitting force-dependent kinetic models to data taken at all forces by maximum likelihood. For most real-world probability distributions, maximum likelihood estimators are guaranteed to be consistent, which is to say that they converge to the correct result as the number of samples increases?. Such a method also allows for the fitting of models in which some transitions are force-dependent others are not, removes the need for binning. In the case of a single, force-dependent transition, such an approach yields an elegant graphical method of estimating the mean zero-force lifetime τ or transition rate k distance to the transition state δ. Consider the simple reaction, bound unbound, describing an irreversible detachment of two molecules, e.g. a single myosin head head an actin filament microtubule as in the experiment of Nishizaka et al.?. Since there are no substeps, dwell times in the bound state are
3 Statistics Case Studies: Single Molecule Dwell Times Math 363 Spring 9 5 A 5 Count Estimated τ Figure : Histogram of mean lifetimes inferred from, data sets simulated as described in the text. These sets were then split into force bins. Histograms A, B, C corresponding to fits to,, 5 force bins, respectively. Mean dwell times Figureforces : were Histogram computed forof eachunloaded bin, ˆτ bin was zero-force then fit to these means using lifetimes the stard weighted inferred leastusing squares method. the ad Dashed hoc bin-fitlinemethod indicates expected frommean, lifetimedata force-bin sets center. simulated Note that the estimated as follows: value of τ, stronglyforces depends on were the bindrawn size from that the optimally sized bins cannot be known in advance. U pn, pn. Exponentially distributed rom lifetimes were then computed for each force using drawn from thea single- Arrhenius-like exponential distribution. rate law Allowing Eq, thek transition B T =4. to bepn force dependent nm, δ = as above, nm, the observed τ = k = dwell time given a measured applied force is drawn from the probability distribution having density 3 arbitrary The log-likelihood units. These functionsets for the were experiment then with split independently into force measured bins. forces Histograms F = F, FA,,.. B,., F N C corresponding corresponding to fits independent to, dwell, times t 5 = t force, t,.. bins,., c N isrespectively. Mean lifetimes forces were computed for each bin, τf =τ e F δ/kbt was then fit to these means using the stard weighted least squares method. log Lτ, δ t, Dashed F = n log vertical τ + δ N n lines indicate thet i e Fiδ/k BT correct value of τ. 3 Maximizing the log likelihood requires that both k B T δ log Lτ, δ t, F = τ Nτ τ log Lτ, δ t, F = F Nτ F i τ B C N t i e Fiδ/kBT = N t i F i e Fiδ/kBT = where F is the mean value of the forces. These equations do not separate to produce closed forms for the maximum likelhood estimators ˆδ ˆτ. However, they can be rearranged to obtain ˆτ = N t i e Fi ˆδ/k B T 4 N ˆτ = N F N t i F i e Fi ˆδ/k B T. 5 This facilitates a simple graphical interpretation for ˆτ ˆδ as the intersections of these two functions of ˆδ. 3
4 Statistics Case Studies: Single Molecule Dwell Times Math 363 Spring " Figure 3: Example of the maximum likelihood estimator using values described in the text. This is found by the intersection of 4 5 in the ˆδ-ˆτ plane. The direction of the axes of the ellipse are derived from the covariance of δ τ as given by the off diagonal elements in the inverse of the Fisher information matrix 6. igure 3: Example of the model-fitting technique presented in Section. The intersection of the In addition to providing estimates ˆτ urves given in Equations 3 4, marked ˆδ, the theory of maximum likelihood estimation also provides stard errors of these estimates in a straightforward manner. with Thea covariance, gives matrix of the asymptotic maximum normallikelihood distribution of estimators these estimators has a simple closed form. The expected Fisher information is given by f the unloaded mean time τ the distance to transition state δ. The ellipse demarcates he region of 95% confidence, the true value /τ Iτ, δ = N of F δ/τ k B T τ is marked with a +. Fit was erformed on, simulated data, with forces F /τ k B T F /k B T between zero. ten piconewtons drawn from Its inverse, the covariance matrix, is uniform distribution dwell times drawn from the force-dependent distribution function iven in Eq., δ= nm, τ =3 arbitrary Iτ, δ units, k τ = l F τ B T =4. k B T pn nm. F. 6 NvarF! τ k B T F k B T Iτ, δ Iτ, δ are the respective stard errors of ˆτ ˆδ. Plotted in Figure?? is an example of the result of this procedure. One thous pairs of independent measurements were simulated by drawing a force uniformly distributed between pn. then drawing the corresponding lifetime from the exponential distribution given in Eq., As above, δ= nm, k B =4. pn nm, τ =3 arbitrary units. The true values of δ τ are marked with a +. The intersection, located numerically marked with a, of the curves given in Equations 4 5 plotted gives the maximum likelihood estimators ˆτ =.99 ±.9, ˆδ =.4 ±.45. The region of 95% confidence is demarcated by the ellipse. 4 Discussion The procedure we present in Section is straightforward computationally rapid. It is easier to apply than the ad hoc method in common use. Moreover, for this ad hoc method, the mean force is highly dependent on the choice of bin size consequently is unreliable as an estimation technique. Locating the intersection of two curves is simpler than the force-binning, fitting single-force likelihoods, then fitting curves to parameters. Likelihood functions for even relatively simple models can become rather complicated. Solutions like 4 5 can sometimes be found. For more complex models, the likelihood given or its logarithm may be maximized using an entirely numerical procedure. In addition to avoiding complications surrounding binning, maximum likelihood estimation of kinetic model parameters allows us to make the most of available data, as estimation is done using all of the information collected in the experiment, at once. One no longer has to perform dwell time measurements repeatedly at very precisely controlled fixed forces, but can use data obtained at widely varying forces as long as the forces are measured or known. In the 3
5 Statistics Case Studies: Single Molecule Dwell Times Math 363 Spring 9 Count 5 5 A 5 5 B Count C D ˆτ ˆδ Figure 4: Histograms of ˆτ A,C ˆδ B,D estimated by finding the intersection of the curves given in Eqs A B were Figure constructed 4: Histograms from, sets of of ˆτ A,C force-lifetime ˆδ pairs B,D generated estimated as for Figure by finding??; C the D were intersection constuctedof fromthe, sets curves of, given such pairs. in Eqs. Dashed 3 vertical 4. lines A indicate B were the true constructed values of τ from δ. Compare,tosets Figure of. force-lifetime pairs generated as for Figure 3; C D were constucted from, sets of, such pairs. simple Dashed one-step vertical caselines treatedindicate Section the, true we see values from the of τcovariance δ. Compare matrix Eq. to 6, Figure which scales. as varf, that we benefit from taking measurements over a wide range of forces, incurring a penalty in estimating the unloaded rate as the mean square force becomes high. We expect the same qualitative behavior a benefit for taking measurements at a broad range of forces but a penalty for speeding up the reaction with high assisting forces when treating more complicated models. Estimating the parameters of force-dependent models by maximum likelihood opens the possibility of using model selection techniques such as information criterion or likelihood ratio tests to complement biochemical biophysical experiments. References [] L. S. Milescu, A. Yildiz, P. R. Selvin, F. Sachs, 6 Biophysical Journal 9, [] B. C. Carter, M. Vershinin, S. P. Gross, 8 Biophysical Journal 94, [3] S. N. Block 7 Biophysical Journal 9, [4] A. Mehta Journal of Cel l Science 4, [5] J.-F. Allem, D. Bensimon, V. Croquette 3 Current Opinion in Structural Biology 3, [6] G. M. Skinner, C. G. Baumann, D. M. Quinn, J. E. Molloy, J. G. Hoggett 4 Journal of Biological Chemistry 79,
6 Statistics Case Studies: Single Molecule Dwell Times Math 363 Spring 9 5 Model STATISTICAL PROCEDURES Let t i, i n be the observed values of n independent exponential rom variables with parameter Here, the model parameters, exp F iδ τ k B T. τ is the characteristic disassociation time with no force, δ is the displacement from equilbrium. In addiition, F i is the force associated to the i-th observation, k B is Boltzmann s constant, T is the absolute temperature Thus, an observed time t based on a force F is drawn from a probability distribution having density ft τ, δ, F = τ ef δ/ exp t τ e F δ/, t In particular, the mean E τ,δt = τ e F δ/k BT. 6 Likelihood Function The likelihood function for the experiment The score function has components Lτ, δ t, F = n log Lτ, δ t, F = n log τ + δ τ log Lτ, δ t, F = n δ log Lτ, δ t, F = = n τ efiδ/ exp t i τ e Fiδ/. 5 F i τ τ nτ F i τ F nτ n t i e Fiδ/. t i e Fiδ/. t i F i e Fiδ/ t i F i e Fiδ/
7 Statistics Case Studies: Single Molecule Dwell Times Math 363 Spring 9 Consequently, the maximum likelihood estimates ˆτ, ˆδ satisfy = ˆτ nˆτ = F nˆτ t i e Fiδ/, or ˆτ = n Fi ˆδ/ t i F i e or ˆτ = n F t i e Fi ˆδ/, t i F i e Fi ˆδ/. Thus, to find the maximum likelihood estimates, graph the values of τ in each of the two relations above as a function of δ. The intersection of these two graphs will give the maximum likelihood estimates. 7 Fisher Information Matrix Write θ = τ, δ let ˆθ n be the maximum likelihood estimator based on n observations. We shall see that in these circumstances that ˆθ n is consistent, that is, if θ is the state of nature, In addition, ˆθ n θ in probability as n. nˆθn θ W in distribution as n where W is a mean zero normal rom vector with covariance matrix the inverse of the Fisher information matrix Iθ. Consequently, the maximum likelihood estimator is asymptotically efficient in the sense that the ratio of any other consistent estimator with the variance of this estimator is at most. To find the Fisher information matrix, note that Next, τ I τ, δ = E τ,δ = n τ log Lτ, δ t, F = n [ τ nτ 3 τ nτ 3 ] log Lτ, δ t, F = n τ t i e Fiδ/. nτ 3 e Fiδ/ E τ,δt i e Fiδ/ τ e Fiδ/ = n δ log Lτ, δ t, F = n nτ t i F i F i efiδ/ [ ] I τ, δ = E τ,δ δ log Lτ, δ t, F = Finally, for the mutual information, F i = n F δτ log Lτ, δ t, F = τ 6 = τ τ nτ = τ = n τ. t i Fi e Fiδ/ Fi e Fiδ/ E τ,δt i t i F i e Fiδ/
8 Statistics Case Studies: Single Molecule Dwell Times Math 363 Spring 9 [ ] I τ, δ = E τ,δ log Lτ, δ t, F δτ = τ F i = n F τ = τ F i e Fiδ/ E τ,δt i Thus, the information matrix /τ F /τ Iτ, δ = n F /τ F / The deteminant is varf /τ the inverse Iτ, δ = τ F / F /τ = τ F τ F nvarf F /τ /τ nvarf τ F Now, we can use stard z-statistics to perform interval estimation hypothesis testing.. 7
Anatoly B. Kolomeisky. Department of Chemistry CAN WE UNDERSTAND THE COMPLEX DYNAMICS OF MOTOR PROTEINS USING SIMPLE STOCHASTIC MODELS?
Anatoly B. Kolomeisky Department of Chemistry CAN WE UNDERSTAND THE COMPLEX DYNAMICS OF MOTOR PROTEINS USING SIMPLE STOCHASTIC MODELS? Motor Proteins Enzymes that convert the chemical energy into mechanical
More informationTransport of single molecules along the periodic parallel lattices with coupling
THE JOURNAL OF CHEMICAL PHYSICS 124 204901 2006 Transport of single molecules along the periodic parallel lattices with coupling Evgeny B. Stukalin The James Franck Institute The University of Chicago
More informationFor slowly varying probabilities, the continuum form of these equations is. = (r + d)p T (x) (u + l)p D (x) ar x p T(x, t) + a2 r
3.2 Molecular Motors A variety of cellular processes requiring mechanical work, such as movement, transport and packaging material, are performed with the aid of protein motors. These molecules consume
More informationNIH Public Access Author Manuscript J Phys Condens Matter. Author manuscript; available in PMC 2014 November 20.
NIH Public Access Author Manuscript Published in final edited form as: J Phys Condens Matter. 2013 November 20; 25(46):. doi:10.1088/0953-8984/25/46/463101. Motor Proteins and Molecular Motors: How to
More informationSupplementary Information
Supplementary Information Switching of myosin-v motion between the lever-arm swing and Brownian search-and-catch Keisuke Fujita 1*, Mitsuhiro Iwaki 2,3*, Atsuko H. Iwane 1, Lorenzo Marcucci 1 & Toshio
More informationLecture 3. G. Cowan. Lecture 3 page 1. Lectures on Statistical Data Analysis
Lecture 3 1 Probability (90 min.) Definition, Bayes theorem, probability densities and their properties, catalogue of pdfs, Monte Carlo 2 Statistical tests (90 min.) general concepts, test statistics,
More informationLecture 7 : Molecular Motors. Dr Eileen Nugent
Lecture 7 : Molecular Motors Dr Eileen Nugent Molecular Motors Energy Sources: Protonmotive Force, ATP Single Molecule Biophysical Techniques : Optical Tweezers, Atomic Force Microscopy, Single Molecule
More informationActo-myosin: from muscles to single molecules. Justin Molloy MRC National Institute for Medical Research LONDON
Acto-myosin: from muscles to single molecules. Justin Molloy MRC National Institute for Medical Research LONDON Energy in Biological systems: 1 Photon = 400 pn.nm 1 ATP = 100 pn.nm 1 Ion moving across
More informationStep-Stress Models and Associated Inference
Department of Mathematics & Statistics Indian Institute of Technology Kanpur August 19, 2014 Outline Accelerated Life Test 1 Accelerated Life Test 2 3 4 5 6 7 Outline Accelerated Life Test 1 Accelerated
More informationSINGLE-MOLECULE PHYSIOLOGY
SINGLE-MOLECULE PHYSIOLOGY Kazuhiko Kinosita, Jr. Center for Integrative Bioscience, Okazaki National Research Institutes Higashiyama 5-1, Myodaiji, Okazaki 444-8585, Japan Single-Molecule Physiology under
More informationStatistical Data Analysis Stat 3: p-values, parameter estimation
Statistical Data Analysis Stat 3: p-values, parameter estimation London Postgraduate Lectures on Particle Physics; University of London MSci course PH4515 Glen Cowan Physics Department Royal Holloway,
More informationSUPPLEMENTARY INFORMATION
doi:10.1038/nature09450 Supplementary Table 1 Summary of kinetic parameters. Kinetic parameters were V = V / 1 K / ATP and obtained using the relationships max ( + m [ ]) V d s /( 1/ k [ ATP] + 1 k ) =,
More informationσ(a) = a N (x; 0, 1 2 ) dx. σ(a) = Φ(a) =
Until now we have always worked with likelihoods and prior distributions that were conjugate to each other, allowing the computation of the posterior distribution to be done in closed form. Unfortunately,
More informationarxiv: v1 [physics.bio-ph] 9 Aug 2011
Phenomenological analysis of ATP dependence of motor protein Yunxin Zhang Laboratory of Mathematics for Nonlinear Science, Centre for Computational System Biology, School of Mathematical Sciences, Fudan
More informationSECOND PUBLIC EXAMINATION. Honour School of Physics Part C: 4 Year Course. Honour School of Physics and Philosophy Part C C7: BIOLOGICAL PHYSICS
2757 SECOND PUBLIC EXAMINATION Honour School of Physics Part C: 4 Year Course Honour School of Physics and Philosophy Part C C7: BIOLOGICAL PHYSICS TRINITY TERM 2011 Monday, 27 June, 9.30 am 12.30 pm Answer
More informationMS&E 226: Small Data
MS&E 226: Small Data Lecture 12: Frequentist properties of estimators (v4) Ramesh Johari ramesh.johari@stanford.edu 1 / 39 Frequentist inference 2 / 39 Thinking like a frequentist Suppose that for some
More informationSec. 2.1 Filaments in the cell 21 PART I - RODS AND ROPES
Sec. 2.1 Filaments in the cell 21 PART I - RODS AND ROPES Sec. 2.1 Filaments in the cell 22 CHAPTER 2 - POLYMERS The structural elements of the cell can be broadly classified as filaments or sheets, where
More informationSTATS 200: Introduction to Statistical Inference. Lecture 29: Course review
STATS 200: Introduction to Statistical Inference Lecture 29: Course review Course review We started in Lecture 1 with a fundamental assumption: Data is a realization of a random process. The goal throughout
More informationTecniche sperimentali: le optical tweezers
Tecniche sperimentali: le optical tweezers Le tecniche di molecola singola rispetto a quelle di insieme ensemble z z z z z z frequency activity activity time z z z Single molecule frequency activity activity
More informationBiochemistry in Singulo:
Biochemistry in Singulo: When Less Means More Carlos Bustamante University of California, Berkeley The Ensemble Approach The study of chemical transforma9ons has been dominated by the ensemble method:
More informationStatistics. Lecture 2 August 7, 2000 Frank Porter Caltech. The Fundamentals; Point Estimation. Maximum Likelihood, Least Squares and All That
Statistics Lecture 2 August 7, 2000 Frank Porter Caltech The plan for these lectures: The Fundamentals; Point Estimation Maximum Likelihood, Least Squares and All That What is a Confidence Interval? Interval
More informationRegularity and synchrony in motor proteins
Regularity and synchrony in motor proteins R. E. Lee DeVille and Eric Vanden-Eijnden Courant Institute of Mathematical Sciences, New York University, New York, NY 10012 1 Corresponding author: R. E. Lee
More informationMaximum-Likelihood Estimation: Basic Ideas
Sociology 740 John Fox Lecture Notes Maximum-Likelihood Estimation: Basic Ideas Copyright 2014 by John Fox Maximum-Likelihood Estimation: Basic Ideas 1 I The method of maximum likelihood provides estimators
More informationc 2006 by Prasanth Sankar. All rights reserved.
c 2006 by Prasanth Sankar. All rights reserved. PHENOMENOLOGICAL MODELS OF MOTOR PROTEINS BY PRASANTH SANKAR M. S., University of Illinois at Urbana-Champaign, 2000 DISSERTATION Submitted in partial fulfillment
More informationNanomotors: Nanoscale machines
Nanomotors: Nanoscale machines October 31, 2016 1 Introduction to nanomotors In this part of the course we will study nanomotors. First we will define what we mean by nanomotor. A motor (of any size) is
More informationPolymerization and force generation
Polymerization and force generation by Eric Cytrynbaum April 8, 2008 Equilibrium polymer in a box An equilibrium polymer is a polymer has no source of extraneous energy available to it. This does not mean
More informationSingle-Molecule Methods I - in vitro
Single-Molecule Methods I - in vitro Bo Huang Macromolecules 2014.03.10 F 1 -ATPase: a case study Membrane ADP ATP Rotation of the axle when hydrolyzing ATP Kinosita group, 1997-2005 Single Molecule Methods
More informationChanges in microtubule overlap length regulate kinesin-14-driven microtubule sliding
Supplementary information Changes in microtubule overlap length regulate kinesin-14-driven microtubule sliding Marcus Braun* 1,2,3, Zdenek Lansky* 1,2,3, Agata Szuba 1,2, Friedrich W. Schwarz 1,2, Aniruddha
More informationModern Methods of Data Analysis - WS 07/08
Modern Methods of Data Analysis Lecture VIc (19.11.07) Contents: Maximum Likelihood Fit Maximum Likelihood (I) Assume N measurements of a random variable Assume them to be independent and distributed according
More informationLikelihood-Based Methods
Likelihood-Based Methods Handbook of Spatial Statistics, Chapter 4 Susheela Singh September 22, 2016 OVERVIEW INTRODUCTION MAXIMUM LIKELIHOOD ESTIMATION (ML) RESTRICTED MAXIMUM LIKELIHOOD ESTIMATION (REML)
More informationAnalysis of the AIC Statistic for Optimal Detection of Small Changes in Dynamic Systems
Analysis of the AIC Statistic for Optimal Detection of Small Changes in Dynamic Systems Jeremy S. Conner and Dale E. Seborg Department of Chemical Engineering University of California, Santa Barbara, CA
More informationAdaptive Response of Actin Bundles under Mechanical Stress
Biophysical Journal, Volume 113 Supplemental Information Adaptive Response of Actin Bundles under Mechanical Stress Florian Rückerl, Martin Lenz, Timo Betz, John Manzi, Jean-Louis Martiel, Mahassine Safouane,
More informationStatistics - Lecture One. Outline. Charlotte Wickham 1. Basic ideas about estimation
Statistics - Lecture One Charlotte Wickham wickham@stat.berkeley.edu http://www.stat.berkeley.edu/~wickham/ Outline 1. Basic ideas about estimation 2. Method of Moments 3. Maximum Likelihood 4. Confidence
More informationFree energy recovery in single molecule experiments
Supplementary Material Free energy recovery in single molecule experiments Single molecule force measurements (experimental setup shown in Fig. S1) can be used to determine free-energy differences between
More informationATP binding controls distinct structural transitions. of Escherichia coli DNA gyrase in complex with DNA
Supplementary Information ATP binding controls distinct structural transitions of Escherichia coli DNA gyrase in complex with DNA Aakash Basu, Allyn J. Schoeffler, James M. Berger, and Zev Bryant Table
More informationBMB November 17, Single Molecule Biophysics (I)
BMB 178 2017 November 17, 2017 14. Single Molecule Biophysics (I) Goals 1. Understand the information SM experiments can provide 2. Be acquainted with different SM approaches 3. Learn to interpret SM results
More informationOperation modes of the molecular motor kinesin
PHYSICAL REVIEW E 79, 011917 2009 Operation modes of the molecular motor kinesin S. Liepelt and R. Lipowsky Max Planck Institute of Colloids and Interfaces, Science Park Golm, 14424 Potsdam, Germany *
More informationHandling uncertainties in background shapes: the discrete profiling method
Journal of Instrumentation OPEN ACCESS Handling uncertainties in background shapes: the discrete profiling method To cite this article: P.D. Dauncey et al 5 JINST P45 View the article online for updates
More informationMulti-Ensemble Markov Models and TRAM. Fabian Paul 21-Feb-2018
Multi-Ensemble Markov Models and TRAM Fabian Paul 21-Feb-2018 Outline Free energies Simulation types Boltzmann reweighting Umbrella sampling multi-temperature simulation accelerated MD Analysis methods
More informationBrandon C. Kelly (Harvard Smithsonian Center for Astrophysics)
Brandon C. Kelly (Harvard Smithsonian Center for Astrophysics) Probability quantifies randomness and uncertainty How do I estimate the normalization and logarithmic slope of a X ray continuum, assuming
More informationHomework #4 Physics 498Bio Spring 2012 Prof. Paul Selvin
Assigned Wednesday Feb. 22, 2012: Due Wednesday February 29, 10:30am. Hand in at start of class. Late homework is not accepted. (Solution sets will be posted shortly after deadline.) Note: Marco will give
More informationBMB Class 17, November 30, Single Molecule Biophysics (II)
BMB 178 2018 Class 17, November 30, 2018 15. Single Molecule Biophysics (II) New Advances in Single Molecule Techniques Atomic Force Microscopy Single Molecule Manipulation - optical traps and tweezers
More informationA graph contains a set of nodes (vertices) connected by links (edges or arcs)
BOLTZMANN MACHINES Generative Models Graphical Models A graph contains a set of nodes (vertices) connected by links (edges or arcs) In a probabilistic graphical model, each node represents a random variable,
More informationDimension Reduction (PCA, ICA, CCA, FLD,
Dimension Reduction (PCA, ICA, CCA, FLD, Topic Models) Yi Zhang 10-701, Machine Learning, Spring 2011 April 6 th, 2011 Parts of the PCA slides are from previous 10-701 lectures 1 Outline Dimension reduction
More informationMeasurements of interaction forces in (biological) model systems
Measurements of interaction forces in (biological) model systems Marina Ruths Department of Chemistry, UMass Lowell What can force measurements tell us about a system? Depending on the technique, we might
More informationWebsite: Selected readings Topics Introduction to Cell Biology Analysis of Cell Mechanics Cell
Session 1 Website: http://faculty.washington.edu/nsniadec/me599/w13/ Selected readings Topics Introduction to Cell Biology Analysis of Cell Mechanics Cell Mechanics Modeling Measuring Cell Forces Mechanotransduction
More informationDiscrete Latent Variable Models
Discrete Latent Variable Models Stefano Ermon, Aditya Grover Stanford University Lecture 14 Stefano Ermon, Aditya Grover (AI Lab) Deep Generative Models Lecture 14 1 / 29 Summary Major themes in the course
More informationParameter estimation! and! forecasting! Cristiano Porciani! AIfA, Uni-Bonn!
Parameter estimation! and! forecasting! Cristiano Porciani! AIfA, Uni-Bonn! Questions?! C. Porciani! Estimation & forecasting! 2! Cosmological parameters! A branch of modern cosmological research focuses
More informationParametric Techniques
Parametric Techniques Jason J. Corso SUNY at Buffalo J. Corso (SUNY at Buffalo) Parametric Techniques 1 / 39 Introduction When covering Bayesian Decision Theory, we assumed the full probabilistic structure
More informationSingle molecule biochemistry using optical tweezers
FEBS 20359 FEBS Letters 430 (1998) 23^27 Minireview Single molecule biochemistry using optical tweezers Amit D. Mehta*, Katherine A. Pullen, James A. Spudich Department of Biochemistry, Stanford University
More informationUntangling the Mechanics of Entangled Biopolymers
Untangling the Mechanics of Entangled Biopolymers Rae M. Robertson-Anderson Physics Department University of San Diego students/postdocs: Cole Chapman, PhD Tobias Falzone, PhD Stephanie Gorczyca, USD 16
More informationIntracellular transport
Transport in cells Intracellular transport introduction: transport in cells, molecular players etc. cooperation of motors, forces good and bad transport regulation, traffic issues, Stefan Klumpp image
More informationThe biological motors
Motor proteins The definition of motor proteins Miklós Nyitrai, November 30, 2016 Molecular machines key to understand biological processes machines in the micro/nano-world (unidirectional steps): nm,
More informationt For l = 1 a monomer cannot be destroyed or created from nothing: = b p(2,t) a p(1,t).
IITS: Statistical Physics in Biology Assignment # 5 KU Leuven 5/31/2013 Drift, Diffusion, and Dynamic Instability 1. Treadmilling Actin: Actin filaments are long, asymmetric, polymers involved in a variety
More informationStatistical Inference on Constant Stress Accelerated Life Tests Under Generalized Gamma Lifetime Distributions
Int. Statistical Inst.: Proc. 58th World Statistical Congress, 2011, Dublin (Session CPS040) p.4828 Statistical Inference on Constant Stress Accelerated Life Tests Under Generalized Gamma Lifetime Distributions
More informationLecture 1: Introduction. Keywords: Mathematics as a language, Need of learning mathematics, Applications of mathematics in BIology
NPTEL Syllabus Biomathematics - Video course COURSE OUTLINE Graphs and functions, Derivative of a function, Techniques of differentiation Differentiation and its application in Biology, Finding maxima,
More informationFlexing Protein muscles: How to Pull with a "Burning Rope"
Flexing Protein muscles: How to Pull with a "Burning Rope" J. P. Keener 215 Joint International Conference on via Guangzhou University and Huaihua University Huaihua 8/15 p.1/28 Eukaryotic Chromosomal
More informationAppendix from L. J. Revell, On the Analysis of Evolutionary Change along Single Branches in a Phylogeny
008 by The University of Chicago. All rights reserved.doi: 10.1086/588078 Appendix from L. J. Revell, On the Analysis of Evolutionary Change along Single Branches in a Phylogeny (Am. Nat., vol. 17, no.
More informationKinesins with Extended Neck Linkers: A Chemomechanical Model for Variable-Length Stepping
Bull Math Biol (2012) 74:1066 1097 DOI 10.1007/s11538-011-9697-6 ORIGINAL ARTICLE Kinesins with Extended Neck Linkers: A Chemomechanical Model for Variable-Length Stepping John Hughes William O. Hancock
More informationScanning Tunneling Microscopy
Scanning Tunneling Microscopy References: 1. G. Binnig, H. Rohrer, C. Gerber, and Weibel, Phys. Rev. Lett. 49, 57 (1982); and ibid 50, 120 (1983). 2. J. Chen, Introduction to Scanning Tunneling Microscopy,
More informationBiophysik der Moleküle!
Biophysik der Moleküle!!"#$%&'()*+,-$./0()'$12$34!4! Molecular Motors:! - linear motors" 6. Dec. 2010! Muscle Motors and Cargo Transporting Motors! There are striking structural similarities but functional
More informationAdvanced sampling. fluids of strongly orientation-dependent interactions (e.g., dipoles, hydrogen bonds)
Advanced sampling ChE210D Today's lecture: methods for facilitating equilibration and sampling in complex, frustrated, or slow-evolving systems Difficult-to-simulate systems Practically speaking, one is
More informationPhysics of Cellular materials: Filaments
Physics of Cellular materials: Filaments Tom Chou Dept. of Biomathematics, UCLA, Los Angeles, CA 995-766 (Dated: December 6, ) The basic filamentary structures in a cell are reviewed. Their basic structures
More informationProductive Cooperation among Processive Motors Depends Inversely on Their Mechanochemical Efficiency
386 Biophysical Journal Volume July 2 386 395 Productive Cooperation among Processive Motors Depends Inversely on Their Mechanochemical Efficiency Jonathan W. Driver, D. Kenneth Jamison, Karthik Uppulury,
More informationMolecular Machines and Enzymes
Molecular Machines and Enzymes Principles of functioning of molecular machines Enzymes and catalysis Molecular motors: kinesin 1 NB Queste diapositive sono state preparate per il corso di Biofisica tenuto
More informationBiophys J BioFAST, published on May 5, 2006 as doi: /biophysj
Biophys J BioFAS, published on May 5, 26 as doi:.529/biophysj.5.7954 his un-edited manuscript has been accepted for publication in Biophysical Journal and is freely available on BioFast at http://www.biophysj.org.
More informationParametric Techniques Lecture 3
Parametric Techniques Lecture 3 Jason Corso SUNY at Buffalo 22 January 2009 J. Corso (SUNY at Buffalo) Parametric Techniques Lecture 3 22 January 2009 1 / 39 Introduction In Lecture 2, we learned how to
More informationOptical Tweezers -working principles and applications
Optical Tweezers -working principles and applications Photo taken from the WWW Star Trek Picture Page Biophysics with optical tweezers Optical tweezers use forces of laser radiation pressure to trap small
More informationHypothesis testing: theory and methods
Statistical Methods Warsaw School of Economics November 3, 2017 Statistical hypothesis is the name of any conjecture about unknown parameters of a population distribution. The hypothesis should be verifiable
More informationAdjustments for the Display of Quantized Ion Channel Dwell Times in Histograms with Logarithmic Bins
662 Biophysical Journal Volume 78 February 2000 662 667 Adjustments for the Display of Quantized Ion Channel Dwell Times in Histograms with Logarithmic Bins J. Alex Stark and Stephen B. Hladky National
More informationReaction time distributions in chemical kinetics: Oscillations and other weird behaviors
Introduction The algorithm Results Summary Reaction time distributions in chemical kinetics: Oscillations and other weird behaviors Ramon Xulvi-Brunet Escuela Politécnica Nacional Outline Introduction
More informationF & B Approaches to a simple model
A6523 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Spring 215 http://www.astro.cornell.edu/~cordes/a6523 Lecture 11 Applications: Model comparison Challenges in large-scale surveys
More informationRecent Advances in Optical Tweezers
ANNUAL REVIEWS Further Click here for quick links to Annual Reviews content online, including: Other articles in this volume Top cited articles Top downloaded articles Our comprehensive search Annu. Rev.
More informationUse separation of variables to solve the following differential equations with given initial conditions. y 1 1 y ). y(y 1) = 1
Chapter 11 Differential Equations 11.1 Use separation of variables to solve the following differential equations with given initial conditions. (a) = 2ty, y(0) = 10 (b) = y(1 y), y(0) = 0.5, (Hint: 1 y(y
More informationDecision theory. 1 We may also consider randomized decision rules, where δ maps observed data D to a probability distribution over
Point estimation Suppose we are interested in the value of a parameter θ, for example the unknown bias of a coin. We have already seen how one may use the Bayesian method to reason about θ; namely, we
More informationOn the errors introduced by the naive Bayes independence assumption
On the errors introduced by the naive Bayes independence assumption Author Matthijs de Wachter 3671100 Utrecht University Master Thesis Artificial Intelligence Supervisor Dr. Silja Renooij Department of
More informationMATH 251 Week 6 Not collected, however you are encouraged to approach all problems to prepare for exam
MATH 51 Week 6 Not collected, however you are encouraged to approach all problems to prepare for exam A collection of previous exams could be found at the coordinator s web: http://www.math.psu.edu/tseng/class/m51samples.html
More informationConventional kinesin is a processive motor protein that walks
Kinesin crouches to sprint but resists pushing Michael E. Fisher* and Young C. Kim Institute for Physical Science and Technology, University of Maryland, College Park, MD 20742 Contributed by Michael E.
More informationBasic Sampling Methods
Basic Sampling Methods Sargur Srihari srihari@cedar.buffalo.edu 1 1. Motivation Topics Intractability in ML How sampling can help 2. Ancestral Sampling Using BNs 3. Transforming a Uniform Distribution
More informationOutline of GLMs. Definitions
Outline of GLMs Definitions This is a short outline of GLM details, adapted from the book Nonparametric Regression and Generalized Linear Models, by Green and Silverman. The responses Y i have density
More informationConfidence Distribution
Confidence Distribution Xie and Singh (2013): Confidence distribution, the frequentist distribution estimator of a parameter: A Review Céline Cunen, 15/09/2014 Outline of Article Introduction The concept
More informationStatistics. Lecture 4 August 9, 2000 Frank Porter Caltech. 1. The Fundamentals; Point Estimation. 2. Maximum Likelihood, Least Squares and All That
Statistics Lecture 4 August 9, 2000 Frank Porter Caltech The plan for these lectures: 1. The Fundamentals; Point Estimation 2. Maximum Likelihood, Least Squares and All That 3. What is a Confidence Interval?
More informationWhen do diffusion-limited trajectories become memoryless?
When do diffusion-limited trajectories become memoryless? Maciej Dobrzyński CWI (Center for Mathematics and Computer Science) Kruislaan 413, 1098 SJ Amsterdam, The Netherlands Abstract Stochastic description
More informationExperimental biophysics: Optical tweezer lab Supervisor: Stefan Holm,
Experimental biophysics: Optical tweezer lab :, stefan.holm@ftf.lth.se Written by Jason Beech & Henrik Persson, March 2009. Modified 2014 Karl Adolfsson, 2016 Experimental Biophysics: FAF010F, FYST23,
More informationAn Introduction to Bioinformatics Algorithms Hidden Markov Models
Hidden Markov Models Outline 1. CG-Islands 2. The Fair Bet Casino 3. Hidden Markov Model 4. Decoding Algorithm 5. Forward-Backward Algorithm 6. Profile HMMs 7. HMM Parameter Estimation 8. Viterbi Training
More informationROBERTO BATTITI, MAURO BRUNATO. The LION Way: Machine Learning plus Intelligent Optimization. LIONlab, University of Trento, Italy, Apr 2015
ROBERTO BATTITI, MAURO BRUNATO. The LION Way: Machine Learning plus Intelligent Optimization. LIONlab, University of Trento, Italy, Apr 2015 http://intelligentoptimization.org/lionbook Roberto Battiti
More informationCOS513 LECTURE 8 STATISTICAL CONCEPTS
COS513 LECTURE 8 STATISTICAL CONCEPTS NIKOLAI SLAVOV AND ANKUR PARIKH 1. MAKING MEANINGFUL STATEMENTS FROM JOINT PROBABILITY DISTRIBUTIONS. A graphical model (GM) represents a family of probability distributions
More informationBiophysics Biological soft matter
Biophysics Biological soft matter!"#$%&'(&)%*+,-.& /"#$%("%*+,-.0."122,13$(%4(5+& Biophysics lectures outline Biological soft matter 1. Biopolymers 2. Molecular motors 3. The cytoskeleton Biophysics 1.
More informationLecture 8: Information Theory and Statistics
Lecture 8: Information Theory and Statistics Part II: Hypothesis Testing and I-Hsiang Wang Department of Electrical Engineering National Taiwan University ihwang@ntu.edu.tw December 23, 2015 1 / 50 I-Hsiang
More informationParameter estimation and prediction using Gaussian Processes
Parameter estimation and prediction using Gaussian Processes Yiannis Andrianakis and Peter G. Challenor August, 9 Abstract This is a report on methods of parameter estimation and prediction for Gaussian
More informationVelocity, processivity and individual steps of single myosin V molecules in live cells
Biophysical Journal, Volume 96 Supporting Material Velocity, processivity and individual steps of single myosin V molecules in live cells Paolo Pierobon, Sarra Achouri, Sébastien Courty, Alexander R. Dunn,
More informationNeurite formation & neuronal polarization. The cytoskeletal components of neurons have characteristic distributions and associations
Mechanisms of neuronal migration & Neurite formation & neuronal polarization Paul Letourneau letou001@umn.edu Chapter 16; The Cytoskeleton; Molecular Biology of the Cell, Alberts et al. 1 The cytoskeletal
More informationA Simple Protein Synthesis Model
A Simple Protein Synthesis Model James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University September 3, 213 Outline A Simple Protein Synthesis Model
More informationHST 583 FUNCTIONAL MAGNETIC RESONANCE IMAGING DATA ANALYSIS AND ACQUISITION A REVIEW OF STATISTICS FOR FMRI DATA ANALYSIS
HST 583 FUNCTIONAL MAGNETIC RESONANCE IMAGING DATA ANALYSIS AND ACQUISITION A REVIEW OF STATISTICS FOR FMRI DATA ANALYSIS EMERY N. BROWN AND CHRIS LONG NEUROSCIENCE STATISTICS RESEARCH LABORATORY DEPARTMENT
More informationMachine Learning. Lecture 9: Learning Theory. Feng Li.
Machine Learning Lecture 9: Learning Theory Feng Li fli@sdu.edu.cn https://funglee.github.io School of Computer Science and Technology Shandong University Fall 2018 Why Learning Theory How can we tell
More informationChapter 2 Inference on Mean Residual Life-Overview
Chapter 2 Inference on Mean Residual Life-Overview Statistical inference based on the remaining lifetimes would be intuitively more appealing than the popular hazard function defined as the risk of immediate
More informationUndirected Graphical Models
Outline Hong Chang Institute of Computing Technology, Chinese Academy of Sciences Machine Learning Methods (Fall 2012) Outline Outline I 1 Introduction 2 Properties Properties 3 Generative vs. Conditional
More informationSMASIS A STOCHASTIC MECHANO-CHEMICAL MODEL FOR COOPERATIVE MOTOR PROTEIN DYNAMICS
Proceedings of SMASIS 8 8 ASME Conference on Smart Materials, Adaptive Structures and Intelligent Systems October 8-3, 8, Ellicott City, Maryland, USA SMASIS8-585 A STOCHASTIC MECHANO-CHEMICAL MODEL FOR
More informationBinding Theory Equations for Affinity and Kinetics Analysis
Technology Note #101 Binding Theory Equations for Affinity and Kinetics Analysis This technology note summarizes important equations underlying the theory of binding of solute analytes to surface-tethered
More informationElastic Lever-Arm Model for Myosin V
3792 Biophysical Journal Volume 88 June 2005 3792 3805 Elastic Lever-Arm Model for Myosin V Andrej Vilfan J. Stefan Institute, Ljubljana, Slovenia ABSTRACT We present a mechanochemical model for myosin
More information