The Effect of Macromolecular Crowding, Ionic Strength and Calcium Binding on Calmodulin Dynamics

The flexibility in the structure of calmodulin (CaM) allows its binding to over 300 target proteins in the cell. To investigate the structure-function relationship of CaM, we combined methods of computer simulation and experiments based on circular dichroism (CD) to investigate the structural characteristics of CaM that influence its target recognition in crowded cell-like conditions. We developed a unique multiscale solution of charges computed from quantum chemistry, together with protein reconstruction, coarse-grained molecular simulations, and statistical physics, to represent the charge distribution in the transition from apoCaM to holoCaM upon calcium binding. Computationally, we found that increased levels of macromolecular crowding, in addition to calcium binding and ionic strength typical of that found inside cells, can impact the conformation, helicity and the EF hand orientation of CaM. Because EF hand orientation impacts the affinity of calcium binding and the specificity of CaM’s target selection, our results may provide unique insight into understanding the promiscuous behavior of calmodulin in target selection inside cells.

💡 Research Summary

Calmodulin (CaM) is a highly acidic, 148‑residue protein that contains four EF‑hand calcium‑binding motifs organized into two lobes connected by a flexible linker. Its ability to bind over 300 target proteins stems from its remarkable conformational flexibility, which is modulated by calcium binding, ionic strength, and the crowded intracellular milieu. In this study, Wang et al. combined circular dichroism (CD) spectroscopy, thermal denaturation assays, and a novel multiscale computational framework to dissect how these three factors influence CaM structure and dynamics.

Experimentally, CD spectra of apo‑CaM displayed characteristic negative bands at 207 nm and 222 nm, indicative of substantial α‑helical content. Upon saturating Ca²⁺, the negative bands intensified, consistent with either increased helicity or a re‑orientation of existing helices—a conclusion supported by NMR data showing similar helix content in apo‑ and holo‑states. Thermal scans at 222 nm revealed that apo‑CaM unfolds completely above 60 °C under low‑salt conditions, whereas holo‑CaM retains a partially folded population even at 90 °C, demonstrating calcium‑induced stabilization.

Increasing KCl concentration from 0 to 250 mM produced modest, systematic enhancements in the CD signal, plateauing at 500 mM. Fitting the thermal denaturation curves to a three‑state model (Native ↔ Intermediate ↔ Unfolded) showed that higher ionic strength raises the melting temperatures of both the N‑I and I‑U transitions, especially for the C‑terminal lobe, indicating that electrostatic screening stabilizes the lobes.

Macromolecular crowding was mimicked with Ficoll 70 (100–400 mg ml⁻¹). Crowding caused a pronounced increase in CD ellipticity and shifted the thermal midpoints upward by up to 24 °C. When crowding and salt were combined, the stabilizing effects were additive but not strictly linear; the presence of 100 mM KCl already conferred substantial stability, so the incremental benefit of Ficoll was reduced at the highest concentrations.

On the computational side, the authors introduced a coarse‑grained (CG) protein model that explicitly incorporates Debye‑Hückel electrostatics, with residue charges derived from quantum‑chemical calculations for both apo‑ and holo‑CaM. The CG simulations were coupled to the MultiSCAAL reconstruction pipeline, allowing generation of all‑atom ensembles from CG trajectories. Free‑energy landscapes were projected onto two order parameters: the overlap function χ (similarity to the crystal structure) and asphericity Δ (shape descriptor). At low ionic strength, two basins emerged: an extended, dumbbell‑like state (M1) resembling the crystal apo‑CaM, and a compact, spherical state (M2). Raising the salt concentration screened charge–charge attractions, shifting the population toward M1 and increasing overall helicity by ~2–3 %.

Temperature‑dependent radius of gyration (R_g) analyses showed that high salt reduces R_g, reflecting a more compact unfolded ensemble. Simulated CD spectra, constructed from the average helical and β‑strand contents, reproduced the experimental trend of increasing negative ellipticity with both salt and crowding. Importantly, the simulations revealed that crowding and ionic strength not only affect global compaction but also reorient the EF‑hand domains, exposing hydrophobic residues such as Met that are critical for target recognition.

Overall, the study demonstrates that intracellular‑like conditions—high macromolecular crowding, physiological ionic strength, and calcium binding—cooperatively reshape CaM’s conformational ensemble. These environmental cues increase helicity, promote a compact yet flexible architecture, and alter EF‑hand orientation, thereby modulating calcium affinity and the promiscuous selection of target proteins. The integrated experimental‑computational approach provides a quantitative framework for understanding how the crowded, charged cellular interior governs the dynamics of intrinsically flexible signaling proteins like calmodulin.