Skip to main content

Main menu

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
    • Front Matter Portal
    • Journal Club
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Home
Home
  • Log in
  • My Cart

Advanced Search

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
    • Front Matter Portal
    • Journal Club
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
Biological Sciences

Calcium-dependent folding of single calmodulin molecules

Johannes Stigler and Matthias Rief
  1. aPhysik Department E22, Technische Universität München, James-Franck-Strasse, 85748 Garching, Germany; and
  2. bMunich Center for Integrated Protein Science, 81377 München, Germany

See allHide authors and affiliations

PNAS October 30, 2012 109 (44) 17814-17819; https://doi.org/10.1073/pnas.1201801109
Johannes Stigler
aPhysik Department E22, Technische Universität München, James-Franck-Strasse, 85748 Garching, Germany; and
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Matthias Rief
aPhysik Department E22, Technische Universität München, James-Franck-Strasse, 85748 Garching, Germany; and
bMunich Center for Integrated Protein Science, 81377 München, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: Matthias.Rief@mytum.de
  1. Edited by Alan R. Fersht, Medical Research Council Laboratory of Molecular Biology, Cambridge, United Kingdom, and approved June 12, 2012 (received for review March 30, 2012)

  • Article
  • Figures & SI
  • Info & Metrics
  • PDF
Loading

Abstract

Calmodulin is the primary calcium binding protein in living cells. Its function and structure depend strongly on calcium concentration. We used single molecule force spectroscopy by optical tweezers to study the folding of calmodulin in the physiologically relevant range. We find that full-length calmodulin switches from a rich and complex folding behavior at high calcium to a simple folding pathway at apo conditions. Using truncation mutants, we studied the individual domains separately. Folding and stability of the individual domains differ significantly at low calcium concentrations. With increasing calcium, the folding rate constants increase while unfolding rate constants decrease. The complete kinetic as well as energetic behavior of both domains could be modeled using a calcium-dependent three-pathway model. We find that the dominant folding pathway at high calcium concentrations proceeds via a transition state capable of binding one calcium ion. The folding of calmodulin seems to be designed to occur fast robustly over a large range of calcium concentrations and hence energetic stabilities.

  • protein folding
  • EF hand
  • ligand

Proteins interact with numerous partners when they function in our body. This interaction naturally also affects their structure and energetics. In extreme cases, proteins can only fold if they interact with a specific binding partner (1). Understanding of protein structure formation hence requires that the folding process be studied in the presence of those binding partners.

Calmodulin (CaM) is the most important calcium binding protein in our body and is involved in the regulation of numerous calcium-dependent pathways (2). Calcium binding as well as peptide binding properties have been studied in detail over the last decades (3⇓⇓⇓–7). Despite its outstanding relevance, however, not much is known about the folding of calmodulin at physiologically relevant calcium concentrations. At very high calcium concentrations, we reported a complex folding pathway of calmodulin in a recent study (8). In this earlier study, we found that while calmodulin folds fast and robustly when its individual domains fold sequentially, its overall folding time is slowed whenever it folds into one of two off-pathway intermediates. Those intermediates involved, in one case, the non-native pairing of the EF hands 2 and 3, and in the other case the collapse of EF hand 3 onto the already folded N-terminal domain in a non-native conformation. We found that all of the on and off pathway states bind one calcium per folded EF-hand (Fig. 1). However, in vivo calcium concentrations are low. At those conditions, calmodulin switches between structurally distinct apo and holo conformations. Calcium-calmodulin dependent regulation of cellular processes critically relies on those conformations. A first characterization of the calcium-dependent folding of calmodulin at high calcium concentrations was performed in earlier experiments using AFM (9). However, technical limitations in resolution and drift stability limited the accessible range in calcium concentration as well as kinetics. In this paper, we use optical tweezers (10⇓–12) to study single molecule folding/unfolding properties of calmodulin around physiologically relevant calcium concentrations.

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

Cartoon of the folding network of full-length calmodulin at high calcium. At high calcium, full-length calmodulin folds in a complex network of states (8). From the unfolded state (U, red), the molecule can proceed to the native state (F1234, purple) via an intermediate with only the C-terminal domain folded (F34, green), or an intermediate where only the N-terminal domain is folded (F12, light blue). From F12, an additional EF hand can fold to form the off-pathway intermediate F123 (dark blue). From U, the protein can populate another off-pathway intermediate where the EF hands 2 and 3 form a pair with non-native contacts (F23, orange).

Results

In a first set of experiments, we investigated the folding/unfolding mechanics of full-length calmodulin at calcium concentrations of 10 mM, 100 μM, and 0 μM Ca2+. We recorded both stretch and relax cycles at a pulling velocity of 500 nm/s (Fig. 2 A–C, upper traces) as well as made measurements where the trap centers were held at a constant separation and hence imposed an average force-bias on the fluctuating molecule (Fig. 2 A–C, lower traces). Assignment of states indicated by different colors in the constant trap separation measurements was performed using a Hidden-Markov-Model (SI Text).

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

