Dynamics of Ca2+-saturated calmodulin D129N mutant studied by multiple molecular dynamics simulations

Role of water-bridged interactions in metal ion coupled protein allostery

Allosteric communication between distant parts of proteins controls many cellular functions, in which metal ions are widely utilized as effectors to trigger the allosteric cascade. Due to the involvement of strong coordination interactions, the energy landscape dictating the metal ion binding is intrinsically rugged. How metal ions achieve fast binding by overcoming the landscape ruggedness and thereby efficiently mediate protein allostery is elusive. By performing molecular dynamics simulations for the Ca²⁺ binding mediated allostery of the calmodulin (CaM) domains, each containing two Ca²⁺ binding helix-loop-helix motifs (EF-hands), we revealed the key role of water-bridged interactions in Ca²⁺ binding and protein allostery. The bridging water molecules between Ca²⁺ and binding residue reduces the ruggedness of ligand exchange landscape by acting as a lubricant, facilitating the Ca²⁺ coupled protein allostery. Calcium-induced rotation of the helices in the EF-hands, with the hydrophobic core serving as the pivot, leads to exposure of hydrophobic sites for target binding. Intriguingly, despite being structurally similar, the response of the two symmetrically arranged EF-hands upon Ca²⁺ binding is asymmetric. Breakage of symmetry is needed for efficient allosteric communication between the EF-hands. The key roles that water molecules play in driving allosteric transitions are likely to be general in other metal ion mediated protein allostery.

Natural proteins often utilize allostery in executing a variety of functions. Metal ions are typical cofactors to trigger the allosteric cascade. In this work, using the Ca²⁺ sensor protein calmodulin as the model system, we revealed crucial roles of water-bridged interactions in the metal ion coupled protein allostery. The coordination of the Ca²⁺ to the binding site involves an intermediate in which the water molecule bridges the Ca²⁺ and the liganding residue. The bridging water reduces the free energy barrier height of ligand exchange, therefore facilitating the ligand exchange and allosteric coupling by acting as a lubricant. We also showed that the response of the two symmetrically arranged EF-hand motifs of CaM domains upon Ca²⁺ binding is asymmetric, which is directly attributed to the differing dehydration process of the Ca²⁺ ions and is needed for efficient allosteric communication.

The mechanism of complex formation between calmodulin and voltage gated calcium channels revealed by molecular dynamics

Calmodulin, a ubiquitous eukaryotic calcium sensor responsible for the regulation of many fundamental cellular processes, is a highly flexible protein and exhibits an unusually wide range of conformations. Furthermore, CaM is known to interact with more than 300 cellular targets. Molecular dynamics (MD) simulation trajectories suggest that EF-hand loops show different magnitudes of flexibility. Therefore, the four EF-hand motifs have different affinities for Ca2+ ions, which enables CaM to function on wide range of Ca2+ ion concentrations. EF-hand loops are 2-3 times more flexible in apo CaM whereas least flexible in Ca2+/CaM-IQ motif complexes. We report a unique intermediate conformation of Ca2+/CaM while transitioning from extended to compact form. We also report the complex formation process between Ca2+/CaM and IQ CaM-binding motifs. Our results showed how IQ motif recognise its binding site on the CaM and how CaM transforms from extended to compact form upon binding to IQ motif.

Inherent flexibility determines the transition mechanisms of the EF-hands of calmodulin

We explore how inherent flexibility of a protein molecule influences the mechanism controlling allosteric transitions by using a variational model inspired from work in protein folding. The striking differences in the predicted transition mechanism for the opening of the two domains of calmodulin (CaM) emphasize that inherent flexibility is key to understanding the complex conformational changes that occur in proteins. In particular, the C-terminal domain of CaM (cCaM), which is inherently less flexible than its N-terminal domain (nCaM), reveals "cracking" or local partial unfolding during the open/closed transition. This result is in harmony with the picture that cracking relieves local stresses caused by conformational deformations of a sufficiently rigid protein. We also compare the conformational transition in a recently studied even-odd paired fragment of CaM. Our results rationalize the different relative binding affinities of the EF-hands in the engineered fragment compared with the intact odd-even paired EF-hands (nCaM and cCaM) in terms of changes in flexibility along the transition route. Aside from elucidating general theoretical ideas about the cracking mechanism, these studies also emphasize how the remarkable intrinsic plasticity of CaM underlies conformational dynamics essential for its diverse functions.

