Skip to main content

Shaping mechanisms of metal specificity in a family of metazoan metallothioneins: evolutionary differentiation of mollusc metallothioneins



The degree of metal binding specificity in metalloproteins such as metallothioneins (MTs) can be crucial for their functional accuracy. Unlike most other animal species, pulmonate molluscs possess homometallic MT isoforms loaded with Cu+ or Cd2+. They have, so far, been obtained as native metal-MT complexes from snail tissues, where they are involved in the metabolism of the metal ion species bound to the respective isoform. However, it has not as yet been discerned if their specific metal occupation is the result of a rigid control of metal availability, or isoform expression programming in the hosting tissues or of structural differences of the respective peptides determining the coordinative options for the different metal ions. In this study, the Roman snail (Helix pomatia) Cu-loaded and Cd-loaded isoforms (HpCuMT and HpCdMT) were used as model molecules in order to elucidate the biochemical and evolutionary mechanisms permitting pulmonate MTs to achieve specificity for their cognate metal ion.


HpCuMT and HpCdMT were recombinantly synthesized in the presence of Cd2+, Zn2+ or Cu2+ and corresponding metal complexes analysed by electrospray mass spectrometry and circular dichroism (CD) and ultra violet-visible (UV-Vis) spectrophotometry. Both MT isoforms were only able to form unique, homometallic and stable complexes (Cd6-HpCdMT and Cu12-HpCuMT) with their cognate metal ions. Yeast complementation assays demonstrated that the two isoforms assumed metal-specific functions, in agreement with their binding preferences, in heterologous eukaryotic environments. In the snail organism, the functional metal specificity of HpCdMT and HpCuMT was contributed by metal-specific transcription programming and cell-specific expression. Sequence elucidation and phylogenetic analysis of MT isoforms from a number of snail species revealed that they possess an unspecific and two metal-specific MT isoforms, whose metal specificity was achieved exclusively by evolutionary modulation of non-cysteine amino acid positions.


The Roman snail HpCdMT and HpCuMT isoforms can thus be regarded as prototypes of isoform families that evolved genuine metal-specificity within pulmonate molluscs. Diversification into these isoforms may have been initiated by gene duplication, followed by speciation and selection towards opposite needs for protecting copper-dominated metabolic pathways from nonessential cadmium. The mechanisms enabling these proteins to be metal-specific could also be relevant for other metalloproteins.


Metallothioneins (MTs) constitute a superfamily of genetically polymorphic cysteine (Cys)-rich polypeptides that bind, with high affinity, closed-shell metal ions such as Zn2+, Cd2+, Cu+ and others [1, 2]. In many organisms they play multiple roles. By serving as principal cellular stores for Zn2+ and Cu+ they warrant the supply of these essential trace elements in growth and rescue processes and, by their high sequestration power, they shield cell components from deleterious bonding by highly reactive metals such as Cd2+ and overabundant Zn2+ and Cu+ [3, 4]. In addition, they are also thought to serve as quenchers of free radicals [5, 6]. In most animal species examined different MT isoforms show poor or no appreciable differentiation and specialization in their functions and metal-binding preferences [7, 8], although variations in metal selectivity between MT domains exist in some cases [911]. As a consequence, the metal composition of most native metal-MT complexes is very often remarkably promiscuous [12, 13]. Occasionally, apparent metal specificity results from a disproportional oversupply of a certain metal ion, due to particular physiological conditions such as metabolic trace element disorders [14] or cellular overload due to metal exposure [15]. In such cases, this metal ion occupies all the binding sites of an MT molecule which would otherwise form heterometallic complexes. However, true metal specificity requires an exclusive binding preference for a certain metal to an MT peptide due to its innate structural configuration.

Understanding how MTs and other metalloproteins achieved metal specificity through evolution is a key question in the study of their structure/function relationship [2, 16]. Among molluscs, pulmonate snails provide an optimal system with which to study the determinants of metal MT specificity in metazoans. Molluscs comprise a range of economically, medically and ecologically significant species and represent one of the most successful animal phyla, having been able to colonize nearly every habitat on earth [17]. Some gastropod molluscs - particularly from the subclass of pulmonate snails - feature MT isoforms that can be isolated from their tissues as homometallic complexes with either Cd2+ or Cu+ [1821]. Hence, one isoform isolated from Cd-exposed Roman snails (Helix pomatia) exhibited an exclusive metal complement of six equivalents of Cd2+ per mol of protein [22], whereas another isoform from the same species contained 12 equivalents of Cu+ [23]. It has been proposed that these two pulmonate MT isoforms serve metal specific tasks related to cadmium detoxification [2426] or homeostatic copper regulation [27]. Only recently, a third MT isoform, recovered as a mixed Cd2+- and Cu+-containing complex has been detected in a terrestrial pulmonate but, due to its low abundance, this isoform is probably less important to the snail's metal metabolism [21].

So far, all pulmonate MT isoforms have been obtained as native metal-MT complexes purified from snail tissues, where they are primarily involved in the metabolism of the metal ion species bound to the respective isoform [18, 19]. Therefore, it cannot yet be discerned if their specific metal occupation is the result of a rigid control of metal availability, an isoform expression programming in the hosting cells and tissues or of structural differences of the respective peptides determining the coordinative options for the different metal ions. The aim of the present investigation was, therefore, to use the metal-specific snail MTs as model molecules and to test which of the above determinants contributes to the metal specificity in the snail MT system. To this end, the two metal-specific H. pomatia (Hp) MT isogenes (HpCdMT and HpCuMT) were expressed in two different heterologous environments - bacteria and yeast - in order to test their metal binding and functional capacities independent of their native environments. HpCdMT and HpCuMT were recombinantly synthesized in Escherichia coli in the presence of zinc, cadmium and copper and the features of the formed metal complexes compared by optical, chiroptical and mass spectrometric analyses; both isoforms were compared for their functional competence in yeast MT-knockout cells. Furthermore, cell- and tissue-specific expression of the two MT isogenes in the snail was scrutinized by in situ hybridization techniques and their expression regulation pattern was assessed by real-time polymerase chain reaction (PCR). Finally, the phylogenetic radiation of the gastropod MT family cluster inside the mollusc phylum was assessed and evaluated. Overall, our data provide an unprecedented analysis of the mechanisms determining, at different levels, the specificity of functions of paralogous MTs, suggesting clues to how these could have been achieved through evolution.


Analysis of recombinant HpCdMT and HpCuMT metal complexes reveals sequence-determined specialization of metal binding

Recombinant expression of HpCdMT and HpCuMT in metal-exposed E. coli cells was expected to reveal their in vivo binding ability for zinc, cadmium and copper, independent of which is the natively coordinated metal ion.

Binding of Zn2+ and Cd2+by HpCdMT

Mass spectrometric analysis documents that recombinant synthesis of HpCdMT in E. coli cultured in Cd2+ or Zn2+-enriched media led to the production of only a single species with a fixed content of either six equivalents of cadmium or zinc (Figure 1a and 1b). The two metal complexes - Cd6-HpCdMT and Zn6-HpCdMT - display optical spectra with steep rises of absorbance below 270 nm and below 240 nm which is typical of tetrahedral bonding of both metal ions to multiple thiolate ligands (Figure 2a and 2c). The metal-to-sulphur linkages also manifest themselves by the intense positive and negative circular dichroism (CD) bands associated with the absorption envelopes (Figure 2b and 2d). These signals arise in part from the dissymmetric excitonic interactions of the sulphur-based transitions in pairs of doubly coordinated metal-connecting cysteine residues (of bridging thiolate ligands), thereby signifying the collective bonding of Cd2+ or Zn2+ in oligonuclear metal thiolate complexes [2830]. The spectral features of the recombinant Cd6-HpCdMT product are indistinguishable from those of the native Cd6-HpCdMT previously isolated from the tissue of Cd-exposed Roman snails [23]. They also reappear when, at neutral pH, the full complement of six equivalents of Cd2+ is added in vitro to the metal-free protein apo-HpCdMT [22] or when Zn2+ is replaced in recombinant Zn6-HpCdMT by exposure to the much more firmly binding Cd2+ (Figure 2b). Thus, in both in vivo and in vitro, the structure of the protein product of the HpCdMT gene is seen to direct bonding of Cd2+ in a single energetically favoured complex fitting its supposed role of shielding the snail tissue from this highly toxic metal ion [20, 24].

Figure 1
figure 1

Deconvoluted electrospray ionization time-of-flight mass spectography spectra of the different metal-metallothionein (MT) complexes recombinantly synthesized: (a) HpCdMT obtained in Cd-enriched medium; (b) HpCdMT produced in Zn-enriched medium; and (c) HpCuMT synthesized in Cu-enriched medium under low aeration conditions. Inductively coupled plasma atomic emission spectroscopy (ICP-AES) analysis of these preparations indicated a respective mean content of 6.2 Cd/MT in (a) and 5.8 Zn/MT in (b) for HpCdMT and of 12.2 Cu/MT for HpCuMT in (c).

Figure 2
figure 2

