# Mechanism of Zn 2+ and Ca 2+ Binding to Human S100A1

Viktoriia E Baksheeva, Andrei Yu Roman, Claude Villard, François Devred, Byrne Deborah, Dahbia Yatoui, Arth O Zalevsky, Alisa A Vologzhanniko, Andrey S Sokolov, Sergei E Permyakov, et al.

► To cite this version:

Viktoriia E Baksheeva, Andrei Yu Roman, Claude Villard, François Devred, Byrne Deborah, et al. Mechanism of Zn 2+ and Ca 2+ Binding to Human S100A1. Biomolecules, 2021, 11 (12), pp.1823. 10.3390/biom11121823. [hal-03464582](https://hal.science/hal-03464582)

HAL Id: hal-03464582

https://hal.science/hal-03464582

Submitted on 3 Dec 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Mechanism of Zn<fnref n="2" />⁺ and Ca<fnref n="2" />⁺ Binding to Human S100A1

Viktoriia E. Baksheeva <fnref n="1" />,<fnref n="2" />, Andrei Yu. Roman <fnref n="1" />,<fnref n="3" />, Claude Villard <fnref n="1" />, François Devred ¹ID, Deborah Byrne <fnref n="4" />,
Dahbia Yatoui <fnref n="1" />, Arthur O. Zalevsky ⁵ID, Alisa A. Vologzhannikova ⁶ID, Andrey S. Sokolov <fnref n="6" />,
Sergei E. Permyakov ⁶ID, Andrey V. Golovin <fnref n="5" />,<fnref n="7" />,⁸ID, Gary S. Shaw ⁹ID, Philipp O. Tsvetkov <fnref n="1" />,*ID and
Evgeni Yu. Zernii <fnref n="2" />,*ID

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Baksheeva, V.E.; Roman, A.Y.; Villard, C.; Devred, F.; Byrne, D.; Yatoui, D.; Zalevsky, A.O.; Vologzhannikova, A.A.; Sokolov, A.S.; Permyakov, S.E.; et al. Mechanism of Zn<fnref n="2" />⁺ and Ca<fnref n="2" />⁺ Binding to Human S100A1. *Biomolecules* **2021**, *11*, 1823. https://doi.org/10.3390/biom11121823

Academic Editor: Alessandro Alaimo

Received: 22 October 2021

Accepted: 30 November 2021

Published: 3 December 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2021 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** S100A1 is a member of the S100 family of small ubiquitous Ca<fnref n="2" />⁺-binding proteins, which participates in the regulation of cell differentiation, motility, and survival. It exists as homo- or heterodimers. S100A1 has also been shown to bind Zn<fnref n="2" />⁺, but the molecular mechanisms of this binding are not yet known. In this work, using ESI-MS and ITC, we demonstrate that S100A1 can coordinate 4 zinc ions per monomer, with two high affinity (K_D~4 and 770 nm) and two low affinity sites. Using competitive binding experiments between Ca<fnref n="2" />⁺ and Zn<fnref n="2" />⁺ and QM/MM molecular modeling we conclude that Zn<fnref n="2" />⁺ high affinity sites are located in the EF-hand motifs of S100A1. In addition, two lower affinity sites can bind Zn<fnref n="2" />⁺ even when the EF-hands are saturated by Ca<fnref n="2" />⁺, resulting in a 2Ca<fnref n="2" />⁺:S100A1:2Zn<fnref n="2" />⁺ conformer. Finally, we show that, in contrast to calcium, an excess of Zn<fnref n="2" />⁺ produces a destabilizing effect on S100A1 structure and leads to its aggregation. We also determined a higher affinity to Ca<fnref n="2" />⁺ (K_D~0.16 and 24 µm) than was previously reported for S100A1, which would allow this protein to function as a Ca<fnref n="2" />⁺/Zn<fnref n="2" />⁺-sensor both inside and outside cells, participating in diverse signaling pathways under normal and pathological conditions.

**Keywords:** S100A1; zinc; calcium; ESI-MS; ITC; nanoDSF

## 1. Introduction

The S100 protein family encompasses small (~10 kDa) ubiquitous two-EF-hand Ca<fnref n="2" />⁺-binding proteins, normally existing as homo- or heterodimers and acting as both secreted and intracellular signaling molecules [1]. In cells, the S100 proteins participate in the regulation of mobility, differentiation and survival through interactions with numerous targets, such as effector enzymes, cytoskeletal proteins, receptors, and transcription factors. When released from cells, they display cytokine-like activity and bind to cell surface receptors (like toll-like receptors and G-protein-coupled receptors) or interact with growth

<footnote n="1">¹ Institut de Neurophysiopathologie, INP, CNRS, Faculté des Sciences Médicales et Paramédicales, Aix-Marseille Université, 13005 Marseille, France; vbaksheeva@belozersky.msu.ru (V.E.B.); andrei.roman@ubc.ca (A.Y.R.); claude.villard@univ-amu.fr (C.V.); francois.devred@univ-amu.fr (F.D.); yatouidahbia19@gmail.com (D.Y.)</footnote>

<footnote n="2">² Belozersky Institute of Physico-Chemical Biology, Lomonosov Moscow State University, 119992 Moscow, Russia</footnote>

<footnote n="3">³ Institute of Physiologically Active Compounds, Russian Academy of Sciences, 142432 Chernogolovka, Russia</footnote>

<footnote n="4">⁴ Institut de Microbiologie de la Méditerranée, CNRS, FR3479, Aix-Marseille Université, 13402 Marseille, France; Deborah.byrne@univ-amu.fr</footnote>

<footnote n="5">⁵ Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry, Russian Academy of Sciences, 117997 Moscow, Russia; aozalevsky@fbb.msu.ru (A.O.Z.); golovin@fbb.msu.ru (A.V.G.)</footnote>

<footnote n="6">⁶ Institute for Biological Instrumentation, Pushchino Scientific Center for Biological Research of the Russian Academy of Sciences, 142290 Pushchino, Russia; lisiks.av@gmail.com (A.A.V.); 212sok@gmail.com (A.S.S.); permyakov.s@gmail.com (S.E.P.)</footnote>

<footnote n="7">⁷ Faculty of Bioengineering and Bioinformatics, Lomonosov Moscow State University, 119991 Moscow, Russia</footnote>

<footnote n="8">⁸ Sirius University of Science and Technology, 354340 Sochi, Russia</footnote>

<footnote n="9">⁹ Department of Biochemistry, Schulich School of Medicine & Dentistry, University of Western Ontario, London, ON N6A 3K7, Canada; gshaw1@uwo.ca</footnote>

<footnote>* Correspondence: philipp.tsvetkov@univ-amu.fr (P.O.T.); zerni@belozersky.msu.ru (E.Y.Z.)</footnote>
factors, thereby regulating the activity of different immune cells, chondrocytes, myoblasts, cardiomyocytes, epitheliocytes, astrocytes, glial cells and neurons (reviewed in [2]). Abnormalities in the expression of S100 proteins are found in a number of cancers [3]. Altered expression of brain-specific S100 proteins is associated with decreased motor and cognitive skills and neurological damage characteristic of Alzheimer's disease (AD) and Parkinson's disease, and several members of the S100 family are found in the amyloid plaques and cerebrospinal fluid of AD patients [4–7]. Secreted S100 proteins are commonly found in the aqueous humor of patients with glaucoma, which is accompanied by overexpression of these proteins in the damaged retina [8,9]. Despite extensive knowledge of the cellular and tissue functions of S100 proteins, understanding of the exact mechanisms underlying their signaling activity remains scarce, especially regarding Ca²⁺-sensitivity and interrelations between multiple isoforms of these proteins.

Typical representatives of S100 proteins are S100A1 and S100B, expressed in the skeletal muscles and heart, as well as in the neurons, astrocytes and microglia, where they exist as homo- and heterodimers [10–12]. In the absence of a target, both of these proteins exhibit very low calcium affinity: K<sub>D</sub> exceeds 10 µM, which is outside the intracellular concentration range for this cation [13]. There are several possible explanations for how these proteins can sense alterations in physiological concentrations of calcium. One common hypothesis is that the affinity of their Ca²⁺-coordinating sites can be adjusted via allosteric regulation provided by target binding [14]. Consistently, in the presence of one such target, the ryanodine receptor, S100A1 exhibits 100-fold increased Ca²⁺ affinity [15]. Another option is the binding of Zn²⁺, which was suggested to increase the Ca²⁺-affinity of certain S100 proteins. In recent years, Zn²⁺ has been recognized as the second major signaling ion in the nervous system (in addition to Ca²⁺) [16]. It regulates the normal functioning of the neurons and participates in the pathogenesis of neurological and neuro-ophthalmological disorders, predominantly AD and glaucoma [17–19]. For instance, Zn²⁺ can act as a signaling ion in the nervous system by binding to Ca²⁺-sensor proteins and modifying their structural and functional properties, such as Ca²⁺-sensitivity, affinity to cellular membranes or target proteins, and a propensity for non-covalent/covalent (disulfide) dimerization and aggregation [20–22]. In the S100B homodimer, the binding of Zn²⁺ greatly improves the Ca²⁺-affinity of the protein, promoting its activation by Ca²⁺ signals and the recognition of regulatory targets [13,23–25]. At the same time, Zn²⁺ has been reported to have no such effect on S100B/S100A1 heterodimers, which leads to the suggestion that S100A1 serves as an inhibitor of the Zn²⁺-dependent function of S100B [13]. S100A1 is also able to function independently from S100B, as it forms homodimers and can directly bind Ca²⁺ and Zn²⁺ [13]. Notably, previous studies have shown that the metal-binding properties of S100A1 drastically differ from those of S100B, pointing to its unique mode of its function.

In this study, using electrospray mass-spectrometry, isothermal titration calorimetry (ITC), dynamic light scattering (DLS), differential scanning fluorimetry (nanoDSF), and QM/MM molecular modeling, we investigated the mechanisms of Ca²⁺ and Zn²⁺ coordination in S100A1, including the structure and interrelations of its metal-binding sites. Our findings suggest that S100A1 can function as a Ca²⁺/Zn²⁺-sensor protein, mediating the fluctuations in both ions as a part of the diverse signaling pathways in the nervous system. Moreover, we demonstrate that an excess of Zn²⁺ produces a destabilizing effect on S100A1 structure and leads to the aggregation of the protein, which could contribute to its neurotoxicity under the pathological conditions of glaucoma and AD.

## 2. Materials and Methods

### 2.1. *Protein Purification*

The recombinant S100A1 protein was purified from *Escherichia coli* culture, as described previously [26]. Briefly, the bacterial extract was subjected to two steps of chromatography: anion exchange chromatography on a HiTrap Q column (GE Healthcare, Chicago, IL, USA), using a linear gradient of 0–1 M NaCl; and Ca²⁺-dependent affinity chromatography on a phenyl-Sepharose column (GE Healthcare), where the S100A1-containing fractions from
the first step were pooled and loaded on the column in 5 mM CaCl₂, and the pure protein was eluted with 5 mM ethylene glycol-bis(β-aminoethyl ether)-N,N,N',N'-tetraacetic acid (EGTA) (Sigma-Aldrich, Saint-Loius, MO, USA). Residual Ca²⁺ was removed from the protein according to the previously described procedure [20,27]. For circular dichroism (CD) studies (see below), recombinant S100A1 was obtained as described previously [28].

## 2.2. Electrospray Mass-Spectrometry (ESI-MS)

Mass spectrometric measurements of 15 µM of S100A1 in the presence of 200 µM of zinc or calcium ions were performed on an LCQ Deca XP ion trap mass spectrometer (Thermo Finnigan, San-Jose, CA, USA) equipped with an electrospray ionization source, as described previously [29]. All experiments were performed in a 10 mM ammonium bicarbonate buffer (pH 7.0). The infusion flow rate was 0.5 µL/min, the spray voltage was 4 kV, and the nitrogen sheath gas setting was 20 au. Full scans of m/z 1500–4000 were acquired using the Excalibur 1.3 data system (Thermo Finnigan, San-Jose, CA, USA).

## 2.3. Isothermal Titration Calorimetry

The binding of ions to the S100A1 protein was studied on an iTC200 isothermal titration calorimeter (MicroCal, Los Angeles, CA, USA) at 37 °C, in 50 mM Tris buffer and 1 mM tris(2-carboxyethyl)phosphine (TCEP) at pH 7.5, as described previously [20,30]. The protein concentration in the calorimetric cell was 50 µM whereas the concentration of zinc and calcium ions in the titrating syringe was 0.5 mM. The volume of the injections was 0.5 µL before the ion concentration reached the equimolar ratio with protein, and then set to 2.5 µL. The dilution heat was measured by titration of the ion-containing buffer into the same buffer without the protein. The obtained curve was analyzed by MicroCal Origin 7.0 software (Los Angeles, CA, USA) using the models of the two types of binding sites. Therefore, the stoichiometry (N), constant (K_a), enthalpy (ΔH), and entropy of binding (ΔS) were calculated from the standard thermodynamic equations.

## 2.4. Dynamic Light Scattering (DLS)

All DLS experiments were carried out at 37 °C using a Zetasizer Nano ZS (Malvern Instruments, Malvern, UK) as described previously [31]. S100A1 size was determined in the presence of 0.5–10 mM of Zn²⁺ or Ca²⁺ or in the absence of divalent metal ions. The protein concentration was 50 µM. The hydrodynamic diameter was determined using the software provided by the manufacturer. The dispersant viscosity was set to 0.6909 cP and the protein refractive index was set to 1.450 (at 37 °C). Several replicate measurements were performed; each one consisting of 10–15 runs of 10 s.

## 2.5. Differential Scanning Fluorimetry (nanoDSF) and Light Scattering

Protein thermostability was estimated using a label-free fluorimetric analysis in a Prometheus NT.Plex instrument (NanoTemper Technologies, Munich, Germany) as described in our previous works [20,31]. A total of 37 µM of S100A1 (the concentration corresponds to the protein concentration in a calorimetric cell at the end of ITC titration) was loaded into NanoDSF grade capillaries (NanoTemper Technologies, Munich, Germany) in the presence of either 1 mM Ca²⁺, 1 mM Zn²⁺, or neither. The concentration of zinc or calcium ions varied from 18.5 to 296 µM (excesses 0.5 to 8). The capillaries were loaded onto the instrument and heated from 20 °C to 110 °C at a 1 K/min heating rate with the excitation laser power set at 100%. The unfolding transition points (T_m) were determined from the first derivative of the changes in the emission wavelengths of ratio of tryptophan fluorescence intensities at 350 nm and 330 nm (I_350/I_330), which were automatically identified by the pre-installed PR.ThermControl software (NanoTemper Technologies, Munich, Germany). In addition, the aggregation temperatures (T_agg) were determined based on first derivative of temperature dependence of light scattering at 350 nm, using the same equipment.
## 2.6. Circular Dichroism

Circular dichroism (CD) studies were carried out with a J-810 spectropolarimeter JASCO Inc., Tokyo, Japan), equipped with a Peltier-controlled cell holder, as previously described [32]. Protein concentration was 34–37 µM in 50 mM Tris-HCl, 1 mM TCEP pH 7.5. A total of 1 mM CaCl₂ (for Ca²⁺-loaded form) or 1 mM ZnCl₂ (for Zn²⁺-loaded form) was added to the metal-depleted protein and pre-incubated for 5 min. A total of 296 µM ZnCl₂ (Zn²⁺:protein ratio 8:1) was added to the Ca²⁺-loaded form or 296 µM CaCl₂ (Ca²⁺:protein ratio 8:1) was added to the Zn²⁺-loaded form. The buffer contribution was subtracted from the experimental spectra. Quantitative estimates of the secondary structure contents were made using the CDPro software package (Colorado State University, Fort Collins, CO, USA) [33]. The finalized secondary structure fractions represent the averaged values from several measurements.

