The Solution Structure of a Plant Calmodulin and the CaM-binding Domain of the Vacuolar Calcium-ATPase BCA1 Reveals a New Binding and Activation Mechanism*

The type IIb class of plant Ca2+-ATPases contains a unique N-terminal extension that encompasses a calmodulin (CaM) binding domain and an auto-inhibitory domain. Binding of Ca2+-CaM to this region can release auto-inhibition and activates the calcium pump. Using multidimensional NMR spectroscopy, we have determined the solution structure of the complex of a plant CaM isoform with the CaM-binding domain of the well characterized Ca2+-ATPase BCA1 from cauliflower. The complex has a rather elongated structure in which the two lobes of CaM do not contact each other. The anchor residues Trp-23 and Ile-40 form a 1-8-18 interaction motif. Binding of Ca2+-CaM gives rise to the induction of two helical parts in this unique target peptide. The two helical portions are connected by a highly positively charged bend region, which represents a relatively fixed angle and positions the two lobes of CaM in an orientation that has not been seen before in any complex structure of calmodulin. The behavior of the complex was further characterized by heteronuclear NMR dynamics measurements of the isotope-labeled protein and peptide. These data suggest a unique calcium-driven activation mechanism for BCA1 and other plant Ca2+-ATPases that may also explain the action of calcium-CaM on some other target enzymes. Moreover, CaM activation of plant Ca2+-ATPases seems to occur in an organelle-specific manner.

its N-terminal region between residues Ala-19 and Leu-43. This region is also part of the auto-inhibitory domain of the membrane pump, and it contains a novel CaM-binding motif (19 and Fig. 1a). Previously, we have characterized the interactions between different SCaM isoforms and a 25-residue synthetic peptide corresponding to the CaMBD of BCA1 (20). Isothermal titration calorimetry studies have revealed that both SCaM1 and SCaM4 can strongly bind to the BCA1 peptide with similar thermodynamic parameters. The results were also consistent with those obtained for well characterized CaM-target peptide interactions from mammalian proteins such as smooth muscle myosin light chain kinase (smMLCK) and CaM kinase I peptides. The CaMBD regions of these proteins bind with a fully ␣-helical structure to Ca 2ϩ -CaM. However, unexpectedly, the solution structure of the BCA1 peptide determined in 30% trifluoroethanol was not a complete ␣-helix, as under these conditions, the BCA1 peptide possessed a disordered region in the middle of the ␣-helical peptide (20).
Here, we describe the solution structure of soybean CaM isoform 4 complexed with the CaMBD of BCA1. The bound BCA1 peptide adopts two separate ␣-helical regions that are connected by a bend rendering a unique angle between the two helical segments. Consequently, the relative orientation of the two lobes of SCaM4 is distinct from those seen in other CaM⅐peptide complexes. We have utilized backbone H-N residual dipolar coupling (RDC) NMR data and rotational diffusion dynamics obtained from 15 N NMR relaxation data to confirm this novel orientation of the two lobes. We have also utilized the isotope-labeled peptide to study the motional and structural properties of the bound BCA1 peptide. This CaM binding mode drastically expands the variety of the possible sequence motifs, which can interact tightly with CaM. We also propose a unique mechanism of CaM target activation that can be elucidated on the basis of the SCaM4⅐BCA1 complex structure.

