The role of Phe-92 in the Ca(2+)-induced conformational transition in the C-terminal domain of calmodulin

Linear Response Path Following: A Molecular Dynamics Method To Simulate Global Conformational Changes of Protein upon Ligand Binding

Molecular functions of proteins are often fulfilled by global conformational changes that couple with local events such as the binding of ligand molecules. High molecular complexity of proteins has, however, been an obstacle to obtain an atomistic view of the global conformational transitions, imposing a limitation on the mechanistic understanding of the functional processes. In this study, we developed a new method of molecular dynamics (MD) simulation called the linear response path following (LRPF) to simulate a protein's global conformational changes upon ligand binding. The method introduces a biasing force based on a linear response theory, which determines a local reaction coordinate in the configuration space that represents linear coupling between local events of ligand binding and global conformational changes and thus provides one with fully atomistic models undergoing large conformational changes without knowledge of a target structure. The overall transition process involving nonlinear conformational changes is simulated through iterative cycles consisting of a biased MD simulation with an updated linear response force and a following unbiased MD simulation for relaxation. We applied the method to the simulation of global conformational changes of the yeast calmodulin N-terminal domain and successfully searched out the end conformation. The atomistically detailed trajectories revealed a sequence of molecular events that properly lead to the global conformational changes and identified key steps of local-global coupling that induce the conformational transitions. The LRPF method provides one with a powerful means to model conformational changes of proteins such as motors and transporters where local-global coupling plays a pivotal role in their functional processes.

Modulation of calmodulin lobes by different targets: an allosteric model with hemiconcerted conformational transitions

Calmodulin is a calcium-binding protein ubiquitous in eukaryotic cells, involved in numerous calcium-regulated biological phenomena, such as synaptic plasticity, muscle contraction, cell cycle, and circadian rhythms. It exibits a characteristic dumbell shape, with two globular domains (N- and C-terminal lobe) joined by a linker region. Each lobe can take alternative conformations, affected by the binding of calcium and target proteins. Calmodulin displays considerable functional flexibility due to its capability to bind different targets, often in a tissue-specific fashion. In various specific physiological environments (e.g. skeletal muscle, neuron dendritic spines) several targets compete for the same calmodulin pool, regulating its availability and affinity for calcium. In this work, we sought to understand the general principles underlying calmodulin modulation by different target proteins, and to account for simultaneous effects of multiple competing targets, thus enabling a more realistic simulation of calmodulin-dependent pathways. We built a mechanistic allosteric model of calmodulin, based on an hemiconcerted framework: each calmodulin lobe can exist in two conformations in thermodynamic equilibrium, with different affinities for calcium and different affinities for each target. Each lobe was allowed to switch conformation on its own. The model was parameterised and validated against experimental data from the literature. In spite of its simplicity, a two-state allosteric model was able to satisfactorily represent several sets of experiments, in particular the binding of calcium on intact and truncated calmodulin and the effect of different skMLCK peptides on calmodulin's saturation curve. The model can also be readily extended to include multiple targets. We show that some targets stabilise the low calcium affinity T state while others stabilise the high affinity R state. Most of the effects produced by calmodulin targets can be explained as modulation of a pre-existing dynamic equilibrium between different conformations of calmodulin's lobes, in agreement with linkage theory and MWC-type models.

X-ray structures of magnesium and manganese complexes with the N-terminal domain of calmodulin: insights into the mechanism and specificity of metal ion binding to an EF-hand