An altered mode of calcium coordination in methionine-oxidized calmodulin

Oxidation of methionine residues in calmodulin (CaM) lowers the affinity for calcium and results in an inability to activate target proteins fully. To evaluate the structural consequences of CaM oxidation, we used infrared difference spectroscopy to identify oxidation-dependent effects on protein conformation and calcium liganding. Oxidation-induced changes include an increase in hydration of alpha-helices, as indicated in the downshift of the amide I' band of both apo-CaM and Ca(2+)-CaM, and a modification of calcium liganding by carboxylate side chains, reflected in antisymmetric carboxylate band shifts. Changes in carboxylate ligands are consistent with the model we propose: an Asp at position 1 of the EF-loop experiences diminished hydrogen bonding with the polypeptide backbone, an Asp at position 3 forms a bidentate coordination of calcium, and an Asp at position 5 forms a pseudobridging coordination with a calcium-bound water molecule. The bidentate coordination of calcium by conserved glutamates is unaffected by oxidation. The observed changes in calcium ligation are discussed in terms of the placement of methionine side chains relative to the calcium-binding sites, suggesting that varying sensitivities of binding sites to oxidation may underlie the loss of CaM function upon oxidation.

Structures and metal-ion-binding properties of the Ca2+-binding helix–loop–helix EF-hand motifs

The ‘EF-hand’ Ca2+-binding motif plays an essential role in eukaryotic cellular signalling, and the proteins containing this motif constitute a large and functionally diverse family. The EF-hand is defined by its helix–loop–helix secondary structure as well as the ligands presented by the loop to bind the Ca2+ ion. The identity of these ligands is semi-conserved in the most common (the ‘canonical’) EF-hand; however, several non-canonical EF-hands exist that bind Ca2+ by a different co-ordination mechanism. EF-hands tend to occur in pairs, which form a discrete domain so that most family members have two, four or six EF-hands. This pairing also enables communication, and many EF-hands display positive co-operativity, thereby minimizing the Ca2+ signal required to reach protein saturation. The conformational effects of Ca2+ binding are varied, function-dependent and, in some cases, minimal, but can lead to the creation of a protein target interaction site or structure formation from a molten-globule apo state. EF-hand proteins exhibit various sensitivities to Ca2+, reflecting the intrinsic binding ability of the EF-hand as well as the degree of co-operativity in Ca2+ binding to paired EF-hands. Two additional factors can influence the ability of an EF-hand to bind Ca2+: selectivity over Mg2+ (a cation with very similar chemical properties to Ca2+ and with a cytoplasmic concentration several orders of magnitude higher) and interaction with a protein target. A structural approach is used in this review to examine the diversity of family members, and a biophysical perspective provides insight into the ability of the EF-hand motif to bind Ca2+ with a wide range of affinities.

An Approximate Method in Using Molecular Mechanics Simulations To Study Slow Protein Conformational Changes

The broad range of characteristic motions in proteins has limited the applicability of molecular dynamics simulations in studying large-scale conformational transitions. We present an approximate method, making use of standard MD simulations and using a much larger integration time step, to obtain the structural changes for slow systematic motions of large complex systems. We show the applicability of this method by simulating the open to closed Calmodulin calcium binding domain conformational changes. Starting with the Ca2+-bound X-ray structure, and after the removal of the Ca2+ ions, our calculation yielded intermediate conformations during the rearrangement of helices in each Ca2+ binding pocket, leading to a structure with a lowest rmsd of 1.56 A compared to the NMR apo-calmodulin structure.

Observation of microsecond time-scale protein dynamics in the presence of Ln3+ ions: application to the N-terminal domain of cardiac troponin C

The microsecond time-scale motions in the N-terminal domain of cardiac troponin C (NcTnC) loaded with lanthanide ions have been investigated by means of a (1)H(N) off-resonance spin-lock experiment. The observed relaxation dispersion effects strongly increase along the series of NcTnC samples containing La(3+), Ce(3+), and Pr(3+) ions. This rise in dispersion effects is due to modulation of long-range pseudocontact shifts by micros time-scale dynamics. Specifically, the motion in the coordination sphere of the lanthanide ion (i.e. in the NcTnC EF-hand motif) causes modulation of the paramagnetic susceptibility tensor which, in turn, causes modulation of pseudocontact shifts. It is also probable that opening/closing dynamics, previously identified in Ca(2+)-NcTnC, contributes to some of the observed dispersions. On the other hand, it is unlikely that monomer-dimer exchange in the solution of NcTnC is directly responsible for the dispersion effects. Finally, on-off exchange of the lanthanide ion does not seem to play any significant role. The amplification of dispersion effects by Ln(3+) ions is a potentially useful tool for studies of micros-ms motions in proteins. This approach makes it possible to observe the dispersions even when the local environment of the reporting spin does not change. This happens, for example, when the motion involves a 'rigid' structural unit such as individual alpha-helix. Even more significantly, the dispersions based on pseudocontact shifts offer better chances for structural characterization of the dynamic species. This method can be generalized for a large class of applications via the use of specially designed lanthanide-binding tags.

