Long-range effects on calcium binding and conformational change in the N-domain of calmodulin

Proteins within the EF-hand protein family exhibit different conformational responses to Ca(2+) binding. Calmodulin and other members of the EF-hand protein family undergo major changes in conformation upon binding Ca(2+). However, some EF-hand proteins, such as calbindin D9k (Clb), bind Ca(2+) without a significant change in conformation. Here, we investigate the effects of replacement of a leucine at position 39 of the N-terminal domain of calmodulin (N-Cam) with a phenylalanine derived from Clb. This variant is studied alone and in the context of other mutations that affect the conformational properties of N-Cam. Strikingly, the introduction of Phe39, which is distant from the calcium binding sites, leads to a significant enhancement of Ca(2+) binding affinity, even in the context of other mutations which trap the protein in the closed form. The results yield novel insights into the evolution of EF-hand proteins as calcium sensors versus calcium buffers.     

Fluorescence investigations of calmodulin hydrophobic sites

Calmodulin activation of target enzymes depends on the interaction between calmodulin hydrophobic regions and some enzyme areas. The Ca2+ induced exposure of calmodulin hydrophobic sites was studied by means of 2-p-toluidinylnaphthalene-6-sulfonate, a fluorescent probe. Scatchard and Job plots showed that the calmodulin-Ca42+ complex bound two molecules of this hydrophobic probe, with KD congruent to 1.4 X 10(-4) M. These sites are not totally exposed until calmodulin has bound four Ca2+ per molecule, so the conformational change is not over before the four specific Ca2+ - binding sites are saturated with Ca2+.     

An interdomain linker increases the thermostability and decreases the calcium affinity of the calmodulin N-domain.

A hydrophobic core is a widely accepted determinant of protein stability. However, regulatory proteins undergoing ligand-induced conformational switching may expose interior residues to solvent and cannot afford to be extremely rigid. Optimizing the energetic balance between stability and binding is challenging. The addition of five interdomain residues to rat and Paramecium calmodulin N-domain fragments (residues 1-75) increased their thermostability by 9 degrees C and lowered their calcium affinity by a factor of 4. This demonstrates that the flexible linker regulates functional properties as well as tethering the neighboring domains and that protein stability may be increased markedly by minor modifications of the C-terminus. The sensitivity of this domain to few and conservative variations in helices A and D (D2E, S17A, T70S and M71L) is demonstrated by the rat CaM fragments having lower stability and higher calcium affinity than fragments of the same length derived from Paramecium CaM.     

Functional significance of the central helix in calmodulin

The 3-A crystal structure of calmodulin indicates that it has a polarized tertiary arrangement in which calcium binding domains I and II are separated from domains III and IV by a long central helix consisting of residues 65-92. To investigate the functional significance of the central helix, mutated calmodulins were engineered with alterations in this region. Using oligonucleotide-primed site-directed mutagenesis, Thr-79 was converted to Pro-79 to generate CaMPM. CaMPM was further mutated by insertion of Pro-Ser-Thr-Asp between Asp-78 and Pro-79 to yield CaMIM. Calmodulin, CaMPM, and CaMIM were indistinguishable in their ability to activate calcineurin and Ca2+-ATPase. All mutated calmodulins would also maximally activate cGMP-phosphodiesterase and myosin light chain kinase, however, the concentrations of CaMPM and CaMIM necessary for half-maximal activation (Kact) were 2- and 9-fold greater, respectively, than CaM23. Conversion of the 2 Pro residues in CaMIM to amino acids that predict retention of helical secondary structure did not restore normal calmodulin activity. To investigate the nature of the interaction between mutated calmodulins and target enzymes, synthetic peptides modeled after the calmodulin binding region of smooth and skeletal muscle myosin light chain kinase were prepared and used as inhibitors of calmodulin-dependent cGMP-phosphodiesterase. The data suggest that the different kinetics of activation of myosin light chain kinase by CaM23 and CaMIM are not due to differences in the ability of the activators to bind to the calmodulin binding site of this enzyme. These observations are consistent with a model in which the length but not composition of the central helix is more important for the activation of certain enzymes. The data also support the hypothesis that calmodulin contains multiple sites for protein-protein interaction that are differentially recognized by its multiple target proteins.     