## 2.7. QM/MM Molecular Modeling

Quantum mechanics/molecular mechanics (QM/MM) simulations using the GROMACS-DFTb/Plumed package (https://www.gromacs.org/) were performed as described in [20,34]. The starting conformation of the S100A1 protein was obtained from the NMR solution structure (PDB 2LP3). Ca²⁺ was replaced with Zn²⁺ in two EF-hand sites of chain A. All residues harboring oxygen atoms in proximity to Zn²⁺ were included in the QM part of the system. In the case of oxygen atoms originating from the sidechain, QM/MM partitioning was performed on the CA-CB bond. For the backbone oxygens, partitioning was performed on the CA-C and CA-N bonds. The simulation system was filled with water molecules (presented by the tip3p model), and Na⁺ or Cl⁻ ions to neutralize the charge. The water and ions were equilibrated around the protein-Zn²⁺ complex by carrying out a 2ns MD simulation with a restrained position of protein. All atoms (including water) closer to Zn²⁺ than 3 Å were also included in the QM system. The total length of each simulation was set to 100 ps with a 0.2 fs timestep.

# 3. Results

## 3.1. Stoichiometry of S100A1-Metal Binding

The stoichiometry of Ca²⁺ and Zn²⁺ binding to S100A1 was first analyzed by electro-spray ionization mass-spectrometry (ESI-MS), as this method allows the direct monitoring of complexes of proteins with different metals. To optimize the detection of S100A1 complexes with cations, mass spectrometry experiments were carried out at pH 7.0 in the presence of 200 µM Ca²⁺ or 200 µM Zn²⁺, or both cations together [13]. In all experimental conditions, the main peak in the MS spectra was observed at 1175.8 m/z, which corresponds to the eighth charge state of the S100A1 species [S100A1⁺+8H⁺]⁹⁺ with a calculated average mass of 10,571.2 ± 1.4 Da (Figure 1). This value differs slightly from the theoretical mass (10,545.5 Da) due to the N-formylmethionine (fMet) modification at the N-terminus of the protein, which is unlikely to affect metal binding and therefore was not taken into account in the subsequent analysis. In the presence of calcium or zinc alone, S100A1 existed as five species corresponding to apo-protein and its conformers with one, two, three or four bound cations, the content of which decreases in the order 1Me²⁺-S100A1 > 2Me²⁺-S100A1 > 3Me²⁺-S100A1 > 4Me²⁺-S100A1 (Figure 1A,B). Although this pattern does not directly reflect the S100A1 state in solution due to the destabilizing ESI-MS-related effects, it points to the existence of high/medium-affinity and low-affinity Ca²⁺/Zn²⁺-binding sites in the protein. The preincubation of S100A1 in the presence of both Ca²⁺ and Zn²⁺ yielded multiple conformers with maximum four metal ions per protein with the predominance of single-ion and two-ion forms, including a substantial amount of a mixed Zn²⁺:S100A1:Ca²⁺ form (Figure 1C). All calcium conformers contained no more than two Ca²⁺, which are most likely coordinated in the EF-hands [35,36]. Notably, the only form of S100A1 with four ions contains two Ca²⁺ and two Zn²⁺ per monomer (2Zn²⁺:S100A1:2Ca²⁺). Given the destabilizing conditions of ESI-MS experiments and the absence of other four ion-bound variants of S100A1, this form seems to represent a relatively stable conformer of the protein.
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Mass-spectra of Ca²⁺- and Zn²⁺-bound S100A1. Samples containing 15 µM of S100A1 preincubated with (A) 200 µM Ca²⁺, (B) 200 µM Zn²⁺ or (C) both, were analyzed with electrospray mass spectrometry. Calcium ions are shown as green circles and zinc ions are shown as purple circles.</figcaption></figure>

## 3.2. Affinity of S100A1-Metal Binding

To characterize Zn²⁺ and Ca²⁺ binding to S100A1, the formation of the complexes was monitored using ITC (Figure 2) by the titration of 50 µM apo-S100A1 with 500 µM Ca²⁺ or 500 µM Zn²⁺. The fitting of the resulting titration curves using a “two sets of sites” model revealed the existence of one high- and one medium-affinity binding site for both ions (Table 1). It should be noted that both sites have two orders of magnitude higher association constants for Zn²⁺ than for Ca²⁺. For both metals, binding in the high-affinity site is enthalpy driven (ΔH < 0 and ΔS < 0), whereas binding in the medium-affinity site is entropy driven (ΔH > 0 and ΔS > 0). This indicates different conformational changes of the protein upon ion binding to either high or medium-affinity sites.
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Isothermal titration calorimetry (ITC) analysis of ions binding to S100A1. A total of 50 µM of S100A1 was titrated by (A) 500 µM Ca²⁺ or (B) 500 µM Zn²⁺. Upper panels show ITC curves and lower panels show binding isotherms. Best fits are shown as solid black curves (see Table 1).</figcaption></figure>

<table><caption>**Table 1.** Thermodynamic parameters of Ca²⁺ and Zn²⁺ binding to S100A1.</caption><thead><tr><th>Ion</th><th>N<sup>1</sup></th><th>K<sup>1</sup><sub>a</sub>, M<sup>-1</sup></th><th>K<sup>1</sup><sub>D</sub>, M</th><th>ΔH<sup>1</sup>, kcal M<sup>-1</sup></th><th>ΔS<sup>1</sup>, cal M<sup>-1</sup> K<sup>-1</sup></th><th>N<sup>2</sup></th><th>K<sup>2</sup><sub>a</sub>, M<sup>-1</sup></th><th>K<sup>1</sup><sub>D</sub>, M</th><th>ΔH<sup>2</sup>, kcal M<sup>-1</sup></th><th>ΔS<sup>2</sup>, cal M<sup>-1</sup> K<sup>-1</sup></th></tr></thead><tbody><tr><td>Ca<sup>2+</sup></td><td>0.9</td><td>6.1 ± 2.1 × 10<sup>6</sup></td><td>1.6 ± 0.6 × 10<sup>-7</sup></td><td>-16.9 ± 0.2</td><td>-23.3</td><td>1.0</td><td>4.2 ± 0.6 × 10<sup>4</sup></td><td>2.4 ± 0.3 × 10<sup>-5</sup></td><td>4.0 ± 0.8</td><td>34.1</td></tr><tr><td>Zn<sup>2+</sup></td><td>0.9</td><td>2.6 ± 0.7 × 10<sup>8</sup></td><td>3.8 ± 1.0 × 10<sup>-9</sup></td><td>-13.2 ± 0.2</td><td>-4.2</td><td>0.9</td><td>1.3 ± 0.4 × 10<sup>6</sup></td><td>7.7 ± 2.4 × 10<sup>-7</sup></td><td>3.3 ± 0.4</td><td>38.5</td></tr><tr><td>Ca<sup>2+</sup> in Zn<sup>2+</sup><br>Zn<sup>2+</sup> in Ca<sup>2+</sup></td><td></td><td></td><td></td><td></td><td></td><td></td><td colspan="3">no binding<br>no binding</td><td></td></tr></tbody></table>

All parameters were determined using the “two sets of sites” model.

Since it is known that the two main calcium-binding sites are located in the two EF-hands of S100A1, ion competition experiments allowed us to establish whether the EF-hands of the protein can also chelate zinc ions. The presence of zinc at saturating concentrations results in the absence of heat exchange upon protein titration by calcium (and vice versa). This observation strongly suggests the existence of competition between Zn²⁺ and Ca²⁺ for both sites at S100A1, thus indicating that the Zn²⁺-binding sites are also located in the EF-hands of the protein.

Moreover, since the binding of Ca²⁺ to its high-affinity site is stronger than the binding of Zn²⁺ to its medium-affinity site (Table 1) and since, in the presence of equimolar concentrations of both ions, the Zn²⁺:S100A1:Ca²⁺ species is much more abundant than the 2Zn²⁺:S100A1 or S100A1:2Ca²⁺ species (Figure 1C), we can conclude that high-affinity sites for these cations are localized in different EF-hands.
### 3.3. Effect of Ca²⁺ and Zn²⁺ on Secondary Structure of S100A1

Possible alterations to secondary structure induced by Ca²⁺ and Zn²⁺ binding to S100A1 were assessed by registering the circular dichroism (CD) spectra of the protein in the excess of these ions (Supplementary Materials, Figure S1). It was found that the addition of either 1 mM Ca²⁺ or 1 mM Zn²⁺ to S100A1 induces a moderate reduction to its α- and β-structures and increased its share of turns and disordered elements (Table 2). As a result, the Ca²⁺- and Zn²⁺-loaded forms of the protein exhibited similar patterns of secondary structure elements. Thus, the filling of S100A1 with calcium or zinc yielded conformers with similar secondary structures, both of which slightly differed from the apo from.

<table><caption>**Table 2.** Secondary structure fractions estimated from CD spectra for apo, Ca²⁺-loaded (1 mM CaCl₂), and Zn²⁺-loaded (1 mM ZnCl₂) forms of S100A1.</caption><thead><tr><th>Protein State</th><th>α-Helices, %</th><th>β-Structure, %</th><th>Turns, %</th><th>Unordered Structure, %</th></tr></thead><tbody><tr><td>apo</td><td>70.58 ± 0.56</td><td>5.17 ± 0.09</td><td>8.01 ± 0.40</td><td>16.02 ± 0.25</td></tr><tr><td>Ca²⁺-loaded</td><td>68.05 ± 0.56</td><td>4.05 ± 0.09</td><td>9.28 ± 0.40</td><td>18.65 ± 0.25</td></tr><tr><td>Zn²⁺-loaded</td><td>65.10 ± 0.56</td><td>4.13 ± 0.09</td><td>11.17 ± 0.40</td><td>19.50 ± 0.25</td></tr></tbody></table>

### 3.4. Putative Zn²⁺-Binding Sites in S100A1

The results of the ESI-MS and ITC studies indicated that (1) Zn²⁺ can bind to the EF-hands of S100A1 and (2) this binding is characterized by nanomolar affinity. To assess the structural aspects underlying this high-affinity interaction, we next predicted the Zn²⁺-binding sites in the EF-hands of S100A1 in silico, based on the available NMR structure of the Ca²⁺-bound protein (PDB 2LP3 [36]), since, as shown by our CD experiments (Table 1), the Ca²⁺-bound form of the protein is similar to its Zn²⁺-bound form in secondary structure. To this end, we performed QM/MM simulations of the molecular dynamics associated with Zn²⁺ binding instead of calcium in both EF-hand sites of chain A. Prior to the simulations, the solvent components, namely water molecules and sodium and chloride ions, were added to the system to neutralize the charge and create a physiological ionic strength of 0.14 M. Although initially the structure lacked water molecules, the balancing of the positions of such molecules around the protein resulted in their embedding in both sites in the vicinity of the zinc ions. The N-terminal and C-terminal EF-hands acquired one and three water molecules, respectively. The residues of both metal-binding sites containing Zn²⁺ and the respective water molecules were included in separate QM systems and subjected to four independent molecular dynamics simulations in 500,000 steps.

As can be seen in Figure 3A, the original coordination of calcium in the medium-affinity N-terminal EF-hand is achieved by backbone oxygens of S19, E22, D24, and K27, as well as single carboxyl oxygen of E32. The sixth coordinator might be a water molecule, which is absent in the NMR structure. In the case of zinc binding (Figure 3B), the oxygen atom of the K27 backbone is replaced by the water molecule from the corresponding coordination sphere, and the position of this molecule is stabilized by the backbone oxygen of K27 and the second carboxyl oxygen of E32. At the same time, the OG atom of S19 enters the coordination sphere, thereby creating an octahedral environment favorable for zinc binding. By contrast, in the C-terminal EF-hand, Ca²⁺ is favorably coordinated by the carboxyl groups of D62, D66, and E73 as well as the backbone oxygens of N64 and E68 (a water molecule is absent in the structure of NMR) (Figure 3C). This optimal coordination does not change during the simulation with zinc, except that the respective water molecule, stabilized by carboxyl atoms of the D66 and E73 residues completes the octahedral coordination of this metal.
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Coordination spheres of Ca²⁺ and Zn²⁺ in N-terminal (A,B) and C-terminal (C,D) EF-hands of S100A1. (A,C) NMR structure of Ca²⁺-bound S100A1 (PDB 2LP3). Calcium ions are colored green. Oxygen atoms in the coordination sphere are colored red (protein) or orange (water). The distances between chelator atoms and the metal ion are denoted in angstroms. (B,D) QM/MM models of Zn²⁺ coordination. Zinc ions are colored purple. The other designations are the same as in (A,C).</figcaption></figure>

In the aggregate, our simulations suggest that the N-terminal EF-hand of S100A1 is better suited for Zn²⁺, whereas the C-terminal EF-hand seems to be a fit for both metals with a slight advantage of Ca²⁺.

### 3.5. Effects of Metal Binding on Dimeric Structure of S100A1

Given the results of the metal-binding experiments, S100A1 may exist as three main states (apo, Ca²⁺-loaded, Zn²⁺-loaded) and several mixed states (mainly Zn²⁺:S100A1:Ca²⁺ and 2Zn²⁺:S100A1:2Ca²⁺), all of which seem physiologically relevant, especially considering both the intracellular and extracellular localizations of the protein. As such, we next analyzed the multimeric structure, conformational properties and stability of these forms. The examination of the hydrodynamic properties of S100A1 by DLS revealed that in the absence of metals, its average volume diameter (Dv50) is 5.63 ± 0.422 d.nm for 6 independent measurements, confirming that the protein exists as the dimer (Figure 4). The presence of increasing concentrations of Ca²⁺ (up to 10 mM) or Zn²⁺ (up to 4 mM) does not significantly affect the average diameter of S100A1, indicating the persistence of the dimer (Figure 4A,B). The notable change in protein size, namely a decrease down to 4.6 ± 0.4 d.nm, was observed only in the presence of 10 mM Zn²⁺, apparently reflecting
complete filling of all low-affinity sites of the protein and suggesting a more compacted spherical form (similar shape was observed for homologous protein S100A11 in its apo-form; PDB: 2LUC) (Figure 4C).

<figure><img src="image_7.png" /></figure>

**Figure 4.** Dynamic light scattering (DLS) analysis of Ca²⁺- and Zn²⁺-loaded S100A1. A total of 50µM of S100A1 was incubated with increasing concentrations (0.5–10 mM) of (A) Ca²⁺ and (B) Zn²⁺ at 37 °C and volume distribution of S100A1 hydrodynamic diameter was determined. Lower concentrations are shown as darker coloured lines and higher concentrations are shown as lighter coloured lines. (C) The influence of increasing divalent ions on the major peak size of triplicate concentrations.

3.6. *Effects of Metal Binding on Conformational Properties of S100A1*

To confirm the existence of the suggested states of S100A1, we analyzed its conforma-
tional changes induced by Zn²⁺ and/or Ca²⁺, by measuring fluorescence of the protein
at 350 and 330 nm (I₃₅₀ and I₃₃₀) and 37 °C using nanoDSF. The ratio of I₃₅₀/I₃₃₀, which
correlates with the environment of the only tryptophan (Trp90) located close to the C-
terminus, was plotted as the function of the molar ratio [Me²⁺]/[S100A1] (Figure 5). The
changes in the conformation of S100A1 upon the titration of its apo-form with Ca²⁺ and
Zn²⁺ up to four-fold excess could be associated with metal coordination by the two EF-
hand sites. Indeed, we can observe the saturation of these sites at a four-fold excess in
the ITC experiments (Figure 2). The filling of the first high-affinity site is accompanied by
pronounced conformational changes (molar ratio 0→1 in Figure 5), while the filling of the
medium-affinity site does not significantly impact tryptophan fluorescence (molar ratio
1→4 in Figure 5). Intriguingly, the binding of Zn²⁺, despite its higher affinity, produces
moderate effects on protein conformation, as seen from the smaller change in the I₃₅₀/I₃₃₀
ratio.

The absence of conformational changes upon the titration of the Zn²⁺-saturated protein by Ca²⁺ was in good agreement with our ITC data, thus confirming that the first two zinc ions fill and block the EF-hands. In contrast, the decrease of the I₃₅₀/I₃₃₀ ratio upon titration of Ca²⁺-loaded protein by Zn²⁺ points on the latter binding to Ca²⁺-loaded S100A1. At a molar ratio 4 it results in formation of a conformer that is similar to S100A1:2Ca²⁺. Apparently, this conformer contains two calcium ions in the EF-hands and 1-4 zinc ions which replaced Ca²⁺ in some nonspecific sites outside the EF-hands. Given the results of the ESI-MS experiments (see Figure 1C), this conformer most likely represents 2Zn²⁺:S100A1:2Ca²⁺. Thus, we conclude that each S100A1 monomer competitively binds two Ca²⁺ or two Zn²⁺ in the EF-hands, and can competitively bind these ions in low-affinity (nonspecific) sites outside the EF-hands, which display a preference for Zn²⁺.
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Effects of Ca²⁺ and Zn²⁺ on conformational properties of S100A1. 37 µM of S100A1 in a buffer containing Ca²⁺, Zn²⁺, or their combinations was loaded into capillaries and intrinsic tryptophan fluorescence of the protein was registered at 37 °C. The figure shows the influence of increasing concentrations of metal ions (37–296 µM) on fluorescence intensities at 350 nm and 330 nm. The titration of apo-S100A1 with Ca²⁺ is shown as green circles and black lines; titration of apo-S100A1 with Zn²⁺ is shown as purple circles with black lines; titration with Zn²⁺-saturated S100A1 with Ca²⁺ is shown as green circles with purple lines; titration of Ca²⁺-saturated S100A1 with Zn²⁺ is shown as purple circles with green lines.</figcaption></figure>

3.7. *Effects of Metal Binding on Thermal Stability and Aggregation of S100A1*

Since S100A1 contains only one tryptophan residue, the changes in the fluorescence
of the protein reflect only local conformational modifications in the protein structure.
Thus, to investigate the impact of zinc and calcium ions on the global protein structure,
we examined thermal stability of S100A1 by monitoring the temperature dependence of
intrinsic tryptophan fluorescence at 350 nm and 330 nm in the presence of increasing
concentrations of Ca²⁺ and Zn²⁺ (Figure 6). Interestingly, although the binding of the first
calcium ion in the high-affinity site induced important conformational changes in the local
environment of Trp90, it did not affect the thermostability of the protein (Figure 6A,D).
In contrast, the binding of the second calcium ion in the medium-affinity site did not
impact the Trp90 environment, but did progressively stabilize the S100A1 structure. As for
Zn²⁺, its binding to both high- and medium-affinity sites progressively destabilizes S100A1
structure, but this effect was more pronounced for the high-affinity site (Figure 6A,D).
The titration of the Zn²⁺-saturated protein by Ca²⁺ does not affect its stability, thereby
confirming the absence of Ca²⁺ binding (Figure 6B,D). By contrast, the titration of the
Ca²⁺-saturated protein by Zn²⁺ reduces its stability, but no more than to the baseline level,
similar to that of the apo-form (~71 °C) (Figure 6C,D). Given that Zn²⁺ binding to high-
affinity sites in the EF-hands would result in a more pronounced drop in the temperatures
of denaturation (Figure 6D), this effect indicates that Zn²⁺ replaces Ca²⁺ only from low-
affinity sites, yielding a 2Zn²⁺:S100A1:2Ca²⁺ form. We also determined the aggregation
temperatures of these forms by monitoring light scattering at 350 nm, which yielded very
similar results (Figure S2).
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Effect of Ca²⁺ and Zn²⁺ on the thermal stability of S100A1. A total of 37 µM of S100A1 in a buffer containing Ca²⁺, Zn²⁺, or their combinations was loaded into the capillaries and the intrinsic tryptophan fluorescence of the protein was monitored at a temperature range of 20–110 °C. (A–C) Temperature dependence of the first derivative of I₃₅₀/I₃₃₀ for (A) decalcified S100A1 in the presence of 37–296 µM Ca²⁺ (green) or Zn²⁺ (purple), (B) Zn²⁺-saturated S100A1 in the presence of 37–296 µM Ca²⁺, (C) Ca²⁺-saturated S100A1 in the presence of 37–296 µM Zn²⁺. (D) Mid-transition temperatures (Tₘ) for S100A1 in presence of different ions, determined from the first derivative of I₃₅₀/I₃₃₀.</figcaption></figure>

**4. Discussion**

This study reports new data regarding the mechanisms of Ca²⁺ and Zn²⁺ binding to S100A1, a crucial multitarget regulatory protein abundant in the heart, nervous system and the retina. Using ESI-MS, we directly demonstrate the existence of the Ca²⁺- and Zn²⁺-containing complexes of S100A1. Indeed, this method is well suited for establishing the stoichiometry of metals binding to proteins [37,38]. Our results indicate that the S100A1 homodimer binds up to 4 calcium ions per monomer, although the most abundant forms contain one or two calcium ions bound to EF-hands, which is in agreement with the available NMR and crystal structures of the protein [35,36]. According to the results of the ITC studies (performed in the presence of TCEP to reduce possible disulfide complexes of the protein [39,40]), S100A1 binds first two calcium ions with submicromolar (K_D = 164 nM) and micromolar (K_D = 24 µM) affinities. The second constant corresponds to the site located in the N-terminal EF-hand, which is sometimes called the pseudo-EF-hand, due to weak calcium coordination [10]. The additional third and fourth calcium ions are captured outside of the EF-hands, and their binding is non-specific, which is supported by the ITC
data. These observations are generally in agreement with the early data obtained using flow dialysis, although in our study the affinity of the best Ca²⁺-binding site is found to be two orders of magnitude higher than in previous assessments [13].

In the case of Zn²⁺, we observed a similar stoichiometry, but with stronger binding: the best two sites in each of S100A1 monomers exhibited nanomolar affinities (K_D ~ 4 and 770 nM). These findings differ drastically from the data obtained in previous work, where the dissociation constants of Zn²⁺-S100A1 complexes were estimated to be > 10 µM [13]. This discrepancy may be related to the thionylation of S100A1 by β-mercaptoethanol, which is known to alter the Ca²⁺-sensitivity of the protein and may affect its affinity to Zn²⁺ [39]. A more recent study employing ITC revealed only one Zn²⁺-binding site in S100A1, with K_D of 3.5 µM [41]. We suppose that the failure in determining the second high-affinity Zn²⁺-binding site was due to the insufficient decalcification of S100A1, since, according to our observations, its incubation with Chelex-100 resin (the method used for the decalcification of S-100 proteins in [41]) does not completely remove calcium from the protein.

Interestingly, our ITC and nanoDSF experiments demonstrate that Ca²⁺ and Zn²⁺ compete for the high-affinity binding sites in S100A1, thereby confirming that both metals bind to EF-hands. Moreover, the highest affinities were registered in the case of the first Zn²⁺-binding site (N-terminal EF-hand) and the first Ca²⁺-binding site (C-terminal EF-hand), pointing to the existence of a mixed Zn²⁺:S100A1:Ca²⁺ form of the protein. Indeed, it dominates in MS spectra, along with Ca²⁺-loaded and Zn²⁺-loaded conformers. The difference in the Ca²⁺/Zn²⁺ affinities of the EF-hand sites can be explained by comparing the coordination spheres of these cations based on NMR data and the results of the QM/MM simulations performed in our study [36]. The key dissimilarity between these coordination spheres is the number of charged acidic residues, namely one carboxyl in the N-terminal EF-hand and three carboxyls in the C-terminal EF-hand. The binding of calcium in the N-terminal EF-hand is less favorable due to suboptimal coordination by backbone oxygen atoms and the low compensation of positive charge of the metal in the absence of negatively charged ligands. In the case of zinc binding, one of the backbone oxygens is replaced with a water molecule (stabilized by this oxygen and the carboxyl group of E32) which, together with the oxygen from the S19 side group, ensures octahedral coordination, making Zn²⁺ tightly bound through the maximal number of chelators (coordination number of Zn²⁺ in proteins is 4–6) [42]. In the C-terminal EF-hand, Zn²⁺ is coordinated in a similar fashion by capturing all five chelators and the water molecule. In this site, Ca²⁺ also binds in optimal coordination by these six chelators, including three carboxyls. Notably, the principle difference between Zn²⁺ and Ca²⁺ is the more hindered dehydration of Ca²⁺ when chelating with amide or carbonyl oxygen atoms (analogs of backbone oxygen atoms) as compared to chelation by carboxyl oxygen atoms [43,44]. Overall, Zn²⁺ binding is by default more favorable in the N-terminal EF-hand of S100A1, whereas Ca²⁺ binding is slightly more favorable in the C-terminal EF-hand, although in the latter case the affinities are comparable.

These new data significantly update our current understanding of the composition and behavior of the S100A1 metal-bound forms under physiological conditions. It is known that the extracellular concentration of Ca²⁺ is 1.5–2 mM, whereas the intracellular levels of this metal range from 50 to 100 nM in the mammalian brain and up to 600 nM in the retina [45,46]. The internal concentration of Zn²⁺ in the neurons is at a low picomolar level, whereas its external concentration in some cases may reach up to 300 µM, such as in the synaptic clefts of glycinergic neurons [47]. In addition, the level of free Zn²⁺ may aberrantly increase outside the synapses and provide neurotoxicity in neurodegenerative diseases, such as AD and glaucoma [18,48]. Given our findings regarding the much higher Ca²⁺-sensitivity of S100A1, one can suspect that the protein could actually be sensitive to Ca²⁺ signals within a cell, even in the absence of its targets. Furthermore, the nanomolar affinity of S100A1 to Zn²⁺ revealed in our study indicates that it could function as a Zn²⁺-signaling protein by exchanging the so-called loosely-bound zinc (K_D ~ 10⁻⁸–10⁻⁶ M⁻¹) with other
proteins in response to stimuli [19,49]. Thus, apo, Ca²⁺-loaded and Zn²⁺-loaded S100A1 seem to be the most probable states of the protein within a cell. In turn, being secreted into the extracellular space, S100A1 might adopt a Ca²⁺-loaded conformation, as the calcium concentration outside the cell is 2 orders of magnitude higher than the K_D for the medium-affinity Ca²⁺-binding site of the protein [50,51]. Yet, even Ca²⁺-loaded S100A1 can respond to increased free Zn²⁺ by forming 2Zn²⁺:S100A1:2Ca²⁺ and can thereby participate in the pathophysiological processes leading to neurodegeneration.

How would high-affinity Ca²⁺ and Zn²⁺ binding reflect on the stability and conformational states of S100A1? To answer this question, we first performed DLS studies, as this method allows the monitoring of the effect of cations on the global structural integrity and multimerization of proteins [52]. Our results demonstrate that, at physiological temperature (37 °C), S100A1 remains as a dimer regardless of the Ca²⁺ or Zn²⁺ concentrations although this dimer adopts a more compact conformation at high (10 mM) Zn²⁺ levels. However, according to intrinsic fluorescence (nanoDSF) studies, Ca²⁺-bound and Zn²⁺-bound S100A1 already adopt different conformations at the binding of the first metal ion per monomer. The most pronounced change was observed upon binding of the first two ions i.e., in the high-affinity EF-hand sites, in agreement with previous estimations [13]. Notably, the binding of Ca²⁺ increases the thermal stability of S100A1, whereas Zn²⁺ significantly decreases the temperatures of its denaturation and aggregation. Thus, the most probable physiological forms of S100A1, namely apo, Ca²⁺-loaded and Zn²⁺-loaded, differ in overall structure, conformational flexibility and environmental resistance, and therefore may display different physiological activities. The validation and structural characterization of the revealed additional S100A1 forms, such as Zn²⁺:S100A1:Ca²⁺ and 2Zn²⁺:S100A1:2Ca²⁺, requires further studies.

The discovery in this study of the new physiologically relevant Zn²⁺-loaded and Zn²⁺(Ca²⁺)-bound conformers of S100A1 suggests that Zn²⁺-dependent regulation could be a novel aspect in its function. Thus, Zn²⁺ can modulate both the Ca²⁺-dependent and Ca²⁺-independent activities of S100A1, which is known to play important roles in cells not only at high Ca²⁺, but also under resting Ca²⁺ conditions [53,54]. Zinc binding was shown to have a direct impact on the activity of S100A1. For instance, it enhances S100A1-dependent activation of twitchin kinase, although it is unclear if this effect is mediated by Zn²⁺ binding to Ca²⁺-loaded S100A1, the target protein, or both [55]. Some S100 proteins could function as metallochaperones by releasing Zn²⁺ to be taken up by their target proteins [10]. S100A1 binds p53, as well as a few other Zn²⁺-binding proteins, including tubulin, tau-protein, phosphoglucomutase, and different S100 proteins [56–58]. For instance, in the case of tau, the dissociation constant for Zn²⁺ is ~10⁻⁷ M, which is comparable to those determined for S100A1 in this work, which could theoretically lead to exchange of zinc ions between the two proteins in both directions. By contrast, in the case of tubulin, the constant is significantly lower (K_D ~ 10⁻⁵ M), which means that S100A1 could act as a Zn²⁺ acceptor rather than a donor when forming a complex with tubulin [56,59]. This suggestion is supported by the stabilizing effect of S100A1 on microtubules in the presence of Zn²⁺, which otherwise induces their disassembly [60].

Notably, our results show that the excessive binding of Zn²⁺ destabilizes S100A1 and increases its susceptibility to aggregation. Zinc-induced protein aggregation is a hallmark of neurodegenerative diseases such as Alzheimer's disease and amyotrophic lateral sclerosis [31,61]. Since S100A1 can be secreted in the extracellular space, where high amplitude Zn²⁺ fluctuations are observed, we propose that it can contribute to pathogenesis of Zn²⁺-dependent neurodegenerative diseases such as AD and glaucoma [50,51]. S100A1 is found in the amyloid plaques of AD patients and is known to contribute to plaque formation, neuroinflammation and cell death in AD models [62]. In addition, S100A1, together with amyloid peptides, was detected in the aqueous humor of glaucoma patients [8,9]. S100A1 expression is dysregulated in type I diabetes, which is characterized by altered Zn²⁺ homeostasis [63–65]. Furthermore, the upregulation of S100A1 is found in tumors such as melanoma (> 100-fold increase), renal cell carcinoma and ovarian cancer [65–68]. S100A1
was reported to control the migration and proliferation of cancer cells and its elevated expression is a known marker of tumor malignancy [65]. As cancer cells are characterized by redox imbalance and oxidative stress is a major source of excess Zn²⁺ in cells (the oxidation of sulfhydryl groups in Zn²⁺-buffer proteins leads to massive Zn²⁺ release), S100A1 can predominantly act as a Zn²⁺-saturated conformer in cancer cells [69]. This should be considered while searching for approaches to suppressing the pathological function of S100A1. Indeed, targeting the Zn²⁺-dependent metabolic and signaling pathways in tumors with specific chelators is regarded as a prospective strategy, because it could help sensitize cancer cells to established chemotherapy drugs [70,71].

In conclusion, this work demonstrates, for the first time, the high affinity binding of Ca²⁺ and Zn²⁺ by a S100A1 homodimer, which results in the formation of physiologically relevant, Ca²⁺/Zn²⁺-bound conformers of this protein that are structurally and, potentially, functionally different. New data acquired by direct analysis of the complex by ESI-MS mass-spectrometry and monitoring of the interactions by ITC shows that two calcium and two zinc ions compete for the same sites within the S100A1 structure. In correspondence with the results of ITC, novel Zn²⁺-binding sites within the EF-hands of the protein were predicted using QM/MM molecular modeling. Furthermore, two additional zinc ions can bind outside of the EF-hands, yielding the 4Zn²⁺:S100A1 or 2Zn²⁺:S100A1:2Ca²⁺ forms of the protein, which possess somewhat compromised structural stability. Although further studies are required for more detailed examination of the functional properties of the newly revealed S100A1 forms, the results of this study generally contribute to the understanding of Ca²⁺/Zn²⁺ interplay in cell signaling in health and disease.

**Supplementary Materials:** The following are available online at https://www.mdpi.com/article/10.3390/biom11121823/s1, Figure S1: Far-UV CD spectra of S100A1; Figure S2: S100A1 aggregation temperatures.

**Author Contributions:** Conceptualization, P.O.T.; methodology, G.S.S., A.Y.R., C.V., D.B., P.O.T., A.V.G., S.E.P.; software, A.V.G., A.O.Z.; validation, V.E.B., P.O.T.; formal analysis, A.Y.R., C.V., D.B., P.O.T., A.V.G.; investigation, A.Y.R., C.V., D.B., D.Y., V.E.B., A.O.Z., A.A.V., A.S.S., A.V.G.; resources, G.S.S.; data curation, V.E.B., D.B., P.O.T., A.V.G., A.O.Z.; writing—original draft preparation, V.E.B., E.Y.Z.; writing—review and editing, E.Y.Z., P.O.T., F.D., G.S.S.; visualization, C.V., P.O.T., A.V.G., A.O.Z.; supervision, F.D., S.E.P., E.Y.Z.; project administration, F.D., E.Y.Z.; funding acquisition, E.Y.Z., S.E.P. All authors have read and agreed to the published version of the manuscript.

**Funding:** Investigation of stoichiometry and affinity of S100A1-metal binding, putative Zn²⁺-binding sites in S100A1, as well as the effects of metal binding on the dimeric structure, conformational properties, thermal stability and aggregation of S100A1 (Figures 1–6; Table 1; Figure S2) was supported by the Russian Science Foundation (Grant No. 21-15-00123 to E.Y.Z.). The secondary structure of S100A1 (Table 2; Figure S1) was investigated with the support of the Russian Science Foundation (Grant No. 19-14-00289 to S.E.P.).

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Acknowledgments:** Microcalorimetry and nanoDSF experiments were performed in Timone Microcalorimetry Platform, PINT, Marseille, France. A Prometheus NT.Plex instrument was acquired with the financial support of Région PACA (APR-EX 2016 Neuro-Plasma).

**Conflicts of Interest:** The authors declare no conflict of interest.

**Abbreviations**

ESI-MS—electrospray mass spectrometry, ITC—isothermal titration calorimetry, DLS—dynamic light scattering, nanoDSF—nano-format of differential scanning fluorimetry, QM/MM—quantum mechanics/molecular mechanics.
References

1. Donato, R. Functional roles of S100 proteins, calcium-binding proteins of the EF-hand type. *Biochim. Biophys. Acta* **1999**, *1450*, 191-231. [CrossRef]

2. Donato, R.; Cannon, B.R.; Sorci, G.; Riazzi, F.; Hsu, K.; Weber, D.J.; Geczy, C.L. Functions of S100 proteins. *Curr. Mol. Med.* **2013**, *13*, 24-57. [CrossRef]

3. Bresnick, A.R.; Weber, D.J.; Zimmer, D.B. S100 proteins in cancer. *Nat. Rev. Cancer* **2015**, *15*, 96-109. [CrossRef]

4. Cristovao, J.S.; Gomes, C.M. S100 Proteins in Alzheimer's Disease. *Front. Neurosci.* **2019**, *13*, 463. [CrossRef] [PubMed]

5. Liu, J.; Wang, H.; Zhang, L.; Xu, Y.; Deng, W.; Zhu, H.; Qin, C. S100B transgenic mice develop features of Parkinson's disease. *Arch. Med. Res.* **2011**, *42*, 1-7. [CrossRef]

6. Mori, T.; Koyama, N.; Arendash, G.W.; Horikoshi-Sakuraba, Y.; Tan, J.; Town, T. Overexpression of human S100B exacerbates cerebral amyloidosis and gliosis in the Tg2576 mouse model of Alzheimer's disease. *Glia* **2010**, *58*, 300-314. [CrossRef]

7. Gerlai, R.; Wojtowicz, J.M.; Marks, A.; Roder, J. Overexpression of a calcium-binding protein, S100 beta, in astrocytes alters synaptic plasticity and impairs spatial learning in transgenic mice. *Learn. Mem.* **1995**, *2*, 26-39. [CrossRef] [PubMed]

8. Kliuchnikova, A.A.; Samokhina, N.I.; Ilina, I.Y.; Karpov, D.S.; Pyatnitskiy, M.A.; Kuznetsova, K.G.; Toropygin, I.Y.; Kochergin, S.A.; Alekseev, I.B.; Zgoda, V.G.; et al. Human aqueous humor proteome in cataract, glaucoma, and pseudoexfoliation syndrome. *Proteomics* **2016**, *16*, 1938-1946. [CrossRef] [PubMed]

9. Jiang, B.; Harper, M.M.; Kecova, H.; Adamus, G.; Kardon, R.H.; Grozdanic, S.D.; Kuehn, M.H. Neuroinflammation in advanced canine glaucoma. *Mol. Vis.* **2010**, *16*, 2092-2108.

10. Baudier, J.; Deloulme, J.C.; Shaw, G.S. The Zn²⁺ and Ca²⁺-binding S100B and S100A1 proteins: Beyond the myths. *Biol. Rev. Camb. Philos. Soc.* **2020**, *95*, 738-758. [CrossRef]

11. Duarte-Costa, S.; Castro-Ferreira, R.; Neves, J.S.; Leite-Moreira, A.F. S100A1: A major player in cardiovascular performance. *Physiol. Res.* **2014**, *63*, 669-681. [CrossRef]

12. Donato, R.; Heizmann, C.W. S100B Protein in the Nervous System and Cardiovascular Apparatus in Normal and Pathological Conditions. *Cardiovasc. Psychiatry Neurol.* **2010**, *2010*, 929712. [CrossRef]

13. Baudier, J.; Glasser, N.; Gerard, D. Ions binding to S100 proteins. I. Calcium- and zinc-binding properties of bovine brain S100 alpha alpha, S100a (alpha beta), and S100b (beta beta) protein: Zn²⁺ regulates Ca²⁺ binding on S100b protein. *J. Biol. Chem.* **1986**, *261*, 8192-8203. [CrossRef]

14. Zimmer, D.B.; Weber, D.J. The Calcium-Dependent Interaction of S100B with Its Protein Targets. *Cardiovasc. Psychiatry Neurol.* **2010**, *2010*, 728052. [CrossRef] [PubMed]

15. Treves, S.; Scutari, E.; Robert, M.; Groh, S.; Ottolia, M.; Prestipino, G.; Ronjat, M.; Zorzato, F. Interaction of S100A1 with the Ca²⁺ release channel (ryanodine receptor) of skeletal muscle. *Biochemistry* **1997**, *36*, 11496-11503. [CrossRef] [PubMed]

16. Frederickson, C.J.; Koh, J.Y.; Bush, A.I. The neurobiology of zinc in health and disease. *Nat. Rev. Neurosci.* **2005**, *6*, 449-462. [CrossRef] [PubMed]

17. Krall, R.F.; Tzounopoulos, T.; Aizenman, E. The Function and Regulation of Zinc in the Brain. *Neuroscience* **2021**, *457*, 235-258. [CrossRef]

18. Li, Y.; Andereggen, L.; Yuki, K.; Omura, K.; Yin, Y.; Gilbert, H.Y.; Erdogan, B.; Asdourian, M.S.; Shrock, C.; de Lima, S.; et al. Mobile zinc increases rapidly in the retina after optic nerve injury and regulates ganglion cell survival and optic nerve regeneration. *Proc. Natl. Acad. Sci. USA* **2017**, *114*, E209-E218. [CrossRef]

19. Ugarte, M.; Osborne, N.N. Recent advances in the understanding of the role of zinc in ocular tissues. *Metallomics* **2014**, *6*, 189-200. [CrossRef]

20. Tsvetkov, P.O.; Roman, A.Y.; Baksheeva, V.E.; Nazipova, A.A.; Shevelyova, M.P.; Vladimirov, V.I.; Buyanova, M.F.; Zinchenko, D.V.; Zamyatnin, A.A., Jr.; Devred, F.; et al. Functional Status of Neuronal Calcium Sensor-1 Is Modulated by Zinc Binding. *Front. Mol. Neurosci.* **2018**, *11*, 459. [CrossRef] [PubMed]

21. Permyakov, S.E.; Cherskaya, A.M.; Wasserman, L.A.; Khokhlova, T.I.; Senin, I.I.; Zargarov, A.A.; Zinchenko, D.V.; Zernii, E.Y.; Lipkin, V.M.; Philippov, P.P.; et al. Recoverin is a zinc-binding protein. *J. Proteome. Res.* **2003**, *2*, 51-57. [CrossRef]

22. Baksheeva, V.E.; Baldin, A.V.; Zalevsky, A.O.; Nazipova, A.A.; Kazakov, A.S.; Vladimirov, V.I.; Gorokhovets, N.V.; Devred, F.; Philippov, P.P.; Bazhin, A.V.; et al. Disulfide Dimerization of Neuronal Calcium Sensor-1: Implications for Zinc and Redox Signaling. *Int. J. Mol. Sci.* **2021**, *22*, 12602. [CrossRef]

23. Ostendorp, T.; Diez, J.; Heizmann, C.W.; Fritz, G. The crystal structures of human S100B in the zinc- and calcium-loaded state at three pH values reveal zinc ligand swapping. *Biochim. Biophys. Acta* **2011**, *1813*, 1083-1091. [CrossRef]

24. Wilder, P.T.; Varney, K.M.; Weiss, M.B.; Gitti, R.K.; Weber, D.J. Solution structure of zinc- and calcium-bound rat S100B as determined by nuclear magnetic resonance spectroscopy. *Biochemistry* **2005**, *44*, 5690-5702. [CrossRef]

25. Barber, K.R.; McClintock, K.A.; Jamieson, G.A., Jr.; Dimlich, R.V.; Shaw, G.S. Specificity and Zn²⁺ enhancement of the S100B binding epitope TRTK-12. *J. Biol. Chem.* **1999**, *274*, 1502-1508. [CrossRef]

26. Marlatt, N.M.; Spratt, D.E.; Shaw, G.S. Codon optimization for enhanced Escherichia coli expression of human S100A11 and S100A1 proteins. *Protein. Expr. Purif.* **2010**, *73*, 58-64. [CrossRef]

27. Wang, B.; Boeckel, G.R.; Huynh, L.; Nguyen, L.; Cao, W.; De La Cruz, E.M.; Kaftan, E.J.; Ehrlich, B.E. Neuronal Calcium Sensor 1 Has Two Variants with Distinct Calcium Binding Characteristics. *PLoS ONE* **2016**, *11*, e0161414. [CrossRef]
28. Kazakov, A.S.; Sofin, A.D.; Avkhacheva, N.V.; Denesyuk, A.I.; Deryusheva, E.I.; Rastrygina, V.A.; Sokolov, A.S.; Permyakova, M.E.; Litus, E.A.; Uversky, V.N.; et al. Interferon Beta Activity Is Modulated via Binding of Specific S100 Proteins. *Int. J. Mol. Sci.* **2020**, *21*, 9473. [CrossRef]

29. Makarov, A.A.; Tsvetkov, P.O.; Villard, C.; Esquieu, D.; Pourroy, B.; Fahy, J.; Braguer, D.; Peyrot, V.; Lafitte, D. Vinflunine, a novel microtubule inhibitor, suppresses calmodulin interaction with the microtubule-associated protein STOP. *Biochemistry* **2007**, *46*, 14899–14906. [CrossRef]

30. Petrushanko, I.Y.; Lobachev, V.M.; Kononikhin, A.S.; Makarov, A.A.; Devred, F.; Kovacic, H.; Kubatiev, A.A.; Tsvetkov, P.O. Oxidation of capital ES, Cyrillicsmall a, Cyrillic²⁺-Binding Domain of NADPH Oxidase 5 (NOX5): Toward Understanding the Mechanism of Inactivation of NOX5 by ROS. *PLoS ONE* **2016**, *11*, e0158726. [CrossRef]

31. Garnier, C.; Devred, F.; Byrne, D.; Puppo, R.; Roman, A.Y.; Malesinski, S.; Golovin, A.V.; Lebrun, R.; Ninkina, N.N.; Tsvetkov, P.O. Zinc binding to RNA recognition motif of TDP-43 induces the formation of amyloid-like aggregates. *Sci. Rep.* **2017**, *7*, 6812. [CrossRef]

32. Vologzhannikova, A.A.; Shevelyova, M.P.; Kazakov, A.S.; Sokolov, A.S.; Borisova, N.I.; Permyakov, E.A.; Kircheva, N.; Nikolova, V.; Dudev, T.; Permyakov, S.E. Strontium Binding to alpha-Parvalbumin, a Canonical Calcium-Binding Protein of the “EF-Hand” Family. *Biomolecules* **2021**, *11*, 1158. [CrossRef]

33. Sreerama, N.; Venyaminov, S.Y.; Woody, R.W. Estimation of protein secondary structure from circular dichroism spectra: Inclusion of denatured proteins with native proteins in the analysis. *Anal. Biochem.* **2000**, *287*, 243–251. [CrossRef]

34. Golovin, A.V.; Devred, F.; Yatoui, D.; Roman, A.Y.; Zalevsky, A.O.; Puppo, R.; Lebrun, R.; Guerlesquin, F.; Tsvetkov, P.O. Zinc Binds to RRM2 Peptide of TDP-43. *Int. J. Mol. Sci.* **2020**, *21*, 9080. [CrossRef] [PubMed]

35. Melville, Z.; Aligholizadeh, E.; McKnight, L.E.; Weber, D.J.; Pozharski, E.; Weber, D.J. X-ray crystal structure of human calcium-bound S100A1. *Acta Crystallogr. F Struct. Biol. Commun.* **2017**, *73*, 215–221. [CrossRef]

36. Nowakowski, M.; Ruszczynska-Bartnik, K.; Budzinska, M.; Jaremko, L.; Jaremko, M.; Zdanowski, K.; Bierzynski, A.; Ejchart, A. Impact of calcium binding and thionylation of S100A1 protein on its nuclear magnetic resonance-derived structure and backbone dynamics. *Biochemistry* **2013**, *52*, 1149–1159. [CrossRef]

37. Deng, L.; Sun, N.; Kitova, E.N.; Klassen, J.S. Direct quantification of protein-metal ion affinities by electrospray ionization mass spectrometry. *Anal. Chem.* **2010**, *82*, 2170–2174. [CrossRef]

38. Potier, N.; Rogniaux, H.; Chevreux, G.; Van Dorsselaer, A. Ligand-metal ion binding to proteins: Investigation by ESI mass spectrometry. *Methods Enzymol.* **2005**, *402*, 361–389. [CrossRef] [PubMed]

39. Zhukov, I.; Ejchart, A.; Bierzynski, A. Structural and motional changes induced in apo-S100A1 protein by the disulfide formation between its Cys 85 residue and beta-mercaptoethanol. *Biochemistry* **2008**, *47*, 640–650. [CrossRef]

40. Goch, G.; Vdovenko, S.; Kozlowska, H.; Bierzynski, A. Affinity of S100A1 protein for calcium increases dramatically upon glutathionylation. *FEBS J.* **2005**, *272*, 2557–2565. [CrossRef] [PubMed]

41. Wheeler, L.C.; Donor, M.T.; Prell, J.S.; Harms, M.J. Multiple Evolutionary Origins of Ubiquitous Cu²⁺ and Zn²⁺ Binding in the S100 Protein Family. *PLoS ONE* **2016**, *11*, e0164740. [CrossRef]

42. Laitaoja, M.; Valjakka, J.; Janis, J. Zinc coordination spheres in protein structures. *Inorg. Chem.* **2013**, *52*, 10983–10991. [CrossRef] [PubMed]

43. Peschke, M.; Blades, A.T.; Kebarle, P. Binding Energies for Doubly-Charged Ions M²⁺ = Mg²⁺, Ca²⁺ and Zn²⁺ with the Ligands L = H₂O, Acetone and N-methylacetamide in Complexes M for n = 1 to 7 from Gas Phase Equilibria Determinations and Theoretical Calculations. *J. Am. Chem. Soc.* **2000**, *122*, 10440–10449. [CrossRef]

44. Dudev, T.; Lim, C. Monodentate versus Bidentate Carboxylate Binding in Magnesium and Calcium Proteins: What Are the Basic Principles? *J. Phys. Chem. B* **2004**, *108*, 4546–4557. [CrossRef]

45. Sitaramayya, A. Calcium-dependent activation of guanylate cyclase by S100b. *Adv. Exp. Med. Biol.* **2002**, *514*, 389–398. [CrossRef]

46. Egelman, D.M.; Montague, P.R. Calcium dynamics in the extracellular space of mammalian neural tissue. *Biophys. J.* **1999**, *76*, 1856–1867. [CrossRef]

47. Assaf, S.Y.; Chung, S.H. Release of endogenous Zn²⁺ from brain tissue during activity. *Nature* **1984**, *308*, 734–736. [CrossRef]

48. Lyubartseva, G.; Lovell, M.A. A potential role for zinc alterations in the pathogenesis of Alzheimer's disease. *Biofactors* **2012**, *38*, 98–106. [CrossRef]

49. Haase, H.; Rink, L. Functional significance of zinc-related signaling pathways in immune cells. *Annu. Rev. Nutr.* **2009**, *29*, 133–152. [CrossRef]

50. Brett, W.; Mandinova, A.; Remppis, A.; Sauder, U.; Ruter, F.; Heizmann, C.W.; Aebi, U.; Zerkowski, H.R. Translocation of S100A1(1) calcium binding protein during heart surgery. *Biochem. Biophys. Res. Commun.* **2001**, *284*, 698–703. [CrossRef]

51. Huttunen, H.J.; Kuja-Panula, J.; Sorci, G.; Agneletti, A.L.; Donato, R.; Rauvala, H. Coregulation of neurite outgrowth and cell survival by amphoterin and S100 proteins through receptor for advanced glycation end products (RAGE) activation. *J. Biol. Chem.* **2000**, *275*, 40096–40105. [CrossRef] [PubMed]

52. Zernii, E.Y.; Grigoriev, I.I.; Nazipova, A.A.; Scholten, A.; Kolpakova, T.V.; Zinchenko, D.V.; Kazakov, A.S.; Senin, I.I.; Permyakov, S.E.; Dell'Orco, D.; et al. Regulatory function of the C-terminal segment of guanylate cyclase-activating protein 2. *Biochim. Biophys. Acta* **2015**, *1854*, 1325–1337. [CrossRef] [PubMed]

53. Landar, A.; Rustandi, R.R.; Weber, D.J.; Zimmer, D.B. S100A1 utilizes different mechanisms for interacting with calcium-dependent and calcium-independent target proteins. *Biochemistry* **1998**, *37*, 17429–17438. [CrossRef] [PubMed]
54. Zimmer, D.B.; Cornwall, E.H.; Reynolds, P.D.; Donald, C.M. S100A1 regulates neurite organization, tubulin levels, and proliferation in PC12 cells. *J. Biol. Chem.* **1998**, *273*, 4705–4711. [CrossRef]

55. Heierhorst, J.; Mann, R.J.; Kemp, B.E. Interaction of the recombinant S100A1 protein with twitchin kinase, and comparison with other Ca²⁺-binding proteins. *Eur. J. Biochem.* **1997**, *249*, 127–133. [CrossRef] [PubMed]

56. Roman, A.Y.; Devred, F.; Byrne, D.; La Rocca, R.; Ninkina, N.N.; Peyrot, V.; Tsvetkov, P.O. Zinc Induces Temperature-Dependent Reversible Self-Assembly of Tau. *J. Mol. Biol.* **2019**, *431*, 687–695. [CrossRef]

57. Magneson, G.R.; Puvathingal, J.M.; Ray, W.J., Jr. The concentrations of free Mg²⁺ and free Zn²⁺ in equine blood plasma. *J. Biol. Chem.* **1987**, *262*, 11140–11148. [CrossRef]

58. Hesketh, J.E. Zinc-stimulated microtubule assembly and evidence for zinc binding to tubulin. *Int. J. Biochem.* **1982**, *14*, 983–990. [CrossRef]

59. Hesketh, J.E. Zinc binding to tubulin. *Int. J. Biochem.* **1983**, *15*, 743–746. [CrossRef]

60. Deinum, J.; Baudier, J.; Briving, C.; Rosengren, L.; Wallin, M.; Gerard, D.; Haglid, K. The effect of S-100a and S-100b proteins and Zn²⁺ on the assembly of brain microtubule proteins in vitro. *FEBS Lett.* **1983**, *163*, 287–291. [CrossRef]

61. Cuajungco, M.P.; Faget, K.Y. Zinc takes the center stage: Its paradoxical role in Alzheimer's disease. *Brain Res. Brain Res. Rev.* **2003**, *41*, 44–56. [CrossRef]

62. Afanador, L.; Roltsch, E.A.; Holcomb, L.; Campbell, K.S.; Keeling, D.A.; Zhang, Y.; Zimmer, D.B. The Ca²⁺ sensor S100A1 modulates neuroinflammation, histopathology and Akt activity in the PSAPP Alzheimer's disease mouse model. *Cell Calcium* **2014**, *56*, 68–80. [CrossRef] [PubMed]

63. Bhatty, A.; Baig, S.; Fawwad, A.; Rubab, Z.E.; Shahid, M.A.; Waris, N. Association of Zinc Transporter-8 Autoantibody (ZnT8A) with Type 1 Diabetes Mellitus. *Cureus* **2020**, *12*, e7263. [CrossRef]

64. Lin, C.C.; Huang, Y.L. Chromium, zinc and magnesium status in type 1 diabetes. *Curr. Opin. Clin. Nutr. Metab. Care* **2015**, *18*, 588–592. [CrossRef] [PubMed]

65. Tian, T.; Li, X.; Hua, Z.; Ma, J.; Liu, Z.; Chen, H.; Cui, Z. S100A1 promotes cell proliferation and migration and is associated with lymph node metastasis in ovarian cancer. *Discov. Med.* **2017**, *23*, 235–245. [PubMed]

66. Morii, K.; Tanaka, R.; Takahashi, Y.; Kuwano, R. Cloning of cDNAs encoding human S-100 alpha and beta subunits and their differential expression in human tumor cell lines. *J. Neurosci. Res.* **1992**, *32*, 27–33. [CrossRef]

67. Kato, K.; Haimoto, H.; Ariyoshi, Y.; Horisawa, M.; Washida, H.; Kimura, S. High levels of S-100a0 (alpha alpha) protein in tumor tissues and in sera of patients with renal cell carcinoma. *JPN J. Cancer Res.* **1985**, *76*, 856–862.

68. Cocchia, D.; Michetti, F.; Donato, R. Immunochemical and immuno-cytochemical localization of S-100 antigen in normal human skin. *Nature* **1981**, *294*, 85–87. [CrossRef]

69. Sheline, C.T.; Cai, A.L.; Zhu, J.; Shi, C. Serum or target deprivation-induced neuronal death causes oxidative neuronal accumulation of Zn²⁺ and loss of NAD+. *Eur. J. Neurosci.* **2010**, *32*, 894–904. [CrossRef]

70. Makhov, P.; Golovine, K.; Uzzo, R.G.; Rothman, J.; Crispen, P.L.; Shaw, T.; Scoll, B.J.; Kolenko, V.M. Zinc chelation induces rapid depletion of the X-linked inhibitor of apoptosis and sensitizes prostate cancer cells to TRAIL-mediated apoptosis. *Cell Death Differ.* **2008**, *15*, 1745–1751. [CrossRef]

71. Kontoghiorghes, G.J.; Efstathiou, A.; Ioannou-Loucaides, S.; Kolnagou, A. Chelators controlling metal metabolism and toxicity pathways: Applications in cancer prevention, diagnosis and treatment. *Hemoglobin* **2008**, *32*, 217–227. [CrossRef]