A statistical approach to the interpretation of molecular dynamics simulations of calmodulin equilibrium dynamics

A sample of 35 independent molecular dynamics (MD) simulations of calmodulin (CaM) equilibrium dynamics was prepared from different but equally plausible initial conditions (20 simulations of the wild-type protein and 15 simulations of the D129N mutant). CaM's radius of gyration and backbone mean-square fluctuations were analyzed for the effect of the D129N mutation, and simulations were compared with experiments. Statistical tests were employed for quantitative comparisons at the desired error level. The computational model predicted statistically significant compaction of CaM relative to the crystal structure, consistent with the results of small-angle X-ray scattering (SAXS) experiments. This effect was not observed in several previously reported studies of (Ca2+)(4)-CaM, which relied on a single MD run. In contrast to radius of gyration, backbone mean-square fluctuations showed a distinctly non-normal and positively skewed distribution for nearly all residues. Furthermore, the D129N mutation affected the backbone dynamics in a complex manner and reduced the mobility of Glu123, Met124, Ile125, Arg126, and Glu127 located in the adjacent alpha-helix G. The implications of these observations for the comparisons of MD simulations with experiments are discussed. The proposed approach may be useful in studies of protein equilibrium dynamics where MD simulations fall short of properly sampling the conformational space, and when the comparison with experiments is affected by the reproducibility of the computational model.

Comparative molecular dynamics-Similar folds and similar motions?

Proteins possessing the same fold may undergo similar motions, particularly if these motions involve large conformational transitions. The increasing amounts of structural data provide a useful starting point with which to test this hypothesis. We have performed a total of 0.29 micros of molecular dynamics across a series of proteins within the same fold family (periplasmic binding proteinlike) in order to address to what extent similarity of motion exists. Analysis of the local conformational space on these timescales (10-20 ns) revealed that the behavior of the proteins could be readily distinguished between an apo-state and a ligand-bound state. Moreover, analysis of the root-mean-square fluctuations reveals that the presence of the ligand exerts a stabilizing effect on the protein, with similar motions occurring, but with reduced magnitude. Furthermore, the conformational space in the presence of the ligand appears to be dictated by sequence but not by the type of ligand present. In contrast, apo-simulations showed considerable overlap of conformational space across the fold as a result of their ability to undergo larger fluctuations. Indeed, we observed several transitions from different simulations between states corresponding to the closed-cleft and open-cleft forms of the fold, with the predominant motions being conserved across the different proteins. Thus, large-scale conformational changes do indeed appear to be conserved across this fold architecture, but smaller conformational motions appear to reflect the differences in sequence and local fold.

Unwinding the helical linker of calcium-loaded calmodulin: A molecular dynamics study

The fold of calmodulin (CaM) consists of two globular domains connected by a helical segment (the linker), whose conformational properties play a crucial role for the protein's molecular recognition processes. Here we investigate the structural properties of the linker by performing a 11.5 ns molecular dynamics (MD) simulation of calcium-loaded human CaM in aqueous solution. The calculations are based on the AMBER force field. The calculated S2 order parameters are in good accord with NMR data: The structure of the linker in our simulations is much more flexible than that emerging from the Homo sapiens X-ray structure, consistently with the helix unwinding observed experimentally in solution. This process occurs spontaneously in a nanosecond timescale, as observed also in a very recent simulation based on the GROMOS force field. A detailed description of the mechanism that determines the linker unwinding is provided, in which electrostatic contacts between the two globular domains play a critical role. The orientation of the domains emerging from our MD calculations is consistent both with former X-ray scattering data and a recent NMR work. Based on our findings, a rationale for the experimentally measured entropy cost associated to binding to the protein's cellular partners is also given.