Functional expression of chicken calmodulin in yeast

The coding region of a chicken calmodulin cDNA was fused to a galactose-inducible GAL1 promoter, and an expression system was constructed in the yeast Saccharomyces cerevisiae. Expression of calmodulin was demonstrated by purifying the heterologously expressed protein and analyzing its biochemical properties. When the expression plasmid was introduced into a calmodulin gene (cmd1)-disrupted strain of yeast, the cells grew in galactose medium, showing that chicken calmodulin could complement the lesion of yeast calmodulin functionally. Repression of chicken calmodulin in the (cmd1)-disrupted strain caused cell cycle arrest with a G2/M nucleus, as observed previously with a conditional-lethal mutant of yeast calmodulin. These results suggest that the essential function of calmodulin for cell proliferation is conserved in cells ranging from yeast to vertebrate cells.     

A fluorescent calmodulin that reports the binding of hydrophobic inhibitory ligands.

Ca2+ binding to calmodulin in the pCa range 5.5-7.0 exposes hydrophobic sites that bind hydrophobic inhibitory ligands, including calmodulin antagonists, some Ca2+-antagonists and calmodulin-binding proteins. The binding of these hydrophobic ligands to calmodulin can be followed by the approx. 80% fluorescence increase they produce in dansylated (5-dimethylaminonaphthalene-1-sulphonylated) calmodulin (CDRDANS). In the presence of Ca2+, calmodulin binds the calmodulin inhibitor, R24571, with an affinity of approx. 2-3 nM and hydrophobic ligands, including trifluoperazine (TFP), W-7 [N-(6-aminohexyl)-5-chloronaphthalene-1-sulphonamide], fendiline, felodipine and prenylamine, with affinities in the micromolar range. This binding is strongly Ca2+-dependent and Mg2+-independent. Calmodulin shows a reasonably high degree of specificity in its binding of these ligands over other ligands tested. CDRDANS, therefore, provides a convenient and simple means of monitoring the interaction of a variety of hydrophobic ligands with the Ca2+-dependent regulatory protein, calmodulin. CDRDANS binds to phospholipid vesicles made of (dimyristoyl)phosphatidylcholine (DMPC) or (dipalmitoyl)phosphatidylcholine (DPPC) and produces fluorescence increases only in the presence of Ca2+ and at temperatures above their gel-to-liquid crystalline phase transition. Although the fluorescence changes in CDRDANS accurately report phase transitions in these liposomes, its binding to these vesicles is weak. Calmodulin probably requires a high-affinity lipid-bound receptor protein for its high-affinity binding to natural membranes.     

Molecular dynamics simulations of calcium-free calmodulin in solution

A 4-ns molecular dynamics simulation of calcium-free calmodulin in solution has been performed, using Ewald summation to treat electrostatic interactions. Our simulation results were mostly consistent with solution experimental studies, including NMR, fluorescence and x-ray scattering. The secondary structures within the N- and C-terminal domains were conserved in the simulation, with trajectory structures similar to the NMR-derived model structure 1CFD. However, the relative orientations of the domains, for which there are no NMR restraints, differed in details between the simulation and the 1CFD model. The most interesting information provided by the simulations is that the dynamics of calcium-free calmodulin in solution is dominated by slow rigid body reorientations of the domains. The interdomain distance fluctuated between 29 and 39 A, and interdomain orientation angle, defined as the pseudo-dihedral formed by the four calcium binding sites, varied between -2 degrees and 108 degrees. Similarly, the domain linker region also exhibited significant fluctuations, with its length varying in the 34-45 A range and its bend angle in the 10-100 degrees range. The simulations are in accord with fluorescence results suggesting that calcium-free calmodulin is more compact and more flexible than the calcium activated form. Surprisingly, quite similar solvent accessibilities of the hydrophobic patches were seen in the calcium-free trajectory described in this work and previously generated calcium-loaded calmodulin simulations. Thus, our simulations suggest a reexamination of the standard model of the structural change of calmodulin upon calcium binding, involving exposure of the hydrophobic patches to solvent.     