Full-length calmodulin at varying calcium concentrations. Stretch-and-relax cycles (Upper) and constant trap separation traces (Lower) of full-length calmodulin at 10 mM Ca2+ (A), 100 μM Ca2+ (B) and 0 mM Ca2+/10 mM EDTA conditions (C). The constant trap separation traces are shown at intermediate biasing forces where all states are populated. Raw data points are colored after classification by a Hidden-Markov model. At 10 mM Ca2+ and 100 μM Ca2+, all six states described in Fig. 1 were identified. At 0 mM Ca2+/10 mM EDTA conditions, only two states remained. The colors in the lower parts of (A) and (B) correspond to the coloring in Fig. 1. The colors in the lower part of (C) are red: unfolded, dark blue: N-domain folded.

At 10 mM Ca2+, we find the complex pattern of various intermediate states reported earlier (8). In brief, calmodulin can populate the six states sketched in Fig. 1: fully folded (purple), N-terminal domain unfolded (green), C-terminal domain unfolded (light blue), one EF-hand of the C-terminal domain unfolded (dark blue), non-native pairing of EF-hands 2 and 3 (orange), fully unfolded (red). At concentrations of 100 μM, the pattern of populated states essentially remains, albeit at lower forces (8) (Fig. 2B). It is important to note that even at a concentration of 100 μM, calcium binding sites are saturated with calcium (13). At 0 mM calcium (Fig. 2C), the observed kinetic pattern changes drastically. The richness of intermediate states, observed under high calcium, can no longer be observed. In stretch-and-relax cycles (Fig. 2C, Upper), a rapid conformational transition close to equilibrium can be clearly observed at forces around 5 pN. This is further supported by constant trap separation measurements (lower traces of Fig. 2C). A two state transition (dark blue and red) with millisecond kinetics can be observed. The contour length change during this transition (24.4 nm, Fig. 2C) suggests that, in full length calmodulin, only one domain exhibits sufficient thermodynamic stability so that its folding under force can be readily observed in our experiments. An additional transition at lower lengths can only be vaguely inferred from contour lengths fit to the very low force regime (< 3 pN) of the force-extension curve (Fig. 2C). To investigate the changes in kinetics and stability of calmodulin at various calcium concentrations in more detail, we studied the N- and C-terminal domains separately using truncation constructs (Figs. 3 and 4 and Table S1). In the following, the N-terminal domain containing EF-hands 1 and 2 is called CaM12; the C-terminal domain with EF-hands 3 and 4 is called CaM34 (SI Text).

Fig. 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 3.

The N-terminal domain at calcium concentrations of 10 mM (Left), 100 μM (Middle) and 0 mM (Right). Shown are data from three individual molecules. (A) Stretch (black) and relax (blue) cycles at 500 nm/s. Dashed lines are fits to the folded and unfolded stretches. (B) Constant-trap separation traces at intermediate forces. The raw data are colored according to the classification by a Hidden-Markov model. Red: unfolded state, light blue: domain folded with calcium bound, dark blue: domain folded with no calcium bound. (C) Force-lifetime scatter plot for the traces shown in (B). Each dot is represented by the average force of one dwell plotted against its lifetime. (D) Force-dependent folding (red) and unfolding (blue) rate constants. Dotted lines are fits to the data.

Fig. 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 4.

The C-terminal domain at calcium concentrations of 10 mM (Left), 100 μM (Middle), and 0 mM (Right). Shown are data from three individual molecules. (A) Stretch (black) and relax (blue) cycles at 500 nm/s. Dashed lines are fits to the folded and unfolded stretches. (B) Constant-trap separation traces at intermediate forces. Arrows indicate occurrences of short-lived folded states incompatible with lifetimes of the light-green folded state. The raw data are colored according to the classification by a Hidden-Markov model. Red: unfolded state, light green: domain folded with calcium bound, dark green: domain folded with no calcium bound. (C) Force-lifetime scatter plot for the traces shown in (B). Each dot is represented by the average force of one dwell plotted against its lifetime. (D) Force-dependent folding (red) and unfolding (green) rate constants. Dotted lines are fits to the data.

Isolated N-Terminal Domain.