Ultra violet- vis absorption (a, c and e) and circular dichroism (CD (b, d and f) spectra of the following single species recombinantly synthesized in the corresponding metal-enriched media (Cd2+, Zn2+, and Cu2+): Cd 6 -HpCdMT (a, b), Zn 6 -HpCdMT (c, d) and Cu 12 -HpCuMT (e, f). Panel (b) also encloses (thin line) the CD spectra of the Cd6-HpCdMT species prepared in vitro by addition of 6 Cd2+ equivalents to Zn6-HpCdMT. MT, metallothionein.

Binding of Cu+by HpCuMT

The recombinant expression of HpCuMT in E. coli cultures grown in Cu2+-enriched medium led to the formation of homometallic Cu12-HpCuMT as an essentially single molecular species (Figure 1c), equivalent to the native complex purified earlier from snail tissue [23]. The absorption spectra display a progressive rise below 350 nm with a broad shoulder centred at 250 nm (Figure 2e) and in CD positive and negative ellipticity bands (Figure 2f). These features match qualitatively and quantitatively those seen in mammalian Cu12-MT prepared in vitro by adding 12 equivalents of the acetonitrile complex of Cu+ to native MT from rabbit and are attributable to the formation of oligonuclear Cu+ thiolate complexes in trigonal coordination geometry [31]. A molecular species with the same composition and spectral properties was attained in vitro by saturating at acidic pH the recombinant, metal-free, apo-HpCuMT peptide with Cu+ using [Cu(CH3CN)4]ClO4 as a titrating agent. As previously observed [23], this homometallic recombinant product is sensitive to atmospheric O2. The recombinant complex was formed as a single product only when the bacteria were grown under low aeration conditions [32]. At normal oxygenation the same culture produced a heterometallic mixture of several Cu,Zn-HpCuMT species ranging from M4 to M12-HpCuMT (where M = Zn + Cu) and with spectroscopic features clearly different from those of Cu12-HpCuMT. In vitro addition of Cu+ to these species in the form of [Cu(CH3CN)4]ClO4 (see above) also failed to transform these products into the homometallic single form.

Complexes of HpCdMT and HpCuMT with non-cognate metal partners

In contrast to the single, well-defined MT species resulting from the recombinant expression of the HpCdMT and HpCuMT genes in cultures enriched with their cognate metal partners, only poorly defined products were obtained when the partners were interchanged (Figure 3). Thus, expression of the HpCdMT gene in the presence of copper resulted in mixtures of heterometallic Zn, Cu-HpCdMT species of varying metal-to-protein stoichiometry, when grown under normal and low aeration conditions. In the first case, their metal complement varied from three to seven equivalents and in the second from eight to 12 equivalents (Figure 3a). The CD spectra of the mixtures (not shown) displayed signals typical of oligonuclear Cu+-thiolate bonding but varied in shape and amplitudes in different preparations.

Figure 3
figure 3

Electrospray ionization (ESI) time-of-flight mass spectroscopy (MS) spectra of the different metal-metallothionein (MT) complexes obtained when recombinantly synthesizing the MT peptides with their 'non-cognate' metal ion. (a) HpCdMT produced in Cu-enriched medium and; HpCuMT in (b) Cd- and (c) Zn-enriched media, where M stands for Zn+Cu, due to the difficulty of ESI-MS to discriminate between these two metal ions. For the Zn and Cd-enriched preparations, the MS spectra show the +4 charge state (intensity vs m/z) obtained for each preparation, while in the protein produced in a Cu-enriched medium, the MS spectrum shows the +5 charge state (intensity vs m/z). inductively coupled plasma atomic emission spectroscopy (ICP-AES) analysis of these preparations indicated a mean content of 0.8 Zn/MT and 7.7 Cu/MT for HpCdMT in (A); 7.8 Cd/MT in (B) and 4.8 Zn/MT in (C) for HpCuMT.

The recombinant synthesis of HpCuMT in Cd2+-enriched cultures led to mixtures of a number of Cd-HpCuMT species with a Cd content ranging from six to more than nine equivalents and also including a sulphide-containing Cd9S-HpCuMT complex (Figure 3b). The CD spectra of these mixtures (not shown) displayed spectropolarimetric features arising at less than 270 nm from oligonuclear Cd-thiolate complexes and close to 280 nm from sulphide bonding to the Cd-thiolate clusters [33]. In a parallel way, the production of HpCuMT in zinc (Zn)-enriched media yielded mixtures of Zn-containing species ranging from Zn3 to Zn7-HpCuMT, with the different forms varying in abundance and in different preparations (Figure 3c). The CD spectra were indicative of Zn-thiolate coordination but differed widely in intensity (data not shown).

Therefore, these two MTs behave in accordance with their high metal specificity when recombinantly synthesized by bacteria grown in cultures enriched with the non-cognate metal. The HpCdMT isoform thus rendered several mixed Zn, Cu-MT complexes of different stoichiometries when biosynthesized in a Cu-rich medium. Following an equivalent behaviour, the HpCuMT isoform produced a mixture of species of different stoichiometry when synthesized by bacteria grown on Zn- or Cd-supplemented media. The classification of MTs according to the metal-binding behaviour shown when recombinantly synthesized in cultures enriched with different metals and the correspondence of this classification with other 'metal-specificity' criteria, have been fully reviewed [34].

Stability of the Cd6-HpCdMT and Cu12-HpCuMT complexes is documented by their metal exchange inertness

In order to study the lability/inertness of the recombinantly synthesized Cd6-HpCdMT and Cu12-HpCuMT species and their propensity to exchange their preferentially bound metal ions, an equimolar mixture of these two complexes was allowed to stand for 20 h at 25°C (Figure 4). The invariant electrospray ionization-mass spectrometry (ESI-MS) spectra recorded just after mixing and 20 h later demonstrate that the integrity and individuality of these two species was maintained for a long period of time, which confirms that both metal-HpMT complexes possess an exceptionally high stability and exhibit a persistence attributable to their metal binding specificity.

Figure 4
figure 4

Electrospray ionization (ESI) time-of-flight mass spectroscopy (MS) spectra of an equimolar mixture of recombinant Cd 6 -HpCdMT and Cu 12 -HpCuMT, recorded at (a) t = 0 and (b) t = 20 h, showing that no metal exchange occurs between the two HpMT isoforms at 25°C. MT, metallothionein.

Transformation of HpCdMT and HpCuMT in yeast MT-knockout cells confirms metal-specific roles

In order to advance from metal-specific folding to metal-specific function, the particular performance of the two snail MT isoforms was compared by complementation studies in another heterologous system. Hence, yeast cells deficient in their endogenous MTs (yeast Cup1 and Crs5 knockout cells) were transformed with complimentary DNAs (cDNAs) coding for HpCdMT, HpCuMT, mouse MT1, yeast Cup1 and yeast Crs5 and their growth was examined for Cu2+ and Cd2+ tolerance. When these cells were grown at increasing Cu2+ concentrations in the medium (Figure 5a), the highest copper tolerance was observed for the strain transformed with yeast Cup1, followed by the strain transformed with the HpCuMT cDNA. The strain transformed with HpCdMT cDNA gave no evidence of tolerating copper at all.

Figure 5
figure 5

(a) Copper and (b) cadmium tolerance evaluated by phenotype rescue on DTY4 [metallothionein (MT) deficient] yeast cells. Metal tolerance of each DTY4 strain transformed with different MTs, as indicated in the side boxes, is shown as a percentage of the growth rate exhibited in a non-metal-supplemented medium. Cup1 and Crs5 are the two yeast MTs, and MT1 stands for the mouse MT1 isoform.

In marked contrast, when the cells were grown in media with increasing concentrations of Cd2+, tolerance was greatest in the strain transformed with the endogenous yeast Crs5 which reflects the known preference of this yeast MT for divalent metal ions [32]. The second best was the strain transformed with the cDNA for HpCdMT showing a Cd2+ detoxification capacity that was also much better than that of the cells transformed with the cDNA coding for the mouse MT1 isoform, which natively binds either Zn2+, Cd2+ or Cu+. The strain carrying cDNA for HpCuMT showed almost the same high sensitivity to Cd2+ as the MT-null knockout cells (Figure 5b).

These results show that the two snail MT isoforms also assume metal-specific roles in a heterologous eukaryotic environment (yeast), in accordance with their metal-specific binding preferences revealed by their synthesis in recombinant prokaryotic systems. Significantly, the total equivalence between the features of the metal-MT complexes synthesized in these two hosts (bacteria and yeast) has recently been demonstrated for both cadmium and copper, using the Cup1 MT as a model system [35].

In pulmonate snails, CdMT and CuMT isoforms are products of cell-specific expression

In the midgut gland of the snail Helix pomatia, the Cd-specific isoform (HpCdMT) is synthesized in all cell types of this organ (Figure 6a) and is also produced in the epithelial cells of foot, gut and kidney [24]. In contrast, the messenger RNA (mRNA) coding for HpCuMT is located only in one cell type, the so-called rhogocytes (Figure 6b), which are present in the midgut gland, and in many other organs, and have been shown to be the sites of hemocyanin synthesis [27]. Consequently, both metal-specific MT isoforms can be recovered natively from the snail midgut gland (Figure 6c).

Figure 6
figure 6

Cell-specific visualizuation of HpCdMT and and HpCuMT mRNA, and isolation of both native expressed isoforms from snail midgut gland tissue. In situ hybridization (dark violet precipitations) of (a) HpCdMT messenger RNA (mRNA) in midgut gland cells and (b) HpCuMT mRNA in rhogocytes from midgut gland of cadmium-exposed Helix pomatia. (c) Reverse-phase high-performance liquid chromatogram of purified homogenate supernatants from midgut gland of cadmium-exposed H. pomatia snails, showing by the arrows the HpCdMT and HpCuMT isoforms originating from the different cell types, as characterized in references [18] and [19]. MT, metallothionein.

Pulmonate MT isoform genes display metal-specific transcription patterns

The pattern of metal-specific transcriptional induction of MT isogenes was examined in two pulmonate species: in H. pomatia, the subject of this work, and in Cantareus aspersus because, in this species, a third and so far unknown MT isogene (here called Cd/CuMT; Figure 7b) has been reported [21]. The product of this gene also seems to occur in other pulmonate snails, as first reported in this work (see below). The effect of metal supplementation in the feed of the snails upon transcription was evaluated by measurement of the mRNA copy number. For both species, the expression of the CdMT genes was highly responsive to cadmium exposure (Figure 7a and 7b). While Cd2+ increased the number of transcripts of the CdMT genes in both species at concentrations as low as 0.45 μmol Cd2+/g dry weight (in the feed), no statistically significant enhancement was observed in H. pomatia for a more than a 10-times higher amount of Zn2+ (6.93 μmol Zn2+/g dry feed weight; Figure 7a). In Cu2+ a significant increase of the mRNA copies of the HpCdMT gene was seen only at an effective concentration of 5.05 μmol/g dry weight in the feed (Figure 7a). In contrast to the CdMT genes, no significant metal-dependent enhancement of mRNA copy number - at least at the metal concentrations assayed - was observed for the CuMT and Cd/CuMT genes of the two species (Figure 7a and 7b). These induction patterns are totally in accordance with the constitutive expression of the CuMT gene in rhogocytes [27], whereas the much higher inducible expression of the CdMT genes in epithelial cells supports the view that the product of this gene plays a role in Cd2+ sequestration and detoxification [2426]. In addition, CdMT may also serve other biological functions, very likely in the form of the Zn2+ complex for H. pomatia, as suggested and discussed elsewhere [36].

Figure 7
figure 7

Real-time detection polymerase chain reaction (copy number/10 ng total RNA) of messenger RNA (mRNA) of Helix pomatia HpCdMT and HpCuMT (a) and Cantareus aspersus CaCdMT, CaCuMT and CaCd/CuMT (b). Respective mRNA concentrations (copy numbers/10 ng total RNA) were measured in midgut gland tissue of control (unexposed) snails at the beginning of the experiments and of controls, as well as metal-exposed snails, after a feeding period of 5 days. For each bar, means and standard deviations are shown (n = 5). Asterisks above bars designate significant deviations (T-test, P ≤ 0.05) from control animals at the beginning of the experiment. For controls, copy numbers are specified above bars. Respective metal concentrations in the feed are shown below each bar, expressed as μmol metal/g dry weight of feed. MT, metallothionein.

HpCdMT and HpCuMT as prototypes of isoform families which have evolved metal specificity in pulmonate snails by modulation of non-cysteine amino acid positions

HpCdMT and HpCuMT can be considered as prototypes of a series of orthologous genes also present, except from H. pomatia, in other pulmonate snails. Within molluscs, gastropods and pulmonate snails, in particular, have evolved three MT gene subfamilies, two of them comprising isoforms with a homometallic composition [20, 22, 23] and distinct metal binding behaviour and functional specificities for either Cd2+ or Cu+ [18 - 21, 23; and this work]. Figure 8 shows an alignment of the, so far identified, MT sequences from pulmonate gastropods, including the secondarily aquatic species Biomphalaria glabrata. Throughout, two of the three MT isoforms are alignable with and can thus be assimilated to the known H. pomatia isoforms HpCdMT and HpCuMT [18]. The third sequence, the Cd/CuMT isoform first observed in the terrestrial garden snail (Cantareus aspersus), has been recovered from native sources as a simultaneously Cu+ and Cd2+-containing complex [21] and has been identified in other species too (see Figure 8). However, this isogene is transcribed at low constitutive levels, as it is not inducible by metals at all (Figure 4b) and thus hardly detectable at the protein level. It may, therefore, be only of circumstantial importance for the metal metabolism in its host. However, its discovery is critical for the understanding of the diversification of MTs in this group of organisms.

Figure 8
figure 8

Amino acid sequence alignments of terrestrial and freshwater pulmonate snail metallothioneins [MTs; CdMT, CuMT and Cd/CuMT]. Protein or nucleotide sequences were obtained from GenBank and are specified in the legend of Figure 9. Conserved cysteine (cys) positions shared by all isoform families are shaded in red. Conserved non-cys amino acid positions are shaded as follows: blue, shared by all three isoform families; yellow, shared by the CdMT and the CuMT families; pink, shared by the CdMT and Cd/CuMT families; green, shared by the Cd/CuMT and the CuMT families. Asterisks between two sequences indicate identity; Cys residues are marked in bold. For abbreviations of species names see legend of Figure 9.

The three isoform types share strictly conserved Cys positions in their primary structure, confirming the fundamental importance of the sulphur atoms provided by these residues for metal complex and metal thiolate cluster formation, irrespective of the metal species involved. Besides Cys, a few other amino acid positions, either on the N-terminal tail of the peptides or in the direct neighbourhood of Cys residues, show conserved positions through the members of the three isoform subfamilies (Figure 8). In contrast, there is significant variability across the three isoform types for the non-cysteine amino acid residues interspersed between the conserved positions. This suggests that the different metal specificities of the isoforms were achieved by gene duplication and subsequent speciation by evolutionary modulation of these non-coordinating amino acid positions. Moreover, the alignment pattern shows that the similarity between the members of the CuMT and the Cd/CuMT isoform families is clearly higher than that observed between those and the CdMT isoforms (Figure 8, Table 1).

Table 1 Comparative protein Blast analysis of pulmonate metallothionein (MT) isoform families*

A nucleotide-based neighbour-joining tree shows that pulmonate MT isoform subfamilies are assembled in three distinct branches and are thus clearly distinguishable from all other mollusc forms represented by the group of Bivalvia (Figure 9). This suggests that the differentiation into these isoforms has been an evolutionary process which, within molluscs, remained restricted essentially to pulmonate snails. The protein distance analysis tree (Figure 10) confirms the close relationship between the CuMT and Cd/CuMT isoforms (Table 1), which apparently evolved from a common ancestor that gave rise to the differentiation of the Cu-specific and the less metal-specific Cd/CuMT gene subfamilies, clearly segregated from the CdMT gene subfamily. On the other hand, the three pulmonate MT subfamilies share a common root with all other gastropod MTs (Figures 9 and 10).

Figure 9
figure 9

Nucleotide-based neighbour-joining tree of mollusc metallothioneins (MTs) comprising mussels and gastropods, with Eisenia foetida MT used as out-group. Pulmonate MTs appear grouped in separate clusters of metal-specific subfamilies (CuMTs, Cd/CuMTs and CuMTs). Accession numbers of GenBank entries used were as follows: Helix pomatia CdMT (HpCdMT), AAK84863 and AF399740; H. pomatia CuMT (HpCuMT), AAK84864 and AF399741. Biomphalaria glabrata CdMT, ACS91928 and GQ205374; B. glabrata CuMT, ACS91927 and GQ205373; B. glabrata Cd/CuMT, ACS91929 and GQ205375; Cantareus aspersus CdMT, ABL73910 and EF152281; C. aspersus CuMT, ABM55268 and EF178297; C. aspersus Cd/CuMT, ABM92276 and EF206312; Nesiohelix samarangae CdMT, ACC17831 and EU437399; Megathura crenulata MT, AY102647; Littorina littorea MT, AY034179; Mytilus edulis MT-10, AJ007506 and EF140765; Crassostrea virginica MT-1A, AY331697; C. virginica MT-1B, AY331699; Meretrix lusoria MT, AY525635; Perna viridis MT-2, F092972; Dreissena polymorpha MT, DPU67347; Unio tumidus MT, EF185127; Corbicula fluminea MT, EF185126; Eisenia foetida MT, AK236886.

Figure 10
figure 10

Amino acid sequence-based neighbour-joining tree of pulmonate metallothionein (MT) subfamilies represented in clusters of Cd/CuMTs, CuMTs and CdMTs, with Megathura crenulata MT as out-group . Accession numbers of GenBank entries are as indicated in the legend of Figure. 9.


In MTs, metal binding and metal exchange reactions are mainly governed by the coordination chemistry of thiolate bonding with closed-shell metal ions such as Zn2+, Cd2+, Cu+ [12]. To be more precise, the relative order of in vitro metal binding affinities of apo-MT peptides, as well as the order of displacement capacity of each heavy metal ion within a metal-MT complex (Hg(II) > Cu(I) ~ Ag(I) > Cd(II) > Pb(II) > Co(II) > Zn(II)), follow the rules established for metal-thiolate and metal-sulphide low-molecular mass complexes [37]. However, the assumption that these 'inorganic chemistry' rules are the unique responsible of the metal-MT complex properties would lead to the conclusion that MT polypeptides sharing the same number and position of Cys residues would exhibit equivalent metal binding behaviour. However, this is essentially untrue, as firmly demonstrated in this work for the snail MT system. MT metal specificity is a subject of vivid, current debate [16], as to a larger or lesser extent, all MTs show a degree of metal specificity in native and/or in vivo environments. In this case, metal specificity is understood to be the set of determinants that eventually leads a given MT peptide to natively discriminate among metal ions, thus allowing the formation of particular metal complexes and the performance of the biological function for which they were selected. To date, major factors claimed to explain metal-MT specificity in live systems are protein sequence optimization, metal ion availability and/or metal-responsive transcription programming [2].

The structural features of the resulting MT complexes with different metal ion species arise from the equilibrium between kinetic and thermodynamic requirements [29], so that they converge to a stoichiometric ratio that reflects their energetically most stable state [38]. Hence, the observation by ESI-MS of metal-MT species synthesized in in-vivo environments, either native or recombinant, allows one to read the propensity of the respective peptides to form metal complexes that are uniquely defined from a stoichiometric and thermodynamic point of view [2, 23, 34]. In combination with spectroscopic studies, this leads to a clear appraisal of the distinctness of a MT metal specificity, as exemplified in the present work. It can, therefore, be concluded that the homometallic and unique composition of the complexes formed by HpCdMT and HpCuMT with their cognate metal species upon isolation from recombinant cultures (Figure 1) reflects the innate metal specificity of the two isoforms, rather than being the result of an occasional association with metal ions governed by their intracellular or environmental availability. This is confirmed by the poor metal-binding behaviour of the same isoforms confronted with their non-cognate metal ions (Figure 3) and by the inertness of the two complexes, Cd6-HpCdMT and Cu12-HpCuMT, to exchange metal ions (Figure 4). Such metal ion exchange processes have repeatedly been reported in MTs [39] and would be reasonable to expect [12].

The complete sequential identity of Cys residues and the high degree of conserved positions for other amino acids shared among the three isoform families (Figure 8) suggest, along with their nearly equal size, that metal-specific differentiation of pulmonate snail MTs must have been initiated by gene duplication events, followed by modulatory speciation of amino acid residues located between the cysteine positions. Gene duplication seems to be a common mechanism driving the evolutionary differentiation of MT isoform in animals and is documented for MTs of Drosophila melanogaster, among others [40], and the mussel Crassostrea gigas [41]. Once duplicated, such genes are free to independently generate mutations, upon which selective forces can then act towards evolution of specific features [42]. The example of pulmonate snail MTs also proves that the evolutionary variation of non-cysteine residues can impose a metal-specific character on to the coordination chemistry of a MT peptide. At present, it is not known how this is achieved at a structural level. However, it must be supposed that, due to their particular position in the sequence and chemical nature of their side-chains, non-cysteine amino acids constrain the sulphur ligands provided by the conserved Cys positions to assume only one of several theoretically possible spatial coordination foldings. Determination of the three dimensional structure of the Cd6-HpCdMT and Cu12-HpCuMT complexes is actually in progress which may unveil the detailed structural basis of the metal specificity of the two HpMT isoforms.

Cell-specific expression may also contribute to enhancing distinct metal-related functionality [2]. Roman snails, for example, synthesize HpCuMT exclusively in rhogocytes (Figure 6b), the modified cells of mesodermal origin found in virtually all connective tissues of mollusc organs [43]. Since they are also the sites of hemocyanin synthesis [44], it was suggested that HpCuMT functions as a Cu+ reservoir/donor for the nascent hemocyanin [27]. The constitutive expression of HpCuMT (Figure 7) and the exclusive preference of this isoform for Cu+ (Figure 6c) support this presumed function. The supposed incorporation of Cu+ into the structure of hemocyanin must occur under reducing conditions, which is also consistent with the high susceptibility of native Cu12-HpCuMT complexes to oxidization [20, 23] and the fact that homometallic Cu12-HpCuMT synthesis is only achieved in low-aerated recombinant cultures (this work). The apparent connection between the tasks of Cu-specific MT isoforms and their presence in organisms with Cu-depending hemocyanins is reminiscent of the situation reported for decapod crustaceans [45, 46]. In these animals, concentrations of Cu-MT complexes fluctuate with the metabolic state and the hemocyanin levels during the moulting cycle [47, 48]. Cu-specific MTs are also observed in organisms of other kingdoms, especially in fungi [4951], where their role may be connected to the synthesis of the Cu-containing enzyme tyrosinase, as in Neurospora crassa [52].

After exposure of Roman snails to Cd2+, virtually all of this metal in the digestive tissues was bound to HpCdMT (20) in a similar manner as for CdMT isoforms of other snail species [53]. Consistently, in the Roman snail, HpCdMT is produced in digestive and excretory tissues [24] (Figure 6a), where the corresponding gene is selectively upregulated by Cd2+ exposure (Figure 7). This suggests that absorption of toxic Cd2+ from the surrounding substrate via the alimentary tract may constitute a particular physiological challenge, exacerbated by the evolutionary transition of gastropods to terrestrial life [36]. Moreover, the sensitivity to Cd2+ of important Cu-dominated metabolic pathways [5456] and Zn-dependent enzymes [57] could have been the basis of the generation of a specific MT isoform devoted to Cd2+ detoxification in these animals. Our data also demonstrate the ready formation of homometallic complexes of HpCdMT with Zn2+, which may be a consequence of the comparable coordination preferences of these two d10 metal ions. The much weaker bonding of Zn2+ [58] to this isoform, however, does not prevent the peptide from functioning as a most effective Cd2+ sequestration agent. In the presence of Zn2+ and the absence of Cd2+, the HpCdMT isoform is expressed only at low basal concentrations (Figure 7) and, as in the case of mammalian MTs [59], is thought to serve other functions [36].


Overall, the present study, together with the extensive body of evidence provided by our previous work, suggests that the pair of the metal-specific H. pomatia MT isoforms (HpCdMT and HpCuMT) can be regarded as the prototype of a series of paralogous forms present in pulmonate gastropods. In these organisms, divergent evolution by gene duplication, with subsequent modulation of non-cysteine amino acid positions and a cell specific occurrence and gene expression regulation, has led to the complete separation of their metal-binding preference, cell-specific occurrence, expression regulation and functionality. This resulted in genuine CdMTs becoming inducible forms specializing in the global protection of the organism from the non-essential toxic element Cd and in genuine CuMTs becoming constitutive forms supplying the essential element Cu. Our findings provide experimental evidence and possible answers to how metallotproteins in general, and MTs specifically, were able to achieve partial or complete specificity in their metal binding behaviour and functionality.


Animals and metal exposure

Roman snails (H. pomatia L.) were obtained from a commercial dealer (Exoterra, Dillingen, Germany). Garden snails (C. aspersus) were provided by the Department of Chrono-Environment, University of Franche-Comté, Besançon, France. All animals were reared under laboratory conditions (20°C, 80% humidity, 12:12 h photoperiod) at the Institute of Zoology in Innsbruck, Austria. Twenty-five snails from each species were split equally into five groups and fed over a period of 5 days on metal-enriched lettuce (Lactuca sativa). Metal loading of feed was achieved by soaking salad leaves in a corresponding metal salt solution (CdCl2 in H2O, with 1 and 3 mg Cd2+ L-1; ZnCl2 in H2O, with 5 and 10 mg Zn2+ L-1; CuCl2, with 10 mg Cu2+ L-1) [60]. Resulting metal ion concentrations in the salad feed were as follows (means ± standard deviation, n = 5): Cd2+, 0.45 ± 0.11 and 1.12 ± 0.23 μmol g-1 dry weight; Cu2+, 3.51 ± 0.73 or 5.05 ± 0.97 μmol g-1 dry weight; Zn2+, 4.97 ± 3.42 and 6.93 ± 0.89 μmol g-1 dry weight). These concentrations range from physiologically to moderately elevated levels and are, therefore, representative for what could be the natural conditions encountered by snails. At days 0 and 5, RNA was extracted from the small midgut gland tissue aliquots (~10 mg fresh weight) of at least three animals and processed for cDNA synthesis as detailed below.