Structure and dynamics of calmodulin in solution.

To characterize the dynamic behavior of calmodulin in solution, we have carried out molecular dynamics (MD) simulations of the Ca2+-loaded structure. The crystal structure of calmodulin was placed in a solvent sphere of radius 44 A, and 6 Cl- and 22 Na+ ions were included to neutralize the system and to model a 150 mM salt concentration. The total number of atoms was 32,867. During the 3-ns simulation, the structure exhibits large conformational changes on the nanosecond time scale. The central alpha-helix, which has been shown to unwind locally upon binding of calmodulin to target proteins, bends and unwinds near residue Arg74. We interpret this result as a preparative step in the more extensive structural transition observed in the "flexible linker" region 74-82 of the central helix upon complex formation. The major structural change is a reorientation of the two Ca2+-binding domains with respect to each other and a rearrangement of alpha-helices in the N-terminus domain that makes the hydrophobic target peptide binding site more accessible. This structural rearrangement brings the domains to a more favorable position for target binding, poised to achieve the orientation observed in the complex of calmodulin with myosin light-chain kinase. Analysis of solvent structure reveals an inhomogeneity in the mobility of water in the vicinity of the protein, which is attributable to the hydrophobic effect exerted by calmodulin's binding sites for target peptides.     

Structural dynamics of calmodulin and troponin C

We present the results of computational simulation studies of the structures of calmodulin (CAM) and troponin C (TNC). Possible differences between the structures of these molecules in the crystal and in solution were suggested by results from some recent experimental studies, which implied that their conformations in solution may be more compacted than the characteristic dumbbell shape observed in the crystal. The molecular dynamics simulations were carried out with the CHARMM system of programs, and the environment was modeled with a distance-dependent dielectric permittivity and discrete water molecules surrounding the proteins at starting positions identified in the crystals of CAM and TNC. Methods of macromolecular structure analysis, including linear distance plots, distance matrices and a matrix representation of hydrogen bonding, were used to analyze the nature, the extent and the source of structural differences between the computed structures of the molecules and their conformations in the crystal. Following the longest simulation, in which intradomain structure was conserved, the crystallographically observed dumbbell structure of the molecule changed due to a kinking or bending in the region of the central tether helix connecting the two Ca(2+)-binding domains which moved into close proximity. The resulting structure correlates with experimental observations of complexes between CAM and peptides such as melittin and mastoparan. Analysis of the corresponding pair distance distribution functions in comparison to experimental results suggests the dynamic existence of a non-negligible fraction of the compacted structure in aqueous solutions of CAM. In this more nearly globular shape, CAM reveals to the environment two interior pockets that contain a number of hydrophobic residues, in agreement with NMR data suggesting involvement of such residues in the binding of inhibitors and proteins to CAM.     

Calcium-induced changes in calmodulin structural dynamics and thermodynamics

The thermodynamics of the interaction between Ca(2+) and calmodulin (CaM) was examined using isothermal titration calorimetry (ITC). The chemical denaturation of calmodulin was monitored spectroscopically to determine the stability of Ca(2+)-free (apo) and Ca(2+)-loaded (holo) CaMs. We explored the conformational and structural dynamics of CaM using amide hydrogen-deuterium (H-D) exchange coupled with Fourier transform infrared (FT-IR) spectroscopy. The results of H-D exchange and FT-IR suggest that CaM activation by Ca(2+) binding involves significant conformational changes. The results have also revealed that while the overall conformation of holo-CaM is more stable than that of the apo-CaM, some part of its α-helix structures, most likely the EF-hand domain region, has more solvent exposure, thus, has a faster H-D exchange rate than that of the apo-CaM. The ITC method provides a new strategy for obtaining site-specific Ca(2+) binding properties and a better estimation of the cooperativity and conformational change contributions of coupled EF-hand proteins.     