Sample traces for stretch and relax curves, as well as constant trap separation measurements obtained at different calcium concentrations on CaM12, are shown in Fig. 3. We observe clear two-state behavior under all calcium conditions. Forces drop while the midpoint transition kinetics get faster with decreasing calcium concentration. To distinguish the holo conformation of calcium-bound calmodulin from the calcium-free apo conformation, we used different coloring for the folded state (light blue vs. dark blue). In the absence of calcium under 10 mM EDTA (Fig. 3 Right), the transitions are indistinguishable both in length and kinetics from the transition observed in full-length calmodulin (Fig. 2C). The dwell-times and forces of the respective states are shown in the scatter plots in Fig. 3C. This analysis corroborates our assumption of a two state process. All dwell-times follow a single exponential distribution (Fig. S1A and SI Text). For a more complete kinetic characterization, we measured the force-dependence of the folding and unfolding rate constants (Fig. 3D). We found that both the folding and the unfolding branches shift as calcium concentrations change. To obtain folding/unfolding rate constants under zero-load conditions, we extrapolated the measured data points to zero force using established methods [(14), see Methods for details, dotted lines in Fig. 3D]. The curvature in the extrapolation of the folding branch is due to the non-linear entropic elasticity of the unfolded polypeptide contracting against load. Since transition state positions for unfolding under force generally lie close to the folded state, we used a linear extrapolation (Bell model) for the unfolding branch.

Isolated C-Terminal Domain.

Similar calcium-dependent measurements were performed for CaM34. In stretch and relax cycles, a clear drop in mechanical stability can be observed with decreasing calcium concentrations. The differences between holo (Fig. 4, Left) and apo states (Fig. 4, Right) are even more pronounced than for CaM12. As with CaM12, the folding/unfolding kinetics become faster with decreasing calcium concentrations. A Hidden-Markov analysis shows that kinetics is two-state at very high calcium (Fig. 4, Left). In the absence of calcium, the midpoint forces of the transition drop to very low values of around 3.5 pN (Fig. S2). The respective Hidden-Markov analysis yields dwell times in the millisecond range, even faster than for CaM12. The respective rate vs. force plot is shown in Fig. 4D, Right. An interesting effect can be observed at intermediate calcium concentrations. The constant trap separation measurements (Fig. 4, Middle) show that in addition to the long dwells marked in light green, also very rapid transient populations of the folded state marked in dark green can be observed. The two populations can be easily discriminated in a scatter plot (Fig. 4C). In addition to the stretched out data cloud at long lifetimes (light green), also, a cloud centered at lifetimes below 1 ms can be observed. This lifetime is identical to the expected lifetime of the apo-state at the respective forces of 5–6 pN. Apparently, for CaM34, at low calcium concentrations a mixture of apo and holo states can be observed online in the single molecule experiments. The rate constants of the whole kinetic network, including the unfolded state (red), the holo state (light green) and the apo state (dark green) can now be extracted from the lifetime measurements (Fig. S1 B and C). The full circles in Fig. 4D indicate folding and unfolding to/from the holo state, while the open circles indicate folding and unfolding to/from the apo state.

Calcium-Dependent Kinetics.

In order to compare the effect of calcium on the folding kinetics, we extrapolated the measured rate constants of folding and unfolding to zero load using the models described above. Fig. 5 A–D shows zero-load folding and unfolding rate constants for CaM12 and CaM34 at calcium concentrations between 0 and 10 mM. Both folding and unfolding rate constants depend on calcium concentration. In the absence of calcium, the unfolding rate constants are approximately 3 s-1 and approximately 50 s-1 for CaM12 and CaM34, respectively (Fig. 5 A and B). Above 10-4 M the unfolding rate constants for CaM12 drop while CaM34 already deviates from apo behavior above 10-5 M. The folding rate constants (Fig. 5 C and D) show an opposite trend. Apo folding rate constants were approximately 13000 s-1 and approximately 1400 s-1 for CaM12 and CaM34. The folding and unfolding rate constants for the apo C-terminal domain are consistent with values obtained in T-jump experiments (4). Folding rate constants for CaM12 rise beyond 10-3 M, while for CaM34 they rise already beyond 10-5 M. The dependence of both folding and unfolding rate constants on calcium can be modeled quantitatively using a kinetic model with a folding transition state able to bind calcium. The kinetic scheme we used is shown in Fig. 6. We assume that the exchange of calcium ions in the folded as well as unfolded state occurs in rapid equilibrium (15, 16). Calcium exchange in the transition state, however, is assumed to be slow compared to its life-time (15, 16). Analytical solutions for the model can be found in the methods section (SI Text).

Fig. 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 5.

Kinetic and equilibrium information for CaM12 (left) and CaM34 (Right) at varying calcium concentrations. (A, B) Zero-force extrapolated unfolding rate constants. (C, D) Zero-force extrapolated folding rate constants. (E, F) Equilibrium free energies of folding. Dashed lines in (C) and (D) are the calculated diffusion-limited on-rates for one calcium ion assuming a rate limit of 108 M-1 s-1. Filled circles represent averages over several molecules at a particular concentration. Empty circles in (B, D, F) are values for the short-lived state as shown in Fig. 4, Middle. Error bars represent the SEM. Colored continuous lines are fits to the data. Vertical grey lines in (A–D) represent the value of KTS1 as obtained from the fit. The shaded areas are the concentration range below Graphic, where the protein is predominantly ligand-free.

