Effect of α-domain substitution on the structure, property and function of human neuronal growth inhibitory factor

Human metallothionein-3 (hMT3), also named human neuronal growth inhibitory factor (hGIF), is attractive due to its distinct neuronal growth inhibitory activity, which is not shown by other human MT isoforms. It has been reported that the neuronal growth inhibitory activity arises from the N-terminal β-domain rather than its C-terminal α-domain. However, previous bioassay results have shown that the single β-domain is less effective at inhibiting the neuron growth than that in intact hMT3 on a molar basis, which suggests that the α-domain is indispensable to the neuronal growth inhibitory activity of hMT3. In order to confirm this assumption, we constructed two domain-hybrid mutants, the β(MT3)–β(MT3) mutant and the β(MT3)–α(MT1) mutant, and investigated their structural and metal binding properties by UV-vis spectroscopy, CD spectroscopy, pH titration, DTNB reaction, EDTA reaction, etc. The results showed that stability of the Cd₃S₉ cluster of the β(MT3)–β(MT3) mutant decreased significantly while the Cd₃S₉ cluster of the β(MT3)–α(MT1) mutant had a similar stability and solvent accessibility to that of hMT3. Interestingly, the bioassay results showed that the neuronal growth inhibitory activity of the β(MT3)–β(MT3) mutant decreased significantly, while the β(MT3)–α(MT1) mutant showed similar inhibitory activity to hMT3. Based on these results, we conclude that the α-domain is indispensable and plays an important role in modulating the stability of the metal cluster in the β-domain by domain–domain interactions, thus influencing the bioactivity of hMT3. Z.-C. Ding and Q. Zheng contributed equally to this work.     

Incorporation of a glycine within the conserved TCPCP motif of human neuronal growth inhibitory factor significantly reduces its bioactivity

It has been reported that the ⁶CPCP⁹ motif near the N-terminus is pivotal to the inhibitory activity of human neuronal growth inhibitory factor (hGIF). In order to better understand the biological significance of this region on the structure, property and function of hGIF, we introduced a highly flexible residue, Gly, either in front of the ⁶CPCP⁹ motif (the IG6 mutant, TGCPCP) or in the middle of it (the IG8 mutant, TCPGCP) and investigated their structural and metal binding properties in detail. The results showed that the overall structure and the stability of the metal-thiolate clusters of the two mutants were comparable to that of hGIF. However, the bioassay results showed that the bioactivity of the IG6 mutant decreased significantly, while the bioactivity of the IG8 mutant was almost abolished. Molecular dynamics simulation results showed that the backbone of the IG6 mutant exhibited high similarity to that of hGIF, and the two prolines could still induce structural constraints on the ⁶CPCP⁹ tetrapeptide and form a similar conformation with that of hGIF, however, the conformation of the first five amino acid residues in the N-terminus was quite different. In hGIF, the five residues are twisted and form a restricted conformation, while in the IG6 mutant this peptide extends more naturally and smoothly, which is similar to that of MT2. As to the IG8 mutant, the Gly insertion broke the ⁶CPCP⁹ motif, thus probably abolishing the interactions with other molecules and eliminating its inhibitory activity. Based on these results, we suggested that although the structure adopted by the ⁶CPCP⁹ motif is the determinant factor of the inhibitory bioactivity of hGIF, other residues within the N-terminal fragment (residue 1-13) may also influence the peptide conformation and contribute to the protein’s bioactivity.     

The role of Thr5 in human neuron growth inhibitory factor