EXPERIMENTAL PROCEDURES
Sample Preparations-Uniformly 15 N-and 15 N, 13 C-labeled SCaM4 were prepared as described previously (21). The unlabeled BCA1 peptide corresponding to the amino acid sequence of residues 19 -43 of BCA1, Ac-ARQRWRSSVSIVKNRAR-RFRMISNL-amide, was synthesized commercially by Anaspec, Inc. (San Jose, CA). Uniformly 15 N and 2 H, 15 N, 13 C-labeled BCA1 peptides corresponding to the amino acid residues 19 -47 of BCA1, PARQRWRSSVSIVKNRARRFRMISNLE-KLA, with one Pro residue at the N terminus as a cloning artifact were prepared using the ketosteroid isomerase-fusion inclusion body approach as described previously (22). The peptides were further purified by reverse-phase HPLC, and their identity was confirmed by mass spectroscopy. All NMR samples contained ϳ0.3 mM 15 N-or 15 N, 13 C-labeled SCaM4, 5 mM CaCl 2 , 100 mM KCl, 10 mM DL-1,4-dithiothreitol-D 10 , 0.03% NaN 3 , and 0.5 mM 2,2-dimethyl-2-silapentane-5-sulfonate in unbuffered 90% H 2 O/10% D 2 O or 99.99% D 2 O solution. The pH/pD values of samples are 7.0 without consideration of isotope effects. The molar ratio of SCaM4 and BCA1 peptide was 1:1.25 to ensure the complete saturation of SCaM4, whereas a ratio of 1:0.85 was used to ensure the absence of free BCA1 peptide in several NMR experiments where the SCaM4-bound BCA1 peptide was monitored. The samples used for RDC measurements also contained 20 mM bis-Tris (pH 7.0), 300 mM KCl, and ϳ16 mg/ml filamentous phage Pf1 (Asla Biotech Ltd.).
NMR Measurements-All NMR experiments were performed at 30°C on Bruker Avance 500 or 700 MHz NMR spectrometers equipped with triple resonance inverse Cryoprobes with a single axis z-gradient. Sequential assignments of HN, N, CO, C␣, and C␤ resonances of BCA1 peptide-bound SCaM4 were achieved using several experiments including CBCANH, CBCA(CO)NH, HNCO, and HN(CA)CO (supplemental Fig.  4). The sequential backbone resonance assignments of the BCA1 peptide bound to SCaM4 were obtained by HNCACB, HN(CO)CACB, HNCO, HNCA, and HN(CO)CA. Side chain assignments of peptide-bound SCaM4 were obtained by C(CCO)NH-TOCSY, H(CCO)NH-TOCSY, HBHA(CBCA-CO)NH, and HCCH-TOCSY experiments. The assignments of the unlabeled SCaM4-bound BCA1 peptide was achieved based on NOE connectivity by analyzing 13 C, 15 N/F 1 , F 2 -filter NOESY, and 13 C/F 1 , F 2 -filter correlation spectroscopy experiments. All NOESY experiments including 15 N-edited NOESY-HSQC, 13 C-edited NOESY-HSQC, and 13 C/F 3 -filter 13 C/F 1edited NOESY-HSQC were measured with a mixing time of 100 ms. H-N RDC measurements were performed using the in-phase/anti-phase ( 1 H, 15 N)-HSQC experiment (23). All NMR experiments for the backbone dynamics studies of the unlabeled BCA1 peptide bound to 15 N-labeled SCaM4 and the 15 N-labeled BCA1 peptide bound to unlabeled SCaM4 were acquired at 50.68 MHz for the 15 N frequency. The 15 N T 2 relaxation data were obtained using a Carr-Purcell-Meiboom-Gill-FIGURE 1. a, amino acid sequence alignment of CaMBDs from various CaM target proteins. The anchoring residues are highlighted in red. The basic and acidic residues are shown in cyan and magenta, respectively, whereas hydrophobic residues are displayed in green. b, a helical wheel representation of the CaMBD of BCA1. The residues are colored in the same manner as in a. NOS, endothelial nitric oxide synthase. type T 2 experiment, where the field strength of the 180°pulse was 5.6 kHz, and the 180°pulses were applied every 1 ms. The 15 N T 1 data were obtained by applying a 2.5 kHz spin-locking pulse during the relaxation delay time. Both T 1 and T 2 experiments were repeated twice to estimate the uncertainty in the peak intensities. { 1 H}-15 N heteronuclear NOE data were obtained using a 5-s train of 120°proton pulses. All spectra were processed using NMRPipe/NMRDraw (24) and analyzed using NMRView software (25).
Structure Calculation-The initial structure of SCaM4⅐BCA1 was calculated with CYANA (version 2.0) (26) using distance restraints derived from the automated NOE assignment protocol implemented in CYANA, hydrogen bond restraints based on the secondary structure from the chemical shift index, and dihedral angle restraints predicted by TALOS (27). A total of 24 Ca 2ϩ -ligand restraints were introduced according to the well known Ca 2ϩ coordination in EF-hand sites (28). Further structure refinement with the addition of backbone H-N RDC restraints were performed by XPLOR-NIH (version 2.19) (29). Initial estimation for the axial component of the molecular alignment tensor (Da) and the rhombicity (R) were obtained from the lowest energy structure calculated by CYANA using PALES (30). Finally, the 25 lowest energy structures from a total of 200 were selected and used for the analysis. All molecular graphics used in this manuscript were created using MolMol (31) or PyMOL software.
Rotational Diffusion Dynamics-The 15 N T 1 and T 2 data were fitted with the program CurveFit (A. G. Palmer, Columbia University). The uncertainties for the peak intensity used for the fitting were estimated from duplicated data sets. The uncertainty for { 1 H}-15 N NOE values were evaluated using the S.D. of the noise in empty spectral regions of the spectra. Those residues that show low NOE values (Ͻ0.65) were excluded from further analysis because these residues have a slow internal motion, which contributes to the T 1 relaxation. Those residues, which are involved in chemical exchange process that affects the T 2 relaxation time, were detected and removed as described by Tjandra et al. (32). For the remaining residues, the rotational correlation time for the global tumbling ( m ) for each residue were estimated from the R 2 /R 1 ratio using the program R 2 R 1 _ m (A. G. Palmer, Columbia University). The m values were globally fit to the N-and C-terminal lobe of SCaM4 in SCaM4⅐BCA1 using the program quadric_diffusion (A. G. Palmer, Columbia University).