Calmodulin (CaM), a member of the EF-hand superfamily, regulates many aspects of cell function by responding specifically to micromolar concentrations of Ca(2+) in the presence of an ~1000-fold higher concentration of cellular Mg(2+). To explain the structural basis of metal ion binding specificity, we have determined the X-ray structures of the N-terminal domain of calmodulin (N-CaM) in complexes with Mg(2+), Mn(2+), and Zn(2+). In contrast to Ca(2+), which induces domain opening in CaM, octahedrally coordinated Mg(2+) and Mn(2+) stabilize the closed-domain, apo-like conformation, while tetrahedrally coordinated Zn(2+) ions bind at the protein surface and do not compete with Ca(2+). The relative positions of bound Mg(2+) and Mn(2+) within the EF-hand loops are similar to those of Ca(2+); however, the Glu side chain at position 12 of the loop, whose bidentate interaction with Ca(2+) is critical for domain opening, does not bind directly to either Mn(2+) or Mg(2+), and the vacant ligand position is occupied by a water molecule. We conclude that this critical interaction is prevented by specific stereochemical constraints imposed on the ligands by the EF-hand β-scaffold. The structures suggest that Mg(2+) contributes to the switching off of calmodulin activity and possibly other EF-hand proteins at the resting levels of Ca(2+). The Mg(2+)-bound N-CaM structure also provides a unique view of a transiently bound hydrated metal ion and suggests a role for the hydration water in the metal-induced conformational change.

Approaches to the assignment of (19)F resonances from 3-fluorophenylalanine labeled calmodulin using solution state NMR

Traditional single site replacement mutations (in this case, phenylalanine to tyrosine) were compared with methods which exclusively employ (15)N and (19)F-edited two- and three-dimensional NMR experiments for purposes of assigning (19)F NMR resonances from calmodulin (CaM), biosynthetically labeled with 3-fluorophenylalanine (3-FPhe). The global substitution of 3-FPhe for native phenylalanine was tolerated in CaM as evidenced by a comparison of (1)H-(15)N HSQC spectra and calcium binding assays in the presence and absence of 3-FPhe. The (19)F NMR spectrum reveals six resolved resonances, one of which integrates to three 3-FPhe species, making for a total of eight fluorophenylalanines. Single phenylalanine to tyrosine mutants of five phenylalanine positions resulted in (19)F NMR spectra with significant chemical shift perturbations of the remaining resonances, and provided only a single definitive assignment. Although (1)H-(19)F heteronucleclear NOEs proved weak, (19)F-edited (1)H-(1)H NOESY connectivities were relatively easy to establish by making use of the (3)J(FH) coupling between the fluorine nucleus and the adjacent fluorophenylalanine delta proton. (19)F-edited NOESY connectivities between the delta protons and alpha and beta nuclei in addition to (15)N-edited (1)H, (1)H NOESY crosspeaks proved sufficient to assign 4 of 8 (19)F resonances. Controlled cleavage of the protein into two fragments using trypsin, and a repetition of the above 2D and 3D techniques resulted in unambiguous assignments of all 8 (19)F NMR resonances. Our studies suggest that (19)F-edited NOESY NMR spectra are generally adequate for complete assignment without the need to resort to mutational analysis.

A Carboxyl-terminal Hydrophobic Interface Is Critical to Sodium Channel Function

Perturbation of sodium channel inactivation, a finely tuned process that critically regulates the flow of sodium ions into excitable cells, is a common functional consequence of inherited mutations associated with epilepsy, skeletal muscle disease, autism, and cardiac arrhythmias. Understanding the structural basis of inactivation is key to understanding these disorders. Here we identify a novel role for a structural motif in the COOH terminus of the heart NaV1.5 sodium channel in determining channel inactivation. Structural modeling predicts an interhelical hydrophobic interface between paired EF hands in the proximal region of the NaV1.5 COOH terminus. The predicted interface is conserved among almost all EF hand-containing proteins and is the locus of a number of disease-associated mutations. Using the structural model as a guide, we provide biochemical and biophysical evidence that the structural integrity of this interface is necessary for proper Na+ channel inactivation gating. We thus demonstrate a novel role of the sodium channel COOH terminus structure in the control of channel inactivation and in pathologies caused by inherited mutations that disrupt it.

Calcium dependence of aequorin bioluminescence dissected by random mutagenesis