GIF, a member of the metallothionein (MT) family (assigned as MT3), is a neuron growth inhibitory factor that inhibits neuron outgrowth in Alzheimer’s disease. The conserved Thr5 is one of the main differences between GIF and other members in the MT family. However, natural sheep GIF has an unusual Ala5, casting doubt on the role of common Thr5. We constructed a series of human GIF mutants at site 5, and characterized their biochemical properties by UV spectroscopy, circular dichroism spectroscopy, EDTA reaction, 5,5′-dithiobis(2-nitrobenzoic acid) (DTNB) reaction, and pH titration. Their inhibitory activity toward neuron survival and neurite extension was also examined. Interestingly, the T5A mutant exhibited distinct metal thiolate activity in the EDTA and DTNB reactions, and also lost its bioactivity. Meanwhile, the T5S mutant had similar biochemical properties and biological activity as wild-type human GIF, indicating the hydroxyl group on the Thr5 was critical to the bioactivity of human GIF. We suggest the hydroxyl group in human GIF may help stabilize the biologically active conformation. On the other hand, lack of the hydroxyl group in sheep GIF may be partially compensated by its abnormal structure.Bin Cai and Qi Zheng contributed equally to this work.     

The cadmium binding domains in the metallothionein isoform Cd7-MT10 from Mytilus galloprovincialis revealed by NMR spectroscopy

The metal–thiolate connectivity of recombinant Cd₇-MT10 metallothionein from the sea mussel Mytilus galloprovincialis has been investigated for the first time by means of multinuclear, multidimensional NMR spectroscopy. The internal backbone dynamics of the protein have been assessed by the analysis of ¹⁵N T ₁ and T ₂ relaxation times and steady state {¹H}–¹⁵N heteronuclear NOEs. The ¹¹³Cd NMR spectrum of mussel MT10 shows unique features, with a remarkably wide dispersion (210 ppm) of ¹¹³Cd NMR signals. The complete assignment of cysteine Hα and Hβ proton resonances and the analysis of 2D ¹¹³Cd–¹¹³Cd COSY and ¹H–¹¹³Cd HMQC type spectra allowed us to identify a four metal–thiolate cluster (α-domain) and a three metal–thiolate cluster (β-domain), located at the N-terminal and the C-terminal, respectively. With respect to vertebrate MTs, the mussel MT10 displays an inversion of the α and β domains inside the chain, similar to what observed in the echinoderm MT-A. Moreover, unlike the MTs characterized so far, the α-domain of mussel Cd₇-MT10 is of the form M₄S₁₂ instead of M₄S₁₁, and has a novel topology. The β-domain has a metal–thiolate binding pattern similar to other vertebrate MTs, but it is conformationally more rigid. This feature is quite unusual for MTs, in which the β-domain displays a more disordered conformation than the α-domain. It is concluded that in mussel Cd₇-MT10, the spacing of cysteine residues and the plasticity of the protein backbone (due to the high number of glycine residues) increase the adaptability of the protein backbone towards enfolding around the metal–thiolate clusters, resulting in minimal alterations of the ideal tetrahedral geometry around the metal centres.     

Study on structure-property-reactivity-function relationship of human neuronal growth inhibitory factor (hGIF)

Human metallothionein-3 (hMT3), also named as human neuronal growth inhibitory factor (hGIF), can inhibit the outgrowth of embryonic cortical neurons in the presence of brain extracts. In order to systematically study the structure-property-reactivity-function relationship of hGIF, our laboratory designed a series of mutants and studied their structure, property, reactivity and functions by a series of chemical and biological tools including UV spectroscopy, CD spectroscopy, NMR, chemical reaction and primary neuronal culture assays. In summary, we concluded that the bioactivity of hGIF was regulated by multiple factors, including the ⁶CPCP⁹ motif, an additional threonine insert at sequence position 5, domain-domain interactions, the structure and stability of the metal-thiolate cluster and the linker. Our studies provide more and more evidences which revealed that the bioactivity of hGIF is mainly related to the essential metal release and its characteristic conformation.     

Studies on amino acid replacement and inhibitory activity of a β-lactamase inhibitory peptide