RESULTS
Structure of SCaM4⅐BCA1-The SCaM4⅐BCA1 structure was calculated based on 1773 distance restraints derived from NOEs, including 188 intermolecular distances (supplemental Table 1). Because there were no interdomain NOEs detected between the N-and the C-terminal lobe of SCaM4, to determine the relative orientation of the two lobes, backbone H-N RDC data were also employed in the structure determination, as these provide information on long range order in multidomain proteins (33)(34)(35). Backbone H-N RDC values were available for 83% of the total residues. Fig. 2a shows the superposition of the 25 lowest energy structures. The backbone r.m.s.d. values in the folded regions are 0.58 Ϯ 0.08 and 0.53 Ϯ 0.10 Å when only the N-terminal unit (the N-terminal lobe of SCaM4 and the C-terminal part of BCA1) and the C-terminal unit (the C-terminal lobe of SCaM4 and the N-terminal part of BCA1) are superposed, respectively. However, the backbone r.m.s.d. value was 1.47 Ϯ 0.59 Å when all of the folded regions of SCaM4⅐BCA1 were used for the analysis (supplemental Table  1). Fig. 2b displays the ribbon representation of the lowest energy structure of SCaM4⅐BCA1. The N-terminal lobe structure of SCaM4 in SCaM4⅐BCA1 is almost identical to that determined without peptide (21), where the backbone r.m.s.d. value in the folded region is 1.06 Ϯ 0.06 Å (supplemental Fig.  1a). On the other hand, the C-terminal lobe structure of SCaM4 in SCaM4⅐BCA1 is somewhat different from that without the BCA1 peptide, in which the structure is altered to a more open conformation by the binding of the peptide (supplemental Fig.  1b). The backbone r.m.s.d. value in this folded region is 1.60 Ϯ 0.06 Å. The structural alteration of SCaM4 upon BCA1 peptide binding could also be detected in the correlation between the experimentally derived RDC values of SCaM4⅐BCA1 and the best-fit RDC values calculated from the free Ca 2ϩ -SCaM4 structure. The obtained correlation factor (R) and the quality factor (Q) are 0.97/0.26 for the N-terminal lobe, and 0.82/0.56 for the C-terminal lobe. The SCaM4-bound BCA1 peptide adopts an ␣-helix in its N-terminal region (residues 20 -33) followed by a bend and another turn of ␣-helix (residues 38 -42). The bend changes the direction of the two separated ␣-helices by ϳ48° (  Fig. 3a). The bend region appears to be highly positively charged (Fig.  3c). The N-and C-terminal lobes of SCaM4 interact with the C-and N-terminal ␣-helices of the BCA1 peptide, respectively, to form the complex in an antiparallel orientation relative to BCA1 peptide as had been estimated previously from an analysis of near-UV CD spectra (20). The aromatic side chain of Trp-23 in the BCA1 peptide provides many intermolecular NOEs to the SCaM4 protein (Fig. 4a), and it is deeply buried in the hydrophobic target binding pocket of the C-terminal lobe of SCaM4 (Figs. 2b and 4b). This indicates that the side chain of Trp-23 provides the N-terminal anchor residue for the BCA1 peptide. The sidechain of residue Val-30 at position 8 also interacts with the hydrophobic patch between helices E and F in the C-terminal lobe of SCaM4. On the other hand, the side chains of both Phe-37 and Ile-40 in the BCA1 peptide show many hydrophobic contacts to the N-terminal lobe of SCaM4 (Fig. 4a). The side chain of Ile-40 is located almost at the center of the hydrophobic pocket of the N-terminal lobe, whereas that of Phe-37 interacts with the hydrophobic patch that is formed between helices A and B (Figs. 2b and 4b). This suggests that Ile-40 acts as the C-terminal anchor residue of the BCA1 peptide rather than the more bulky hydrophobic sidechain of Phe-37. Therefore, the two anchor residues of the BCA1 peptide are separated by 16 residues, and together with the Val residue, they represent a "1-8-18" motif.
Structure Comparison to Other Complexes-When only the N-or C-terminal lobe of the SCaM4⅐BCA1 complex are compared with those of CaM in previously reported CaM⅐peptide complexes, the structures closely resemble one another, where the backbone r.m.s.d. values in the well folded regions are only ϳ1 Å for the majority of available complexes (data not shown). However, the relative orientations and the positions of the two lobes of SCaM4 in the SCaM4⅐BCA1 complex are distinct from any other CaM⅐peptide complexes. In Fig. 5, only the N-terminal lobe of SCaM4 is superposed on other CaM⅐peptide complexes, including CaM⅐smMLCK (36), CaM⅐CaM kinase kinase (37), the CaM⅐ryanodine receptor Ca 2ϩ release channel (RYR1) (38), and the CaM⅐plasma membrane Ca 2ϩ -ATPases isoform 4b (C28W) (39). The difference in the relative orientations/ positions of the C-terminal lobes relative to the N-terminal lobes is emphasized in Fig. 5. It is clear that the C-terminal lobe of SCaM4⅐BCA1 is positioned differently from the others. It is worthwhile to note that although the bound C28W CaMBD also possesses a 1-8-18 motif, the relative lobe orientation of  (20). c, the surface electrostatic properties of the BCA1 peptide. The residues that form a basic cluster in the middle of BCA1 peptide are labeled. C, C-terminal; N, N-terminal.   Table 2. In this analysis, the structure with the closest relative orientation was determined as CaM⅐CaMKK. However, the C-terminal lobe of the CaM⅐CaMKK complex is located much closer to the N-terminal lobe and forms a more compact overall structure than SCaM4⅐BCA1 (Fig. 5). Because of the larger spacing between the two anchor residues (1-18 motif), the SCaM4⅐BCA1 structure seems more ellipsoidal and elongated than the other complexes where the overall dimension of the protein complex is ϳ55 ϫ 35 ϫ 30 Å, whereas the reported dimension for CaM⅐smMLCK is ϳ50 ϫ 30 ϫ 25 Å.