Aequorin bioluminescence is emitted as a rapidly decaying flash upon calcium binding. Random mutagenesis and functional screening were used to isolate aequorin mutants showing slow decay rate of luminescence. Calcium sensitivity curves were shifted in all mutants, and an intrinsic link between calcium sensitivity and decay rate was suggested by the position of all mutations in or near EF-hand calcium-binding sites. From these results, a low calcium affinity was assigned to the N-terminal EF hand and a high affinity to the C-terminal EF-hand pair. In WT aequorin, the increase of the decay rate with calcium occurred at constant total photon yield and thus determined a corresponding increase of light intensity. Increase of the decay rate was underlain by variations of a fast and a slow component and required the contribution of all three EF hands. Conversely, analyses of double EF-hand mutants suggested that single EF hands are sufficient to trigger luminescence at a slow rate. Finally, a model postulating that proportions of a fast and a slow light-emitting state depend on calcium concentration adequately described the calcium dependence of aequorin bioluminescence. Our results suggest that variations of luminescence kinetics, which depend on three EF hands endowed with different calcium affinities, critically determine the amplitude of aequorin responses to biological calcium signals.

Structure of a Trapped Intermediate of Calmodulin: Calcium Regulation of EF-hand Proteins from a New Perspective

Calmodulin (CaM) is a multifunctional Ca2+-binding protein that regulates the activity of many enzymes in response to changes in the intracellular Ca2+ concentration. There are two globular domains in CaM, each containing a pair of helix-loop-helix Ca2+-binding motifs called EF-hands. Ca2+-binding induces the opening of both domains thereby exposing hydrophobic pockets that provide binding sites for the target enzymes. Here, I present a 2.4 A resolution structure of a calmodulin mutant (CaM41/75) in which the N-terminal domain is locked in the closed conformation by a disulfide bond. CaM41/75 crystallized in a tetragonal lattice with the Ca2+ bound in all four EF-hands. In the closed N-terminal domain Ca ions are coordinated by the four protein ligands in positions 1, 3, 5 and 7 of the loop, and by two solvent ligands. The glutamate side-chain in the 12th position of the loop (Glu31 in site I and Glu67 in site II), which in the wild-type protein provides a bidentate Ca2+ ligand, remains in a distal position. Based on a comparison of CaM41/75 with other CaM and troponin C structures a detailed two-step mechanism of the Ca2+-binding process is proposed. Initially, the Ca2+ binds to the N-terminal part of the loop, thus generating a rigid link between the incoming helix (helix A, or helix C) and the central beta structure involving the residues in the sixth, seventh and eighth position of the loop. Then, the exiting helix (helix B or helix D) rotates causing the glutamate ligand in the 12th position to move into the vicinity of the immobilized Ca2+. An adjustment of the phi, psi backbone dihedral angles of the Ile residue in the eighth position is necessary and sufficient for the helix rotation and functions as a hinge. The model allows for a significant independence of the Ca2+-binding sites in a two-EF-hand domain.

Ca2+binding sites in calmodulin and troponin C alter interhelical angle movements

Molecular dynamics analyses were performed to examine conformational changes in the C-domain of calmodulin and the N-domain of troponin C induced by binding of Ca(2+) ions. Analyses of conformational changes in calmodulin and troponin C indicated that the shortening of the distance between Ca(2+) ions and Ca(2+) binding sites of helices caused widening of the distance between Ca(2+) binding sites of helices on opposite sides, while the hydrophobic side chains in the center of helices hardly moved due to their steric hindrance. This conformational change acts as the clothespin mechanism.

Structural preference for changes in the direction of the Ca2+-induced transition: a study of the regulatory domain of skeletal troponin-C