An SHV β-lactamase gene was amplified from a β-lactam resistant Klebsiella pneumoniae K-71 genomic DNA. After expression and purification, we demonstrated that peptide P1 could inhibit the hydrolysis activity of both TEM-1 and SHV β-lactamase in vitro. Three mutations were introduced into P1 in which the first residue S was replaced by F, the 18th residue V was mutated to Y, and the 15th residue Y was substituted with A, C, G, and R to obtain the mutants of P1-A, P1- C, P1-G, and P1-R, respectively. The mutant peptides were purified and their inhibitory constants against TEM-1 and SHV β-lactamase were determined. All these β-lactamase inhibitory peptides could inhibit the activity of both β-lactamases, while the mutant peptides showed stronger inhibitory activities against TEM-1 β-lactamase than against SHV β-lactamase. Inhibition data suggested that P1-A improved the β-lactamase inhibitory activity by over 3-fold compare to P1. When P1-A was incubated with K. pneumoniae K-71 in Luria-Bertani medium containing ampicillin, it showed a much stronger growth of inhibition ratio over P1. This study gives us a good candidate for development of novel β-lactamase inhibitors.     

Isolation of GIF from porcine brain and studies of its zinc transfer kinetics with apo-carbonic anhydrase

Neuronal growth inhibitory factor (GIF) of porcine brain, was isolated and purified by a similar procedure which was used on the isolation of human and bovine GIF. The native porcine protein with stoichiometry of 4Cu⁺, 3Zn²⁺ was obtained for the first time. The kinetics of zinc transfer from Cu₄Zn₃MT-3 to apo-carbonic anhydrase were studied, and zinc transfer rate constants and thermodynamic parameters were obtained. It is found that like other MTs, porcine Cu₄Zn₃MT-3 can also transfer its zinc atom to apoCA, even much easier than other MTs. A possible association mechanism has been proposed, the formation of Cu₄Zn₃MT3-apoCA complex may be the rate-determining step. The obtained data indicate besides its neuronal growth inhibitory function, GIF might play a role in cellular Zn homeostasis in brain.     

The effect of nitric oxide on metal release from metallothionein-3: gradual unfolding of the protein

Metallothionein-3 (MT-3), or neuronal growth inhibitory factor, which exhibits growth inhibitory activity, is a brain-specific metallothionein. In this study, the effect of nitric oxide (NO) on metal release (using Cd²⁺ as a probe) from MT-3 was examined by ¹¹³Cd and 2D [¹H–¹⁵N] heteronuclear single-quantum coherence NMR spectroscopy. The exposure of human MT-3 to NO leads to a nonselective release of the three metals from the β-domain. In contrast to metallothionein-1 and metallothionein-2, two of the bound metals in the α-domain were also partially released, with the domain structure remaining almost unchanged. Further addition of NO resulted in the complete release of metals and concomitant unfolding of the protein. The preference of release of the two metals in the α-domain was attributed to the presence of two slightly different coordination environments for the four cadmium/zinc atoms.     

Different protective roles in vitro of α- and β-domains of growth inhibitory factor (GIF) on neuron injuries caused by oxygen free radicals

It was well known that β-amyloid (Aβ) and tau protein play an important role in pathological procedure of Alzheimer’s disease (AD), a senile dementia. The growth inhibitory factor (GIF, also named metallothionein-3, MT-3) had been demonstrated to inhibit the outgrowth of cortex neurons in the medium with extract of the AD patient brain. In our experiments, it was found that the neurons of cortex and the PC12 (pheochromocytoma) cells could be protected from the cytotoxicity of β-amyloid 25–35 in presence of GIF and its domains. Additionally, GIF can scavenge the hydroxyl radical efficiently in CytC–VitC radical producing system and its α-domain shown more effective potentials than its β-domain. The electron paramagnetic resonance spectra also show that the α-domain has more potential ability for eliminating reactive oxygen free radicals than its β-domain. The results suggest that GIF could act as an efficient scavenger against free radicals in vitro and the α-domain in GIF molecule shows more potential in protecting against reactive oxygen species injury than the β-domain.     