Structure of SCaM4 Complexed with Vacuolar Calcium-ATPase
Rotational Diffusion Dynamics-To confirm the unique relative orientation determined in the SCaM4⅐BCA1 structure, we have performed a rotational diffusion analysis using the available backbone 15 N NMR relaxation data. When the two domains tumble or orient together in the protein structure, the experimentally derived diffusion tensors for each domain should be colinear, and therefore, it can be used to determine the relative domain orientation (40). Due to several overlapped resonances, 81.8% of the total residues (except for one Pro) were available, and m values were estimated from the R 2 /R 1 ratio for 66.9% of the total residues after filtration of the data (supplemental Fig. 2). The m values ranged between 9.31 and 12.30 ns. After fitting the data to all three diffusion models including isotropic, axially symmetric, and fully anisotropic models, the F-static value for the entire SCaM4 protein was 26.21 (p Ͻ 1.0 ϫ 10 Ϫ5 ) in the axially symmetric model that represents a drastic improvement in the fitting over the isotropic model, whereas that in the fully anisotropic model against the axially symmetric model was only 0.51 (p ϭ 0.60), which indicates that the axially symmetric model is the best describing the rotational diffusion properties of SCaM4⅐BCA1. The correlation time for the global molecular tumbling was obtained as 11.20 and 11.10 ns for the N-and C-terminal lobe of SCaM4, respectively. These values are almost identical to the correlation time of 11.12 ns for the entire SCaM4 protein (supplemental Table 3). The principle axes of the diffusion tensors for the N-and C-terminal lobe of SCaM4 are defined separately and visualized on the five lowest energy structures (Fig. 6). The two axes agree very well in all the structures analyzed and are reasonably aligned with the long axis of the SCaM4⅐BCA1 molecule, supporting the relative lobe orientation determined from the RDC data.
Structural and Dynamic Properties of Isotope-labeled BCA1 Peptide in Complex-Because of the unique structure of the BCA1 peptide in the SCaM4⅐BCA1 complex, we sought to investigate the detailed structural and motional properties of the bound BCA1 peptide. Therefore, we have generated an isotope-labeded BCA1 peptide for multinuclear NMR studies. The isotope-labeled BCA1 peptide was bacterially synthesized and designed to be slightly longer than the chemically synthesized amidated BCA1 peptide to ensure that the C-terminal charge does not influence the SCaM4 binding. Supplemental Fig. 3a shows the 1 H, 15 N-HSQC spectrum of the free BCA1 peptide. It is clear that the BCA1 peptide in aqueous solution does not possess any specific conformation, whereas the HSQC spectrum upon addition of Ca 2ϩ -SCaM4 represents characteristic folding of the BCA1 peptide (supplemental Fig. 3b). The backbone resonance assignments were obtained unambiguously except for some residues that experience severe broadening by a chemical exchange process. Subsequently, backbone dynamics experiments including { 1 H}-15 N heteronuclear NOE, T 2 , and T 1 measurements could be carried out for most residues (Fig. 7). These results showed that the BCA1 peptide retains a relatively rigid conformation for the entire SCaM4 bound region (residues 20 -42). Also, the bend region does not appear to be flexible. We also analyzed the ratio between the T 1 and T 2 values to detect possible motions on a timescale of several hundred microseconds (41, Fig. 7). The T 1 /T 2 values except for Ile-40 and several C-terminal residues are very close to one, suggesting the absence of such motions in the bend region of BCA1. The isotope-labeled BCA1 peptide also allowed us to acquire H-N backbone RDC data for the bound BCA1 peptide. The RDC data were available for 27 residues out of 29 and used for structural validation of the bound BCA1 peptide. Although the RDC data were not included in the structure calculation, the R and Q factors for the best-fit calculation against the BCA1 peptide structure using the two helical regions (residues 20 -32 and 39 -42) are 0.93 and 0.15, respectively, supporting the presence of the unique fixed angle between the two helical portions of the bound BCA1 peptide. We note that the residues in the immediate bend region (residues 33-38) were excluded from the analysis as the structure of the bend region was not well determined in our structure calculation due to a paucity of structural restraints in this region. On the other hand, the bestfit calculation against the straight 1-18 C28W peptide using the same RDC data set produced rather poor correlations (R and Q factors are 0.77 and 0.26, respectively).