Fig. 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 6.

Model for folding and unfolding of an isolated domain. The folded state can exist in three conformations: unligated, with one calcium ion bound, and with two calcium ions bound. Similarly, the unfolded state can exist in all three conformations. Folding and unfolding can occur via a transition state that also can bind up to two calcium ions. We assume that the barrier crossing times are faster than the equilibrium of the transition state. Using this scheme, we can calculate the effective folding and unfolding rates. Of the two possible conformations for N1 and U1, only one is shown. They are indistinguishable in our assay.

A global fit to both the unfolding and folding rate data in Fig. 5 A–D using literature values for the calcium affinities KN1, KN2 [log 10KN1 = -4.62, log 10KN2 = -5.17 in molar units for CaM12, log 10KN1 = -5.32, log 10KN2 = -6.21, for CaM34 (13)] of the folded state was applied to the data. The unfolded state was assumed to have low affinity KU1, KU2 for calcium in excess of 1 M. The only free parameters were the folding and unfolding rate constants at 0 M calcium and the calcium affinities KTS1, KTS2 of the transition state for folding. In molar units, log 10 KTS1 was -3.6 for CaM12 and -5.4 for CaM34. Consistently, for both domains, the affinity for binding of a second calcium ion to the transitions state KTS2 was very low. Hence, in the probed range of calcium concentrations, the transition state for folding/unfolding of calmodulin can only bind one calcium ion with a lower affinity of the N-terminal domain as compared to the C-terminal domain. As pointed out earlier, CaM34 at low calcium concentrations exhibits a mixture of apo and holo states. The open circles in Fig. 5 show the calcium-dependence of the short apo states. As expected, the folding/unfolding rate constants for apo calmodulin do not depend on calcium concentration and exhibit values consistent with the measured values under apo conditions (solid dark green circles in Fig. 5).

Equilibrium Free Energies.

Beyond kinetic information, our assay allows us to extract equilibrium folding free energies from the data taken at various loads. Using probability distributions for the population of the folded and the unfolded state as a function of force (Fig. S3), the free energy difference between those two states can be calculated by taking into account the known elastic properties of the DNA and peptide linkers as well as the trap spring constant (12) (SI Text). The apo values we find for CaM12 and CaM34 are 7.7 and 2.5 kBT, respectively. The calcium dependence (Fig. 5 E and F) can be modeled using an equilibrium version of the model above (5). Beyond the calcium affinities of CaM12 and CaM34 (approximately 10-5 M and approximately 10-6 M, respectively), the stability increases with a slope corresponding to the energy gain of 2· ln(10)kBT per decade of calcium concentration expected from the solution chemical potential when two ions bind (see lines in Fig. 5 E and F).

Discussion

Holo to Apo Transition in Full-Length Calmodulin.

At high calcium concentrations, calmodulin folding and unfolding proceeds through a network of six states (see Fig. 1). This network has been described in detail in an earlier study (8). The high calcium concentrations in the millimolar range used in the earlier study are extreme considering normal cellular conditions where calcium concentrations vary in a range from hundreds of nanomolar to micromolar (17). The goal of this study was observing the changes in calmodulin folding as calcium concentrations fall close to the physiologically relevant range and holo to apo conformational changes occur. As long as calcium concentrations remain so high that all intermediate states within this network remain saturated, a qualitative change in the folding pathways is not to be expected. This is confirmed by the data in Fig. 2B taken at 100 μM. Even though the overall stability and hence folding/unfolding forces of all intermediate states drop as compared to 10 mM calcium, the connectivity as well as the number of states remain. If, however, the ligand is not present anymore, more drastic changes in folding pathways are expected. Such a drastic change indeed occurs when calcium is absent under EDTA conditions (Fig. 2C). Stretch and relax cycles now only reveal a single dominant two-state transition while another transition can only be inferred at very low forces (see dashed fits in Fig. 2C, Upper). From bulk experiments, it is known that under apo conditions the N-terminal domain is more stable than the C-terminal domain. Hence, the dominant transition in Fig. 2C represents unfolding of the N-domain while the less pronounced transition corresponds to unfolding of the C-domain. Interestingly, in full-length calmodulin the apo C-domain appears even weaker than for the isolated construct CaM34 (Fig. 4, Right). This is in accord with an earlier report that in full-length calmodulin under apo conditions, the C-domain is weakened at the expense of the N-domain (5).