The molecular mechanism for human metallothionein-3 to protect against the neuronal cytotoxicity of Aβ1–42 with Cu ions

Aggregation and cytotoxicity of Aβ with redox-active metals in neuronal cells have been implicated in the progression of Alzheimer disease. Human metallothionein (MT) 3 is highly expressed in the normal human brain and is downregulated in Alzheimer disease. Zn₇MT3 can protect against the neuronal toxicity of Aβ by preventing copper-mediated Aβ aggregation, abolishing the production of reactive oxygen species (ROS) and the related cellular toxicity. In this study, we intended to decipher the roles of single-domain proteins (α/β) and the α–β domain–domain interaction of Zn₇MT3 to determine the molecular mechanism for protection against the neuronal cytotoxicity of Aβ_1–42 with copper ions. With this in mind, the α and β single-domain proteins, heterozygous β(MT3)–α(MT1), and a linker-truncated mutant ?31–34 were prepared and characterized. In the presence/absence of various Zn₇MT3 proteins, the Aβ_1–42–Cu²⁺-mediated aggregation, the production of ROS, and the cellular toxicity were investigated by transmission electron microscopy, ROS assay by means of a fluorescent probe, and SH-SY5Y cell viability, respectively. The β domain cannot abolish Aβ_1–42–Cu²⁺-induced aggregation, and neither the β domain nor the α domain can quench the production of ROS because of the redox cycling of Aβ–Cu²⁺. Similarly to wild-type Zn₇MT3, the heterozygous β(MT3)–α(MT1) possesses the characteristic of alleviating Aβ_1–42 aggregation and oxidative stress to neuronal cells. Therefore, the two domains through the linker Lys-Lys-Ser form a cooperative unit, and each of them is indispensable in conducting its bioactivity. The α domain plays an important role in modulating the stability of the metal–thiolate cluster, and the α–β domain–domain interaction through the linker is critical for its protective role in the brain.     

Characterization of glycerol dehydratase expressed by fusing its α- and β-subunits

The gdh and gdhr genes, encoding B₁₂-dependent glycerol dehydratase (GDH) and glycerol dehydratase reactivase (GDHR), respectively, in Klebsiella pneumoniae, were cloned and expressed in E. coli. Part of the β-subunit was lost during GDH purification when co-expressing α, β and γ subunit. This was overcome by fusing the β-subunit to α- or γ-subunit with/without the insertion of a linker peptide between the fusion moieties. The kinetic properties of the fusion enzymes were characterized and compared with wild type enzyme. The results demonstrated that the fusion protein GDHALB/C, constructed by linking the N-terminal of β-subunit to the C-terminal of α subunit through a (Gly₄Ser)₄ linker peptide, had the greatest catalytic activity. Similar to the wild-type enzyme, GDHALB/C underwent mechanism-based inactivation by glycerol during catalysis and could be reactivated by GDHR. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Kinetic and structural studies on roles of the serine ligand and a strictly conserved tyrosine residue in nitrile hydratase

Nitrile hydratases (NHase), which catalyze the hydration of nitriles to amides, have an unusual Fe³⁺ or Co³⁺ center with two modified Cys ligands: cysteine sulfininate (Cys-SO₂ ⁻) and either cysteine sulfenic acid or cysteine sulfenate [Cys-SO(H)]. Two catalytic mechanisms have been proposed. One is that the sulfenyl oxygen activates a water molecule, enabling nucleophilic attack on the nitrile carbon. The other is that the Ser ligand ionizes the strictly conserved Tyr, activating a water molecule. Here, we characterized mutants of Fe-type NHase from Rhodococcus erythropolis N771, replacing the Ser and Tyr residues, αS113A and βY72F. The αS113A mutation partially affected catalytic activity and did not change the pH profiles of the kinetic parameters. UV–vis absorption spectra indicated that the electronic state of the Fe center was altered by the αS113A mutation, but the changes could be prevented by a competitive inhibitor, n-butyric acid. The overall structure of the αS113A mutant was similar to that of the wild type, but significant changes were observed around the catalytic cavity. Like the UV–vis spectra, the changes were compensated by the substrate or product. The Ser ligand is important for the structure around the catalytic cavity, but is not essential for catalysis. The βY72F mutant exhibited no activity. The structure of the βY72F mutant was highly conserved but was found to be the inactivated state, with αCys114-SO(H) oxidized to Cys-SO₂ ⁻, suggesting that βTyr72 affected the electronic state of the Fe center. The catalytic mechanism is discussed on the basis of the results obtained.     