DISCUSSION
The structures of the CaMBDs in previously determined CaM⅐peptide complexes have often been classified into "1-5-10", "1-8-14", or "1-16" motifs, although there are some CaMBDs, such as the myristoylated alanine-rich C kinase substrate peptide, that do not conform to this pattern (42)(43)(44)(45). The solution structure of the SCaM4⅐BCA1 complex revealed a novel 1-8-18 CaM-binding mode that was dictated by the unique structure of the bound BCA1 peptide. A 1-8-18 CaMbinding motif has been predicted previously for the CaMBD of the plasma membrane Ca 2ϩ -ATPase isoform 4b (46) and most recently, an NMR solution structure has been determined with the peptide corresponding to its CaMBD (C28W) (39). Unlike the CaM-bound C28W peptide, which forms a nearly perfect ␣-helix for the entire region, the structure of the bound BCA1 peptide in SCaM4⅐BCA1 resembles the reported structure of the BCA1 peptide determined in a 30% trifluoroethanol solution, where the BCA1 peptide consists of two short ␣-helices that are separated by a disordered region around the basic "RARR" cluster ( Fig. 3b) (20). It is known that CaMBD peptides are unstructured in aqueous solution and form a helical structure in trifluoroethanol solution, which resembles that seen in the CaM⅐peptide complexes (47). This earlier result suggested that the basic cluster portion of the BCA1 peptide was intrinsically flexible. A remarkable feature of the bound BCA1 peptide structure in SCaM4⅐BCA1 is that the basic cluster forms a well defined relatively fitted bend and the angle between the two ␣-helices can be defined as 48°Ϯ 5.3° (Fig. 3a). Consequently, the relative orientation of the two lobes of SCaM4 is distinct from any previously reported CaM⅐peptide complexes (Fig. 5). Although the bending region of the BCA1 peptide is relatively poorly converged in the calculated SCaM4⅐BCA1 structures (Figs. 2a and 3a), the correlation times obtained for the N-and C-terminal lobe of SCaM4 are almost identical to that of the entire SCaM4 protein (supplemental Table 3). This indicates that the bend region of the BCA1 peptide is not flexible and that it does not allow the N-and C-terminal units of SCaM4⅐BCA1 to behave independently as is seen for the free CaM protein (48). The two principle axes determined for the diffusion tensors for the N-and C-terminal lobe of SCaM4 are also consistent with the orientation of the two lobes in the structure of SCaM4⅐BCA1 (Fig. 6). However, it is possible that the bend region undergoes a slow internal motion, which could influence our analysis. It has been discussed that RDC data can be influenced by long time averaged information, whereas 15 N relaxation experiments are only sensitive to nanosecond motions (49). The recently determined crystal structure of CaM⅐RYR1 also demonstrated a large "1-17" spacing of the two anchor residues without any interlobe interactions of the two domains of CaM, and the authors obtained a discrepancy between the resulting structure and the H-N RDC data (38). This might also be caused by a slow motion in the middle of the bound RYR1 peptide, and indeed, the presence of such a dynamic region in the CaM-bound RYR1 peptide was suggested by fluorescence resonance energy transfer measurements (50). Importantly, the CaMBD of RYR1 also contains a basic cluster in the middle of the sequence (Fig. 1a). To address these issues for the vacuolar ATPase complex, we have studied the dynamic properties of the bound isotope-labeled BCA1 peptide. The { 1 H}-15 N NOE and T 2 relaxation data and the T 1 /T 2 analysis clearly demonstrated that the bend region has no flexibility compared with the helical regions on the time scales of either nanoseconds or a few hundred microseconds (Fig. 7). We have therefore concluded that the relatively poorly determined structure around the bend is solely due to a lack of structural information rather than dynamic properties of the bound BCA1 peptide. As the hairpin-like loop that was seen in the CaM-bound CaMKK peptide also involved a highly basic region (Fig. 1a) (37,51), it seems that a basic cluster in the middle of a CaMBD peptide can induce a bound peptide conformation that deviates from a regular continuous ␣-helix. Although there is no direct answer to the question why the BCA1 peptide has a bend at a particular angle, we can offer two explanations. First, should BCA1 form a continuous ␣-helical structure, the basic residues located around the bend region, including Arg-33, Arg-35, Arg-36, and Arg-38 will be equally distributed around the surface of the ␣-helix, which is obvious in Fig. 1b. This would cause direct contacts between these positive charges and the hydrophobic target binding patch of the N-terminal lobe of CaM. Second, and more importantly, in our structure, the basic cluster of the bound BCA1 peptide is positioned between helices A and E of SCaM4 and seemingly forms favorable electrostatic contacts with several negatively charged SCaM4 surface residues that include Asp-11, Glu-14, Glu-84, and Glu-87 ( Fig. 3c and Fig. 8). These ion-pair interactions are probably essential to keep the two lobes of SCaM4 in a specific orientation leading to a relatively well defined angle between the two ␣-helices of the bound BCA1 peptide.
The first anchor residue of the BCA1 peptide Trp-32 shows almost identical hydrophobic contacts to the hydrophobic pocket of the C-terminal lobe of SCaM4 to those seen in many other CaM⅐target complexes (e.g. CaM⅐smMLCK, CaM⅐RYR1, and CaM⅐CaM kinase I) (Fig. 4). On the other hand, the second anchor residue Ile-40 is somewhat differently positioned in the hydrophobic pocket of the N-terminal lobe of SCaM4. In Fig.  9a, the interaction of the N-terminal lobe with the peptide is compared with that seen in the CaM⅐smMLCK complex. It appears that Ile-40 is located closer to the center of the hydrophobic pocket than the second anchor residue Leu-813 of smMLCK. The Phe-37 side chain of the BCA1 peptide also interacts with the hydrophobic patch of the N-terminal lobe of SCaM4, unlike the corresponding residue Ile-810 of smMLCK, which points to the other side of the peptide. Consequently, the orientation of the BCA1 peptide relative to the N-terminal lobe of SCaM4 is clearly different from that found in CaM⅐smMLCK (Fig. 9a). On the other hand, the recently reported CaM complex with the NMDA receptor peptide revealed a polar residue, Thr-886, as an anchor residue for the N-terminal lobe that is positioned almost the same as Ile-40 in SCaM4⅐BCA1 (Fig. 9b) (52). In addition, Ile-883 in the CaM⅐NMDA receptor complex also occupies a similar position as Phe-37 in SCaM4⅐BCA1.
Because several different type IIB Ca 2ϩ -ATPases have been identified in different plant species, we have compared the amino acid sequence of the BCA1 peptide with those of the CaMBDs from the other type IIB Ca 2ϩ -ATPases (Fig. 9c). In A. thaliana and rice, these have been classified in four subfam- FIGURE 8. The surface structure of SCaM4 is displayed together with a ribbon structure of SCaM4⅐BCA1. The side chains of the two Ser residues that have been reported to become phosphorylated by protein kinase C are shown (53). The acidic side chains of SCaM4 and the basic side chains of the BCA1 peptide that can form electrostatic interactions are also highlighted. ilies based on genetic analysis (14). It should be noted that no CaMBD could be detected in the N-terminal region in one subfamily that includes ACA12 and ACA13 (19). Therefore, the known CaMBDs of type IIB Ca 2ϩ -ATPases can be classified into three different groups according to their amino acid sequences as shown in Fig. 9d. Group 1 includes BCA1, ACA4, and ACA11, and these calcium pumps are all localized on the prevacuolar or vacuolar membrane, and all the key hydrophobic residues are conserved. Consequently, all the members in this group are expected to have a CaM-binding mode that is identical to that determined here for the SCaM4⅐BCA1 complex. For the group2 ATPases, which include ACA1, ACA2, ACA7, and soybean type IIB Ca 2ϩ -ATPase (SCA1), the Ile residue at position 18 is replaced by a Thr, and instead, an extra Phe residue is located at position 17. This would correspond to the 1-17 CaMBD motif reported for the CaM⅐RYR1 complex. Group 3, which includes ACA8, ACA9, and ACA10, contains a typical 1-8-14 motif because of a deletion at position 5. The latter calcium pumps are therefore more likely to form a CaM complex that resembles CaM⅐smMLCK and CaM⅐CaM kinase I (53). Group 2 and 3 Ca 2ϩ -ATPases appear to be located on the chloroplast/endoplasmic reticulum and plasma membranes, respectively. It is interesting that the Ca 2ϩ -CaM induced structures for CaMBDs in Ca 2ϩ -ATPases located in physically distinct locations in the cell all have a different structure. Consequently, they must have a somewhat different activation mechanism.
Similar to many other CaM-target proteins, the CaMBDs of the type IIB Ca 2ϩ -ATPases can become phosphorylated in vivo. For BCA1, it has been shown that Ser-26 and/or -28 are potential target sites for protein kinase C and that the phosphorylation of Ser-28 had no effect on CaM binding (54). In the SCaM4⅐BCA1 structure, due to the unusual relative positioning of the two lobes of SCaM4, residue Ser-28 is not covered by the N-terminal lobe and is fully exposed to the solvent, which is consistent with the protein phosphorylation results (Fig. 8). On the other hand, Ser-26 is completely covered by the hydrophobic patch of the C-terminal lobe of SCaM4, and hence, phosphorylation of Ser-26 should markedly influence and probably prevent the binding of Ca 2ϩ -CaM to BCA1.
Here, we propose a possible mechanism in which the bend CaMBD structure that is caused by SCaM4 binding is impor-tant for the activation of the calcium pump. BCA1 consists of seven transmembrane domains, and the catalytic domain is located between the second and the third transmembrane domain (residues 251-315) (Fig. 10) (54). The N-terminal auto inhibitory domain partially overlaps with the CaMBD of BCA1. In resting cells with a low Ca 2ϩ concentration, SCaM4 would be either free from BCA1 or it interacts with its CaMBD through the C-terminal lobe because it is well known that the Ca 2ϩ affinity of the C-terminal lobe of CaM is very high in the presence of a target (Fig. 10b). This would facilitate an immediate response to a Ca 2ϩ influx. However, Ca 2ϩ binding to the N-terminal lobe of SCaM4 is required for full activation of BCA1. If the CaM-CaMBD⅐AID complex remains at the catalytic domain after their binding, there could still be steric hindrance to prevent substrate binding. Therefore, it has been discussed that additional mechanisms must be required to dislocate the AID from the catalytic domain. In the case of the activation of smMLCK and NAD kinase by CaM, the importance of additional interactions between the surface of the CaM-CaMBD⅐AID complex and other region of these enzymes have been described as a requirement for full activation (55)(56)(57). Such interactions are not required for BCA1, as removal of the AID from the catalytic domain can be directly achieved by lifting up the AID and the neighboring region by CaM binding (Fig. 10c). The unique bend complex structure would create enough room to allow access for a relatively small ATP molecule to the active site of BCA1.
As the vacuole is the main Ca 2ϩ storage site in plant cells, vacuolar Ca 2ϩ -ATPases play a central role in Ca 2ϩ homeostasis in plants. In this work, we have described the novel solution structure of SCaM4⅐BCA1, a plant CaM in complex with a target protein from a plant species, and proposed a novel CaMtarget activation mechanism. Because of the highly conserved amino acid sequence and tertiary structure of CaM in various eukaryotes, a similar mechanism of action can be expected for CaM target protein activation in other organisms including mammals. Importantly, the observation that CaMBDs can have a break in the middle of an ␣-helical sequence drastically expands the variety of potential CaM-binding motifs. The standard helical wheel prediction that is commonly used to predict anchor residues for CaMBDs obviously does not apply in this case. In fact, a recent protein array study to screen for CaM-binding targets among 37,200 human brain proteins revealed 72 novel CaM-binding proteins, most of which did not contain a typical CaMBD motif (58). Finally, although the CaMBDs identified in the other isoforms of type IIB Ca 2ϩ -ATPases with a different subcellular localization appear to belong to well known motifs (1-14 or 1-17), they all possess the same cluster of basic residues that forces a bend in the middle of the BCA1 CaMBD sequence (Fig. 9c). Therefore, we cannot rule out the possibility that the bound CaMBDs from the other FIGURE 10. Model for a possible target activation mechanism from SCaM4⅐BCA1. a, only the N-terminal cytosolic domain, which contains the CaMBD and AID, and the second cytosolic domain, which contains the catalytic domain, are shown. b, at a low level of Ca 2ϩ , SCaM4 would be either free from the CaMBD or bound to the CaMBD only through its C-terminal lobe. c, at a high level of Ca 2ϩ , the N-terminal lobe binds Ca 2ϩ and interacts with the CaMBD of BCA1 and induces a bend causing the release of the AID from the active site region leading to the full activation of BCA1. DECEMBER 3, 2010 • VOLUME 285 • NUMBER 49

JOURNAL OF BIOLOGICAL CHEMISTRY 38509
Ca 2ϩ -ATPase isoforms will also adopt a bend and have a similar activation mechanism for full enzyme activation. Further structural investigations of the other isoforms is required to fully understand the isoform-and organelle-specific CaM regulation of plant type IIB Ca 2ϩ -ATPases.