For in-situ-hybridization of HpCdMT isoform mRNAs, five individuals of H. pomatia were exposed over 14 days to a concentration of 14.97 μmol Cd g-1 dry weight. At the end of the exposure period, animals were sacrificed and their organs (midgut gland, midgut, kidney, mantle and foot) used for in-situ-hybridization analysis as described below.

Metal analyses

Metal-enriched salad samples were oven-dried at 60°C for several days. Dried samples (snail tissues: 50-100 mg dry weight; salad samples: 100 - 500 mg dry weight) were wet-digested in screw-capped polypropylene tubes (Greiner, Kremsmünster, Austria) with a mixture of HNO3 (suprapure; Merck, Darmstadt, Germany) and distilled water (1:1) by heating at 70°C for several days. At the end of digestion, a few drops of H2O2 were added to the heated samples. The remaining solutions were diluted to a known volume with distilled water and analysed for metal concentrations (Cd, Zn, Cu) either by flame (model 2380 instrument, Perkin Elmer, Massachusetts, USA) or graphite furnace atomic absorption spectrophotometry (Hitachi Z-8200) with polarized Zeeman background correction (Hitachi, Tokyo, Japan).

Real-time detection PCR

RNA sampling for real-time detection PCR was done in control snails and animals exposed to metals over a 5 day period (see above). This time range was chosen because, in pulmonate snails, maximal induction of the CdMT gene by Cd2+ is reached only after several days. Total RNA was isolated from the homogenized midgut gland tissue of H. pomatia and C. aspersus individuals (Ultra Turrax T25, IKA Maschinenbau, Staufen, Germany) and quantified after DNaseI digestion (Fermentas, St Leon-Rot, Germany) by means of RiboGreen® RNA Quantitation Kit (Molecular Probes, OR, USA) with calibration curves derived from RNA standards using a fluorescence plate reader (Molecular Devices, CA, USA). Of total RNA, 450 ng was applied for cDNA synthesis using RevertAid™ H Minus M-MuLV Reverse Transcriptase (Fermentas) with hexamer primers in a 50 μL approach. Quantification of the RNA copy number was performed on a 7500 real-time PCR (RT-PCR) instrument from Applied Biosystems (CA, USA) using the Power SYBR Green approach (Applied Biosystems). Calibration curves from amplicon plasmids were used for copy number analysis for each of the MT isoforms involved, using primers designed with the Primer Express 3.0 software (Applied Biosystems) based on the known cDNA sequences for MT isoforms published in GenBank (H. pomatia HpCdMT and HpCuMT, accession numbers: AF399740 and AF399741, respectively; C. aspersus CaCdMT, CaCuMT and CaCd/CuMT, accession numbers: EF152281, EF178297, and EF206312). PCR primers used were as follows: HpCdMT: sense primer, 5'AAAGTGCACCTCAGCTTGCA 3'; antisense primer: 5' GCAGGCGGCACA TGTACAG 3'; amplicon length, 85 bp. HpCuMT: sense primer, 5' CCTTGCAGCTGTGGT AACGA 3'; antisense primer, 5' CAAGAACTGCATCGGTCACAA 3'; amplicon length, 65 bp; CaCdMT: sense primer, 5' GCCGCCTGTAAGACTTGCA 3'; antisense primer: 5' CACG CCTTGCCACACTTG 3'; amplicon length, 56 bp. CaCuMT: sense primer, 5' AACAGCAACCCTTGCAACTGT 3'; antisense primer, 5' CGAGCACTGCATTGATCACAA 3'; amplicon length, 74 bp. CaCd/CuMT: sense primer, 5' TGTGGAGCCGGCTGTTCT 3'; antisense primer, 5' CAGGTGTCATTGTTGCATTGG 3'; amplicon length, 59 bp. Optimal primer concentrations were determined by means of dissociation curves established for each primer pair. Two microlitres of cDNA were applied for RT detection PCR in a 20-μL approach (1x Power SYBR Green PCR Mastermix, 1x U-BSA, 900 mM sense primer, 300 mM antisense primer for HpCdMT and HpCuMT; 300 mM sense and 900 mM antisense primer for CaCdMT; 900 mM for sense and antisense primer for CaCuMT; 99 mM sense and 300 mM antisense primer for CaCd/CuMT). The PCR conditions were as follows: 50°C, 2 min; 95°C, 10 min; 40 repeats of 95°C, 15 s; and 60°C, 1 min.