The determinants for specificity in the Ca(2+)-dependent response of the regulatory N-terminal domain of skeletal troponin-C are a combination of intrinsic and induced properties. We characterized computationally the intrinsic propensity of this domain for structural changes similar to those observed experimentally in the Ca(2+)-induced transition. The preference for such changes was assessed by comparing the structural effect of the harmonic and quasiharmonic vibrations specific for each Ca(2+) occupancy with crystallographic data. Results show that only the Ca(2+)-saturated form of the protein features a slow vibrational motion preparatory for the transition. From the characteristics of this mode, we identified a molecular mechanism for transition, by which residues 42-51 of helix B and of the adjacent linker move toward helices (A, D), and bind to the surface used by the protein to interact with troponin-I. By obstructing the access of the target to hydrophobic residues important in the formation of the complex, helix B and the adjacent linker act as an autoinhibitory structural element. Specific properties of the methionines at the interaction surface were found to favor the binding of the autoinhibitory region. Located over hydrophobic residues critical for binding, the methionines are easily displaceable to increase the accessibility of these residues to molecular encounter.

Phenylalanine fluorescence studies of calcium binding to N-domain fragments of Paramecium calmodulin mutants show increased calcium affinity correlates with increased disorder

Calmodulin (CaM) is a ubiquitous, essential calcium-binding protein that regulates diverse protein targets in response to physiological calcium fluctuations. Most high-resolution structures of CaM-target complexes indicate that the two homologous domains of CaM are equivalent partners in target recognition. However, mutations between calcium-binding sites I and II in the N-domain of Paramecium calmodulin (PCaM) selectively affect calcium-dependent sodium currents. To understand these domain-specific effects, N-domain fragments (PCaM(1-75)) of six of these mutants were examined to determine whether energetics of calcium binding to sites I and II or conformational properties had been perturbed. These PCaM((1-75)) sequences naturally contain 5 Phe residues but no Tyr or Trp; calcium binding was monitored by observing the reduction in intrinsic phenylalanine fluorescence at 280 nm. To assess mutation-induced conformational changes, thermal denaturation of the apo PCaM((1-75)) sequences, and calcium-dependent changes in Stokes radii were determined. The free energy of calcium binding to each mutant was within 1 kcal/mole of the value for wild type and calcium reduced the R(s) of all of them. A striking trend was observed whereby mutants showing an increase in calcium affinity and R(s) had a concomitant decrease in thermal stability (by as much as 18 degrees C). Thus, mutations between the binding sites that increased disorder and reduced tertiary constraints in the apo state promoted calcium coordination. This finding underscores the complexity of the linkage between calcium binding and conformational change and the difficulty in predicting mutational effects.

A novel Ca2+ binding protein associated with caldesmon in Ca2+-regulated smooth muscle thin filaments: evidence for a structurally altered form of calmodulin

Smooth muscle thin filaments are made up of actin, tropomyosin, the inhibitory protein caldesmon and a Ca2+-binding protein. Thin filament activation of myosin MgATPase is Ca2+-regulated but thin filaments assembled from smooth muscle actin, tropomyosin and caldesmon plus brain or aorta calmodulin are not Ca2+-regulated at 25 degrees C/50 mM KCl. We isolated the Ca2+-binding protein (CaBP) from smooth muscle thin filaments by DEAE fast-flow chromatography in 6 M urea and phenyl sepharose chromatography using sheep aorta as our starting material. CaBP combines with smooth muscle actin, tropomyosin and caldesmon to reconstitute a normally regulated thin filament at 25 degrees C/50 mM KCl. It reverses caldesmon inhibition at pCa5 under conditions where CaM is largely inactive, it binds to caldesmon when complexed with actin and tropomyosin rather than displacing it and it binds to caldesmon independently of [Ca2+]. Amino acid sequencing, and electrospray mass spectrometry show the CaBP is identical to CaM. Structural probes indicate it is different: calmodulin increases caldesmon tryptophan fluorescence but CaBP does not. The distribution of charged species in electrospray mass spectrometry and nozzle skimmer fragmentation patterns are different indicating a less stable N-terminal lobe for CaBP. Brief heating abolishes these special properties of the CaBP. Mass spectrometry in aqueous buffer showed no evidence for the presence of any covalent or non-covalently bound adduct. The only remaining conclusion is that CaBP is calmodulin locked in a metastable altered state.