Molecular Evolution and Characterization of Fungal Indoleamine 2,3-Dioxygenases

Indoleamine 2,3-dioxygenase (IDO) and tryptophan 2,3-dioxygenase (TDO) are tryptophan-degrading enzymes. Mammalian IDO expression is induced by cytokines and has antimicrobial and immunomodulatory effects. A major role of mammalian TDO is to supply nicotinamide adenine dinucleotide (NAD⁺). In fungi, the IDO homologue is thought to be expressed constitutively and supply NAD⁺, as TDO is absent from their genomes. Here, we reveal the distribution of IDO genes among fungal species and characterize their enzymatic activity. The yeast, Saccharomyces cerevisiae has only one IDO gene, whereas the koji-mold, Aspergillus oryzae has two genes, IDOα and IDOβ. The A. oryzae IDOα showed more similar enzymatic properties to those of S. cerevisiae IDO than IDOβ, suggesting that the A. oryzae IDOα is a functional homologue of the S. cerevisiae IDO. From the IDOβ gene, two isoforms, IDOβ and IDOβ⁺ could be generated by alternative splicing. The latter contained a 17 amino acids insertion which were encoded by the first intron of IDOβ gene. In comparison to IDOβ⁺, bacterially expressed IDOβ showed much lower K _m value and more than five-times faster V _max value, resulting in 85 times higher catalytic efficiency; i.e., the removal of the domain encoded by the first intron from IDOβ⁺ increases its enzymatic activity drastically. This might be a unique regulation mechanism of the l-Trp metabolism in the A. oryzae. The levo-1-methyl tryptophan (l-1MT) is a good inhibitor of both IDO1 and IDO2. However, the activity of fungal IDOs tested was not inhibited at all by l-1MT.     

Unusual hydrophobic linker region of β-glucosidase (BGLII) from Thermoascus aurantiacus is required for hyper-activation by organic solvents

A gene encoding a putative β-glucosidase was isolated from Thermoascus aurantiacus IFO9748 and designated as bgl2. The recombinant enzyme showed β-glucosidase activity when p-nitrophenyl-β-glucose (pNP-Glc) was used as substrate. We also found that the enzyme activity was increased in the presence of organic solvents. An addition of 20 % (v/v) 1-octanol resulted in 54-fold higher activity of pNP-Glc hydrolysis, and transglycosylation activity was also found to be activated. The results of tryptophan fluorescence spectral analysis revealed the changes in the tertiary structure of the enzyme in the presence of 1-hexanol that may cause increased enzyme activity. BGLII has a distinctive hydrophobic linker region between N- and C-terminal domains. A chimeric enzyme in which the linker region was substituted by the corresponding region of another β-glucosidase failed to be activated by organic solvents, suggesting that the hydrophobic linker region may act as a molecular switch in BGLII.     

Transferase and hydrolytic activities of the laminarinase from rhodothermus marinus and its M133A, M133C, and M133W mutants