In situhybridization techniques

Cell- and tissue-specific expression of both HpMT isoforms was demonstrated by in situ hybridization (ISH). Construction of digoxigenin-11-UTP-labelled sense and antisense RNA probes for ISH of both MT isoform mRNAs, as well as ISH, antibody exposure and staining of parafomaldehyde-phosphate buffered saline (PBS)-fixed paraffin sections (5 μm) from tissues (midgut gland, midgut, kidney, mantle and foot) of control and metal-exposed animals were performed exactly as described previously [24]. Control sections (exposed to either hybridization antisense or sense probes) were treated and incubated in the same way as samples but without anti-digoxigenin-alkaline phosphatase antibodies. For microscopy, all sections were embedded in Entellan (Merck) [24].

Construction of the recombinant expression vectors for wild-type Roman snail MT isoforms

The H. pomatia coding regions for both MT isoforms were amplified by PCR using the respective cDNAs synthesized during a previous investigation as a template [24]. In order to facilitate their in frame cloning into the pGEX-4T1 expression vector (Amersham GE Healthcare Bio-Sciences AB, Uppsala, Sweden), BamHI and SalI restriction sites were generated just before the anti-thymocyte globulin (ATG) and after the stop codon. The oligonucleotides used for these PCR amplifications were: 5' ACAGGATCCGGACGAGGAAAGAACTGC 3' and 5' ATTGGATCCGGGAAAG GAAAAGGAGAAAAGTG 3' as HpCuMT and HpCdMT upstream primers, and 5' AGGCGTCGACTTGTCGTTTATTTGCAG 3' and 5' ATGCGTCGACTTGTCCTGC GGTTACT 3' as the HpCuMT and HpCdMT downstream primers. 35-cycle PCR reactions were performed under the following conditions: 94°C 30 s, 55°C 30 s and 72°C 30 s, using Deep Vent (New England Biolabs, Massachusetts, USA) thermostable DNA polymerase. PCR products were isolated from 2% agarose gels, digested with BamHI-SalI restriction enzymes (New England Biolabs) and cloned into the corresponding sites of pGEX-4T-1, for glutathione-S-transferase (GST)-MT fusion protein synthesis. The product resulting after purification was used in the second PCR as reverse megaprimer together with the forward primer mentioned earlier. In the final amplification product, the desired mutation had been introduced and the flanking restriction sites (BamHI and SalI) allowed the cloning in frame in pGEX-4T-1. Prior to the protein synthesis assays, all the DNA constructs were confirmed by automatic DNA sequencing (ABI 370, Perkin Elmer Life Sciences), using BigDye Terminator (Applied Biosystems). DH5α was the E. coli host strain used for cloning and sequencing purposes and, thereafter, the expression plasmids were transformed into the E. coli protease-deficient strain BL21 for recombinant protein overexpression.

