Subtle pH differences trigger single residue motions for moderating conformations of calmodulin

This study reveals the essence of ligand recognition mechanisms by which calmodulin (CaM) controls a variety of Ca2+ signaling processes. We study eight forms of calcium-loaded CaM each with distinct conformational states. Reducing the structure to two degrees of freedom conveniently describes main features of conformational changes of CaM via simultaneous twist-bend motions of the two lobes. We utilize perturbation-response scanning (PRS) technique, coupled with molecular dynamics simulations to analyze conformational preferences of calcium-loaded CaM, initially in extended form. PRS is comprised of sequential application of directed forces on residues followed by recording the resulting coordinates. We show that manipulation of a single residue, E31 located in one of the EF hand motifs, reproduces structural changes to compact forms, and the flexible linker acts as a transducer of binding information to distant parts of the protein. Independently, using four different pKa calculation strategies, we find E31 to be the charged residue (out of 52), whose ionization state is most sensitive to subtle pH variations in the physiological range. It is proposed that at relatively low pH, CaM structure is less flexible. By gaining charged states at specific sites at a pH value around 7, local conformational changes in the protein will lead to shifts in the energy landscape, paving the way to other conformational states. These findings are in accordance with FRET measured shifts in conformational distributions towards more compact forms with decreased pH. They also corroborate mutational studies and proteolysis results which point to the significant role of E31 in CaM dynamics.

💡 Research Summary

This paper investigates how subtle variations in pH can trigger large‑scale conformational changes in calmodulin (CaM), a ubiquitous calcium‑binding messenger protein. The authors examine eight calcium‑loaded CaM structures: one extended apo‑like form and seven ligand‑bound forms that differ in ligand size (from a few atoms to a 26‑residue peptide) and binding location. By aligning these structures and computing overall and domain‑specific RMSDs, they find that the primary structural differences arise from the relative orientation of the N‑ and C‑lobes rather than from internal rearrangements within each lobe. Most pairs show overall RMSDs of 15–16 Å, while intra‑lobe RMSDs remain below 1.5 Å, indicating a global “twist‑bend” motion of the two lobes.

To uncover the mechanistic basis of these motions, the authors apply Perturbation‑Response Scanning (PRS), a linear‑response technique that sequentially applies a small directed force to each residue in a molecular dynamics (MD) trajectory and records the resulting displacement of all atoms. The covariance matrix derived from the last 90 ns of a 100 ns MD simulation serves as the response operator. For each residue, the overlap (Oi) between the PRS‑generated response vector (ΔRi) and the experimentally observed conformational change vector (ΔS) is calculated. Random perturbations on single residues yield modest average overlaps (~0.56), but a striking exception is found for residue Glu‑31 (E31) located in EF‑hand I. When a force is applied to E31 in a specific direction, the overlap rises to ~0.70 ± 0.03 for five of the seven target structures, effectively reproducing the experimentally observed transition from the extended to the compact conformations. Neighboring residues such as L69 contribute modestly in a few cases, but combined perturbations of E31 with any other residue do not improve the overlap, confirming the singular importance of E31.

The authors then explore why E31 is uniquely effective. Using four independent pKa prediction methods (PROPKA, H++, MCCE, DelPhi), they calculate the pKa values of all 52 ionizable residues in CaM. Only E31 shows a pronounced shift in its protonation state within the physiological pH range (≈6.5–7.5). At lower pH, E31 is largely neutral, reducing electrostatic repulsion and rendering the protein more rigid; the extended conformation is favored. Near neutral pH, E31 becomes negatively charged, altering the electrostatic landscape of the flexible linker that connects the two lobes. This charge change promotes attractive interactions that pull the lobes together, facilitating the twist‑bend motion toward the compact state. The computational findings align with fluorescence resonance energy transfer (FRET) experiments that report a shift toward more compact CaM conformations as pH decreases, as well as with mutagenesis and proteolysis studies that highlight the functional relevance of E31.

To further simplify the picture, the authors reduce CaM to a two‑degree‑of‑freedom model representing the twist and bend of the lobes. Even in this coarse‑grained representation, a single perturbation at the E31 position reproduces the global conformational shift, underscoring that a localized electrostatic perturbation can be transduced through the flexible linker to drive a large‑scale structural rearrangement.

The study concludes that (1) CaM’s conformational landscape is dominated by a low‑dimensional twist‑bend motion of its two lobes; (2) PRS combined with MD provides a powerful framework for identifying single‑residue “allosteric hotspots” that can control such motions; and (3) the pH‑sensitive protonation of E31 acts as a molecular switch that modulates CaM’s flexibility and its ability to adopt multiple functional conformations. These insights deepen our molecular understanding of how subtle environmental cues, such as pH fluctuations, can regulate calcium signaling pathways through calmodulin.