Calcium-dependent hydrophobic chromatography of calmodulin, S-100 protein and troponin-C

We have demonstrated calcium-dependent hydrophobic interactions among calmodulin, S-100 protein and troponin-C and a homologous series of omega-aminoalkyl-agaroses. The three Ca2+-binding proteins were retained on the column of agarose substituted with omega- aminooctyl or even longer with alkylamine, in the presence of Ca2+ and 0.15 M NaCl. As these proteins were not retained on the column with shorter alkylamine 'arms' (N = 2, 4), they are probably successively absorbed with a higher affinity to the hydrophobic agarose column. Calmodulin and S-100 protein were eluted from the aminoocytl -agarose column with 1 mM EGTA in the presence of 0.15 M NaCl and the elution of troponin-C was Ca2+-independently carried out with 0.3 M NaCl. On the other hand, S-100 and troponin-C were eluted Ca2+-dependently from aminodecyl -agarose in the presence of 1 M NaCl and half the amount of the calmodulin applied was eluted with 1 M NaCl. As there are obvious differences among the three Ca2+-binding proteins with regard to chromatographic behavior on omega-aminoalkyl-agarose columns, our results suggest that these three proteins expose different hydrophobic regions following Ca2+-induced conformational changes and, if so, such would explain the interaction with aminoalkyl-agaroses.     

Structure and dynamics of calcium-activated calmodulin in solution.

Two 4-ns molecular dynamics simulations of calcium loaded calmodulin in solution have been performed, using both standard nonbonded cutoffs and Ewald summation to treat electrostatic interactions. Our simulation results are generally consistent with solution experimental studies of calmodulin structure and dynamics, including NMR, cross-linking, fluorescence and x-ray scattering. The most interesting result of the molecular dynamics simulations is the detection of large-scale structural fluctuations of calmodulin in solution. The globular N- and C-terminal domains tend to move approximately like rigid bodies, with fluctuations of interdomain distances within a 7 A range and of interdomain angles by up to 60 deg. Essential dynamics analysis indicates that the three dominant types of motion involve bending of the central helix in two perpendicular planes and a twist in which the domains rotate in opposite directions around the central helix. In the more realistic Ewald trajectory the protein backbone remains mostly within a 2-3 A root-mean-square distance from the crystal structure, the secondary structure within the domains is conserved and middle part of the central helix becomes disordered. The central helix itself exhibits limited fluctuations, with its bend angle exploring the 0-50 degrees range and the end-to-end distance falling in 39-43 A. The results of the two simulations were similar in many respects. However, the cutoff trajectory exhibited a larger deviation from the crystal, loss of several helical hydrogen bonds in the N-terminal domain and lack of structural disorder in the central helix.     

Ca2+-induced hydrophobic site on calmodulin: application for purification of calmodulin by phenyl-Sepharose affinity chromatography

Calmodulin binds quantitatively to phenyl-Sepharose and octyl-Sepharose affinity columns in the presence of micromolar concentrations of Ca²⁺. In addition to EGTA, calmodulin also can be eluted from these affinity columns with low ionic strength buffer, non-ionic detergent (i.e., 1% Triton X-100), or ethylene glycol (50%), suggesting hydrophobic interaction. Using hydrophobic interaction chromatography calmodulin can be purified to homogeneity from bovine brain homogenate in a single step. For large-scale purification the protein fraction containing calmodulin was concentrated by isoelectric precipitation prior to application to the affinity column. The yield obtained by this procedure (160–180 mg calmodulin per kg brain) is significantly greater, and the time required (～ 5 hr) is substantially less, than that of previously described procedures for calmodulin purification. It is apparent that phenyl-Sepharose offers several advantages over phenothiazine-Sepharose for affinity purification of calmodulin.     

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.     