Recombinant synthesis and purification of the metal-HpMT complexes

All HpMT metal complexes analysed in this work were biosynthesized in 2-L Erlenmeyer cultures of the corresponding transformed E. coli cells grown in LB medium with 100 mg mL-1 ampicillin and the following metal supplements: 300 μM ZnCl2 or CdCl2 for the zinc- or cadmium-rich media, or 500 μM CuSO4 for the copper-rich medium. Copper cultures were performed under two aeration conditions (high and low aeration) as previously described [32]. GST-MT synthesis was induced with isopropyl-1-thio-β-D-galactopyranoside at a final concentration of 100 mM 30 min before the addition of the metal solution. After a 2.5 h-induction, cells were harvested by centrifugation. In order to prevent oxidation of the metal-HpMT complexes, argon was bubbled in all the steps of the purification following cell disruption.

For protein purification, cells were re-suspended in ice-cold PBS (1.4 M NaCl, 27 mM KCl, 101 mM Na2HPO4, 18 mM KH2PO4)-0.5% v/v β-mercaptoethanol, disrupted by sonication and centrifuged at 12,000 g for 30 min. The recovered supernatant was used to purify the GST-HpMT polypeptides by batch affinity chromatography with glutathione sepharose 4B (GE Healthcare, Buckinghamshire, UK) incubating the mixture with gentle agitation for 60 min at room temperature. After three washes in PBS and, since the GST-HpMT fusions include a thrombin recognition site, this protease was added (10 μ per mg of fusion protein) and digestion was carried out overnight at 23°-25°C. This allowed separation of the GST fragment of the fusion proteins, which remained bound to the gel matrix from the metal-HpMT portions that were eluted together with thrombin. Therefore, the eluate was concentrated using Centriprep Concentrators (Amicon; Millipore, MA, USA) with a cut-off of 3 kDa and subsequently fractionated using fast protein liquid chromatography (FPLC), through a Superdex-75 column (GE Healthcare) equilibrated with 50 mM Tris-HCl, pH 7.0, and run at 1 mL min-1. Fractions were collected and analysed for protein content by their absorbance at 254 nm. Aliquots of the protein-containing FPLC fractions were analysed by 15% SDS-PAGE and stained by Coomassie Blue. HpMT-containing samples were pooled and stored at -70°C until further use. Due to the specific recombinant expression conditions, the three synthesized snail MT isoforms contained one additional amino acid residue (G) at their N-termini in relation to the native isoforms previously isolated [20]. These modifications do not interfere with the metal-binding capacity, as previously shown for both vertebrate [61] and invertebrate [8] MT isoforms.

Analysis of recombinantly expressed and in vitroprepared metal-HpMT complexes

The recombinantly expressed metal-MT complexes were analysed for element composition (S, Zn, Cd and Cu) by inductively coupled plasma atomic emission spectroscopy (ICP-AES) on a Polyscan 61E spectrometer (Thermo Jarrell Ash Corporation, MA, USA) at appropriate wavelengths (S, 182.040 nm; Zn, 213.856 nm; Cd, 228.802 nm; Cu, 324.803 nm), either under 'conventional' (dilution with 2% HNO3 (v/v)) or under 'acidic' (incubation in 1 M HCl at 65°C for 5 min) conditions [62]. MT concentration in the recombinant preparations was calculated from the acidic ICP sulphur measurements, thus assuming the only contribution to their S content was that made by the HpCuMT and HpCdMT peptides. Protein concentrations were confirmed by standard amino acid analysis performed on an Alpha Plus Amino acid Autoanalyzer (Pharmacia LKB Biotechnology, Uppsala, Sweden) after sample hydrolysis in 6 M HCl (22 h at 110°C). Ser, Lys and Gly contents were used to extrapolate sample concentrations.