Even though, under apo conditions, from our data alone we cannot completely exclude a presence of non-native intermediates like F123 or F23 (see structural cartoons in Fig. 1), the absence of calcium will likely prevent their formation, thus simplifying the overall folding pathway as compared to higher calcium concentrations. Supporting this conclusion, solution experiments have shown that F23 cannot form in the absence of calcium (18). Even though the individual domains of calmodulin fold more slowly at low calcium, an absence of intermediates and kinetic traps can increase its overall folding speed. While the average folding time at 10 mM calcium lies in the range of seconds (8), under apo conditions, in the absence of folding traps the average folding time will be in the millisecond range. Hence at low calcium conditions, which typically are found in a living cell, calmodulin will fold fast and robustly.

Calcium-Dependence of Folding Kinetics and Free Energies of the Isolated Domains.

Our experiments on the isolated N-terminal domain (CaM12) readily reveal this domain as the stable domain also in full-length calmodulin under EDTA conditions. Both the kinetics as well as the midpoint forces of unfolding are identical to the dominant transition observed in wild-type (compare Fig. 2C to Fig. 3 A and B, Right). The isolated C-terminal apo domain (CaM34) (Fig. 4, Right) is considerably less stable than apo CaM12 (2.5 kBT vs. 7.7 kBT). The stability values we extract from our single molecule measurements for both domains are consistent with bulk values (13).

For a discussion of the calcium-dependence of folding and unfolding rate constants, it is important to note that we also observe changes of the transition state positions with calcium concentration. This can be seen in the rate vs. force plots of Fig. 3D where the slopes of the folding as well as unfolding branches differ among different calcium conditions. Hence, to compensate for this effect and simplify modeling of calcium-dependence, we extrapolated all folding and unfolding rate constants to zero force conditions (Graphic and Graphic in Fig. 5 A–D). Interestingly, both folding and unfolding rate constants depend on calcium. In an earlier study, an effect of calcium was only reported on the calmodulin folding rate constant (9). However, in this earlier study, measurements were performed in the presence of a calmodulin binding peptide. Moreover, a separation of folding and unfolding kinetics was not directly possible due to limited resolution in those measurements. The overall dependence of folding and unfolding of CaM12 and CaM34 could be successfully modeled by the kinetic scheme in Fig. 6 (SI Methods). At calcium concentrations below the dissociation constant Graphic for calcium binding, the folding and unfolding rate constants for both CaM12 and CaM34 remain constant (shaded areas in Fig. 5 A–D). Above KD, folding rate constants increase while unfolding rate constants decrease with calcium. The slope of the increase of the folding rate constants towards high calcium concentrations directly reveals that, as soon as KTS1 is exceeded, one calcium ion binds to the transition state during folding even at the highest calcium concentrations measured. A detailed calculation of the number of calcium ions bound to the native state as well as the transition state for both domains can be found in Fig. S4. Hence the dominant folding pathway for both domains at high calcium concentrations (Fig. 6) involves rapid and weak binding of one calcium ion to the unfolded state (U0 → U1) and subsequent fast folding over the one calcium bound transition state TS1 into the native state N1. Binding of a second calcium ion will occur rapidly in a diffusion-controlled manner (N1 → N2). Unfolding occurs in the reverse order. For the N-terminal domain, the transition state apparently exhibits a somewhat lower affinity to calcium (approximately 250 μM) such that in a range between 10 and 250 μM calcium, the domain folds with no calcium bound and two calcium ions then bind later to the folded state. It is important to note that all the structural interpretation we give here do not come from direct time-resolved structural evidence but are a consequence of the transition state model of protein folding (Fig. 6).

The dashed lines in Fig. 5 C and D represent a diffusion limited binding of calcium ions with a rate constant of 108 M-1 s-1 (6, 19). Apparently, at high calcium, the folding rate constants are limited by the binding of calcium to calmodulin. This observation also explains why folding rate constants for the two domains differ under apo conditions but converge at high calcium.

We note that the kinetic scheme of Fig. 6 does not distinguish between a calcium-free protein that adopts a holo-like fold and calmodulin in the apo conformation. A rapid equilibrium of such conformations has been postulated in NMR experiments (20) as well as from MD simulations (21). Since the energetic difference between such conformations is low, we cannot distinguish them in our experiments. Likewise, we do not distinguish between holo calmodulin and apo calmodulin with bound calcium ions (22, 23).

Direct Observation of the Coexistence of Holo and Apo Conformations.

