Original Research ARTICLE
Arabidopsis Calmodulin-Like Proteins, CML15 and CML16 Possess Biochemical Properties Distinct from Calmodulin and Show Non-overlapping Tissue Expression Patterns
- 1Department of Biology, Queen's University, Kingston, ON, Canada
- 2Protein Function Discovery Laboratory, Department of Biomedical and Medical Sciences, Queen's University, Kingston, ON, Canada
Calcium ions are used as ubiquitous, key second messengers in cells across eukaryotic taxa. In plants, calcium signal transduction is involved in a wide range of cellular processes from abiotic and biotic stress responses to development and growth. Calcium signals are detected by calcium sensor proteins, of which calmodulin (CaM), is the most evolutionarily conserved and well-studied. These sensors regulate downstream targets to propagate the information in signaling pathways. Plants possess a large family of calcium sensors related to CaM, termed CaM-like (CMLs), that are not found in animals and remain largely unstudied at the structural and functional level. Here, we investigated the biochemical properties and gene promoter activity of two closely related members of the Arabidopsis CML family, CML15 and CML16. Biochemical characterization of recombinant CML15 and CML16 indicated that they possess properties consistent with their predicted roles as calcium sensors. In the absence of calcium, CML15 and CML16 display greater intrinsic hydrophobicity than CaM. Both CMLs displayed calcium-dependent and magnesium-independent conformational changes that expose hydrophobic residues, but the degree of hydrophobic exposure was markedly less than that observed for CaM. Isothermal titration calorimetry indicated two and three calcium-binding sites for CML15 and CML16, respectively, with affinities expected to be within a physiological range. Both CML15 and CML16 bound calcium with high affinity in the presence of excess magnesium. Promoter-reporter analysis demonstrated that the CML16 promoter is active across a range of Arabidopsis tissues and developmental stages, whereas the CML15 promoter activity is very restricted and was observed only in floral tissues, specifically anthers and pollen. Collectively, our data indicate that these CMLs behave biochemically like calcium sensors but with properties distinct from CaM and likely have non-overlapping roles in floral development. We discuss our findings in the broader context of calcium sensors and signaling in Arabidopsis.
Cellular signaling has evolved in eukaryotes as a means of detecting and deciphering information and organizing appropriate and timely responses to internal and environmental stimuli. The use of small molecules and ions as second messengers is a highly conserved aspect of signal transduction and information processing in eukaryotes. Although most cells employ a range of second messengers, the divalent metal Ca2+ is considered the most universal and ubiquitous as it functions across a diverse range of cellular signaling events (Edel and Kudla, 2015; Marchadier et al., 2016; Edel et al., 2017).
In plants, Ca2+ signaling plays important roles in responses to stimuli such as abiotic stress (e.g., drought, salinity, temperature stress), biotic stress (e.g., microbial attack, viral attack), as well as many developmental cues (Ranty et al., 2006; DeFalco et al., 2010). In the standard Ca2+ signal transduction model, resting cytosolic Ca2+ [Ca2+]cyt is maintained at nM (~100 nM) levels through active transport into Ca2+ stores (e.g., central vacuole) or out of the cell into the apoplast, thereby creating Ca2+ gradients several orders of magnitude (Steinhorst and Kudla, 2014). Complex patterns of cytosolic Ca2+ influx and efflux, controlled via channel and pump activities, are evoked in response to various stimuli. The Ca2+-signature hypothesis posits that these spatial-temporal patterns of Ca2+ signals encode information about the nature of the evoking stimulus (Berridge et al., 2003; Bender and Snedden, 2013). A corollary of this model is that different stimuli lead to distinct Ca2+ signatures that can be decoded by mechanisms within the cell. Others have suggested that Ca2+ signals act less like coded signatures and more like a switch or lock-and-key in combination with additional messengers (Scrase-Field and Knight, 2003; Plieth, 2016). These models are not mutually exclusive and recognize that, regardless of the nature of the signal, specific Ca2+-binding proteins must serve as sensors to detect changes in Ca2+ levels and function to activate appropriate signaling cascades. Calmodulin (CaM) is the canonical Ca2+ sensor and represents one of the most evolutionarily-conserved eukaryotic proteins. CaM is a sensor relay protein in that it lacks inherent catalytic activity and functions instead through the direct binding and regulation of downstream targets including ion-pumps and channels, transcription factors, metabolic enzymes, cytoskeletal proteins, and a wide range of other proteins (Ikura and Ames, 2006; Bender and Snedden, 2013). CaM is a small (~17 kDa) protein that has been well-characterized at the structural level, consisting of N- and C-terminal globular regions separated by a flexible central helix, giving rise to what has been described as a dumbbell-shaped molecule. Each globular region of CaM contains a pair of Ca2+-binding EF-hand motifs. A classic EF-hand consists of a 29-residue helix-loop-helix structure (Strynadka and James, 1989). Seven oxygen-containing groups, five in the EF loop (at positions 1, 3, 5, 7, and 9) and two donated by the glutamate outside of the loop (at position 12), are responsible for coordinating a Ca2+ ion in a pentagonal bipyramidal arrangement (Strynadka and James, 1989; Falke et al., 1994; Gifford et al., 2007; Grabarek, 2011). Residues 1, 3, 5, 7, 9, and 12 are referred to as X, Y, Z, –Y, –X, and –Z, respectively. Ca2+-binding to the EF-hands of sensors triggers reversible conformational changes in tertiary structure that lead to the exposure of hydrophobic residues. These hydrophobic regions facilitate interactions with downstream effector proteins, which in turn elicit a physiological response to the initial Ca2+ signal (Yamniuk et al., 2004). In CaM, as in most EF-hand Ca2+ sensors, EF-hands are paired and exhibit positive cooperativity with each hand displaying an increased affinity for Ca2+ after the partner has bound Ca2+ (Akke et al., 1991; Gifford et al., 2007).
In the genetic model plant, Arabidopsis, there are about 250 EF-hand proteins encoded in the genome (Day et al., 2002). Similarly, 262 EF-hand proteins are predicted in soybean, suggesting such complexity is found throughout plant taxa (Zeng et al., 2017). In addition to CaMs, plants possess three additional large families of Ca2+ sensors: Ca2+-dependent protein kinases (CPKs), calcineurin B-like proteins (CBLs), and CMLs (DeFalco et al., 2010; Ranty et al., 2016). CaMs, CBLs, and CMLs fall into the class of sensor-relay proteins, whereas CPKs are considered sensor responders as they possess both Ca2+ binding and kinase activities. Remarkably, although CaM is highly conserved throughout all eukaryotes, CBLs, CPKs, and CMLs are found only in plants and some protists (DeFalco et al., 2010; Beckmann et al., 2016). The Arabidopsis genome encodes 7 CaM genes (AtCaM1-AtCaM7), with identical isoforms and 4 distinct but highly conserved isoforms that differ by up to 4 amino acids (McCormack and Braam, 2003). In addition, Arabidopsis possesses 50 CML proteins (AtCMLs) that range from 16 to 75% sequence identity with conserved CaM (e.g., AtCaM2) and are arranged into 9 phylogenetic subgroups (McCormack and Braam, 2003). This large diversity of CMLs is seen across plant taxa (Boonburapong and Buaboocha, 2007; Zhu et al., 2015). Like CaM, these CMLs are thought to localize predominantly to cytosolic and nuclear compartments where Ca2+ signals are important (McCormack and Braam, 2003; DeFalco et al., 2010), although CML30 and CML3 have been observed in mitochondria and peroxisomes, respectively (Chigri et al., 2012).
Structurally, plant CMLs resemble CaM and are predicted to possess EF-hands, and no other known functional domains. Arabidopsis CaMs have 149 amino acids and possess 4 EF-hands, whereas AtCMLs range from 80 to 330 amino acids and possess from 2 to 4 predicted EF-hands with CML12 (6 EF-hands) deviating from this paradigm (McCormack and Braam, 2003). The evolutionary divergence of CMLs from CaM, their conservation across plant taxa, and the fact that CMLs are restricted to plants and some protists, suggests that CMLs play unique roles in plant Ca2+ signaling pathways. Although a number of recent reports have implicated plant CMLs in processes ranging from biotic and abiotic stress response to development, the vast majority of CMLs in Arabidopsis, and in plants in general, remain unstudied.
In order to understand how plants coordinate physiological responses from such a common second messenger, a comprehensive understanding of the Ca2+ sensors unique to plants is required. To this end, we have been systematically studying the biochemical and physiological properties of Arabidopsis CMLs including CML42, CML43, CML37, CML38, CML39 (Chiasson et al., 2005; Vanderbeld and Snedden, 2007; Dobney et al., 2009; Bender et al., 2013, 2014; Scholz et al., 2014). Here, were describe an investigation into two putative Ca2+ sensors of unknown function, CML15 and CML16, which are two closely-related paralogs (73.9% identity) from Arabidopsis CML subfamily four. We discuss their distinct biochemical properties and gene promoter activity patterns in light of Ca2+ signaling in different tissues and developmental contexts.
Materials and Methods
Plant Material and Growth Conditions
All experiments were performed with Columbia (Col-0) accession Arabidopsis thaliana, referred to as Arabidopsis. All Arabidopsis seeds sown to soil were stratified at 4°C for 3 days in the dark, before transfer to a growth incubator set at a 16-h light (50 μmol/m2s)/8-h dark photoperiod at 22°C and at 70% relative humidity. Plants were supplemented with 1 g/L 20-20-20 (N-P-K) fertilizer every 3 weeks until senescence. In some experiments, seeds were surface sterilized as described (Bender et al., 2013) prior to plating and stratification on sterile petri dishes with 0.5X Murashige and Skoog (MS) medium (Murashige and Skoog, 1962; Caisson Laboratories, Inc.) supplemented with 0.8% (w/v) agar followed by cold-stratification as described above. Transgenic lines for CML promoter-β-glucuronidase (GUS) reporter activity analysis (CMLpro::GUS plants) were selected on media containing kanamycin as described (Bender et al., 2013) supplemented with 1% sucrose. Transgenic Arabidopsis CML15pro::GUS and CML16pro::GUS transgenic lines were grown for two subsequent generations to obtain homozygous lines for analyses of CMLpro::GUS reporter activity. Preliminary analysis showed comparable and reliable qualitative patterns of expression among several lines tested and CML15pro::GUS lines 1A and 3C were used to assay CML15pro::GUS reporter activity, and lines 12A and 14B were used to assay CML16pro::GUS reporter activity.
Plasmid Construct Design and cDNA Cloning
The protein-coding regions of CML15 and CML16, both intronless genes, were cloned into the pET21b vector (Novagen) for recombinant expression in Escherichia coli (E. coli). The genomic region from the ATG translation start site to the TGA translation stop site of CML15 (AT1G18530), and CML16 (ATG3G25600), were amplified by PCR and independently subcloned into the pET21b expression vector (Novagen) using the CML-F and CML-R primers described in Supplementary Table 1. The evolutionarily-conserved CaM, CaM81 from petunia, which is 100% identical to Arabidopsis CaM2, and cloned into a pET5a vector (Novagen), was a gift from Hillel Fromm (Tel Aviv University). All constructs were transformed into E. coli strain BL21 (DE3) RPRIL or PLysS (Novagen), for the expression of recombinant CML15 and CML16, or CaM81, respectively.
Cloning the CML15 Promoter and CML16 Promoter into pBI101
The putative promoter regions of CML15 and CML16 were amplified from genomic DNA using CMLPRO-F and CMLPR-R primers (Supplementary Table 1) and cloned into the binary vector pBI101 to control the expression of the gene-reporter, β-glucuronidase (GUS) (Jefferson et al., 1987). The regions of genomic DNA upstream of CML15 (CML15pro) and CML16 (CML16pro) from the stop site of the nearest upstream gene to the ATG start translation site of each respective CML were subcloned into the binary vector pBI101 to yield CML15pro::GUS and CML16pro::GUS constructs. CML15pro (591 bp) and CML16pro (700 bp) are referred to as the putative promoter regions of CML15 and CML16, respectively. Each construct was confirmed by sequencing (SickKids, Toronto) and chemically transformed into Agrobacterium tumefaciens strain GV3101 for stable transformation into Arabidopsis via the floral dip method (Clough and Bent, 1998).
Histochemical GUS Assays
Transgenic Arabidopsis for CMLpro::GUS were stained for GUS activity, as a proxy for CML-promoter activity, using the protocol essentially described (Jefferson et al., 1987; Vanderbeld and Snedden, 2007) with the following exceptions. Plant tissue was fixed in 90% acetone at 4°C for 1 h, then rinsed in double-deionized water 3 times before incubation in the GUS-staining solution (100 mM NaH2PO4, 100 mM Na2HPO4, pH 7.0, containing 10 mM EGTA and 0.1% (v/v) TritonX-100) with the substrate 5-bromo-4-chloro-3-indolyl-β-D-glucuronide (X-gluc) added to 0.1 mg/ml. Tissues were incubated for ~18 h at 37°C with shaking, and were subsequently rinsed several times with 70% ethanol. Stained tissue was then stored in 70% ethanol. Developmental histochemical GUS assays were performed on 1-, 7-, 10-, 22-, 25-day old, and fully mature Arabidopsis tissues. The representative developmental stages were selected based on Boyes et al. (2001).
Fluorometric GUS Assays
Transgenic Arabidopsis for CMLpro::GUS analyses were quantitatively assayed for GUS activity as described (Vanderbeld and Snedden, 2007). Seedling, mature leaf, and mature floral tissues were harvested in 1.5-mL microfuge tubes and snap frozen with liquid nitrogen. Frozen tissue was ground with a micropestle in 500 μL of QB buffer (0.1 M KPO4, 1 mM EDTA, 0.1% (v/v) Triton X-100, 10% (v/v) glycerol, and 0.1% (w/v) N-lauryl sarcosine sodium salt) in the presence of about 100 mL sterile sand and 100 mL 1.5 mm metal beads to facilitate grinding. The ground tissue was further homogenized using the Next Advanced Bullet Blender® at high speed for 5 min, centrifuged at 10,000 × g for 10 min at 4°C, and the resulting supernatant was assayed for protein concentration using the Bradford reagent (BioShop) as described (Bradford, 1976). Protein extract (50 μL) was added to 450 μL QB buffer supplemented with 1 mM DTT and 0.5 mg/mL 4-Methylumbelliferyl-β-D-glucuronide (MUG; BioShop) and samples were incubated at 37°C for 2 h after which time reactions were quenched with 400 μL of 0.2 M Na2CO3. For quantitative analysis of GUS activity, triplicate samples were assayed in a black Fluotrac-200 96-well microplate and the emission fluorescence of 4-MU (4-methylumbelliferone) at 455 nm was measured after excitation with wavelength of 365 nm using a Spectramax Gemini XS spectrophotometer. Average fluorescence values were subsequently converted to units of specific activity (pmoles 4-MU/μg protein/min).
Recombinant Protein Expression and Purification
Preparation of recombinant CaM and CMLs was essentially as described (Zielinski, 2002). E. coli cultures (10 mL) expressing CaM81, CML15, or CML16, respectively, were grown overnight (~10 h) at 37°C. The 10-mL cultures were used to inoculate 1–1.25 L of LB media, until an OD600 = 0.6–0.8 was reached, at which point isopropyl-β-D-thiogalactopyranoside (IPTG) was added (0.5 mM) to induce recombinant protein expression. After 4–6 h E. coli cells were isolated by centrifugation and suspended in a lysis buffer (50 mM Tris-Cl, pH 7.5, 1.5 mM EDTA, 1 mM DTT, 1 mM benzamidine, 1 mM PMSF) at 1/20 culture volume and were lysed using a French® Pressure Cell Press (Thermo Scientific). Recombinant proteins were purified by Ca2+-dependent phenyl-sepharose column chromatography, essentially as described (Zielinski, 2002) with the following exceptions. The binding of CaM81, CML15, and CML16 (in CaM-binding buffer: 50 mM Tris-Cl, pH7.5) to phenyl-sepharose was performed in the presence of CaM-binding buffer containing 2 mM CaCl2, followed by 2–5 column-volume washes of CaM-binding buffer with 200 mM NaCl, 2–5 column-volume washes of CaM-binding buffer with 400 mM NaCl, prior to elution of adsorbed protein with CaM-binding buffer containing 2 mM EGTA. Protein was eluted over 8 fractions in volumes of elution buffer equivalent to the bed volume of the resin (~1–2 mL). The eluted fractions were resolved on 17% SDS-PAGE gels which were then stained with Coomassie Brilliant Blue (R250) and subsequently destained [50% ddH20 (v/v), 40% MeOH (v/v), and 10% glacial acetic acid (v/v)] to evaluate purity. Insufficiently pure fractions were pooled, and re-purified over phenyl-sepharose 1–2 more times using the method described above. In order to achieve near-homogeneity of recombinant CML16, following phenyl-sepharose chromatography, pooled CML16 samples were resolved by fast protein liquid chromatography (FPLC) using a Superdex-75 gel-filtration column to eliminate contaminating E. coli proteins of ~37 and ~65 kD. Prior to FPLC, the gel filtration column was pre-equilibrated with 2 mM EGTA (in 50 mM Tris-Cl, pH 7.5). Fractions (2 mL) were collected and analyzed for purity by SDS-PAGE and Commassie staining. Pure recombinant CMLs were concentrated using 4,000- or 10,000-kD cutoff Amicon® Ultra-4 spin-column centrifugal filters (Millipore). Protein samples were then dialyzed overnight into various buffers suited for specific biochemical analyses as described below. In an effort to minimize Ca2+ contamination of samples, for all biochemical analyses, milli-Q (Millipore Sigma) purified water was used for all solutions and for rinsing of glassware. Dialysis tubing was treated with 1 mM EDTA and 1 mM EGTA, followed by extensive rinsing with milli-Q purified water.
ANS Fluorescence Spectroscopy
The hydrophobic properties of recombinant CaM, CML15 and CML16 were monitored using ANS (8-anilinonapthalene-1-sulfonic acid) fluorescence emission spectroscopy using an excitation wavelength of 380 nm, and scanning emission spectra 430–600 nm, performed essentially as described (Dobney et al., 2009). CaM and CML samples were dialyzed overnight against 1 L of ANS buffer (10 mM Tris-CL, pH 7.5, containing 1 mM DTT, and 100 mM KCl). The fluorescence emission of 15 μM CML15 and 15 μM CML16 was monitored using 250 μM ANS in ANS buffer. Spectra were recorded at room temperature (22°C) under various conditions of CaCl2, MgCl2, EGTA, or Na3C6H6O7 in ANS buffer, as described in Figure Legends. Ca2+-CaM81 (15 μM) was monitored as a positive control for Ca2+-induced hydrophobicity. ANS background fluorescence was monitored and subtracted from the fluorescence readings of the protein-ANS buffer conditions. Average fold-increases in peak florescence emission of ANS in the presence of 15 μM CML15, 15 μM CML16 or 15 μM CaM81, in the presence of 1 mM CaCl2 + 5 mM MgCl2 vs. 1 mM MgCl2, were calculated as the fold-increase in fluorescence due to Ca2+-binding. Fold-induction was averaged over four independent biological replicates with three technical replicates per condition.
Phenyl-Sepharose Chromatography of Recombinant CMLs under Various Conditions
In order to assess the dependence of hydrophobic exposure of CMLs on Ca2+ binding, the ability of CML15 and CML16 to bind to, and elute from, phenyl-sepharose under various ionic conditions was tested. Protein samples and the associated wash buffers were adjusted to final concentrations of either 1 mM CaCl2, or 1 mM MgCl2, or 1 mM CaCl2 + 1 mM MgCl2, or 1 mM EGTA, in 50 mM Tris-Cl, pH7.5, respectively, prior to phenyl-sepharose binding. Proteins were subjected to hydrophobic interaction column chromatography, essentially as described above. Post phenyl-sepharose binding, columns were washed (10 mL) using the same buffer and conditions as noted above but with the addition 200 mM KCl, followed a second wash (10 mL) which included 400 mM KCl. All column elutions were performed using 50 mM Tris-Cl, pH7.5 containing either 2 mM EGTA (for samples loaded with 1 mM CaCl2), or 3 mM EDTA (for samples loaded with 1 mM MgCl2), or 2 mM EGTA + 3 mM EDTA (for samples loaded with 1 mM CaCl2 and 1 mM MgCl2), or 1 mM EGTA (for samples loaded with 1 mM EGTA). Samples of the unbound protein (flow-through), washes, elutions, and the pre-column protein samples were resolved by SDS-PAGE and analyzed by staining with Coomassie Brilliant blue (R250).
Circular Dichroism (CD) Spectroscopy
Far-UV CD spectroscopy was performed to evaluate the secondary structure of recombinant CML15 and CML16 as described (Bender et al., 2013). Spectra were collected and analyzed at the Queen's Protein Function Discovery facility using a Chirascan CD spectrophotometer. Samples of CML15 and CML16 were independently dialysed into 2 mM Tris-Cl (pH 7.5), quantitated (45 μM and 30 μM, respectively), and analyzed by CD spectroscopy, using a cylindrical quartz cuvette with a pathlength of 0.1 mm. Spectra for both CML15 and CML16 were initially obtained in the presence of 1 mM EGTA, 1 mM CaCl2, or 1 mM MgCl2. To confirm that the spectral changes under these conditions represented either complete Ca2+ chelation or Ca2+ saturation, samples were also analyzed in the presence of 2 mM EGTA or 2 mM CaCl2. For each experimental condition, spectra were obtained as an average of 6 raw scans which were subsequently reference-corrected vs. buffer-only scans and adjusted for molar concentration. Average, corrected CD spectra obtained for each experimental condition were deconvoluted using both the OLIS GlobalWorks and CDNN deconvolution software to estimate protein secondary structural composition.
Isothermal Titration Calorimetry (ITC)
ITC experiments were performed using a MicroCal VP-ITC microcalorimeter to determine the affinities of Ca2+ for the EF-hand domains of CML15 and CML16, and to examine the effect of Mg2+ on these Ca2+-binding events. Prior to ITC titrations, CML15 and CML16 samples were dialysed into 20 mM HEPES (pH 7.5) containing 100 mM KCl, either with or without 5 mM MgCl2. Samples of CML15 or CML16 were independently titrated at 30°C with twenty-nine 10 μL injections of either 600 μM CaCl2 or 600 μM MgCl2 at 6-min intervals. For each experimental condition, titrations were performed in triplicate. Concentrations of CML15 and CML16 used for data analysis of isotherms (see Figure Legends) were determined by amino-acid analysis (SickKids Hospital, Toronto). Origin 7.0 software (MicroCal) was used to determine the enthalpies (ΔH), association constants (Ka), binding entropies (ΔS), and stoichiometries of CML15 and CML16 with respect to Ca2+ and Mg2+ binding. Kd values were determined using 1/Ka.
Molecular Modeling of CMLs
The protein homology/analogy recognition engine (Phyre2.0) website (http://www.sbg.bio.ic.ac.uk/phyre2) was used to predict structural models CML15 and CML16 based on the highest-ranked homologous model in the database using default parameters (Kelley et al., 2015). Model images were generated using the PyMOL molecular graphics system, version 1.9. Schrodinger, LLC (http://www.pymol.org). Grand average of hydropathy (GRAVY) scores were determined using the ProtParam tool from ExPasy online bioinformatics portal (www.expasy.org, Artimo et al., 2012).
Sequence Comparison of CML15, CML16, and Conserved CaM
CML15 and CML16 are closely related paralogs showing 74% sequence identity with each other and 39.5% identity with the evolutionarily-conserved CaM, AtCaM2 (Figure 1). CML15 and CML16, along with CML17 and CML18 comprise the four members of subgroup IV in the Arabidopsis CaM/CML phylogeny (McCormack et al., 2005). A recent reorganization of the Arabidopsis phylogeny places CML15 and CML16 in subgroup II (Zhu et al., 2015). Orthologs of both these CMLs are found across plant taxa (Zhu et al., 2015; Zeng et al., 2017). Although CML15 and CML16 are predicted to possess 4 EF-hand motifs, variations within the Ca2+-binding loops and adjacent regions are present when compared to CaM2. The backbone glycine in position 6 (Gly6), and Glu12, are conserved across the four EF-loops of CML15 and CML16. Conversely, across all four EF-loops of the two CMLs and CaM, the residues found in positions 7 and 9 are variable. This holds true for the residue at position 8 as well, which though variable, is consistently hydrophobic. EF-loop I of CML15 and CML16 do not differ substantially from EF-loop I of AtCaM2 with respect to the residues typically involved in forming the Ca2+-coordination sphere. Residues 1, 3, 5, 9, and 12 of EF-loop I are identical to those of AtCaM2. The third residue of EF-hand II and IV, however, for both CML15 and CML16, has an asparagine (Asn) substitution rather than an aspartate typically found in this position (Gifford et al., 2007). Finally, position 9 of EF-loop III of CML15, CML16, and GmCaM4, and position 9 of EF-loop IV of CML15 and CML16, is occupied by a serine (Ser9). This is unlike AtCaM2, which has a threonine (Thr9) in EF-loop III and an Asp9 in EF-loop IV. This range of substitutions in the EF hands of these CMLs suggests their Ca2+-binding properties are likely distinct from those of CaM.
Figure 1. Amino-acid sequence comparison of Arabidopsis CML15 (accession Q9FZ75) and CML16 (accession Q9LI84) to the evolutionarily-conserved CaM, CaM2 (accession NP_850344). Identical and conserved residues are shaded in black and gray, respectively, and alignment gaps are indicated by dashes. Positions of Ca2+-binding EF-hands for CaM2, and the corresponding regions of the CMLs, are indicated with a line. Residue positions are presented to the left and right of each sequence. Sequences were aligned using Clustal Omega (Sievers et al., 2011).
Calcium-Binding Properties of CML15 and CML16
The increased electrophoretic mobility of CaM in the presence of Ca2+ is a hallmark property that has been used as a comparison when evaluating various putative Ca2+ sensors for their ability to bind Ca2+ (Garrigos et al., 1991; Dobney et al., 2009). Although both CML15 and CML16 displayed shifts in electrophoretic mobility in the presence of Ca2+, the relative mobility shift of CML16 was more pronounced than that observed for CML15 which was barely perceptible (Figure 2).
Figure 2. Ca2+-dependent electrophoretic mobility analysis of recombinant CML15 and CML16 using SDS-PAGE in the presence of (A) 5 mM EGTA or (B) 5 mM CaCl2. Recombinant CaM served as a positive control for electrophoretic mobility shifts in the presence of Ca2+ and purified recombinant glutathione S-transferase (GST) served as the negative control. SDS-PAGE gels were stained with Coomassie Brilliant Blue R250. Molecular weight markers (MW, kDa) are indicated on the left. Representative images from a minimum of three experimental replicates are presented.
An additional property of CaM and related Ca2+ sensors are reversible conformation changes that expose hydrophobic regions upon Ca2+ binding (Yamniuk et al., 2004; Ikura and Ames, 2006). We investigated CML15 and CML16 for Ca2+-mediated increases in hydrophobic exposure using ANS-based fluorescence and phenyl-sepharose chromatography. ANS-fluorescence spectroscopy was conducted to assess the nature of exposed hydrophobicity for CML15 and CML16 under various in vitro conditions. A blue shift and associated increase in the intensity of the fluorescence emission of ANS was observed in the presence of Ca2+-CML15 and Ca2+-CML16 (Figure 3) compared to protein in ANS buffer alone. The greatest relative blue-shifts in emission spectra were observed for CML15 and CML16 in the presence of 1 mM CaCl2 or in the presence of both 1 mM CaCl2 and 5 mM MgCl2. Only a very weak increase in ANS fluorescence was observed for CML15 or CML16 in the presence of 1 mM MgCl2 alone. There was a 4.05 ± 0.31-fold increase in ANS fluorescence for Ca2+-CML15 (in a 5 mM MgCl2 background) compared to Mg2+-CML15. Likewise, there was a 2.44 ± 0.23-fold increase in ANS fluorescence for Ca2+-CML16 (in a 5 mM MgCl2 background) compared to Mg2+-CML16. Based on their amino acid compositions, we note that both CML15 and CML16 have grand average of hydrophobicity (GRAVY) scores of 0.009 and −0.171, respectively, which are higher than that for CaM at −0.619 (Supplementary Figure 3). However, CaM showed a much greater increase in exposed hydrophobicity upon binding Ca2+ than either CML15 or CML16 (Supplementary Figure 1). Interestingly, there were slight increases in ANS fluorescence, suggesting increased hydrophobicity of CML15 and CML16, in the presence of 1 mM EGTA, even in the absence of Ca2+ or with 1 mM EDTA (in the absence of Mg2+) relative to protein in ANS buffer alone or in the presence of 1 mM MgCl2 (Figure 3). However, ANS fluorescence of CML15 and CML16 in the presence of 1 mM Na3C6H5O7 (trisodium citrate), another Ca2+-chelator, was similar to that of CML15 and CML16 in ANS buffer (Figure 3). Collectively, these data indicate that CML15 and CML16 possess intrinsic hydrophobicity in the absence of Ca2+ or Mg2+, and that they display a marked increase in hydrophobic exposure in the presence of Ca2+. In addition, CML16 displays greater intrinsic hydrophobicity than CML15 at peak fluorescence emission. These data also show that EGTA and/or EDTA have a positive effect on the hydrophobic exposure of CML15 and CML16 that is less than that induced by Ca2+, but greater than their intrinsic hydrophobic exposure. For both CML15 and CML16, the increase in exposed hydrophobicity in response to Ca2+ binding is reminiscent of CaM but the magnitude of change and the high level of intrinsic hydrophobicity in the absence of Ca2+ are quite different than observed for CaM (Supplementary Figure 1), indicating unique structural properties for these CMLs.
Figure 3. Ca2+- and Mg2+-induced changes in exposed hydrophobicity of recombinant CML15 and CML16 as demonstrated by ANS fluorescence. Fluorometric scans (430-600 nm) were recorded following the addition of 15 μM (A) CML15 or (B) CML16 to 250 μM ANS where all sample conditions used ANS buffer (10 mM Tris-Cl pH 7.5, 100 mM KCl, and 1 mM DTT). The y-axis depicts fluorescence in international arbitrary units (IAU). Fluorescence values were analyzed relative to a protein-free ANS control sample that was subtracted from the data as background. Apo-CML (CML + ANS alone) scans represent baseline CML hydrophobicity. Scans were recorded under the conditions noted in the figure panels by adding various compounds to apo-CML. Each data set is representative of five experimental replicates.
For comparison with our ANS-based data, we further explored the Ca2+- and Mg2+-mediated hydrophobic properties of CML15 and CML16 using phenyl-sepharose chromatography under various in vitro conditions (Supplementary Figure 2). Both CMLs bound to phenyl-sepharose in the presence of 1 mM CaCl2 and in the presence of 1 mM CaCl2 + 1 mM MgCl2, however, neither CML bound in the presence of 1 mM MgCl2 alone. Interestingly, there was weak binding of CML15 and CML16 to phenyl-sepharose in the presence of 1 mM EGTA (i.e., absence of Ca2+), and a slow elution during column washing with 0.2 and 0.4 M KCl prior to elution with EGTA. Collectively, these data suggest that most of the hydrophobic exposure of these CMLs is triggered by Ca2+ binding but there is also a notable degree of intrinsic hydrophobicity present in the absence of Ca2+, observations consistent with the ANS fluorescent data (Figure 3).
In order to evaluate the secondary structural characteristics of CML15 and CML16 in the apo-form, and in the presence of Mg2+ and Ca2+, we performed far-UV CD spectroscopy (Figure 4). CD spectra of CML15 and CML16 show that both proteins are structured, each presenting predominantly α-helical content in the absence of any divalent cation, as indicated by the peaks of molar ellipticity at 190 nm and troughs at 208 and 222 nm, respectively. CML15 showed a slight (3–6%) decrease in α-helical content in the presence of saturating Ca2+ (1 mM CaCl2), whereas CML16 showed ~8–15% increase in α-helical content in the presence of 1 mM CaCl2. Interestingly, Mg2+-CML16 displayed similar CD spectra to that of Ca2+-CML16, suggesting that both Ca2+ and Mg2+ binding elicits structural increases in the α-helical content of CML16.
Figure 4. Far-UV circular-dichroism (CD) spectroscopy of recombinant (A) CML15 and (B) CML16. Proteins were analyzed in 2 mM Tris-Cl (pH 7.5) in the presence of 1 mM EGTA, 1 mM CaCl2, or 1 mM MgCl2 as noted in figure panels. Far-UV CD spectra (180–260 nm) were collected at room temperature. Each spectrum is representative of 6 scans and is presented in units of molar ellipticity. CML15 and CML16 were analyzed at final concentrations of 45 μM and 30 μM, respectively.
The energetics of Ca2+ (and Mg2+) binding to recombinant CML15 and CML16 were studied by ITC (Figure 5). Raw ITC data representing the endothermic and/or exothermic nature of divalent cation binding events were used to generate Wiseman ITC plots which in turn allowed for estimation of the stochiometries and affinities for Ca2+ or Mg2+ binding, as well as the relative entropic and enthalpic contributions of each binding event (Table 1). Amino-acid sequence data indicates that both CMLs are predicted to have four EF-hands (Figure 1; McCormack and Braam, 2003), however, ITC analysis suggests both CMLs may only have a subset of sites that are active. Data were fitted using the Origin 7.0 software to give stochiometries of two Ca2+-binding sites for CML15 and three Ca2+-binding sites for CML16. Analysis of ITC data predicted Ca2+-binding affinities for CML15 of 0.22 and 1.23 μM for site 1 and site 2, respectively, in the absence of Mg2+. These affinities are reduced to 1.56 and 7.25 μM when excess Mg2+ is present (Table 1). In either the presence or absence of Mg2+, thermodynamic data (Table 1) suggests that favorable entropy drives site-1 Ca2+ binding whereas site-2 binding is driven predominantly by favorable enthalpy (Figures 5A–C). Although it is not possible to determine stochiometries from the Mg2+ titration, a macroscopic binding constant of 6.95 μM for Mg2+ was estimated (Table 1). Unlike the Ca2+-binding event, which is exothermic, Mg2+ binding to CML15 (Figure 5C) is endothermic; consequently, entropy is the predominant driving force for Mg2+ binding.
Figure 5. ITC analysis of Ca2+ and Mg2+ binding to CML15 (A–C) and CML16 (D–F). Representative isotherms are presented for Ca2+ titration into CML samples in either the absence (A,D) or presence (B,E) of 5 mM Mg2+. Also presented are Mg2+ titrations in a Ca2+-free background (C,F) for CML15 and CML16, respectively. Binding was monitored at 30°C in 25 mM HEPES, 100 mM KCl, pH 7.5. The upper panels show raw data from calorimetric titrations, using 5 μL injections, of recombinant CMLs with 600 μM CaCl2 or 600 μM MgCl2 as noted above. The lower panels present the corresponding integrated binding isotherm modeled for two Ca2+-binding events for CML15 and three Ca2+-binding events for CML16. Each isotherm is representative of a minimum of three technical replicates. Protein was used at the following concentrations (A = 26 μM, B = 26 μM, C = 31 μM, D = 23 μM, E = 19 μM, F = 23 μM).
For CML16, in a Mg2+-free background, we observed three Ca2+-binding sites ranging in affinity from 0.03 μM up to 3.76 μM. For site-1, both entropic and enthalpic values were favorable, whereas the remaining sites were driven either by favorable entropy (site-2) or enthalpy (site-3). In contrast, in the presence of Mg2+, only two Ca2+-binding sites were observed with markedly reduced affinities of 0.91 and 7.09 μM, respectively. When Mg2+ is present, favorable entropy and enthalpy drive Ca2+ binding to the higher-affinity site, whereas the weaker site is enthalpically driven. In the absence of Ca2+, CML16 displays a macroscopic binding constant for Mg2+ of 54 μM (Figure 5F). In contrast to CML15, Mg2+ binding appears to be both entropically and enthapically driven.
CML15- and CML16-Promoter Activity in Arabidopsis Tissues
To gain insight into the tissues and cells types that these CMLs function in, qualitative histochemical CML-promoter::GUS reporter analysis was performed to elucidate the spatial patterns of CML15 and CML16 promoter activity across a representative set of developmental stages from radical emergence to floral maturation. CML15pro::GUS analysis indicates that CML15 promoter activity is very restricted and was detectable only in mature pollen and anther tissue (Figure 6). In marked contrast, CML16pro::GUS analysis revealed strong CML16 promoter activity across most Arabidopsis tissues (Figure 7). CML16pro::GUS activity was detected in the root tip of mature embryos, whole 7- and 10-day-old seedlings including guard cells, and the vasculature of adult leaves and inflorescences (Figures 7A–E). In addition, CML16pro::GUS activity was observed in vasculature tissue of petals, the stigma, and in the vasculature of filaments in mature (~stage 10–13) floral organs (Figure 7F). Quantitative fluorometric CML-pro::GUS analysis using soluble protein extracts from whole seedlings (stage 3), mature leaves (stage 5) and mature floral (stage 7) tissues from CML15pro::GUS and CML16pro::GUS plants corroborated the qualitative promoter activity data (Figure 8). As demonstrated by the histochemical GUS analysis, CML15 promoter activity was only found within stage 10–13 floral tissue extract (Figure 8). Conversely, CML16 promoter activity was detected in all stages assayed for GUS enzyme activity, with the greatest specific activity detected in whole seedling tissue extracts. Taken together, these data suggest a very restricted expression for CML15 in specific male reproductive tissues compared to broad expression of CML16 across various tissues and throughout development. Notably, CML16 is largely, if not exclusively, absent from the floral tissue where CML15 was observed suggesting these CMLs have non-overlapping patterns of expression.
Figure 6. Promoter analysis of CML15 reveals very specific floral activity. Representative images of histochemical CML15pro::GUS reporter activity in Arabidopsis tissues, (A) 7-day-old seedling, (B) 22-day-old rosette, (C) mature silique, (D,F) stage 6–13 flowers, (E) anther at anthesis, and (G) 25 day old roots. Images are representative from a minimum of two independent transgenic lines. Bars = 1 mm (A,C,E,G), 0.1 mm (D), 5 mm (B), 2 mm (F).
Figure 7. Promoter analysis of CML16 reveals broad tissue and developmental-stage activity. Representative images of histochemical CML16pro::GUS reporter activity in Arabidopsis tissues: (A) 1-day-old seedling, (B) 7-day-old seedling, inset: closeup of guard cells, (C) 10-day-old seedling, (D) 22-day-old rosette, (E) primary inflorescence at 25 days, (F) stage ~13 flower, (G) mature roots, (H) mature silique. Images are representative from a minimum of two independent transgenic lines. Bars = 0.2 mm (A,B inset), 1 mm (B,C), 5 mm (D,E), 0.5 mm (F), 2 mm (G,H).
Figure 8. Quantitative fluorometric analysis of CML15pro::GUS and CML16pro::GUS reporter activity across whole seedling, mature leaf, and mature floral Arabidopsis tissues. Specific CML15pro::GUS and CML16pro::GUS activities were assayed as described in the section Materials and Methods. Specific GUS activity (mean ± SE) is shown for each line. Extracts from WT Arabidopsis were assayed as negative controls for background fluorescence and yielded no detectable signal. Annotations 1C, 3C, 12A, 14B, refer to the names of the independent transgenic lines used for analyses.
Despite representing the largest class of Ca2+ sensors in plants, most CMLs remain unstudied. Among the 50 CMLs in Arabidopsis, only a few have been examined to date at any biochemical or physiological level. Given the importance of Ca2+ signaling in plants, elucidating the properties and roles of plant-specific Ca2+ sensors is a necessary step toward a broader understanding of the mechanisms underlying the actions of such a universal second messenger. It has long been hypothesized that the remarkable diversity of Ca2+ sensors in plants, and the presence of plant-specific families such as CMLs, reflect the need for plants as sessile organisms to coordinate cellular responses to various stimuli. While such a hypothesis is difficult to test directly, analysis of CMLs provides insight into how they differ from CaM and other Ca2+ sensors and therefore contribute to cell signaling in distinct ways.
Our examination of two Arabidopsis paralogs, CML15 and CML16, from subgroup 4 indicates that they possess unique biochemical properties in comparison to CaM and other CMLs examined to date. At 74% identity to each other, CML15 and CML16 are closely related and carry the hallmark predicted EF-hands of Ca2+ sensors but the variation within these domains (Figure 1), and our biochemical analyses (Figures 2–5), suggest they likely respond to Ca2+ signals differently than CaM. Multiple lines of evidence indicate that both CML15 and CML16 possess the general biochemical properties expected of Ca2+ sensors. Far-UV CD revealed that both CMLs display notable α-helical character in the absence of Ca2+. CML16 shows a modest increase in α-helical character in the presence of Ca2+, reminiscent of CaM (Martin and Bayley, 1986), whereas helical content for CML15 undergoes only a slight change upon Ca2+ binding. From this data, we cannot conclude that these minor changes in helical content are critical structural responses in these CMLs to Ca2+ binding. Importantly, in the case of CaM, the changes in response to Ca2+ binding reflect predominantly a reorientation of the helices, as opposed to strictly a change in α-helical content (Finn et al., 1995; Zhang et al., 1995). Likewise, Ca2+-CML15 and Ca2+-CML16 CD spectra may indicate similar changes in helical orientation to accommodate Ca2+-binding and expose regions of hydrophobicity that are likely important for target interaction as is the case for CaM. As downstream targets of CML15 and CML16 have not yet been identified, it is unknown at this point whether these differences in structure contribute to properties such as target specificity. The CD data for CML15, where Ca2+ binding appears to have little impact on secondary structure, is reminiscent of CML43 (Dobney et al., 2009) and CML42 (Bender et al., 2014). Interestingly, the conserved CaM, AtCaM1, shows similar behavior when the two N-terminal EF-hands are rendered incapable of binding Ca2+ (Astegno et al., 2016). Conversely, the changes in secondary structure of CML16 suggested by CD analysis are similar to previously reported data for CML37 (Scholz et al., 2014), CML39 (Bender et al., 2013) and CML36 (Astegno et al., 2017) where Ca2+-binding increases helical content. It is not yet clear how these subtle structural distinctions contribute to the functions of different CMLs but it underscores the biochemical complexity of this large family of Ca2+ sensors.
The reversible exposure of hydrophobic regions in response to Ca2+ binding is a defining biochemical feature of CaM and is important for target interaction (Ikura and Ames, 2006; Gifford et al., 2007). One of the most striking differences we observed for CML15 and CML16 in comparison to CaM is that both CMLs display substantial hydrophobic exposure even in the absence of Ca2+ (Figure 3, Supplementary Figure 1). As a consequence, in the presence of Ca2+ and Mg2+, CML15 and CML16 display fold-inductions in ANS fluorescence that are markedly lower than that of Ca2+-CaM (Supplementary Figure 1). Although both CML15 and CML16 show intrinsic hydrophobicity in the absence of Ca2+, importantly, neither undergo a notable increase in hydrophobic exposure in the presence of excess MgCl2 alone. This observation is consistent with the hypothesis that these CMLs function as Ca2+ sensors in a cellular environment where background Mg2+ levels are expected to be much higher than Ca2+. It is interesting that while both CML15 and CML16 show greater intrinsic hydrophobicity in the absence of Ca2+, the magnitude of change in exposed hydrophobicity upon Ca2+ binding is much greater for CaM (Supplementary Figure 1). This suggests that the tertiary structures of the Ca2+-bound CMLs have less exposed hydrophobicity relative to CaM. Despite having higher GRAVY scores, the reduced hydrophobic exposure of Ca2+-bound forms of CML15 and CML16 compared to CaM might in part be accounted for by their lower relative methionine (Met) content given that Met residues account for 46% of the solvent accessible hydrophobic regions of CaM (O'Neil and DeGrado, 1990; Yuan et al., 1999). Determining the structure of the holo- and apo-forms of these CMLs via x-ray crystallography or NMR is needed to conclusively address these questions. Interestingly, CML15 and CML16 show slightly increased hydrophobic exposure in the presence of EGTA/EDTA (Figure 3). While the mechanism underlying this phenomenon is unclear, it is reminiscent of that observed for the C-domain of rabbit skeletal Ca2+ sensor troponin-c which exhibits greater ANS fluorescence in the presence of 0.5 mM EGTA than in 10 mM MgCl2 (Grabarek, 2011). An alternatve Ca2+-chelator (trisodium citrate) did not evoke this response from CML15 or CML16, suggesting that hydrophobic exposure in the presence of EGTA/EDTA might be an artifact of direct interaction with the CMLs, perhaps mimicking in some respect CML-target interaction, although this remains speculative. Regardless, both CML15 and CML16 show greater intrinsic hydrophobicity than CaM in the absence of Ca2+ but retain the ability to respond to Ca2+ signals with conformational changes that increase exposed hydrophobic regions, features expected of Ca2+ sensors. The levels of exposed hydrophobicity, both in apo- and holo-forms of CMLs vary among family members. For example, CML36 resembles CMl15 and CML16 in that only a modest increase in exposed hydrophobicity is observed upon Ca2+ binding (Astegno et al., 2017). In contrast, CML37, CML42, CML43, behave more like CaM, displaying low intrinsic hydrophobic exposure in the apo form and much greater relative increases in response to Ca2+ (Dobney et al., 2009; Bender et al., 2014; Scholz et al., 2014). Interestingly, CML14, binds a single Ca2+ atom and shows no change in exposed hydrophobicity at all (Vallone et al., 2016). This structural diversity among CMLs is consistent with the hypothesis that they play distinct though possibly overlapping roles, predominantly as Ca2+ sensors. The challenge in future work will be to assess the importance of these different structural properties in target binding or other functional contexts.
ITC analysis revealed the Ca2+-binding affinities and thermodynamic parameters of CML15 and CML16 and also shed light on the Mg2+-binding properties of these CMLs. Importantly, both CMLs are able to bind Ca2+ with Kd values in the high nanomolar to low μM range in a mM Mg2+ background (Table 1), consistent with what would be considered physiologically relevant Ca2+ levels (Rizzuto and Pozzan, 2006; Dodd et al., 2010; Steinhorst and Kudla, 2014), thereby supporting their predicted roles as Ca2+ sensors. Various environmental stimuli induce cytosolic Ca2+ fluxes in plant cells that reach into the μM range and microdomains near Ca2+ channels are expected to spike at levels in the tens of μM (Rizzuto and Pozzan, 2006; Steinhorst and Kudla, 2014). Having a repertoire of Ca2+ sensors tuned to a range of Ca2+ signals is likely a key aspect of successful information processing in plants. As might be expected for a large family of Ca2+ sensors, the Ca2+ affinities among CMLs tested to date vary but generally fall within what would be considered a physiological range (nM–μM). The Ca2+ dissociation values for CML15 and CML16 are comparable to those observed for CML42 (Dobney et al., 2009), CML43 (Bender et al., 2014), CML14 (Vallone et al., 2016), and CML36 (Astegno et al., 2017). The presence of both high- and low-affinity Ca2+ sites in CMLs appears to be an emerging pattern (Bender et al., 2014; Astegno et al., 2017) but how these features relate to in vivo function remains unclear. A recent study on a tobacco CML, rgsCaM, revealed a low affinity for Ca2+ in vitro, calling into question whether rgsCaM could function as a Ca2+ sensor in vivo (Makiyama et al., 2016). These authors speculate that target binding might enhance the affinity of rgsCaM for Ca2+ as has been observed for CaM (Gifford et al., 2007).
It is noteworthy that while sequence analysis predicts four EF hands for both CMLs, ITC indicates two and three Ca2+ binding events for CML15 and CML16, respectively. This is reminiscent of the differences between predicted and observed binding sites noted for other CMLs (Dobney et al., 2009; Bender et al., 2014; Vallone et al., 2016). However, we cannot exclude the possibility that, despite extensive dialysis against EGTA, some Ca2+ remained bound to these CMLs throughout the purification process and led to underestimations in binding site numbers. Nevertheless, as discussed above, the variations in the EF-hand loops and nearby regions of these CMLs indicates they likely possess Ca2+ binding properties that differ from those of CaM. Synthetic CaM (SynCaM), a hybrid of mammalian and plant CaM, exhibits two exothermic (enthalpically-driven) binding events followed by two endothermic (entropically-driven) Ca2+-binding events (Gilli et al., 1998). Interestingly, ITC data demonstrates that for CML15, the Ca2+-binding event that occurs first is primarily entropically driven (site 1), and the second Ca2+-binding event is enthalpically driven (Table 2). In the absence of a strong entropic contribution to the first binding event we would expect the interaction with the more exothermic enthalpy to occur first. However, the entropic contribution to site-1 Ca2+ binding is sufficiently strong in the presence of Mg2+ to generate a more negative free energy change than for the second binding event, thereby driving site-1 Ca2+-binding to occur first. Interestingly, the titration of CML15 with Mg2+ generates an endothermic Mg2+-binding event, indicating this event is entropically driven. This phenomenon is unusual, as electrostatic attraction (favorable enthalpy) typically drives ion binding (Linse et al., 1991; Gifford et al., 2013). As a result, positive entropy is responsible for driving Mg2+ binding to CML15. Interestingly, the N-terminal region of the Ca2+-binding domain of NADPH oxidase 5 (NCaBD of NOX5) generated a similar ITC isotherm to that of CML15 titrated with Ca2+ (Linse and Chazin, 1995; Wei et al., 2012). ITC analysis of EF-hand I and EF-hand II of NCaBD suggests that EF-hand II has the higher affinity for Ca2+, triggering conformational changes that promote the binding of Ca2+ to EF-hand I. Like CaM, these Ca2+-binding sites in NCaBD display positive cooperativity (Wei et al., 2012). As proteins with considerable structural similarity to CaM, it is reasonable to speculate that Ca2+ binding to CML15 and CML16 occurs with positive cooperativity, but a more detailed structural analysis is needed to empirically assess this possibility.
Table 2. Thermodynamic parameters and dissociation valuesa,b for Mg2+-binding to CML15 and CML16 in a Ca2+-free background.
In the case of CML16, one of the Ca2+-binding sites was only detectable in the absence of Mg2+. Collectively, the ITC data for CML16 suggests that Mg2+ binding might entirely preclude Ca2+ binding to one of the EF-hands. Without additional structural information, it is unclear how Mg2+ would affect Ca2+ binding to CML16. In vivo, [Mg2+]cyt is likely about 3 orders of magnitude greater than [Ca2+]cyt which is estimated to be in the 100–200 nM range (Waters, 2011), making possible competition or allosteric effects of Mg2+ potentially more relevant. It should also be noted that there may be multiple weak Mg2+ binding sites on these CMLs that led to our underestimating Mg2+ Kd values using a single-site binding model. Weaker Ca2+ affinities in the presence of Mg2+ are common for EF-hand containing proteins based on suggestions that Mg2+ might be a direct competitor for Ca2+ binding to EF-hands (Ohki et al., 1997; Malmendal et al., 1999; Clapham, 2007; Gifford et al., 2007). The isotherms for the Ca2+ titration of CML16 are remarkably similar to those observed for GmCaM4, a soybean CML with 77% identity to mammalian CaM (Gifford et al., 2013). In the case of GmCaM4, Mg2+ stabilizes the closed conformation and acts as a competitive antagonist, decreasing the affinity for Ca2+ by ~6.3-fold. We speculate that Mg2+ fulfills a similar role for CML16 resulting in a lag in saturation for the titration of CML16 with Ca2+, and a decrease in the affinity of CML16 for Ca2+. Interestingly, Mg2+ binding to the N-terminal domain of GmCaM4 has been suggested to increase the affinity of the EF-hands for Ca2+ by folding the domains, and thus paying the energetic costs of restructuring the EF-hands (Gifford et al., 2013). This means that Mg2+ serves as both a competitive antagonist and an allosteric activator for GmCaM4. As mentioned, quantitatively similar changes in α-helical character occur at the level of secondary structure for CML16 in the presence of Ca2+ and in the presence of Mg2+. CD data shows that CML16 exhibits comparable increases in α-helical character in the presence of Ca2+ or Mg2+. Taken together, this data indicates that Mg2+ may play the role of an allosteric activator for CML16 but, as noted above, Mg2+ does not induce conformational changes that lead to hydrophobic exposure for either CML15 or CML16 (Figure 3, Supplementary Figure 2) and thus would not be expected to substitute for Ca2+ in target interaction.
Aside from EF-loop I, which is the most conserved among the predicted EF hands in comparison to CaM (Figure 1), it is difficult to speculate which residue changes in the CMLs are responsible for the differences in metal-binding properties relative to CaM. Moreover, it remains possible that regions outside of the EF hands themselves contribute to structural differences that impact cation binding. Regardless, the picture that emerges is that although both CML15 and CML16 are biochemically distinct from CaM, these CMLs possess high-affinity Ca2+ binding sites and respond to Ca2+ binding, but not Mg2+ binding, with substantial conformational changes that expose hydrophobic regions. The latter property is likely critical in considering these CMLs as bona fide Ca2+ sensors whose postulated roles are to regulate downstream targets. How the biophysical differences from CaM contribute to CML physiological function remains an open question for further study. In future, identification of targets for these CMLs will be important to assess how CML-target interaction is affected by cation binding and vice versa. Moreover, caution is always merited in trying to extrapolate the in vitro biochemical properties of Ca2+ sensors to their behavior under cellular conditions.
The very different patterns of CML15 and CML16 promoter activity we observed suggest that spatial expression of these genes likely overlap little, if at all, in terms of their tissue distribution. The broad pattern of CML16pro::GUS activity across a range of tissues and developmental stages is in marked contrast to that of CML15pro::GUS which was observed exclusively in floral tissue, specifically mature anthers and developing pollen grains. While CML16pro::GUS activity was also detected in floral tissue, it was excluded from the anthers where CML15pro::GUS activity was quite strong. These data suggest that CML15 and CML16 likely play distinct roles in floral reproductive structures. The expression of CML15 in pollen is consistent with the fact that Ca2+ signaling is very important in pollen germination and tube growth and a number of Ca2+ sensors have been implicated in these processes (Steinhorst and Kudla, 2013). The restriction of CML15 promoter activity to anther and pollen tissue is reminiscent of CML39 (Vanderbeld and Snedden, 2007). However, CML39 expression is strongly inducible across a range of tissues by abiotic stress and the hormone methyl jasmonate (Vanderbeld and Snedden, 2007). In contrast, we did not observe any induction of CML15 promoter activity when plants were exposed to different abiotic stress conditions (data not shown) suggesting the role of CML15 is primarily associated with male floral organ development.
It is interesting to note that CML16 promoter activity was strong in guard cells in seedlings but not mature plants, suggesting a role for CML16 in early guard cell development. This pattern is similar to that reported previously for CML38 (Vanderbeld and Snedden, 2007). As the gatekeepers of CO2 entry and water vapor exit, guard cells play key roles in both water relations and carbon fixation by adjusting their turgidity and thereby controlling stomatal pore aperture. Ca2+ is well-established as a key second messenger in guard cell signaling and the strong activity of CML16pro::GUS in these cells suggests that CML16 may be serving as a Ca2+ sensor in these cells during stomatal pore opening and/or closing, both of which involve Ca2+ signals (Murata et al., 2015; Agurla and Raghavendra, 2016). CML16pro::GUS activity was also strong in vascular tissues at all growth stages examined, indicating that CML16 may participate in vascular tissue development or function.
To date, information on CML expression patterns comes predominantly from global transcriptomic analyses although there have been a number of promoter-reporter studies. For Arabidopsis, transcriptome database mining predicted that CaMs are broadly expressed across tissues and developmental stages whereas CMLs cluster into five main expression groups that range from highly restrictive to widespread like the CaMs (McCormack et al., 2005). Likewise, transcriptome analysis for soybean CMLs indicates similar expression-based clustering with many displaying broad tissue distribution and at least 7 CMLs being floral specific (Zeng et al., 2017). Promoter-reporter analysis has corroborated or expanded upon these data for a few Arabidopsis CMLs. For example, CML42 shows broad expression across a wide range of tissues and developmental stages whereas its close paralog, CML43, is restricted to root tips under normal growth conditions (Dobney et al., 2009; Bender et al., 2014). This contrast in expression patterns between paralogs has been observed for other CMLs. CML39 is expressed predominantly in male floral tissue whereas CML38 shows much broader tissue expression (Vanderbeld and Snedden, 2007). These examples are thus reminiscent of the patterns observed in the present study for CML15 and CML16. While our GUS reporter data agree well with the expression patterns for CML15 and CML16 as suggested by public transcriptome databases (e.g., bar.utoronto.ca (Toufighi et al., 2005), several caveats should be kept in mind when interpreting pro::GUS expression data. Although both CML15 and CML16 are intronless genes, we cannot exclude the possibility that some regulatory elements may be absent from the genomic regions we used for analysis of CML promoter activity. In addition, GUS enzyme activity is not a direct indicator of transcription or translation rates or protein stability. Nevertheless, our GUS reporter analysis has provided new insight into CML15 and CML16 tissue- and cell-specific promoter activity and emphasized the sharp differences in spatial expression for these two genes. Although reverse genetic approaches, such as gene-knockout analysis, often provide insight into gene/protein function, we did not observe any phenotypic differences using T-DNA insertional knockout lines for CML15, CML16, or with CML15/16 double-knockout lines (data not shown). Functional redundancy among the large CML family is the most likely explanation for the absence of detectable phenotypes in these mutants. However, it remains likely that under some specific environmental condition the loss of these CMLs would manifest itself in a fitness cost or aberrant phenotype; we simply don't yet know what such conditions would be. Generation of higher order mutants of multiple CML knockouts in future studies may help unravel the role of CML15 and CML16 in plant development.
It is necessary to also point out that the subcellular localizations of CML15 and CML16 have not been empirically demonstrated. However, it seems highly likely that both of these proteins are cytosolic given that (i) neither possesses an N- or C-terminal extension relative to CaM that might function as a localization sorting sequence, and (ii) prediction algorithms such as found at Suba4 (suba.live, Hooper et al., 2017) and Aramemnon (aramemnon.uni-koeln.de/multi-prediction, Schwacke et al., 2003) present consensus predictions for cytosolic localization with high probability.
The breadth of biochemical properties, tissue distributions, and target interactions among Ca2+ sensors, including CMLs, almost certainly reflects both the importance and complexity of Ca2+ signaling in plants. Going forward, future analysis aimed at identifying downstream targets of CML15 and CML16 is needed in order to delineate the Ca2+ signaling pathways and cellular events that these proteins participate in. Target identification will also allow for a comparison with CaM of the structural mechanisms through which these Ca2+ sensors interact with their effectors. It is important to continue investigations into the properties and functions of CMLs as a unique and highly diverse family of Ca2+ sensors as it brings us closer to a more complete picture of Ca2+ signaling in plants and of information processing in general.
WS: designed the experiments and supervised the study; AO, AD, MU, and KM performed the experiments with assistance from WS. All authors analyzed, interpreted, and discussed the results and contributed to writing of the manuscript.
This research was financially supported by Natural Sciences and Engineering Research Council (NSERC) through a Discovery Grant (23888-2013) to WS and a graduate scholarship to AO.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We thank Dr. Hillel Fromm (Univ. Tel Aviv) for the kind gift of pET5a:CaM81. We also thank Dr. Robert Campbell (Queen's) for molecular modeling assistance.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2017.02175/full#supplementary-material
Supplementary Figure 1. Ca2+- and Mg2+-induced changes in exposed hydrophobicity of an evolutionarily-conserved recombinant CaM (CaM81) as demonstrated by ANS fluorescence. Fluorometric scans (430–600 nm) were recorded following the addition of 15 μM CaM to 250 μM ANS where all conditions used ANS buffer (10 mM Tris-Cl, pH 7.5, 100 mM KCl and 1 mM DTT). The y-axis depicts fluorescence in international arbitrary units (IAU). Fluorescence values were analyzed relative to a protein-free ANS control sample that was subtracted from the data as background. Each data set is representative of five experimental replicates.
Supplementary Figure 2. Hydrophobic (phenyl-sepharose) column chromatography of recombinant CML15 (A) and CML16 (B) as visualized by Coomassie-stained SDS-PAGE. Chromatography was performed in the presence of either 1 mM CaCl2, 1 mM CaCl2 and 1 mM MgCl2, 1 mM MgCl2, or 1 mM EGTA, as indicated. All column-binding events were followed with a subsequent wash step using 10 bed-volumes of 0.2 M KCl (“W1-W10”), and a wash step using 10 bed-volumes with 0.4 M KCl (“W11-W20”). “UB” denotes the starting protein sample that was loaded onto the column, “FT” denotes the unbound protein flow-through, and “E” denotes eluted fractions.
Supplementary Figure 3. Homology models of CML15 and CML16 based on a Ca2+-CaM crystal structure. Models were generated using Phyre2.0 software and PyMol molecular graphics system as described in the section Materials and Methods. The structure of Paramecium tetraurelia CaM (UNIPROT identifier P07463) was selected by Phyre2.0 as the CaM with the greatest sequence identity in the structural database for comparison to CML15 (UNIPROT identifier Q9FZ75) and CML16 (UNIPROT identifier Q9LI84). Hydrophobic residues are shaded gray, polar but uncharged residues are pale cyan, charged residues are blue, prolines are magenta, cysteines are green, and Ca2+ atoms are presented as green spheres in the CaM model. The grand average of hydropathy (GRAVY) score was determined as described in the section Materials and Methods and is presented below each model.
Supplementary Table 1. Oligonucleotide primers used for PCR.
Agurla, S., and Raghavendra, A. S. (2016). Convergence and divergence of signaling events in guard cells during stomatal closure by plant hormones or microbial elicitors. Front. Plant Sci. 7:1332. doi: 10.3389/fpls.2016.01332
Akke, M., Forsen, S., and Chazin, W. J. (1991). Molecular basis for co-operativity in Ca2+ binding to calbindin D9k. 1H nuclear magnetic resonance studies of (Cd2+)1-bovine calbindin D9k. J. Mol. Biol. 220, 173–189. doi: 10.1016/0022-2836(91)90389-N
Astegno, A., Bonza, M. C., Vallone, R., La Verde, V., D'Onofrio, M., Luoni, L., et al. (2017). Arabidopsis calmodulin-like protein CML36 is a calcium (Ca2+) sensor that interacts with the plasma membrane Ca2+-ATPase isoform ACA8 and stimulates its activity. J. Biol. Chem. 292, 15049–15061. doi: 10.1074/jbc.M117.787796
Astegno, A., La Verde, V., Marino, V., Dell'Orco, D., and Dominici, P. (2016). Biochemical and biophysical characterization of a plant calmodulin: role of the N- and C-lobes in calcium binding, conformational change, and target interaction. Biochim. Biophys. Acta 1864, 297–307. doi: 10.1016/j.bbapap.2015.12.003
Beckmann, L., Edel, K. H., Batistic, O., and Kudla, J. (2016). A calcium sensor - protein kinase signaling module diversified in plants and is retained in all lineages of Bikonta species. Sci. Rep. 6:31645. doi: 10.1038/srep31645
Bender, K. W., Dobney, S., Ogunrinde, A., Chiasson, D., Mullen, R. T., Teresinski, H. J., et al. (2014). The calmodulin-like protein CML43 functions as a salicylic-acid-inducible root-specific Ca2+ sensor in Arabidopsis. Biochem. J. 457, 127–136. doi: 10.1042/BJ20131080
Bender, K. W., Rosenbaum, D. M., Vanderbeld, B., Ubaid, M., and Snedden, W. A. (2013). The Arabidopsis calmodulin-like protein, CML39, functions during early seedling establishment. Plant J. 76, 634–647. doi: 10.1111/tpj.12323
Boonburapong, B., and Buaboocha, T. (2007). Genome-wide identification and analyses of the rice calmodulin and related potential calcium sensor proteins. BMC Plant Biol. 7:4. doi: 10.1186/1471-2229-7-4
Boyes, D. C., Zayed, A. M., Ascenzi, R., McCaskill, A. J., Hoffman, N. E., Davis, K. R., et al. (2001). Growth stage-based phenotypic analysis of Arabidopsis: a model for high throughput functional genomics in plants. Plant Cell 13, 1499–1510. doi: 10.1105/TPC.010011
Bradford, M. M. (1976). A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal. Biochem. 72, 248–254. doi: 10.1016/0003-2697(76)90527-3
Chiasson, D., Ekengren, S. K., Martin, G. B., Dobney, S. L., and Snedden, W. A. (2005). Calmodulin-like proteins from Arabidopsis and tomato are involved in host defense against Pseudomonas syringae pv. tomato. Plant Mol. Biol. 58, 887–897. doi: 10.1007/s11103-005-8395-x
Chigri, F., Flosdorff, S., Pilz, S., Kölle, E., Dolze, E., Gietl, C., et al. (2012). The Arabidopsis calmodulin-like proteins AtCML30 and AtCML3 are targeted to mitochondria and peroxisomes, respectively. Plant Mol. Biol. 78, 211–222. doi: 10.1007/s11103-011-9856-z
Dobney, S., Chiasson, D., Lam, P., Smith, S. P., and Snedden, W. A. (2009). The calmodulin-related calcium sensor CML42 plays a role in trichome branching. J. Biol. Chem. 284, 31647–31657. doi: 10.1074/jbc.M109.056770
Edel, K. H., and Kudla, J. (2015). Increasing complexity and versatility: how the calcium signaling toolkit was shaped during plant land colonization. Cell Calcium 57, 231–246. doi: 10.1016/j.ceca.2014.10.013
Finn, B. E., Evenäs, J., Drakenberg, T., Waltho, J. P., Thulin, E., and Forsén, S. (1995). Calcium-induced structural changes and domain autonomy in calmodulin. Nat. Struct. Biol. 2, 777–783. doi: 10.1038/nsb0995-777
Garrigos, M., Deschamps, S., Viel, A., Lund, S., Champeil, P., Moller, J. V., et al. (1991). Detection of Ca2+-binding proteins by electrophoretic migration in the presence of Ca2+ combined with 45Ca2+ overlay of protein blots. Anal. Biochem. 194, 82–88. doi: 10.1016/0003-2697(91)90154-L
Gifford, J. L., Jamshidiha, M., Mo, J., Ishida, H., and Vogel, H. J. (2013). Comparing the calcium binding abilities of two soybean calmodulins: towards understanding the divergent nature of plant calmodulins. Plant Cell 25, 4512–4524. doi: 10.1105/tpc.113.113183
Gifford, J. L., Walsh, M. P., and Vogel, H. J. (2007). Structure and metal-ion-binding properties of the Ca2+-binding helix-loop-helix EF-hand motifs. Biochem. J. 405, 199–221. doi: 10.1042/BJ20070255
Gilli, R., Lafitte, D., Lopez, C., Kilhoffer, M., Makarov, A., Briand, C., et al. (1998). Thermodynamic analysis of calcium and magnesium binding to calmodulin. Biochemistry 37, 5450–5456. doi: 10.1021/bi972083a
Grabarek, Z. (2011). Insights into modulation of calcium signaling by magnesium in calmodulin, troponin C and related EF-hand proteins. Biochim. Biophys. Acta 1813, 913–921. doi: 10.1016/j.bbamcr.2011.01.017
Hooper, C. M., Castleden, I., Tanz, S. K., Aryamanesh, N., and Millar, A. H. (2017). SUBA4: the interactive data analysis centre for Arabidopsis subcellular protein locations. Nucleic Acids Res. 45, D1064–D1074. doi: 10.1093/nar/gkw1041
Ikura, M., and Ames, J. B. (2006). Genetic polymorphism and protein conformational plasticity in the calmodulin superfamily: two ways to promote multifunctionality. Proc. Natl. Acad. Sci. U. S. A. 103, 1159–1164. doi: 10.1073/pnas.0508640103
Kelley, L. A., Mezulis, S., Yates, C. M., Wass, M. N., and Sternberg, M. J. (2015). The Phyre2 web portal for protein modeling, prediction and analysis. Nat. Protoc. 10, 845–858. doi: 10.1038/nprot.2015.053
Linse, S., Johansson, C., Brodin, P., Grundstrom, T., Drakenberg, T., and Forsen, S. (1991). Electrostatic contributions to the binding of Ca2+ in calbindin D9k. Biochemistry 30, 154–162. doi: 10.1021/bi00215a023
Makiyama, R. K., Fernandes, C. A., Dreyer, T. R., Moda, B. S., Matioli, F. F., Fontes, M. R., et al. (2016). Structural and thermodynamic studies of the tobacco calmodulin-like rgs-CaM protein. Int. J. Biol. Macromol. 92, 1288–1297. doi: 10.1016/j.ijbiomac.2016.08.016
Malmendal, A., Linse, S., Evenäs, J., Forsén, S., and Drakenberg, T. (1999). Battle for the EF-hands: magnesium-calcium interference in calmodulin. Biochemistry 38, 11844–11850. doi: 10.1021/bi9909288
Marchadier, E., Oates, M. E., Fang, H., Donoghue, P. C., Hetherington, A. M., and Gough, J. (2016). Evolution of the calcium-based intracellular signaling system. Genome Biol. Evol. 8, 2118–2132. doi: 10.1093/gbe/evw139
Martin, S. R., and Bayley, P. M. (1986). The effects of Ca2+ and Cd2+ on the secondary and tertiary structure of bovine testis calmodulin. A circular-dichroism study. Biochem. J. 238, 485–490. doi: 10.1042/bj2380485
O'Neil, K. T., and DeGrado, W. F. (1990). How calmodulin binds its targets: sequence independent recognition of amphiphilic alpha-helices. Trends Biochem. Sci. 15, 59–64. doi: 10.1016/0968-0004(90)90177-D
Ranty, B., Aldon, D., Cotelle, V., Galaud, J. P., Thuleau, P., and Mazars, C. (2016). Calcium sensors as key hubs in plant responses to biotic and abiotic stresses. Front. Plant Sci. 7:327. doi: 10.3389/fpls.2016.00327
Scholz, S. S., Vadassery, J., Heyer, M., Reichelt, M., Bender, K. W., Snedden, W. A., et al. (2014). Mutation of the Arabidopsis calmodulin-like protein CML37 deregulates the jasmonate pathway and enhances susceptibility to herbivory. Mol. Plant 7, 1712–1726. doi: 10.1093/mp/ssu102
Schwacke, R., Schneider, A., Van Der Graaff, E., Fischer, K., Catoni, E., Desimone, M., et al. (2003). ARAMEMNON, a novel database for arabidopsis integral membrane proteins. Plant Physiol. 131, 16–26. doi: 10.1104/pp.011577
Sievers, F., Wilm, A., Dineen, D., Gibson, T. J., Karplus, K., Li, W., et al. (2011). Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol. Syst. Biol. 7:539. doi: 10.1038/msb.2011.75
Toufighi, K., Brady, S. M., Austin, R., Ly, E., and Provart, N. J. (2005). The botany array resource: e-northerns, expression angling, and promoter analyses. Plant J. 43, 153–163. doi: 10.1111/j.1365-313X.2005.02437.x
Vallone, R., La Verde, V., D'Onofrio, M., Giorgetti, A., Dominici, P., and Astegno, A. (2016). Metal binding affinity and structural properties of calmodulin-like protein 14 from Arabidopsis thaliana. Protein Sci. 25, 1461–1471. doi: 10.1002/pro.2942
Vanderbeld, B., and Snedden, W. (2007). Developmental and stimulus-induced expression patterns of Arabidopsis calmodulin-like genes CML37, CML38, and CML39. Plant Mol. Biol. 64, 683–697. doi: 10.1007/s11103-007-9189-0
Wei, C. C., Reynolds, N., Palka, C., Wetherell, K., Boyle, T., Yang, Y. P., et al. (2012). Characterization of the 1st and 2nd EF-hands of NADPH oxidase 5 by fluorescence, isothermal titration calorimetry, and circular dichroism. Chem. Cent. J. 6:29. doi: 10.1186/1752-153X-6-29
Yamniuk, A. P., Nguyen, L. T., Hoang, T. T., and Vogel, H. J. (2004). Metal ion binding properties and conformational states of calcium- and integrin-binding protein. Biochemistry 43, 2558–2568. doi: 10.1021/bi035432b
Yuan, T., Ouyang, H., and Vogel, H. J. (1999). Surface exposure of the methionine side chains of calmodulin in solution. A nitroxide spin label and two-dimensional NMR study. J. Biol. Chem. 274, 8411–8420. doi: 10.1074/jbc.274.13.8411
Zeng, H., Zhang, Y., Zhang, X., Pi, E., and Zhu, Y. (2017). Analysis of EF-Hand proteins in soybean genome suggests their potential roles in environmental and nutritional stress signaling. Front. Plant Sci. 8:877. doi: 10.3389/fpls.2017.00877
Keywords: plant, Arabidopsis, calcium signaling, calmodulin, CML, gene expression, protein biochemistry
Citation: Ogunrinde A, Munro K, Davidson A, Ubaid M and Snedden WA (2017) Arabidopsis Calmodulin-Like Proteins, CML15 and CML16 Possess Biochemical Properties Distinct from Calmodulin and Show Non-overlapping Tissue Expression Patterns. Front. Plant Sci. 8:2175. doi: 10.3389/fpls.2017.02175
Received: 08 August 2017; Accepted: 12 December 2017;
Published: 22 December 2017.
Edited by:Ute Vothknecht, University of Bonn, Germany
Copyright © 2017 Ogunrinde, Munro, Davidson, Ubaid and Snedden. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Wayne A. Snedden, email@example.com