CD spectroscopy was performed using a model J-715 spectropolarimeter (JASCO, Gross-Umstadt, Germany) equipped with a computer (J-700 software, JASCO). Measurements were carried out at a constant temperature of 25°C maintained by a Peltier PTC-351 S apparatus (TE Technology Inc, MI, USA). Electronic absorption was measured on an HP-8453 diode-array ultra violet (UV)-vis spectrophotometer (GMI Inc, MN, USA), using 1-cm capped quartz cuvettes, and correcting for the dilution effects by means of the GRAMS 32 software (Thermo Fisher Scientific Inc, MA, USA).

Molecular mass determination was performed by electrospray ionization mass spectrometry equipped with a time-of-flight analyser (ESI-TOF MS) using a Micro Tof-Q Instrument (Brucker Daltonics GmbH, Bremen, Germany) calibrated with NaI (200 ppm NaI in a 1:1 H2O: isopropanol mixture), interfaced with a Series 1100 HPLC pump (Agilent Technologies, CA, USA) equipped with an autosampler, both controlled by the Compass Software. The experimental conditions for analysing proteins with divalent metals (Zn, Cd) were: 20 μL of the sample were injected through a PEEK long tube (1.5 m × 0.18 mm i.d.) at 40 μL/min under the following conditions: capillary-counterelectrode voltage, 5.0 kV; desolvation temperature, 90-110°C; dry gas 6 L/min. Spectra were collected throughout an m/z range from 800 to 2000. The proteins that contain copper were analysed injecting 20 μL of the sample at 30 μL/min; capillary-counterelectrode voltage, 4.0 kV; desolvation temperature, 80°C; m/z range from 800 to 2000. The liquid carrier was a 90:10 mixture of 15 mM ammonium acetate and acetonitrile, pH 7.0. For the analysis at acidic pH the conditions used were the same as those used in the analysis of the case for divalent metals, except in the composition of the carrier liquid which, in this case, was a 95:5 mixture of formic acid and acetonitrile at pH 2.4. All samples were injected at least in duplicate to ensure reproducibility. In all cases, molecular masses were calculated according to the reported method [63].

Metal replacement titrations were performed by adding the corresponding metal ions (Cd2+ or Cu+) at equivalent molar ratios to the recombinant Zn-HpMT complexes. Titrations were carried out following previously described procedures [64, 65]. The resulting in vitro complexes were analysed by UV-Vis and CD spectroscopy as well as mass spectrometry. All assays were carried out in an Ar atmosphere and the pH for all experiments remained constant throughout, without the addition of any extra buffers.

Metal tolerance complementation assays in transformed yeast MT-knockout cells

The Saccharomyces cerevisiae DTY4 strain (MATα, leu2-3, 112his3Δ1, trp1-1, ura3-50, gal1, cup1::URA3) was used for metal tolerance complementation assays. This strain is characterized by a total MT deficiency due to cup1 disruption and Crs5 truncation [66].

The cDNAs coding for the different MTs assayed - the two snail MT isoforms (snail HpCdMT and snail HpCuMT), the two yeast MTs (Cup1 and Crs5) and the mouse MT1 - were ligated into the BamHI/PstI sites of the yeast vector p424, which contains TRP1 for selection, the constitutive GPD (glyceraldehyde-3-phosphate dehydrogenase) promoter for heterologous gene expression, and the CYC1 (cytochrome-c-oxidase) transcriptional terminator [67]. The recombinant p424 vectors were introduced into the DTY4 cells using the lithium acetate procedure [68]. Transformed cells were selected according to their capacity to grow in synthetic complete medium (SC) without Trp and Ura.

For metal tolerance tests, transformed yeast cells were initially grown in selective SC-Trp-Ura medium at 30°C and 220 rpm until saturation. These cells were then diluted to OD600 0.01 and used to re-inoculate tubes with 3 mL of fresh medium supplemented with CuSO4 added at 15, 30, 45, 60, 75, 90 and 105 μM final concentrations or CdCl2 at 1, 2, 4, 6 and 8 μM final concentrations. These cultures were allowed to grow for 18 h and the final OD600 was recorded and plotted as a percentage of the OD600 reached by the culture grown without metal supplement. For each concentration, and each kind of transformation, two replicas were run.

MT sequence alignment and phylogenetic analyses

MT amino acid sequences used for pulmonate MT alignments were derived mostly from amino acid sequence files or translated cDNA open reading frame sequences published in GenBank (; see accession numbers in the legend to Figure 9). The editing and alignment were done manually in combination with ClustalX software Version 2.0.9 [69]. For phylogenetic analyses, nucleotide sequences of the coding region of mollusc MT cDNAs or MT genes as well as protein primary sequences were used as published in GenBank (for accession numbers see legend of Figure 9). Phylogenetic reconstructions we performed with neighbour-joining [70] using the computer program PAUP* (version 4). The robustness of the phylogenetic hypothesis was tested by bootstrapping (1000 replicates) [71].



anti-thymocyte globulin


circular dichroism


complimentary DNA


electrospray ionization


fast protein liquid chromatography


glyceraldehyde-3-phosphate dehydrogenase




Helix pomatia


inductively coupled plasma atomic emission spectroscopy


in situ hybridization


glyceraldehyde-3-phosphate dehydrogenase


messenger RNA


mass spectrometry




phosphate buffered saline


real time polymerase chain reaction