Oxidatively modified calmodulin binds to the plasma membrane Ca-ATPase in a nonproductive and conformationally disordered complex

Oxidation of either Met(145) or Met(146) in wheat germ calmodulin (CaM) to methionine sulfoxide prevents the CaM-dependent activation of the plasma membrane (PM) Ca-ATPase (D. Yin, K. Kuczera, and T. C. Squier, 2000, Chem. Res. Toxicol. 13:103-110). To investigate the structural basis for the inhibition of the PM-Ca-ATPase by oxidized CaM (CaM(ox)), we have used circular dichroism (CD) and fluorescence spectroscopy to resolve conformational differences within the complex between CaM and the PM-Ca-ATPase. The similar excited-state lifetime and solvent accessibility of the fluorophore N-1-pyrenyl-maleimide covalently bound to Cys(26) in unoxidized CaM and CaM(ox) indicates that the globular domains within CaM(ox) assume a native-like structure following association with the PM-Ca-ATPase. However, in comparison with oxidized CaM there are increases in the 1) molar ellipticity in the CD spectrum and 2) conformational heterogeneity between the opposing globular domains for CaM(ox) bound to the CaM-binding sequence of the PM-Ca-ATPase. Furthermore, CaM(ox) binds to the PM-Ca-ATPase with high affinity at a distinct, but overlapping, site to that normally occupied by unoxidized CaM. These results suggest that alterations in binding interactions between CaM(ox) and the PM-Ca-ATPase block important structural transitions within the CaM-binding sequence of the PM-Ca-ATPase that are normally associated with enzyme activation.

Solvation energetics and conformational change in EF-hand proteins

Calmodulin and other members of the EF-hand protein family are known to undergo major changes in conformation upon binding Ca(2+). However, some EF-hand proteins, such as calbindin D9k, bind Ca(2+) without a significant change in conformation. Here, we show the importance of a precise balance of solvation energetics to conformational change, using mutational analysis of partially buried polar groups in the N-terminal domain of calmodulin (N-cam). Several variants were characterized using fluorescence, circular dichroism, and NMR spectroscopy. Strikingly, the replacement of polar side chains glutamine and lysine at positions 41 and 75 with nonpolar side chains leads to dramatic enhancement of the stability of the Ca(2+)-free state, a corresponding decrease in Ca(2+)-binding affinity, and an apparent loss of ability to change conformation to the open form. The results suggest a paradigm for conformational change in which energetic strain is accumulated in one state in order to modulate the energetics of change to the alternative state.

Importance of phenylalanine residues of yeast calmodulin for target binding and activation

Recent genetic studies of yeast calmodulin (yCaM) have shown that alterations of different sets of Phe residues result in distinct functional defects (Ohya, Y., and Botstein, D. (1994) Science 263, 963-966). To examine the importance of Phe residues for target binding and activation, we purified mutant yCaMs containing single or double Phe to Ala substitutions and determined their ability to bind and activate two target proteins, calcineurin and CaM-dependent protein kinase (CaMK). Binding assays using the gel overlay technique and quantitative analyses using surface plasmon resonance measurements indicated that the binding of yCaM to calcineurin is impaired by either double mutations of F16A/F19A or a single mutation of F140A, while binding to CaMK is impaired by F89A, F92A, or F140A. These same mutant yCaMs fail to activate calcineurin and CaMK, respectively, in vitro. In addition, F19A exhibited a severe defect in activation of both enzymes. F12A activated calcineurin to only 50% of the level achieved by wild-type calmodulin but fully activated CaMK. These results suggest that each target protein requires a specific and distinct subset of Phe residues in yCaM for target binding and activation.