A great advantage of single molecule measurements over bulk studies is their ability to directly observe sub-populations without ensemble averaging. For slowly exchanging ligands with strong affinity, single molecule mechanical measurements have shown that a stable ligand-bound form can coexist with a mechanically weaker apo form (24⇓–26). However, in the case of calmodulin, the exchange is rapid and so far, only average properties could be measured in single molecule assays (9). At calcium concentrations that lie in the transition regime from apo to holo conformations, one would expect to observe a simultaneous coexistence of apo and holo conformations. We could indeed observe such a mixture for CaM34 (Fig. 4, Middle, and Fig. 5, open symbols); however, not for CaM12. At first sight, this comes as a surprise because both domains should exhibit such a coexistence regime albeit at different calcium concentrations. The important difference between the N-terminal domain CaM12 and the C-terminal domain CaM34 lies in their calcium exchange kinetics. At calcium concentrations of 100 μM, the diffusion-limited on-rate constant for calcium binding is 104 s-1 (3, 6). To observe unfolding from an apo state, it is necessary that the apo state lifetime lie below the binding rate constant for calcium, because otherwise, calmodulin would have bound calcium and subsequently switched to the holo conformation. The sample trace in Fig. 4B (Middle) for CaM34 at 100 μM calcium is recorded at forces of 5–6 pN. At these forces, apo CaM34 will unfold at a rate constant of ca. 2000 s-1 (Fig. 4D, Right). This value is comparable to the rate constant of calcium binding and hence, in a significant number of events (approximately 17%), the apo conformation will unfold before it is transformed to the holo conformation and we will record a short lived event (arrows in Fig. 4B, Middle). For CaM12, the forces where folding/unfolding is observed in equilibrium lie around 6.5–7.5 pN (see sample trace in Fig. 3B, Middle). Apo CaM12 will unfold at a rate constant of 400 s-1 at those forces (Fig. 3D, Right). Hence, in most cases (> 96%), calcium will bind before the apo domain unfolds. This rapid binding therefore likely prevents the observation of pure apo states and lifetimes will be distributed single exponentially.

Conclusion

In summary, we were able to study the calcium-dependent folding of a single molecule of calmodulin in a physiologically relevant regime around the apo to holo transition. We were able to identify the kinetic folding pathway in detail. It has been reported that in the calcium free state, calmodulin undergoes rapid transitions between apo and holo-like conformations (20, 21). Distinguishing such populations in our single molecule assay will be an important challenge for the future. So far, this has not been possible, owing to limitations in time and force resolution. Moreover, it will be interesting to see how the kinetic folding pathway couples to the simultaneous interaction with calmodulin binding targets, which is the primary function of this calcium-sensing molecule.

Methods

Experimental Procedures.

CaM12 (CaM34) comprised residues 2–78 (76–149) of human calmodulin (P62158). The construct was genetically inserted between two ubiquitins with terminal cysteines that served as spacers (8). Using the cysteines as attachment points for thiol-modified DNA handles with biotin/digoxigenin functionalized ends of a length of approximately 180 nm each, we attached the construct to micron-sized neutravidin/anti-digoxigenin silica beads in a dumbbell configuration. Data were collected on a custom-built dual beam optical tweezers setup. Measurements were performed in 50 mM Tris/150 mM KCl, pH 8 with varying concentrations of CaCl2 or EDTA (SI Text).

Data Analysis.

Transitions between states were detected using a Hidden-Markov model on the unfiltered raw data of the difference signal of the two traps [(27, 28), SI Text].

Acknowledgments

We thank Fabian Ziegler for technical assistance, and Lorenz Rognoni, Fabian Ziegler, Gabriel Zoldak, and Ziad Ganim for helpful discussions. This work was supported by the SFB 863/A2 project of Deutsche Forschungsgemeinschaft. JS was supported by the Elitenetzwerk Bayern in the framework of the doctorate program Material Science of Complex Interfaces. MR acknowledges funding through a DIP grant of Deutsche Forschungsgemeinschaft.

Footnotes

  • ↵1To whom correspondence should be addressed. E-mail: Matthias.Rief{at}mytum.de.
  • Author contributions: J.S. and M.R. designed research; J.S. performed research; J.S. analyzed data; and J.S. and M.R. wrote the paper.

  • The authors declare no conflict of interest.

  • This article is a PNAS Direct Submission.

  • This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1201801109/-/DCSupplemental.

References

  1. ↵
    1. Turjanski AG,
    2. Gutkind JS,
    3. Best RB,
    4. Hummer G
    (2008) Binding-induced folding of a natively unstructured transcription factor. PLoS Comput Biol 4:e1000060.
    OpenUrlCrossRefPubMed
  2. ↵
    1. Yamniuk AP,
    2. Vogel HJ
    (2004) Calmodulin’s flexibility allows for promiscuity in its interactions with target proteins and peptides. Mol Biotechnol 27:33–57.
    OpenUrlCrossRefPubMed
  3. ↵
    1. Linse S,
    2. Helmersson A,
    3. Forsén S
    (1991) Calcium binding to calmodulin and its globular domains. J Biol Chem 266:8050–8054.
    OpenUrlAbstract/FREE Full Text
  4. ↵
    1. Rabl C-R,
    2. Martin SR,
    3. Neumann E,
    4. Bayley PM
    (2002) Temperature jump kinetic study of the stability of apo-calmodulin. Biophys Chem 101–102:553–564.
    OpenUrl
  5. ↵
    1. Masino L,
    2. Martin SR,
    3. Bayley PM
    (2000) Ligand binding and thermodynamic stability of a multidomain protein, calmodulin. Protein Sci 9:1519–1529.
    OpenUrlPubMed
  6. ↵
    1. Martin SR,
    2. Andersson Teleman A,
    3. Bayley PM,
    4. Drakenberg T,
    5. Forsén S
    (1985) Kinetics of calcium dissociation from calmodulin and its tryptic fragments. A stopped-flow fluorescence study using Quin 2 reveals a two-domain structure. Eur J Biochem 151:543–550.
    OpenUrlPubMed
  7. ↵
    1. Bayley P,
    2. Ahlström P,
    3. Martin SR,
    4. Forsen S
    (1984) The kinetics of calcium binding to calmodulin: Quin 2 and ANS stopped-flow flurescence studies. Biochem Biophys Res Commun 120:185–191.
    OpenUrlCrossRefPubMed
  8. ↵
    1. Stigler J,
    2. Ziegler F,
    3. Gieseke A,
    4. Gebhardt JCM,
    5. Rief M
    (2011) The complex folding network of single calmodulin molecules. Science 334:512–516.
    OpenUrlAbstract/FREE Full Text
  9. ↵
    1. Junker JP,
    2. Ziegler F,
    3. Rief M
    (2009) Ligand-dependent equilibrium fluctuations of single calmodulin molecules. Science 323:633–637.
    OpenUrlAbstract/FREE Full Text
  10. ↵
    1. Cecconi C,
    2. Shank EA,
    3. Bustamante C,
    4. Marqusee S
    (2005) Direct observation of the three-state folding of a single protein molecule. Science 309:2057–2060.
    OpenUrlAbstract/FREE Full Text
  11. ↵
    1. Woodside MT,
    2. et al.
    (2006) Nanomechanical measurements of the sequence-dependent folding landscapes of single nucleic acid hairpins. Proc Natl Acad Sci USA 103:6190–6195.
    OpenUrlAbstract/FREE Full Text
  12. ↵
    1. Gebhardt JCM,
    2. Bornschlögl T,
    3. Rief M
    (2010) Full distance-resolved folding energy landscape of one single protein molecule. Proc Natl Acad Sci USA 107:2013–2018.
    OpenUrlAbstract/FREE Full Text
  13. ↵
    1. Bayley PM,
    2. Findlay WA,
    3. Martin SR
    (1996) Target recognition by calmodulin: Dissecting the kinetics and affinity of interaction using short peptide sequences. Protein Sci 5:1215–1228.
    OpenUrlPubMed
  14. ↵
    1. Schlierf M,
    2. Berkemeier F,
    3. Rief M
    (2007) Direct observation of active protein folding using lock-in force spectroscopy. Biophys J 93:3989–3998.
    OpenUrlCrossRefPubMed
  15. ↵
    1. Fersht AR
    (2004) Phi value versus psi analysis. Proc Natl Acad Sci USA 101:17327–17328.
    OpenUrlFREE Full Text
  16. ↵
    1. Bodenreider C,
    2. Kiefhaber T
    (2005) Interpretation of protein folding psi values. J Mol Biol 351:393–401.
    OpenUrlCrossRefPubMed
  17. ↵
    1. Clapham DE
    (2007) Calcium signaling. Cell 131:1047–1058.
    OpenUrlCrossRefPubMed
  18. ↵
    1. Lakowski T,
    2. Lee G,
    3. Okon M,
    4. Reid R,
    5. McIntosh L
    (2007) Calcium-induced folding of a fragment of calmodulin composed of EF-hands 2 and 3. Protein Sci 16:1119–1132.
    OpenUrlCrossRefPubMed
  19. ↵
    1. Linse S,
    2. et al.
    (1991) Electrostatic contributions to the binding of Ca2+ in calbindin D9k. Biochemistry 30:154–162.
    OpenUrlCrossRefPubMed
  20. ↵
    1. Malmendal A,
    2. Evenäs J,
    3. Forsén S,
    4. Akke M
    (1999) Structural dynamics in the C-terminal domain of calmodulin at low calcium levels. J Mol Biol 293:883–899.
    OpenUrlCrossRefPubMed
  21. ↵
    1. Chen Y-G,
    2. Hummer G
    (2007) Slow conformational dynamics and unfolding of the calmodulin C-terminal domain. J Am Chem Soc 129:2414–2415.
    OpenUrlCrossRefPubMed
  22. ↵
    1. Park HY,
    2. et al.
    (2008) Conformational changes of calmodulin upon Ca2+ binding studied with a microfluidic mixer. Proc Natl Acad Sci USA 105:542–547.
    OpenUrlAbstract/FREE Full Text
  23. ↵
    1. Grabarek Z
    (2005) Structure of a trapped intermediate of calmodulin: Calcium regulation of EF-hand proteins from a new perspective. J Mol Biol 346:1351–1366.
    OpenUrlCrossRefPubMed
  24. ↵
    1. Cao Y,
    2. Er KS,
    3. Parhar R,
    4. Li H
    (2009) A force-spectroscopy-based single-molecule metal-binding assay. Chemphyschem 10:1450–1454.
    OpenUrlCrossRefPubMed
  25. ↵
    1. Cao Y,
    2. Li H
    (2011) Dynamics of protein folding and cofactor binding monitored by single-molecule force spectroscopy. Biophys J 101:2009–2017.
    OpenUrlCrossRefPubMed
  26. ↵
    1. Puchner EM,
    2. Gaub HE
    (2010) Exploring the conformation-regulated function of titin kinase by mechanical pump and probe experiments with single molecules. Angew Chem Int Ed Engl 49:1147–1150.
    OpenUrlCrossRefPubMed
  27. ↵
    1. Rabiner LR
    (1989) A tutorial on Hidden Markov models and selected applications in speech recognition. Proc IEEE 77:257–286.
    OpenUrlCrossRef
  28. ↵
    1. Stigler J,
    2. Rief M
    (2012) Hidden Markov analysis of trajectories in single-molecule experiments and the effects of missed events. Chemphyschem 13:1079–1086.
    OpenUrlCrossRefPubMed
PreviousNext
Back to top
Article Alerts
Email Article

Thank you for your interest in spreading the word on PNAS.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Calcium-dependent folding of single calmodulin molecules
(Your Name) has sent you a message from PNAS
(Your Name) thought you would like to see the PNAS web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Folding at low calcium concentrations
Johannes Stigler, Matthias Rief
Proceedings of the National Academy of Sciences Oct 2012, 109 (44) 17814-17819; DOI: 10.1073/pnas.1201801109

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Request Permissions
Share
Folding at low calcium concentrations
Johannes Stigler, Matthias Rief
Proceedings of the National Academy of Sciences Oct 2012, 109 (44) 17814-17819; DOI: 10.1073/pnas.1201801109
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Mendeley logo Mendeley

Article Classifications

  • Biological Sciences
  • Biophysics and Computational Biology
Proceedings of the National Academy of Sciences: 109 (44)
Table of Contents

Submit

Sign up for Article Alerts

Jump to section

  • Article
    • Abstract
    • Results
    • Discussion
    • Conclusion
    • Methods
    • Acknowledgments
    • Footnotes
    • References
  • Figures & SI
  • Info & Metrics
  • PDF

You May Also be Interested in

Setting sun over a sun-baked dirt landscape
Core Concept: Popular integrated assessment climate policy models have key caveats
Better explicating the strengths and shortcomings of these models will help refine projections and improve transparency in the years ahead.
Image credit: Witsawat.S.
Model of the Amazon forest
News Feature: A sea in the Amazon
Did the Caribbean sweep into the western Amazon millions of years ago, shaping the region’s rich biodiversity?
Image credit: Tacio Cordeiro Bicudo (University of São Paulo, São Paulo, Brazil), Victor Sacek (University of São Paulo, São Paulo, Brazil), and Lucy Reading-Ikkanda (artist).
Syrian archaeological site
Journal Club: In Mesopotamia, early cities may have faltered before climate-driven collapse
Settlements 4,200 years ago may have suffered from overpopulation before drought and lower temperatures ultimately made them unsustainable.
Image credit: Andrea Ricci.
Steamboat Geyser eruption.
Eruption of Steamboat Geyser
Mara Reed and Michael Manga explore why Yellowstone's Steamboat Geyser resumed erupting in 2018.
Listen
Past PodcastsSubscribe
Birds nestling on tree branches
Parent–offspring conflict in songbird fledging
Some songbird parents might improve their own fitness by manipulating their offspring into leaving the nest early, at the cost of fledgling survival, a study finds.
Image credit: Gil Eckrich (photographer).

Similar Articles

Site Logo
Powered by HighWire
  • Submit Manuscript
  • Twitter
  • Facebook
  • RSS Feeds
  • Email Alerts

Articles

  • Current Issue
  • Special Feature Articles – Most Recent
  • List of Issues

PNAS Portals

  • Anthropology
  • Chemistry
  • Classics
  • Front Matter
  • Physics
  • Sustainability Science
  • Teaching Resources

Information

  • Authors
  • Editorial Board
  • Reviewers
  • Subscribers
  • Librarians
  • Press
  • Site Map
  • PNAS Updates
  • FAQs
  • Accessibility Statement
  • Rights & Permissions
  • About
  • Contact

Feedback    Privacy/Legal

Copyright © 2021 National Academy of Sciences. Online ISSN 1091-6490