synthetic complete medium.


  1. Sigel A, Sigel H, Sigel RKO: Metal ions in Life Sciences; v 5: Metallothioneins and Related Chelators. 2009, Cambridge: Royal Society of Chemistry, 1-514. full_text.

    Book  Google Scholar 

  2. Blindauer CA, Leszczyszyn OI: Metallothioneins: unparalleled diversity in structures and functions for metal ion homeostasis and more. Nat Prod Rep. 2010, 27: 720-741. 10.1039/b906685n.

    Article  CAS  PubMed  Google Scholar 

  3. Klaassen CD, Liu J, Choudhuri S: Metallothionein: An intracellular protein to protect against cadmium toxicity. Ann Rev Pharmacol Toxicol. 1999, 39: 267-294. 10.1146/annurev.pharmtox.39.1.267.

    Article  CAS  Google Scholar 

  4. Egli D, Yepiskoposyan H, Selvaraj A, Balamurugan K, Rajaram R, Simons A, Multhaup G, Mettler S, Vardanyan A, Georgiev O, Schaffner W: A family knockout of all four Drosophila metallothioneins reveals a central role in copper homeostasis and detoxification. Mol Cell Biol. 2006, 26: 2286-2296. 10.1128/MCB.26.6.2286-2296.2006.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  5. Palmiter RD: The elusive function of metallothioneins. Proc Natl Acad Sci USA. 1998, 95: 8428-8430. 10.1073/pnas.95.15.8428.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  6. Davis SR, Cousins RJ: Metallothionein expression in animals: a physiological perspective on function. J Nutr. 2000, 130: 1085-1088.

    CAS  PubMed  Google Scholar 

  7. Chen P, Muñoz A, Nettesheim D, Shaw CF-III, Petering DH: Stoichiometry and cluster specificity of copper binding to metallothionein: homogeneous metal clusters. Biochem J. 1996, 317: 395-402.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Egli D, Domenech J, Selvaraj A, Balamurugan K, Hua H, Capdevila M, Georgiev O, Schaffner W, Atrian S: The four members of the Drosophila metallothionein family exhibit distinct yet overlapping roles in heavy metal homeostasis and detoxification. Genes to Cells. 2006, 11: 647-658. 10.1111/j.1365-2443.2006.00971.x.

    Article  CAS  PubMed  Google Scholar 

  9. Stillman MJ, Zelazowski AJ: Domain specificity of Cd2+ and Zn2+ binding to rabbit liver metallothionein2. Biochem J. 1989, 262: 181-188.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  10. Li H, Otvos JD: 111Cd NMR studies of the domain specificity of Ag+ and Cu+ binding to metallothionein. Biochemistry. 1996, 35: 13929-13936. 10.1021/bi961401n.

    Article  CAS  PubMed  Google Scholar 

  11. Tio L, Villarreal L, Atrian S, Capdevila M: Functional differentiation in the mammalian metallothionein gene family. Metal binding features of Mouse MT4 and comparison with its paralog MT1. J Biol Chem. 2004, 279: 24403-24413. 10.1074/jbc.M401346200.

    Article  CAS  PubMed  Google Scholar 

  12. Kägi JHR, Kojima Y: Chemistry and biochemistry of metallothionein. Metallothionein II. Edited by: Kägi JHR,Kojima Y. 1987, Basel: Birkhäuser Verlag, 25-61. Exp Suppl 52.

    Chapter  Google Scholar 

  13. Vasák M: Standard isolation procedure for metallothionein. Meth Enzymol. 1991, 205: 41-44.

    Article  PubMed  Google Scholar 

  14. Lerch K, Johnsen GF, Grushoff PS, Sternlieb I: Canine hepatic lysosomal copper protein: identification as metallothionein. Arch Biochem Biophys. 1985, 243: 108-114. 10.1016/0003-9861(85)90778-7.

    Article  CAS  PubMed  Google Scholar 

  15. Overnell J: Occurrence of cadmium in crabs (Cancer pagurus) and the isolation and properties of cadmium metallothionein. Environ Hlth Perspect. 1986, 65: 101-105. 10.2307/3430169.

    CAS  Google Scholar 

  16. Waldron KJ, Robinson NJ: How do bacterial cells ensure that metalloproteins get the correct metal?. Nature Rev Microb. 2010, 6: 25-35.

    Google Scholar 

  17. Ponder W, Lindberg DR: Phylogeny and Evolution of the Mollusca. 2008, CA: University of California Press, 1-488.

    Book  Google Scholar 

  18. Dallinger R, Berger B, Hunziker PE, Birchler N, Hauer C, Kägi JHR: Purification and primary structure of snail metallothionein. Similarity of the N-terminal sequence with histones H4 and H2A. Eur J Biochem. 1993, 216: 739-746. 10.1111/j.1432-1033.1993.tb18193.x.

    Article  CAS  PubMed  Google Scholar 

  19. Berger B, Dallinger R, Gehrig P, Hunziker PE: Primary structures of a copper-binding metallothionein from mantle tissue of the terrestrial gastropod Helix pomatia. L Biochem J. 1997, 328: 219-224.

    Article  CAS  PubMed  Google Scholar 

  20. Dallinger R, Berger B, Hunziker PE, Kägi JHR: Metallothionein in snail Cd and Cu metabolism. Nature. 1997, 388: 237-238. 10.1038/40785.

    Article  CAS  PubMed  Google Scholar 

  21. Hispard F, Schuler D, deVaufleury A, Scheifler R, Badot PM, Dallinger R: Metal distribution and metallothionein induction after cadmium exposure in the terrestrial snail Helix aspersa (gastropoda, pulmonata). Environ Toxicol Chem. 27: 1533-1542. 10.1897/07-232.1. 208

  22. Dallinger R, Wang Y, Berger B, Mackay EA, Kägi JHR: Spectroscopic characterization of metallothionein from the terrestrial snail, Helix pomatia. Eur J Biochem. 2001, 268: 4126-4133. 10.1046/j.1432-1327.2001.02318.x.

    Article  CAS  PubMed  Google Scholar 

  23. Gehrig PM, You C, Dallinger R, Gruber C, Brouwer M, Kägi JHR, Hunziker PE: Electrospray ionization mass spectrometry of Zn, Cd and Cu metallothioneins: evidence for metal-binding cooperativity. Prot Sci. 2000, 9: 395-402. 10.1110/ps.9.2.395.

    Article  CAS  Google Scholar 

  24. Chabicovsky M, Niederstaetter H, Thaler R, Hödl E, Parson W, Rossmanith W, Dallinger R: Localization and quantification of Cd- and Cu-specific metallothionein isoform mRNA in cells and organs of the terrestrial gastropod Helix pomatia. Toxicol Appl Pharmacol. 2003, 190: 25-36. 10.1016/S0041-008X(03)00148-0.

    Article  CAS  PubMed  Google Scholar 

  25. Chabicovsky M, Klepal W, Dallinger R: Mechanisms of cadmium toxicity in terrestrial pulmonates: programmed cell death and metallothionein overload. Environ Toxicol Chem. 2004, 23: 648-655. 10.1897/02-617.

    Article  CAS  PubMed  Google Scholar 

  26. Manzl C, Krumschnabel G, Schwarzbaum PJ, Dallinger R: Acute toxicity of cadmium and copper in hepatopancreas cells from the Roman snail (Helix pomatia). Comp Biochem Physiol C. 2004, 138: 45-52. 10.1016/j.cbpb.2004.02.019.

    Article  Google Scholar 

  27. Dallinger R, Chabicovsky M, Hödl E, Prem C, Hunziker P, Manzl C: Copper in Helix pomatia (Gastropoda) is regulated by one single cell type: differently responsive metal pools in rhogocytes. Am J Physiol. 189: R1185-R1195. 205

  28. Vasák M, Kägi JHR, Hill HAO: Zinc(II), cadmium(II) and mercury(II) thiolate transitions in metallothionein. Biochemistry. 1981, 20: 2852-2856.

    Article  PubMed  Google Scholar 

  29. Willner H, Vasák M, Kägi JHR: Cadmium-thiolate clusters in metallothionein: spectrometric and spectropolarimetric features. Biochemistry. 1987, 26: 6287-6292. 10.1021/bi00393a049.

    Article  CAS  PubMed  Google Scholar 

  30. Willner H, Bernhard WR, Kägi JHR: Optical properties of metallothioneins. Origin of the optical properties of mammalian metallothioneins. Metallothioneins, Synthesis, Structure and Properties of Metallothioneins, Phytochelatins and Metal-Thiolate Complexes. Edited by: Stillman MJ, Shaw CF, Suzuki KT. 1992, New York: VCH Publishers, 128-143.

    Google Scholar 

  31. Presta A, Rae Green A, Zelazowski A, Stillman MJ: Copper binding to rabbit liver metallothionein: formation of a continuum of Copper(I)-thiolate stoichiometric species. Eur J Biochem. 1885, 227: 226-240. 10.1111/j.1432-1033.1995.tb20380.x.

    Article  Google Scholar 

  32. Pagani A, Villarreal L, Capdevila M, Atrian S: The Saccharomyces cerevisiae Crs5 metallothionein metal-binding abilities and its role in the response to zinc overload. Mol Microbiol. 2007, 63: 256-269. 10.1111/j.1365-2958.2006.05510.x.

    Article  CAS  PubMed  Google Scholar 

  33. Capdevila M, Domènech J, Pagani A, Tío L, Villarreal L, Atrian S: Zn- and Cd-metallothionein recombinant species from the most diverse phyla may contain sulfide (S2-) ligands. Angew Chem Int Ed Engl. 2005, 44: 4618-4622. 10.1002/anie.200501183.

    Article  CAS  PubMed  Google Scholar 

  34. Bofill R, Capdevila M, Atrian S: Independent metal-binding features of recombinant metallothioneins convergently draw a step gradation between Zn- and Cu-thioneins. Metallomics. 2009, 1: 229-234. 10.1039/b904953c.

    Article  CAS  PubMed  Google Scholar 

  35. Orihuela R, Monteiro F, Pagani A, Capdevila M, Atrian S: Evidence of native metal-(S2-)-MT complexes confirmed by the analysis of Cup1 divalent metal ion binding properties. Chem A Eur J. 16: 12363-12372. 10.1002/chem.201001125. 210

  36. Egg M, Höckner M, Brandstätter A, Schuler D, Dallinger R: Structural and bioinformatic analysis of the Roman snail Cd-Metallothionein gene uncovers molecular adaptation towards plasticity in coping with multifarious environmental stress. Mol Ecol. 2009, 18: 2426-2443. 10.1111/j.1365-294X.2009.04191.x.

    Article  CAS  PubMed  Google Scholar 

  37. Nielson KB, Atkin CL, Winge DR: Distinc Metal-binding Configurations in Metallothioneins. J Biol Chem. 1985, 260: 5342-5350.

    CAS  PubMed  Google Scholar 

  38. Stillman MJ, Rae Green A, Gui Z, Fowle D, Presta PA: Circular dichroism, emission, and exafs studies of Ag(I), Cd(II), Cu(I), and Hg(II) binding to metallothioneins and modelling the metal binding site. Metallothionein. Edited by: Klaassen C. 1999, Basel: Birkhäuser Verlag, 23-35. Vol IV

    Chapter  Google Scholar 

  39. Otvos JD, Engeseth HR, Nettesheim DG, Hilt CR: Interprotein metal exchange reactions of metallothionein. Experientia Suppl. 1987, 52: 171-178.

    Article  CAS  PubMed  Google Scholar 

  40. Maroni G, Wise J, Young JE, Otto E: Metallothionein gene duplications and metal tolerance in natural populations of Drosophila melanogaster. Genetics. 1987, 117: 739-744.

    PubMed Central  CAS  PubMed  Google Scholar 

  41. Tanguy A, Moraga D: Cloning and characterization of a gene coding for a novel metallothionein in the Pacific Oyster Crassostrea gigas (CgMT2): a case of adaptive response to metal-induced stress?. Gene. 2001, 273: 123-130. 10.1016/S0378-1119(01)00577-7.

    Article  CAS  PubMed  Google Scholar 

  42. Kondrashov FA, Kondrashov AS: Role of selection in fixation of gene duplications. J Theoret Biol. 2006, 239: 141-151. 10.1016/j.jtbi.2005.08.033.

    Article  CAS  Google Scholar 

  43. Haszprunar G: The molluscan rhogocyte (pore-cell, Blasenzelle, cellule nucale), and its significance for ideas on nephridial evolution. J Mollus Stud. 1996, 62: 185-211. 10.1093/mollus/62.2.185.

    Article  Google Scholar 

  44. Albrecht U, Keller H, Gebauer W, Markl J: Rhogocytes (pore cells) as the site of hemocyanin biosynthesis in the marine gastropod Haliotis tuberculata. Cell Tissue Res. 2001, 304: 455-462. 10.1007/s004410100368.

    Article  CAS  PubMed  Google Scholar 

  45. Brouwer M, Whaling P, Engel DW: Copper-Metallothioneins in the American Lobster, Homarus americanus: Potential role as Cu(I) donors to apohemocyanin. Environ Hlth Persp. 1986, 65: 93-100. 10.2307/3430168.

    CAS  Google Scholar 

  46. Narula SS, Brouwer M, Hua Y, Armitage IM: Three-dimensional solution structure of Callinectes sapidus metallothionein-1 determined by homonuclear and heteronuclear magnetic resonance spectroscopy. Biochemistry. 1995, 34: 620-631. 10.1021/bi00002a029.

    Article  CAS  PubMed  Google Scholar 

  47. Engel DW, Brouwer M: Metal regulation and molting in the blue crab, Callinectes sapidus: Metallothionein function in metal metabolism. Biol Bull. 1987, 173: 239-251. 10.2307/1541876.

    Article  Google Scholar 

  48. Engel DW, Brouwer M: Short-term metallothionein and copper changes in blue crabs at ecdysis. Biol Bull. 1991, 180: 447-452. 10.2307/1542345.

    Article  CAS  Google Scholar 

  49. Butt TR, Sternberg EJ, Gorman JA, Clark P, Hamer D, Rosenberg M, Crooke ST: Copper Metallothionein of yeast, structure of the gene, and regulation of expression. Proc Natl Acad Sci USA. 1984, 81: 3332-3336. 10.1073/pnas.81.11.3332.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  50. Münger K, Lerch K: Copper Metallothionein from the fungus Agaricus bisporus: chemical and spectroscopic properties. Biochemistry. 1985, 24: 6751-6756.

    Article  Google Scholar 

  51. Beltramini M, Lerch K: Primary structure and spectroscopic studies of Neurospora copper Metallothionein. Environ Hlth Persp. 1986, 65: 21-27. 10.2307/3430157.

    CAS  Google Scholar 

  52. Lerch K: Copper metallothionein, a copper-binding protein from Neurospora crassa. Nature. 1980, 284: 368-370. 10.1038/284368a0.

    Article  CAS  PubMed  Google Scholar 

  53. Dallinger R, Berger B, Bauer-Hilty A: Purification of cadmium-binding proteins from related species of terrestrial helicidae (Gastropoda, Mollusca): a comparative study. Mol Cell Biochem. 1989, 85: 135-145. 10.1007/BF00577109.

    Article  CAS  PubMed  Google Scholar 

  54. Münger K, Lerch K, Tschierpe HJ: Metal accumulation in Agaricus bisporus: Influence of Cd and Cu on growth and tyrosinase activity. Cell Mol Life Sci. 1982, 38: 1039-1041.

    Article  Google Scholar 

  55. Brouwer M, Bonaventura C, Bonaventura J: Metal ion interactions with Limulus polyphemus and Callinectes sapidus hemocyanins: stoichiometry and structural and functional consequences of calcium(II), cadmium(II), zinc(II), and mercury(II) binding. Biochemistry. 1983, 22: 4713-4723. 10.1021/bi00289a016.

    Article  CAS  PubMed  Google Scholar 

  56. Bubacco L, Rocco GP, Salvato B, Beltramini M: The binding of Cd(II) to the hemocyanin of the Mediterranean crab Carcinus maenas. Arch Biochem Biphys. 1993, 302: 78-84. 10.1006/abbi.1993.1183.

    Article  CAS  Google Scholar 

  57. Ngo HTT, Gerstmann S, Frank H: Subchronic effects of environment-like cadmium levels on the bivalve Anodonta anatina (Linnaeus 1758): III. Effects on carbonic anhydrase activity in relation to calcium metabolism. Toxicol Environ Chem. 2010, 92 (10): 1029-0486.

    Google Scholar 

  58. Vasak M, Kägi JHR: Spectroscopic properties of metallothionein. Metal Ions in Biological Systems. Edited by: Sigel H. 1983, New York: Marcel Dekker, 213-273. Vol 15

    Google Scholar 

  59. Haq F, Mahoney M, Koropatnick J: Signalling events for metallothionein induction. Mutat Res. 2003, 533: 211-226.

    Article  CAS  PubMed  Google Scholar 

  60. Dallinger R, Chabicovsky M, Berger B: Isoform-specific quantification of metallothionein in the terrestrial gastropod Helix pomatia I. Molecular, biochemical, and methodical background. Environ Toxicol Chem. 2004, 23: 890-901. 10.1897/03-100.

    Article  CAS  PubMed  Google Scholar 

  61. Cols N, Romero-Isart N, Capdevila M, Oliva B, Gonzàlez-Duarte P, Gonzàlez-Duarte R, Atrian S: Binding of excess cadmium(II) to Cd7-metallothionein from recombinant mouse Zn7-metallothionein 1. UV-VIS absorption and circular dichroism studies and theoretical location approach by surface accessibility analysis. J Inorg Biochem. 1997, 68: 157-166. 10.1016/S0162-0134(97)00085-8.

    Article  CAS  PubMed  Google Scholar 

  62. Bongers J, Walton CD, Richardson DE, Bell JU: Micromolar protein concentrations and metalloprotein stoichiometries obtained by inductively coupled plasma. Atomic emission spectrometric determination of sulphur. Anal Chem. 1988, 60: 2683-2686. 10.1021/ac00175a008.

    Article  CAS  PubMed  Google Scholar 

  63. Fabris D, Zaia J, Hathout Y, Fenselau C: Retention of thiol protons in two classes of protein zinc ion coordination centers. J Am Chem Soc. 1996, 118: 12242-12243. 10.1021/ja962863r.

    Article  CAS  Google Scholar 

  64. Capdevila M, Cols N, Romero-Isart N, González-Duarte R, Atrian S, González-Duarte P: Recombinant synthesis of mouse Zn3-α and Zn4-α metallothionein 1 domains and characterization of their cadmium(II) binding capacity. Cell Mol Life Sci. 1997, 53: 681-688. 10.1007/s000180050088.

    Article  CAS  PubMed  Google Scholar 

  65. Bofill R, Palacios O, Capdevila M, Cols N, González-Duarte R, Atrian S, González-Duarte P: A new insight into the Ag+ and Cu+ binding sites in the metallothionein β domain. J Inorg Biochem. 1999, 73: 57-64. 10.1016/S0162-0134(98)10091-0.

    Article  CAS  PubMed  Google Scholar 

  66. Longo VD, Gralla EB, Valentine JS: Superoxide dismutase activity is essential for stationary phase survival in Saccharomyces cerevisiae. J Biol Chem. 1996, 271: 12275-12280. 10.1074/jbc.271.21.12275.

    Article  CAS  PubMed  Google Scholar 

  67. Mumberg D, Müller R, Funk M: Yeast vectors for the controlled expression of heterologous proteins in different genetic backgrounds. Gene. 1995, 156: 119-122. 10.1016/0378-1119(95)00037-7.

    Article  CAS  PubMed  Google Scholar 

  68. Stearns T, Ma H, Botstein D: Manipulating yeast genome using plasmid vectors. Methods Enzymol. 1990, 185: 280-297. full_text.

    Article  CAS  PubMed  Google Scholar 

  69. Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, et al: Clustal W and Clustal X version 2.0. Bioinformatics. 2007, 23: 2947-2948. 10.1093/bioinformatics/btm404.

    Article  CAS  PubMed  Google Scholar 

  70. Saitou N, Nei M: The neighbour-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4: 406-425.

    CAS  PubMed  Google Scholar 

  71. Felsenstein J: Confidence limits on phylogenies: an approach using the bootstrap. Evolution. 1985, 39: 783-791. 10.2307/2408678.

    Article  Google Scholar 