Aluminum-induced conformational changes in calmodulin alter the dynamics of interaction with melittin

Studies were undertaken to examine the impact of aluminum-induced structural changes in bovine brain calmodulin on the protein's interface region with melittin, a model for calmodulin's target enzymes. Both steady-state and time-dependent fluorescence characteristics of the single tryptophanyl residue of melittin were employed to derive information on aluminum-related changes in the fluorophore's microenvironment. In the presence of stoichiometric amounts of aluminum ions, calmodulin's target region with melittin appears to be more polar than that with aluminum absent. As a result, upon association of melittin with aluminum-calmodulin, the enhancement of helical arrays is less pronounced. The fluorophore's average microenvironment also is modified such that its apparent lifetime is shortened when aluminum is present. In the presence of aluminum ions, the solvation structure of calmodulin is possibly changed, which may be unfavorable for a proper fit between calmodulin and target proteins.     

Effect of hydrophobic residue substitutions with glutamine on Ca(2+) binding and exchange with the N-domain of troponin C

Troponin C (TnC) is an EF-hand Ca(2+) binding protein that regulates skeletal muscle contraction. The mechanisms that control the Ca(2+) binding properties of TnC and other EF-hand proteins are not completely understood. We individually substituted 27 Phe, Ile, Leu, Val, and Met residues with polar Gln to examine the role of hydrophobic residues in Ca(2+) binding and exchange with the N-domain of a fluorescent TnC(F29W). The global N-terminal Ca(2+) affinities of the TnC(F29W) mutants varied approximately 2340-fold, while Ca(2+) association and dissociation rates varied less than 70-fold and more than 45-fold, respectively. Greater than 2-fold increases in Ca(2+) affinities were obtained primarily by slowing of Ca(2+) dissociation rates, while greater than 2-fold decreases in Ca(2+) affinities were obtained by slowing of Ca(2+) association rates and speeding of Ca(2+) dissociation rates. No correlation was found between the Ca(2+) binding properties of the TnC(F29W) mutants and the solvent accessibility of the hydrophobic amino acids in the apo state, Ca(2+) bound state, or the difference between the two states. However, the effects of these hydrophobic mutations on Ca(2+) binding were contextual possibly because of side chain interactions within the apo and Ca(2+) bound states of the N-domain. These results demonstrate that a single hydrophobic residue, which does not directly ligate Ca(2+), can play a crucial role in controlling Ca(2+) binding and exchange within a coupled and functional EF-hand system.     

Quantitative aspects of the development of a hydrophobic binding site on calmodulin by calcium binding

The interactive binding by calmodulin of Ca²⁺ and 1-anilinonaphthalene-8-sulfonate (1,8-ANS) has been examined. In the presence of saturating levels of Ca²⁺, calmodulin develops one moderately strong binding site for 1,8-ANS, plus one or more weaker sites. The binding of 1,8-ANS by unliganded, or singly liganded, calmodulin is slight; the development of a strong binding site, as well as the characteristic fluorescence enhancement and CD spectrum, requires the binding of two Ca²⁺ ions. Little further change occurs on binding additional Ca²⁺ ions.     

Functional relevance of the internal hydrophobic cavity of urate oxidase

Urate oxidase from Aspergillus flavus is a 135 kDa homo-tetramer which has a hydrophobic cavity buried within each monomer and located close to its active site. Crystallographic studies under moderate gas pressure and high hydrostatic pressure have shown that both gas presence and high pressure would rigidify the cavity leading to an inhibition of the catalytic activity. Analysis of the cavity volume variations and functional modifications suggest that the flexibility of the cavity would be an essential parameter for the active site efficiency. This cavity would act as a connecting vessel to give flexibility to the neighboring active site, and its expansion under pure oxygen pressure reveals that it might serve as a transient reservoir on its pathway to the active site.