Comparative studies of the transglycosylation and hydrolytic activities have been performed on the Rhodothermus marinus β-1,3-glucanase (laminarinase) and its M133A, M133C, and M133W mutants. The M133C mutant demonstrated near 20% greater rate of transglycosylation activity in comparison with the M133A and M133W mutants that was measured by NMR quantitation of nascent β(1-4) and β(1-6) linkages. To obtain kinetic probes for the wild-type enzyme and Met-133 mutants, p-nitrophenyl β-laminarin oligosaccharides of degree of polymerisation 2–8 were synthesized enzymatically. Catalytic efficiency values, k _cat/K _m, of the laminarinase catalysed hydrolysis of these oligosaccharides suggested possibility of four negative and at least three positive binding subsites in the active site. Comparison of action patterns of the wild-type and M133C mutant in the hydrolysis of the p-nitrophenyl-β-D-oligosac- charides indicated that the increased transglycosylation activity of the M133C mutant did not result from altered subsite affinities. The stereospecificity of the transglycosylation reaction also was unchanged in all mutants; the major transglycosylation products in hydrolysis of p-nitrophenyl laminaribioside were β-glucopyranosyl-β-1,3-D-glucopy- ranosyl-β-1,3-D-glucopyranose and β-glucopyranosyl-β-1, 3-D-glucopyranosyl-β-1,3-D-glucpyranosyl-β-1,3-D- glucopyranoxside. In a memoriam of Dr. Kirill N. Neustroev. All we, his friends and colleagues, mourn for his sudden death. He was a bright and talented scientist, brilliant manager and good friend.     

Cadmium-induced conformational changes in type 2 metallothionein of medicinal plant Coptis japonica: insights from molecular dynamics studies of apo,partially and fully metalated forms

Plants play an important role in the removal of excess heavy metals from soil and water. Medicinal plants can also have non-traditional use in phytoremediation technologies. Among the heavy metals, Cadmium (Cd) is the most abundant and readily taken up by the crop plants. Plant metallothioneins (MTs) are small proteins having cysteine-rich residues and appear to play key roles in metal homoeostasis. Plant metallothionein 2 (MT 2) from Coptis japonica (Gold-thread; CjMT 2) is a typical member of this subfamily and features two cysteine-rich regions containing eight and six cysteine residues, respectively, separated by 42 amino acids long linker region. In-silico analysis of MT 2 protein sequences of C. japonica was performed. In this study, ab initio methods were utilised for the prediction of three-dimensional structure of CjMT 2. After structure validation, heavy metal-binding sites were predicted for the selected modelled structures of CjMT 2. To obtain Cd_i-CjMT 2 (i = 1–7), metalated complex individual docking experiments were performed. The stability of the metalated docked structures was assessed by molecular dynamics (MD) simulation studies. Our study showed that CjMT 2 binds up to 4 Cd²⁺ ions in two distinct domains: a N-terminal β-domain that binds to 2 Cd²⁺ ions and a C-terminal α-domain that binds with 2 Cd²⁺ ions. Our analysis revealed that Cys residues of alpha and beta domain and some residues of spacer region of CjMT 2 protein might be important for the cadmium interaction. MD simulation studies provided insight into metal-induced conformational changes and mechanism of metalation of CjMT 2, an intrinsically disordered protein. This study provides useful insights into mechanism of cadmium-type 2 metallothionein interaction.     

Autocrine secretion of TGF-β1 and TGF-β2 by pre-adipocytes and adipocytes: A potent negative regulator of adipocyte differentiation and proliferation of mammary carcinoma cells

Summary We have developed an in vitro system to examine the influence of adipocytes, a major mammary stromal cell type, on the growth of a murine mammary carcinoma, SP1. Previously, we have shown that 3T3-L1 adipocytes release a mitogenic factor, hepatocyte growth factor, which strongly stimulates proliferation of SP1 cells. We now show that 3T3-L1 pre-adipocytes secrete active inhibitory molecules which inhibit DNA synthesis in SP1 cells. In addition, latent inhibitory activity is present in conditioned media (CM) from both pre-adipocytes and adipocytes, and is activated following acid treatment. CM also inhibited DNA synthesis in Mv1Lu wild type epithelial cells, but not DR27 mutant epithelial cells which lack TGF-β type II receptor. Inhibitory activity of CMs was partially abrogated by neutralizing anti-TGF-β1 and anti-TGF-β2 antibodies, and was removed following ultrafiltration through membranes of 10 000 M_r but not 30 000 M_r pore size. These results show that the inhibitory effect on DNA synthesis is mediated by TGF-β1-like and TGF-β2-like molecules. In addition, acid-treated CM as well as purified TGF-β inhibited differentiation of pre-adipocytes. Untreated pre-adipocyte CM, but not mature adipocyte CM, spontaneously inhibited adipocyte differentiation. Together, these findings indicate that pre-adipocytes spontaneously activate their own secreted TGF-β, whereas mature adipocytes do not, and suggest that activation of TGF-β has a potent negative regulatory effect on adipocyte differentiation and tumor growth. Thus, TGF-β may be an important modulator of tumor growth and adipocyte differentiation via both paracrine and autocrine mechanisms. These findings emphasize the importance of adipocyte-tumor interactions in the regulation of tumor microenvironment.     

Computational Analysis of Cysteine Substitutions Modelled on the α- and β-domains of Cd5,Zn2-Metallothionein 2

Several mutant forms of rat liver Cd₅,Zn₂-metallothionein 2 (Cd₅,Zn₂-MT 2) [1] have been computationally modelled and analysed. All terminal cysteines (5, 13, 19, 21, 26, 29, 33, 36, 41, 48, 57 and 59, Figure 1) have been independently substituted by three other co-ordinating amino-acids (aspartate, glutamate and histidine), and the side-chains of the mutated residues have been modelled to co-ordinate the seven metal ions while minimizing the conformational variations with respect to the wild type protein. We have compared the ability of the putative mutant forms to maintain the MT binding properties. Substitution by aspartate residue best preserves the 3D MT structure. In addition, the mutations C5H plus C21H/E/D show neighbouring impairments that prevent their simultaneous substitution. Although replacement of cysteine by aspartate is feasible in all cases, to our knowledge there is no example of aspartate and cysteine residues co-ordinating to the same zinc atom. Accordingly, the use of histidine or glutamate instead of aspartate cannot be ruled out. The mutant forms in the -domain of Cd₅,Zn₂-MT 2 have yielded more neighbouring contacts than those in the -domain, which is corroborated by the accessible surface areas [2] of the sulfur atoms [3] in the native form.Abbreviations MT metallothionein - CD₅,Zn₂-MT Cadmium, Zinc-metallothionein - RMSD Root Mean Square Deviation - PDB Protein Data Bank - FEP Free Energy Perturbation - CnX mutant form of cysteine n (n = residue number) substituted by X (X = H, E or D, with H = histidine, E = glutamate, D = aspartate) - CnX/Y mutant forms CnX and CnY     

Influence of {\text{NH - }}{{\text{S}}^\gamma } bonding interactions on the structure and dynamics of metallothioneins

Mammalian metallothioneins ( \textM₇^\textIIMTs {\text{M}}_7^{\text{IIMTs}} ) show a clustered arrangement of the metal ions and a nonregular protein structure. The solution structures of Cd₃-thiolate cluster containing β-domain of mouse β-MT-1 and rat β-MT-2 show high structural similarities, but widely differing structure dynamics. Molecular dynamics simulations revealed a substantially increased number of \textNH - \textS^g {\text{NH - }}{{\text{S}}^\gamma } hydrogen bonds in β-MT-2, features likely responsible for the increased stability of the Cd₃-thiolate cluster and the enfolding protein domain. Alterations in the \textNH - \textS^g {\text{NH - }}{{\text{S}}^\gamma } hydrogen-bonding network may provide a rationale for the differences in dynamic properties encountered in the β-domains of MT-1, -2, and -3 isoforms, believed to be essential for their different biological function.