Download references


This work was supported by the Spanish Ministerio de Ciencia y Tecnología grants BIO2009-12513-C02-01 to SA, BIO2009-12513-C02-02 to MC and project No. P19782-B02 from the Austrian Science Foundation to RD. Collaboration between the Spanish and Austrian research groups was financed by the 'Acciones Integradas' grants HU2006-0027 (Spain) and ES 02/2007 (Austria). SPR received a pre-doctoral fellowship from the Department de Química, Universitat Autònoma de Barcelona. Authors from the University of Innsbruck are members of the 'Centre of Molecular Biosciences Innsbruck'. The authors from the Universitat Autònoma de Barcelona (UAB) and Universitat de Barcelona (UB) are members of a 'Grups de Recerca de la Generalitat de Catalunya' (2009SGR-1457). We thank the Serveis Cientificotècnics of UB (GC-FPD, ICP-AES, DNA sequencing) and the Servei d'Anàlisi Química of UAB (CD, UV-vis, ESI-MS) for allocating instrument time to us.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Reinhard Dallinger.

Additional information

Authors' contributions

ÒP, AP, and SPR carried out the recombinant expression studies, participating in the construction of the expression vectors, the synthesis and characterization of the corresponding metal complexes, and also performed the yeast complementation experiments. ME and MH performed quantitative Real-Time PCR experiments and in situ hybridization, and carried out most of the sequencing work of MT cDNAs from different pulmonate snail species. AB calculated and established the phylogenetic trees. MC and SA designed the study together with RD, supervised recombinant DNA and recombinant protein experiments, yeast complementation studies, as well as analytical, spectrometric and spectrophotometric work, and contributed to the drafting of the manuscript. RD supervised and participated in molecular sequencing and Real-Time PCR as well as in situ hybridization, and performed chromatography and Reversed-Phase HPLC. He was responsible for the protein alignments and assisted in calculation of phylogenetic trees. He designed the study together with MC and SA, and drafted the manuscript. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Palacios, Ò., Pagani, A., Pérez-Rafael, S. et al. Shaping mechanisms of metal specificity in a family of metazoan metallothioneins: evolutionary differentiation of mollusc metallothioneins. BMC Biol 9, 4 (2011).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: