Mingxi
Yao
a,
Hu
Chen
ab and
Jie
Yan
*acd
aMechanobiology Institute, National University of Singapore, Singapore 117411. E-mail: phyyj@nus.edu.sg
bDepartment of Physics, Xiamen University, Xiamen, Fujian, 361005 China
cDepartment of Physics, National University of Singapore, Singapore 117542
dCentre for Bioimaging Sciences, National University of Singapore, Singapore 117546
First published on 12th March 2015
In this paper, we outline the theoretical framework for understanding the equilibrium force-dependent folding and unfolding transitions of protein domains and small nucleic acid structures, both having small rigid folded structures and highly flexible unfolded polymeric chain conformations. A complete statistical description of the state described by the probability function ρξ(n,x), is obtained, where n is an index denoting the structural state, and x is the extension of the molecule. ξ denotes an external constraint applied to the molecule, which is either a constant force or a harmonic spring attached to one end of the molecule. The extension probability distribution regardless of the structural state: , the free energy landscape: −kBTln(ρξ(x)), and the probability of the states regardless of the extension: , are analyzed using the force-dependent structural transitions of the classic titin I27 domain as an example. The impact of different external constraints is also discussed.
Insight, innovation, integrationSingle molecule manipulation techniques have been widely used to study the force induced structural transitions of biomolecules such as DNA and proteins. With recent advances in technology, increasing measurements are carried out under equilibrium conditions for a long time scale of several hours, in contrast to previous experiments carried out under non-equilibrium conditions over short experimental time scales of seconds to minutes. However, most of the current theoretical frameworks were originally developed to understand experiments carried out under non-equilibrium conditions with a focus on the kinetics of transitions. In this paper, we outlined a systematic theoretical framework to understand the force dependent unfolding and refolding of macromolecules under equilibrium conditions, which is linked to popular experimental systems such as magnetic tweezers, optical tweezers and AFM by introducing different external constraints into the molecules. We believe that such a theoretical description will be useful for both designing experiments and interpreting experimental results. |
The force response of a biomolecule in a particular structural state and the force-dependent transitions between different structural states can be studied at the single-molecule level using single-molecule manipulation techniques. Currently, there are three major types of such techniques available, namely atomic force microscopy (AFM), optical tweezers and magnetic tweezers. They can mechanically stretch a single biomolecule by applying an external mechanical constraint to the molecule. In the cases of AFM and optical tweezers, it can be considered that the molecule is linked to an external harmonic spring at the end of the molecule, where the equilibrium position of the spring and the spring constant can be controlled. In the case of magnetic tweezers, constant forces are applied on the molecule using magnets through a paramagnetic bead attached to one end of the molecule.
In principle, the equilibrium properties of the force-dependent conformations and structural transitions can be obtained under a given external mechanical constraint by observing the extension fluctuation of the molecule. Such measurements require the maintenance of the external mechanical constraint (e.g., a constant force or a spring placed at a fixed position) for a long enough time to allow the molecule to sample all possible conformations and structures by thermal fluctuation. However, equilibrium fluctuations between different structures often involve overcoming high kinetic barriers and therefore have slow transition kinetics. Due to the rapid mechanical drift of the pulling devices, it has been challenging to study the force-dependent structural transitions under equilibrium conditions for many mechanically stable molecules. What has been typically done is forcing the molecule to unfold in a short time by rapidly moving the spring away from the molecule or by rapidly increasing the external force. Data obtained from such non-equilibrium measurements have been interpreted based on kinetic theories6–12 or non-equilibrium statistical mechanics.11
The most recently developed anti-drift technology, magnetic tweezers, has been reported to be able to maintain a constant force for several hours with negligible force drift and spatial drift.13 Such technical improvement has made it possible to probe the equilibrium mechanical folding and unfolding transitions of stable protein domains14–17 and small nucleic acid structures such as G-quadruplexes.18,19 However, a systematic theoretical framework based on equilibrium statistical physics dedicated to the understanding of the mechanical unfolding and folding of small protein domains and nucleic acid structures is still lacking. Motivated by the need for such a systematic theoretical framework for extracting information from equilibrium experimental measurements, we provide here a complete statistical description of the probability function ρξ(n,x), where n is an index denoting the structural state, x is the extension of the molecule, and ξ denotes an external constraint. We demonstrate the application of this probability function by analyzing the force-dependent structural transitions of the classic titin I27 domain, and discuss the impact of different external constraints on experimental observations.
The folded structure is typically a rigid body, which can be modeled as a rigid rod with a rod length of l0 in the absence of mechanical stretching. Its extension is simply x = l0·. In the unfolded state, the molecule is typically an inextensible flexible polymer chain, with an extension of , where L is the polymer chain contour length.
(1) |
In the unfolded state, the conformation of the molecule can be described as an inextensible polymer chain with a weak bending rigidity. The force–extension curve of such polymers is described using the so-called worm-like chain model, which can be approximated analytically using the Marko–Siggia formula:20
(2) |
The above force–extension curves describe the entropic elasticity of the molecule in a particular structural state n. The resulting elastic energies at an extension x are:
(3) |
(4) |
The derivative is the difference between the tension in the molecule and the force applied to the end from the external constraint, which is the mean force applied to the end of the molecule (mean force applied to the link point between the molecule and the constraint). At the equilibrium extension x = xeq,n, the mean force should be zero. Therefore, xeq,n can be solved from the equation Φξ′n(xeq,n) = 0. Considering Φξn(x) as a free energy landscape, xeq,n corresponds to the free energy minimum extension of the structural state n.
In this work, we focus our discussions on two commonly used constraints. One is the constraint of a constant external force of F applied to one end of the molecule, and the other is an external harmonic spring connected to one end of the molecule, which can be described by the equilibrium position of the spring and the spring constant (κ,R). The two different constraints are illustrated in the schematic shown in Fig. 1. wF(x) is the potential energy of the end of the molecule in the constant force field, while w(κ,R)(x) is the spring energy of the spring:
wF(x) = −Fx, | (5a) |
w(κ,R)(x) = ½κ(x − R)2. | (5b) |
(6) |
(7) |
gξ(n,x) = δn,1(Φξn(x) − μ(x)) + δn,2Φξn(x), | (8) |
(9) |
Fig. 2 An example of the extension-dependent folding energy plotted at μ0 = 10 kBT, λ = 0.4 nm, and L1 = 4 nm. |
At equilibrium, ρξ(n,x) should follow the Boltzmann distribution:
(10) |
(11) |
Gξtotal = −kBTlnZξ. | (12) |
(13) |
The so-called free energy landscape Gξ(x) defined as −kBTln(Zξ(x)) is shown as follows:
(14) |
(15) |
(16) |
An I27 domain has 89 amino acids (a.a.) and is folded into two β-sheets with seven β-strands. Treating it as a rigid rod with an l0 of ∼4 nm, its force–extension curve should follow eqn (1). A large value of f0 = 1000 pN is chosen to restrict its strain change within 10% under force <100 pN. After unfolding, it becomes a flexible peptide chain with a contour length of ∼33 nm and its force–extension curve follows the WLC model with a persistence length of A ∼ 0.8 nm,24 which can be described using the Marko–Siggia formula (eqn (2)).20Fig. 3A shows the force–extension curves of the folded and unfolded I27.
Under the constant force constraint described by eqn (5a), the free energy landscape is calculated based on eqn (14) and plotted in Fig. 3B for several different values of F. In the calculations, an extension-dependent folding energy described by eqn (9) with a native folding energy of μ0 = 10 kBT (ref. 17) and a decay constant of λ = 0.4 nm is used. The resulting free energy landscape GF(x) shows two energy minima which are separated by an energy barrier. The regions to the left and right of the energy barrier, which are a few kBT below the barrier, are dominated by the extension fluctuations of the folded and unfolded structures, respectively. This can be seen by plotting GF(x), ΦF1(x) − μ(x) and ΦF2(x) together (Fig. 3C). As a result, the energy minima of GF(x) are basically the equilibrium extensions of the folded state (n = 1) and the unfolded state (n = 2) determined by Φξ′n(xFn,eq) = 0. Due to the rigid body nature, xF1,eq has a weak force dependence and is close to the l0 at force pN; however, due to the flexible chain elasticity of the unfolded peptide, xF2,eq has a strong force dependence and shifts to the right at increasing forces.
The energy barrier corresponds to the maximum located at of GF(x). The “unfolding transition distance” can be defined as: , which is very small and has a length scale close to λ. This is the general outcome of the unfolding of a rigid body which quickly loses its stability upon slight deformation. In contrast, the “folding transition distance” defined as is dominated by the extension of the unfolded state and shows a strong dependence on force. Further, the “unfolding energy barrier” defined as is dominated by the folding energy, which has a weak force dependence and is ∼μ0; while the “folding energy barrier” defined by is dominated by the force response of the unfolded peptide chain with a strong force dependence. The critical force Fc at which the folded and unfolded states have equal probability is estimated by eqn (16) as Fc ∼ 5 pN.
It should be noted that in such a two-state model, the physical meaning of the energy barrier is an extension state at which the mean force applied to the bead from the folded state is counter balanced by that from the unfolded state, weighted by their energies exponentially, , respectively. It differs from the real “transition state” that corresponds to an ensemble of intermediate structural states. Furthermore, in such a two-state description of the unfolding and folding transitions, the transition rates purely depend on the bead diffusion in the free energy landscape, without taking into consideration the intrinsic time scale involved in the structural transitions that often dominates the rate of transitions. In spite of these limitations, the two-state model captures several generic features of the mechanical unfolding and folding transitions of rigid bodies, including a force-insensitive short unfolding transition distance, a force-sensitive folding transition distance, the relative probability of the folded and unfolded states under a given force constraint, etc.
Fig. 4A shows the average (black, left axis) and the variance (red, right axis) of the extension as a function of force calculated by eqn (15). Increasing from 4 pN to 6 pN, the extension increases quickly from that of the folded state to that of the unfolded state. The variance shows a sharp rising profile during the transition with a peak located at Fc ± 1 pN, indicating a large extension fluctuation in this force range. Such a large extension fluctuation is caused by the equilibrium fluctuation between the two structural states. Fig. 4B shows the probabilities of the folded and unfolded states as a function of F calculated by eqn (16), indicating a switch from the folded state to the unfolded state when force is increased over a narrow range from 4 pN to 6 pN.
Fig. 4 (A) The average and the variance of I27 under a constant force constraint of F. (B) The probabilities of the folded and unfolded states of I27 as a function of F. |
Similar calculations can be done for the spring constraint of a spring constant (κ,R). Fig. 5A shows the free energy landscape Gκ,R(x) obtained at R = 10 nm and different values of κ from 0.5 pN nm−1 to 50 pN nm−1. For κ = 0.5 pN nm−1, two energy minima are observed, with the folded state being more stable, as indicated by its lower level of energy compared to that of the unfolded state. At κ = 5 pN nm−1, two energy minima are still observed, but the folded state becomes less stable than the unfolded state. At κ = 50 pN nm−1, the spring completely dominates the fluctuation of the bead. The extension of the molecule is nearly stretched to the equilibrium position of the spring, and it is in an unfolded state since the extension is much greater than the contour length of the folded state. Fig. 5B shows the probabilities of structural states of the molecule as a function of the spring position (R) for three different values of κ. In general, it shows that, compared to a stiffer spring, a softer spring needs to move a longer distance to unfold the molecule, with a less steep transition profile.
In typical AFM or optical tweezer experiments, the trap/cantilever position is controlled. The extension of the molecule indicated by the bead/cantilever-tip position is measured, and the force is a readout based on the calibrated spring constant of the trap/cantilever stiffness. Fig. 6A shows the average extension of I27 as a function of R calculated by eqn (15) for three different values of κ. It shows that during unfolding, a larger extension jump can be observed under a softer spring constraint than under a stiffer spring constraint. Fig. 6B shows the force change during unfolding as a function of the average extension, revealing that a larger force drop under a stiffer spring constraint than under a softer spring constraint. Importantly, both the magnitude of the unfolding force and the overall force–extension profile show a significant dependence on the spring stiffness. These results indicate that when a molecule is stretched by a spring, both the molecule and the spring should be considered as a combined thermodynamic system. The interpretations of the observed force–extension relationship during transitions must be based on the context of the spring constant used in the experiments.
Fig. 6 (A) Extension of I27 as a function of the spring position R, and (B) force as a function of the average extension, calculated at different spring constants κ. |
Two commonly used external mechanical constraints are discussed: one is the constant force constraint typically implemented in magnetic tweezer experiments, and the other is the external spring constraint typically implemented in AFM/optical tweezer experiments. We have shown that different constraints lead to different experimental observations during the mechanical unfolding and folding of the molecules. In particular, in the case of the external spring constraint, interpretation of experimental data should consider the molecule and the spring as a combined thermodynamic system.
All equations are derived based on equilibrium statistics. In the derived free energy landscape, the difference in the free energies of the folded and unfolded states determines their relative probabilities under a constraint. The energy barrier in the free energy landscape is mainly contributed by the folding energy of the folded state, which does not provide a realistic description of the actual transition state that typically corresponds to an ensemble of intermediate structural states, which are different for different molecules. With a slight modification of the free energy landscape described in the paper, it is possible to describe the kinetics of real experiments more realistically, for example by changing the height of the energy barrier treating it as a free parameter to match the barrier height estimated in experiments.
This journal is © The Royal Society of Chemistry 2015 |