From: W. Trevor King Date: Thu, 11 Mar 2010 17:43:00 +0000 (-0500) Subject: Addressed Mom's concerns about WLC indexing and fixed some sawsim typos X-Git-Tag: v1.0~430 X-Git-Url: http://git.tremily.us/?a=commitdiff_plain;h=672ecfe30ae5ab7cf1a174523e0b377c23ed2cde;p=thesis.git Addressed Mom's concerns about WLC indexing and fixed some sawsim typos --- diff --git a/tex/src/cantilever/methods.tex b/tex/src/cantilever/methods.tex index a1bcdb9..b60d1b5 100644 --- a/tex/src/cantilever/methods.tex +++ b/tex/src/cantilever/methods.tex @@ -29,10 +29,10 @@ Mechanical unfolding experiments were carried out on I27 octomers (AthenaES) in PBS on gold-coated coverslips. We used both cantilevers on Olympus's OMCL-TR400-PSA-1, which are nominally $80$ and $20\U{pN/nm}$. Promising sawtooth curves were selected by eye and fit -to WLCs to identify I27 unfolding events. The results were sorted -into two bins according to cantilever stiffness, and then averaged -across each cantilever-stiffness/pulling-speed group to produce the -following graph. +to WLCs\index{WLC} to identify I27 unfolding events. The results were +sorted into two bins according to cantilever stiffness, and then +averaged across each cantilever-stiffness/pulling-speed group to +produce the following graph. \begin{figure} \includegraphics[width=4in]{figures/cantilever-data/plot_splits} @@ -44,7 +44,7 @@ unfolding parameters $k$ or $\Delta x$. Note that the increase in mean unfolding force is not entirely due to the increased loading rate of the stiffer cantilever, because the difference is still present in the loading rate dependence. The loading rates were extracted from -the data by taking the slope of the fit WLC at unfolding. +the data by taking the slope of the fit WLC\index{WLC} at unfolding. \begin{figure} \asyfig{figures/cantilever-data/loading-rate} diff --git a/tex/src/cantilever/simulations.tex b/tex/src/cantilever/simulations.tex index 85b13e2..7b14616 100644 --- a/tex/src/cantilever/simulations.tex +++ b/tex/src/cantilever/simulations.tex @@ -2,9 +2,9 @@ We can model the I27 multimers as an array of Bell-model unfolders in series with a cantilever. Any unfolded domains also contribute to the -tension according to the WLC tension formula. Completing 1000 -simulated pulls for each cantilever/pulling-speed/multimer-number -combination with our \sawsim Monte Carlo simulator yielded the +tension according to the WLC\index{WLC} tension formula. Completing +1000 simulated pulls for each cantilever/pulling-speed/multimer-number +combination with our \sawsim\ Monte Carlo simulator yielded the following results (\cref{fig:cant:sim:v-dep,fig:cant:sim:load-dep,fig:cant:sim:i-dep}). diff --git a/tex/src/cantilever/theory.tex b/tex/src/cantilever/theory.tex index ec74f92..e6ff413 100644 --- a/tex/src/cantilever/theory.tex +++ b/tex/src/cantilever/theory.tex @@ -25,21 +25,15 @@ tension. The Bell-model unfolding rate is thus \end{align*} and stiffer linkers will increase the mean unfolding force. -Unfolded I27 domains can be well-modeled as wormlike-chains -(WLCs)\cite{carrion-vazquez99b}, producing a tension under extension -approximated by Bustamante's interpolation -formula\cite{bustamante94}.%,marko95}. -\begin{align*} - F_\text{WLC} &= \frac{k_B T}{p}\p[{\frac{1}{4}\p({\frac{1}{(1-x/L)^2}-1})+\frac{x}{L}}] \;, -\end{align*} -where $p \approx 4\U{\AA}$ is the persistence length, and $L \approx -28\U{nm}$ is the contour length of the unfolded domain. Obviously -effective stiffness of an unfolded I27 domain is highly dependent on -the unfolding force, and for tensions $\sim 280\U{pN}$ is $\sim -190\U{pN/nm}$. This is within a factor of four of common cantilever -spring constants, so cantilever stiffness drives the effective spring -constant for the first four domains, after which point I27 stiffness -takes the lead. +Unfolded I27 domains can be well-modeled as wormlike chains (WLCs, +\cref{sec:tension:wlc})\cite{carrion-vazquez99b}, where $p \approx +4\U{\AA}$ is the persistence length, and $L \approx 28\U{nm}$ is the +contour length of the unfolded domain. Obviously effective stiffness +of an unfolded I27 domain is highly dependent on the unfolding force, +and for tensions $\sim 280\U{pN}$ is $\sim 190\U{pN/nm}$. This is +within a factor of four of common cantilever spring constants, so +cantilever stiffness drives the effective spring constant for the +first four domains, after which point I27 stiffness takes the lead. \begin{figure} \begin{tikzpicture} % Inspired by Florian Hollandt's RNA codons diff --git a/tex/src/introduction/main.tex b/tex/src/introduction/main.tex index 722992f..94c26e7 100644 --- a/tex/src/introduction/main.tex +++ b/tex/src/introduction/main.tex @@ -190,13 +190,14 @@ time, facilitating single molecule studies. deflection, $x_u$ is the extension of the unfolded polymer, and $x_f=x_{f1}+x_{f2}$ is the extension of the folded polymer. (b) An experimental force curve from stretching a ubiquitin polymer - with the rising parts of the peaks fitted to the WLC - model\citep{chyan04}. The pulling speed used was $1\U{$\mu$m/s}$. - The irregular features at the beginning of the curve are due to - nonspecific interactions between the tip and the substrate - surface, and the last high force peak is caused by the detachment - of the polymer from the tip or the substrate surface. Note that - the abscissa is the extension of the protein chain $x_t-x_c$.} + with the rising parts of the peaks fitted to the WLC\index{WLC} + model (\cref{sec:tension:wlc})\citep{chyan04}. The pulling speed + used was $1\U{$\mu$m/s}$. The irregular features at the beginning + of the curve are due to nonspecific interactions between the tip + and the substrate surface, and the last high force peak is caused + by the detachment of the polymer from the tip or the substrate + surface. Note that the abscissa is the extension of the protein + chain $x_t-x_c$.} \end{center} \end{figure} diff --git a/tex/src/sawsim/discussion.tex b/tex/src/sawsim/discussion.tex index 151aa41..fc0519c 100644 --- a/tex/src/sawsim/discussion.tex +++ b/tex/src/sawsim/discussion.tex @@ -116,7 +116,7 @@ and and average \end{equation} where $N_f$ and $\kappa$ depend on the domain index $i=N_u$. Curves based on this formula fit the simulated data remarkably well considering the -effective WLC stiffness $\kappa_\text{WLC}$ is the only fitted +effective WLC\index{WLC} stiffness $\kappa_\text{WLC}$ is the only fitted parameter, and that the actual WLC stiffness is not constant, as we have assumed here, but a non-linear function of $F$. @@ -174,29 +174,6 @@ unfolding order and polymer length. \end{center} \end{figure} -\subsection{The effect of polymer inhomogeneity} -\label{sec:sawsim:results-folded-tension} - -The unfolded polypeptide chain has been shown to follow the WLC model -quite well, though other polymer models, such as the Freely-Jointed -Chain (FJC)\citep{verdier70}, can be used to fit the force-extension -relationship\citep{janshoff00}. A short chain of folded proteins, -however, cannot be described well by polymer models. Several studies -have used WLC and FJC models to fit the elastic properties of the -modular protein titin\citep{granzier97,linke98a}, -% TODO: check it really is folded domains \& bulk titin -but native titin contains hundreds of folded and unfolded domains -domains. For the short protein polymers common in mechanical -unfolding experiments, the cantilever dominates the elasticity of the -polymer-cantilever system before any protein molecules unfold. After -the first unfolding event occurs, the unfolded portion of the chain is -already longer and softer than the sum of all the remaining folded -domains, and dominates the elasticity of the whole chain. Therefore, -the details of the tension model chosen for the folded domains has -negligible effect on the unfolding forces, which was also suggested by -\citet{staple08}. Force curves simulated using different models to -describe the folded domains yielded almost identical unfolding force -distributions (data not shown). \subsection{The effect of cantilever force constant} \label{sec:sawsim:cantilever} @@ -233,7 +210,7 @@ It should also be mentioned that the contour length increment from each unfolding event is not equal to the distance between adjacent peaks in the force curve because the chain is never fully stretched. This contour length increase can only be obtained by fitting the curve -to WLC or other polymer models (\cref{fig:expt-sawtooth}). +to WLC\index{WLC} or other polymer models (\cref{fig:expt-sawtooth}). \begin{figure} \begin{center} diff --git a/tex/src/sawsim/methods.tex b/tex/src/sawsim/methods.tex index a71516e..a28eac3 100644 --- a/tex/src/sawsim/methods.tex +++ b/tex/src/sawsim/methods.tex @@ -46,20 +46,11 @@ by Hooke's law, which gives F = \kappa_c x_c \;, \label{eq:sawsim:hooke} \end{equation} where $\kappa_c$ is the bending spring constant and $x_c$ is the -deflection of the cantilever (\cref{fig:unfolding-schematic}). Unfolded domains -are modeled as a Worm-Like Chain (WLC)\citep{marko95,bustamante94}, in -which the tension $F$ is related to its extension (end-to-end -distance) $x_u$ by -\begin{equation} - F = \frac{k_B T}{p_u} - \p[{ \frac{1}{4}\p({ \frac{1}{(1-x_u/L_u)^2} - 1 }) - + \frac{x_u}{L_u} }] \;, - \label{eq:sawsim:wlc} -\end{equation} -where $p_u$ is the persistence length and $L_u=N_uL_{u1}$ is the total -contour length of the unfolded chain. The chain of $N_f$ folded -domains is modeled as a string free to assume any extension up to some -fixed contour length $L_f=N_fL_{f1}$ +deflection of the cantilever (\cref{fig:unfolding-schematic}). +Unfolded domains are modeled as WLCs (\cref{sec:tension:wlc}). + +The chain of $N_f$ folded domains is modeled as a string free to +assume any extension up to some fixed contour length $L_f=N_fL_{f1}$ \begin{equation} F = \begin{cases} 0 & \text{if $x_f