# Functional Status of Neuronal Calcium Sensor-1 Is Modulated by Zinc Binding

Philipp O Tsvetkov, Andrei Yu Roman, Viktoriia E Baksheeva, Aliya A Nazipova, Marina P Shevelyova, Vasiliy I Vladimirov, Michelle F Buyanova, Dmitry V Zinchenko, Andrey A Zamyatnin, François Devred, et al.

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

► To cite this version:

Philipp O Tsvetkov, Andrei Yu Roman, Viktoriia E Baksheeva, Aliya A Nazipova, Marina P Shevelyova, et al. Functional Status of Neuronal Calcium Sensor-1 Is Modulated by Zinc Binding. Frontiers in Molecular Neuroscience, 2018, 11, 10.3389/fnmol.2018.00459. [https://doi.org/10.3389/fnmol.2018.00459](https://doi.org/10.3389/fnmol.2018.00459). [hal-01955226](https://HAL.inhal.science/hal-01955226)

HAL Id: hal-01955226

https://amu.hal.science/hal-01955226

Submitted on 14 Dec 2018

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_2.png" /></figure>
<figure><img src="image_3.png" /></figure>

Check for
updates

# Functional Status of Neuronal Calcium Sensor-1 Is Modulated by Zinc Binding

Philipp O. Tsvetkov<fnref n="1" />, Andrei Yu. Roman², Viktoriia E. Baksheeva³, Aliya A. Nazipova⁴,
Marina P. Shevelyova⁴, Vasiliy I. Vladimirov⁵, Michelle F. Buyanova⁶,
Dmitry V. Zinchenko⁵, Andrey A. Zamyatnin Jr.³,⁷, François Devred<fnref n="1" />,
Andrey V. Golovin³,⁶,⁷,⁸, Sergei E. Permyakov⁴ and Evgeni Yu. Zernii³,⁷*

OPEN ACCESS

**Edited by:**
Jose R. Naranjo,
Spanish National Research Council
(CSIC), Spain

**Reviewed by:**
Yogendra Sharma,
Centre for Cellular and Molecular
Biology (CSI/R), India
Alexander Dizhoor,
Salus University, United States

*Correspondence:
Evgeni Yu. Zernii
zerni@belozersky.msu.ru

**Received:** 21 July 2018
**Accepted:** 28 November 2018
**Published:** 14 December 2018

**Citation:**
Tsvetkov PO, Roman Ayu, Baksheeva VE, Nazipova AA, Shevelyova MP, Vladimirov VI, Buyanova MF, Zinchenko DV, Zamyatnin AA Jr, Devred F, Golovin AV, Permyakov SE and Zernii EYu (2018) Functional Status of Neuronal Calcium Sensor-1 Is Modulated by Zinc Binding. *Front. Mol. Neurosci.* 11:459. doi: 10.3389/fnmol.2018.00459

Neuronal calcium sensor-1 (NCS-1) protein is abundantly expressed in the central nervous system and retinal neurons, where it regulates many vital processes such as synaptic transmission. It coordinates three calcium ions by EF-hands 2-4, thereby transducing $Ca^{2+}$ signals to a wide range of protein targets, including G protein-coupled receptors and their kinases. Here, we demonstrate that NCS-1 also has $Zn^{2+}$-binding sites, which affect its structural and functional properties upon filling. Fluorescence and circular dichroism experiments reveal the impact of $Zn^{2+}$ binding on NCS-1 secondary and tertiary structure. According to atomic absorption spectroscopy and isothermal titration calorimetry studies, apo-NCS-1 has two high-affinity ($4 \times 10^6 M^{-1}$) and one low-affinity ($2 \times 10^5 M^{-1}$) $Zn^{2+}$-binding sites, whereas $Mg^{2+}$-loaded and $Ca^{2+}$-loaded forms (which dominate under physiological conditions) bind two zinc ions with submicromolar affinity. Metal competition analysis and circular dichroism studies suggest that $Zn^{2+}$-binding sites of apo- and $Mg^{2+}$-loaded NCS-1 overlap with functional EF-hands of the protein. Consistently, high $Zn^{2+}$ concentrations displace $Mg^{2+}$ from the EF-hands and decrease the stoichiometry of $Ca^{2+}$ binding. Meanwhile, one of the EF-hands of $Zn^{2+}$-saturated NCS-1 exhibits a 14-fold higher calcium affinity, which increases the overall calcium sensitivity of the protein. Based on QM/MM molecular dynamics simulations, $Zn^{2+}$ binding to $Ca^{2+}$-loaded NCS-1 could occur at EF-hands 2 and 4. The high-affinity zinc binding increases the thermal stability of $Ca^{2+}$-free NCS-1 and favours the interaction of its $Ca^{2+}$-loaded form with target proteins, such as dopamine receptor D2R and GRK1. In contrast, low-affinity zinc binding

<footnote n="1">¹ Aix-Marseille University, CNRS, INP, Institute of Neurophysiopathology, Faculty of Pharmacy, Marseille, France, ² Institute of Physiologically Active Compounds (RAS), Chernogolovka, Russia, ³ Belozersky Institute of Physico-Chemical Biology, Lomonosov Moscow State University, Moscow, Russia, ⁴ Institute for Biological Instrumentation of the Russian Academy of Sciences, Pushchino, Russia, ⁵ Shemyakin and Ovchinnikov Institute of Bioorganic Chemistry, Russian Academy of Sciences, Pushchino, Russia, ⁶ Faculty of Bioengineering and Bioinformatics, Lomonosov Moscow State University, Moscow, Russia, ⁷ Institute of Molecular Medicine, Sechenov First Moscow State Medical University, Moscow, Russia, ⁸ Faculty of Computer Science, Higher School of Economics, Moscow, Russia</footnote>
promotes NCS-1 aggregation accompanied by the formation of twisted rope-like structures. Altogether, our findings suggest a complex interplay between magnesium, calcium and zinc binding to NCS-1, leading to the appearance of multiple conformations of the protein, in turn modulating its functional status.

**Keywords: neuronal calcium sensor-1, zinc, calcium, magnesium, EF-hand motif, dopamine receptor D2R, GRK1, protein aggregation**

INTRODUCTION

Divalent metal ions play a vital role in the vast majority of cellular processes. Among them, two alkaline earth metals, magnesium and calcium, as well as the transition metal zinc, are particularly important, since they are the most abundant ones in the human body. These ions are significantly different in intracellular levels and play different physiological roles. While intracellular concentration of free magnesium is high (0.5–2 mM), the free concentration of calcium and zinc is much lower (40–100 nM and 100 pM, respectively) (Romani and Scarpa, 1992; Sabatini et al., 2002; Krezel and Maret, 2006). In addition, free calcium concentration drastically increases in signaling waves, reaching > 100 µM in magnitude (Augustine et al., 2003). Consistently, magnesium binds to proteins with much lower affinity (equilibrium association constant, $K_a$, below $10^5 M^{-1}$) than calcium ($K_a$ > $10^5 M^{-1}$) and zinc ($K_a$ > $10^7 M^{-1}$) (for review, see (Dudev and Lim, 2003)). Magnesium binding to proteins is less specific compared to $Ca^{2+}$ and $Zn^{2+}$ binding; and, although the common structural motif for $Mg^{2+}$ binding has not been described yet, magnesium is able to occupy calcium and zinc-binding sites. Yet, at high concentrations, zinc can compete with other metals or bind to nonspecific sites in proteins, which may interfere with their structural integrity and normal function, thus contributing to pathology (Choi and Koh, 1998; Barwinska-Sendra and Waldron, 2017).

The great majority of calcium-binding proteins (CaBPs) contain the same helix-loop-helix $Ca^{2+}$-binding motif (referred to as an “EF-hand”), wherein calcium ion is coordinated by six oxygen atoms. Generally, CaBPs comprise several EF-hands, which bind calcium with micromolar dissociation constant. Commonly, this interaction induces a significant rearrangement in the CaBP structure, inducing the exposure of its hydrophobic surface, which is responsible for interaction with target(s), resulting in their activation/deactivation. All CaBPs have specific tissue, cellular and subcellular distribution profiles (for review, see (Yanez et al., 2012)). Along with other tissues, CaBPs are abundantly expressed in the nervous system where they play an essential role in decoding calcium signals and regulating many processes crucial for the viability and functioning of neurons. The calcium signals, emanating from photoreceptor neurons in response to light stimuli, regulate rhodopsin desensitization and cGMP synthesis through interaction with recoverin and guanylate cyclase activating proteins (GCAPs), two members of the neuronal calcium sensor (NCS) family [reviewed in (Burgoyne, 2007; Burgoyne and Haynes, 2012; Koch and Dell'Orco, 2015)]. Neuronal calcium sensor-1 (NCS-1) is another NCS protein, which was found in photoreceptors and

other retinal neurons (De Raad et al., 1995; Baksheeva et al., 2015). In contrast to recoverin and GCAPs, NCS-1 is expressed throughout the nervous system. Its N-terminus contains a myristoyl group, which participates in the interaction of the protein with cellular membranes. According to *in vitro* and *in vivo* data, NCS-1 regulates more than 20 target proteins, including G protein-coupled receptors and their kinases (GRKs). As such, NCS-1 participates in neuronal growth and survival, reception, neurotransmission, synaptic plasticity and other cellular mechanisms [for review, see (Burgoyne and Haynes, 2012)]. NCS-1 contains four EF-hand motifs, but only three of them (EF2, EF3 and EF4) are able to bind calcium with nanomolar to micromolar dissociation constants (Jeromin et al., 2004). In the absence of calcium, two EF-hands (EF2 and EF3) could be occupied by magnesium ions (Warren et al., 2007; Aravind et al., 2008).

Zinc is long known to be a key element in neuronal growth and activity necessary for the normal development and functioning of the brain (Frederickson et al., 2005). Zinc deficiency results, for instance, in lowered glutamate receptor expression and decreased cognitive and motor performance in children (Penland et al., 1997; Gardner et al., 2005). It is also critical for the development, viability and specific function of the retinal neurons (Ugarte and Osborne, 2014). The levels of $Zn^{2+}$ in the retina and retinal pigment epithelium are decreased in the elderly, which may contribute to the pathogenesis of age-related macular degeneration (Wills et al., 2008; Lyubartseva and Lovell, 2012). Expression of retinal proteins involved in zinc homeostasis also becomes downregulated with age (Leung et al., 2012). In general, zinc serves to maintain the structure and function of hundreds of proteins, including enzymes of all known classes, transcription factors, receptors and signaling proteins. While, in many cases, zinc is tightly bound to proteins, it could also interact with their $Zn^{2+}$-binding sites transiently in order to conduct biochemical stimuli. The most abundant $Zn^{2+}$-binding motif in proteins is the Z-finger, which chelates zinc ion with nanomolar to picomolar affinity [for review, see (Maret and Li, 2009)]. Nevertheless, zinc ions can bind to EF-hands [as reported for calmodulin (Warren et al., 2007)] or between two EF-hand motifs, or even in-between two protein subunits in dimer, as observed in some members of the S100 family (Tsvetkov et al., 2010; Moroz et al., 2011). Considering that some members of the NCS family, such as recoverin (Permyakov et al., 2003), are able to bind zinc ions, we hypothesized that NCS-1 is also sensitive to $Zn^{2+}$.

In this study, we demonstrate that apo-NCS-1 has two high-affinity zinc-specific sites and one low-affinity zinc-binding site. Zinc binding to NCS-1 reduces stoichiometry and increases
the affinity of Ca²⁺ binding to the protein. In contrast, physiologically relevant Mg²⁺- and Ca²⁺-loaded NCS-1 forms only bind two zinc ions with high affinity, which stabilizes their structure and may be required for maintaining the functional status of these forms. In addition, our findings suggest that the elevated concentration of free zinc, characteristic of some neurodegenerative and neuro-ophthalmological disorders, may lead to the formation of unstable, prone-to-aggregation pathological NCS-1 forms.

# MATERIALS AND METHODS

## Purification of Proteins and Membranes

NCS-1 was obtained according to the protocol previously developed for recoverin with some modifications. To obtain recombinant myristoylated protein, NCS-1 gene was co-expressed in *Escherichia coli* strain BL21(DE3) Codon Plus RP with N-myristoyl transferase 1 from *Saccharomyces cerevisiae* (4 h, 37°C) in the presence of 200 µg/ml myristic acid. The cells were harvested by centrifugation and lysed by freezing/thawing in extraction buffer (50 mM Tris-HCl (pH 8.0), 100 mM NaCl, 1 mM EDTA, 1 mM PMSF, 1 mM DTT) and subsequent incubation in the presence of 50 µg/ml of egg white lysozyme in the same buffer for 30 min. The lysate was clarified by centrifugation, loaded onto Phenyl Sepharose column (GE Lifesciences) equilibrated with 20 mM Tris-HCl buffer (pH 8.0), 2 mM CaCl₂, 1 mM DTT and NCS-1 was eluted using the same buffer containing 2 mM EGTA. The obtained protein was loaded on HiTrap Q FF anion exchange column (GE Lifesciences) equilibrated with 20 mM Tris-HCl buffer (pH 8.0), 1 mM DTT and eluted by linear gradient of 0-1 M NaCl in the same buffer. NCS-1 (> 90% purity) was present in the fractions eluted at 380-500 mM NaCl. The obtained protein was dialysed overnight against 20 mM Tris (pH 8.0), 1 mM DTT and stored at -70°C. Alternatively, to remove residual calcium NCS-1 sample obtained immediately after anion exchange chromatography was subjected to dialysis against 50 mM Tris-HCl (pH 8.0), 5 mM EDTA (3 h), followed by dialysis against deionized water (3 h), and dialysis against 20 mM Tris-HCl (pH 8.0), 1 mM DTT (Blachford et al., 2009). The degree of NCS-1 myristoylation was determined by analytical HPLC using a reversed-phase column [Phenomenex Luna C18(2)] and was more than 97%. NCS-1 concentration was measured with bicinchoninic acid assay kit (Thermo Fisher Scientific) or spectrophotometrically using previously determined molar extinction coefficient at 280 nm of 21,430 M⁻¹ (Kazakov et al., 2017).

N-terminal fragment of GRK1 (M1-G183) was obtained as GST-fusion construct (GST-N-GRK1) following the previously developed procedure (Komolov et al., 2009).

Dopamine receptor D2 (D2R) peptide (N430-R443) was produced using Fmoc/But solid-phase peptide synthesis.

Polyclonal (monospecific) antibodies against NCS-1 were generated by rabbit immunization and purified from immune serum on a column with immobilized antigen according to the previously published procedure (Zernii et al., 2003).

Photoreceptor membranes were prepared from frozen bovine retinas following the standard protocol with some modifications described in (Grigoriev et al., 2012).

## Equilibrium Dialysis Experiments

Ca²⁺/Zn²⁺ binding to NCS-1 was studied by equilibrium dialysis method using 96-well micro-equilibrium dialysis system (HTDialysis, LLC) (Banker et al., 2003; Waters et al., 2008). Each well (500 µL) of the teflon block was separated by dialysis membrane (regenerated cellulose, 3.5 kDa MWCO). One half of each well was filled with 180 µL of (4.4–6.2) µM solution of NCS-1 in a buffer (10 mM Hepes-KOH, pH 7.6), whereas the other half contained 180 µL of the same buffer with (2–150) µM Zn(NO₃)₂ or (2–50) µM Ca(NO₃)₂ without NCS-1. The wells were tightly sealed and equilibrated by continuous shaking (130 rpm) of the block at (25.0 ± 0.5)°C for 17–20 h. Total concentrations of Ca²⁺/Zn²⁺ in the equilibrated solutions were measured by electrothermal atomization atomic absorption spectrometer iCE 3000 (Thermo Scientific), using argon as an inert gas. Zinc content was evaluated using the absorption bands at 213.9 or 307.6 nm and deuterium background correction. For calcium content estimates, the band at 422.7 nm and Zeeman background correction were used. The analytical signal was calibrated using AAS standard solutions for Ca²⁺ (Sigma-Aldrich #69349) and Zn²⁺ (Sigma-Aldrich #18827). Concentration of Ca²⁺/Zn²⁺ bound to NCS-1 was estimated for each well as a difference between the total metal concentrations measured for both halves of the well, assuming that free Ca²⁺/Zn²⁺ concentrations do not differ between the two halves of the well.

## Analytical Gel-Filtration

Analytical gel-filtration of NCS-1 forms was carried out using fast protein liquid chromatography instrument as described for GCAP2 in (Olshevskaya et al., 1999) with modifications. NCS-1 (180 µM) was pre-incubated for 30 min at 37°C in 50 mM Tris-HCl buffer (pH 8.0), 150 mM NaCl, 1 mM DTT containing either 1 mM EGTA or 1 mM Ca²⁺ or 100 µM Zn²⁺. The obtained protein sample (200 µl) was loaded onto Superdex 200 10/300 GL column (GE Lifesciences) pre-equilibrated with the same buffer and eluted at 0.5 ml/min. Alternatively, gel-filtration was performed using high performance liquid chromatography instrument on Ultropack TSK G 2000 SW column (Pharmacia) at 1 ml/min.

## Isothermal Titration Calorimetry (ITC)

Binding of divalent ions (Zn²⁺, Ca²⁺, and Mg²⁺) to NCS-1 was analyzed by ITC using MicroCal iTC200 instrument as described previously (Tsvetkov et al., 2012, 2013). Experiments were performed at 25°C in 50 mM Tris-HCl buffer (pH 7.5) in the presence of 1 mM TCEP. Protein concentration in the calorimetric cell was 25 µM, whereas the concentration of ions in the syringe varied from 375 to 750 µM. In competition experiments the concentration of competitive ions (Zn²⁺, Ca²⁺, and Mg²⁺) in cell and syringe were 250 µM, 1 mM and 5 mM, respectively. NCS-1 was titrated by repeated injections of 2 µL aliquots of ions solution. If necessarily syringe was refilled with the same solution without cell refilling and the titration was
continued. Each resulting titration peak was integrated and plotted as a function of the NCS-1/ion molar ratio. The baseline was measured by injecting titrant into the protein-free buffer solution. Data were analyzed using Origin software and were fitted with "sequential binding," "one set of sites" or "two set of sites" models via a non-linear least squares minimization method and led to the determination of affinity constants (K<sub>a</sub>, enthalpy changes (ΔH) and stoichiometry. Thermodynamic values are an average of at least three different experiments.

Binding of NCS-1 to D2R peptide was registered using MicroCal VP-ITC instrument, according to previously developed protocol (Pandalaneni et al., 2015) with modifications described in (Vladimirov et al., 2018). Experiments were performed at 25°C in 20 mM Tris-HCl buffer (pH 8.0), 150 mM NaCl, 1 mM EDTA. Alternatively, the buffer contained 5 mM CaCl<sub>2</sub> instead of EDTA, with or without addition of 100 µM ZnCl<sub>2</sub>. Recombinant NCS-1 was dialyzed against the same buffer, and protein concentration of the stock solution was adjusted to 1 mM. Calorimetric cell contained 50 µM peptide, which was titrated by thirty 10 µl injections of NCS-1. Each injection was followed by 5 min stabilization phase. The resulting titration peaks were integrated and plotted as a function of the NCS-1/ion molar ratio. The baseline was measured by injecting the protein into the working buffer solution. Data were analyzed using Origin software and were fitted with "one set of sites" model. Thermodynamic parameters were determined as an average of at least three different experiments.

## Fluorimetry and Light Scattering (LS)

Fluorescence emission spectra of NCS-1 and bis-ANS were measured using Cary Eclipse spectrofluorimeter (Varian Inc.), equipped with a Peltier-controlled cell holder essentially as previously described (Baksheeva et al., 2015; Zernii et al., 2015). Tryptophan fluorescence of NCS-1 (14 µM) was excited at 280 nm and measured at 25°C in 10 mM Hepes-KOH, 100 mM KCl, pH 7.6 buffer under various content of metal ions: either metal-free conditions (1 mM EDTA) or in the presence of Mg<sup>2+</sup> (1 mM MgCl<sub>2</sub>), Ca<sup>2+</sup> (100 µM CaCl<sub>2</sub>) or Zn<sup>2+</sup> (100 µM ZnCl<sub>2</sub>), or their combinations. Fluorescence of bis-ANS (1 µM) complexes with NCS-1 (5 µM) in the same buffer at 20°C was excited at 385 nm. All spectra were corrected for spectral sensitivity of the instrument and fitted to log-normal curves (Burstein and Emelyanenko, 1996) using LogNormal software (IBI RAS, Pushchino, Russia). Spectrofluorimetric temperature scans were performed at the average heating rate of 0.5°C/min. The mid-transition temperatures for conversion from the native to the thermally denatured protein state were estimated from fits of the temperature dependencies of λ<sub>max</sub> by Boltzmann function using OriginPro 9.0 software (OriginLab Corporation, United States).

Alternatively, tryptophan fluorescence and thermal stability of NCS-1 in the presence of different concentrations of divalent ions were measured in 50 mM Tris-HCl, 1 mM TCEP, pH 7.5 buffer using differential scanning fluorimetry (DSF) instrument Prometheus NT.Plex (NanoTemper Technologies) equipped with LS module. NanoDSF grade capillaries were filled with 25 µM NCS-1 solution. Concentrations of Zn<sup>2+</sup>, Ca<sup>2+</sup> or Mg<sup>2+</sup>

varied from 25 to 500 µM. The capillaries were loaded into the Prometheus NT.Plex instrument and the ratio of NCS-1 fluorescence emission intensities at 330 nm (I<sub>330</sub>) and 350 nm (I<sub>350</sub>) was registered at 25°C at low detector sensitivity and excitation power of 10% (excitation wavelength of 280 nm). Then capillaries were heated from 15°C to 95-110°C at rate of 1 K/min. The unfolding mid-transition temperature (T<sub>m</sub>) was determined from first derivative of the temperature dependence of the ratio, as implemented in Prometheus NT.Plex software. The temperatures of protein aggregation (T<sub>agg</sub>) were determined from temperature dependences of the LS at 350 nm.

## Circular Dichroism (CD)

Circular dichroism measurements were carried out with a JASCO J-810 spectropolarimeter (JASCO Inc., Japan), equipped with a Peltier-controlled cell holder as described in ref. (Permyakov et al., 2012). Briefly, CD spectra of NCS-1 (8 µM) were recorded at 25°C in 10 mM Hepes-KOH, 100 mM KCl, pH 7.6 buffer, either under metal-free conditions (1 mM EDTA) or in the presence of Mg<sup>2+</sup> (1 mM MgCl<sub>2</sub>), Ca<sup>2+</sup> (100 µM CaCl<sub>2</sub>) or Zn<sup>2+</sup> (100 µM ZnCl<sub>2</sub>), or their combinations. The secondary structure contents were estimated using CDPro software package (Sreerama and Woody, 2000).

## Modeling and QM/MM Molecular Dynamics

To predict zinc binding sites in NCS-1, the major parameters for Zn<sup>2+</sup> coordination (distance and angle between the cation, coordinator and one of following atoms) were analyzed in 6327 structures available in PDB. The possible range of coordinators was defined as a list of the following types of atoms: SG, ND1, NE2, OD1, OD2, OE1, OE2, OG, OG1, OH and backbone oxygen (O). The maximal distance from the cation to chelator was limited at 3Å (Laitaoja et al., 2013). Based on these data, a search for possible Zn<sup>2+</sup>-binding sites was performed in X-ray structure of NCS-1 (PDB 5AEQ, Pandalaneni et al., 2015) starting from identification of tightly interconnected (can be defined as cliques in undirected graph) zinc coordinators, namely at least 3 atoms at distances less than 6Å. For every found coordinator, all possible positions of zinc were predicted yielding local density areas. The positions that fall within VdW radii of neighboring atoms were subtracted. The resulting putative Zn<sup>2+</sup>-binding areas were ranged according to maximum density, which was visualized as volumetric data in PyMol (DeLano, 2002).

Putative coordination of zinc in the identified areas of Ca<sup>2+</sup>-occupied EF-hands was assessed by evaluating cation coordination stability using QM/MM molecular dynamics simulations. Two best score metal binding sites at distance of more than 3Å from each other were selected in Ca<sup>2+</sup>-binding loops of all three EF-hands. The best site was loaded with Ca<sup>2+</sup>, whereas the second site was loaded with Zn<sup>2+</sup>. The resulting three systems were filled with TIP3P water with 0.1 M NaCl and the net charge was neutralized with additional ions. The solvated protein, Ca<sup>2+</sup> and Zn<sup>2+</sup> were positionally restrained and water together with sodium and chloride ions were equilibrated with molecular dynamics simulation for 100 ps. On the next step, each
simulation system was divided in molecular mechanics (MM) and quantum mechanics (QM) subsystems. MM subsystem was described with parameters from the parm99sb-ildn force field with corrections (Lindorff-Larsen et al., 2010). The QM subsystem was described utilizing DFTB approach (Grimme et al., 2010; Gaus et al., 2012) and defined as any atom including waters at distance less than 5Å from Zn²⁺ or Ca²⁺. The coupling between MM and QM subsystems was performed using ONIOM approach (Dapprich et al., 1999). Linking atoms were introduced in single C-C bonds to preserve unsaturated structures in QM systems. The accordingly prepared systems were subjected to the QM/MM molecular dynamics simulation in NVT ensemble. The time step was set to 0.2 fs. Temperature coupling was performed with Velocity Rescale scheme (Bussi et al., 2007) allowing observing behavior of the systems at 300K. The total length of each simulation was set to 30 ps. All simulations were performed with GROMACS/DFTB package (Abraham et al., 2015; Kubař et al., 2015).

## Equilibrium Centrifugation Assay

The binding of NCS-1 to urea-washed photoreceptor membranes was performed according to the previously described procedure (Weiergräber et al., 2006; Senin et al., 2011) with some modifications. Briefly, 25 µM NCS-1 in buffer containing 20 mM Tris (pH 8.0), 150 mM NaCl and saturating concentration of MgCl₂ (20 mM), was mixed with the membranes in the absence and in the presence of 1 mM CaCl₂ and 0–100 µM ZnCl₂, agitated in a thermostatic shaker for 15 min (37°C, 1000 rpm) and centrifuged (24000 × g, 15 min). The pellet was dissolved in SDS-PAGE buffer and the rate of NCS-1 binding to membranes was measured by densitometric analysis of bands in polyacrylamide gel, using GelAnalyzer software<fnref n="1" />.

## Pull-Down Assay

Interaction of Zn²⁺-bound NCS-1 with GST-tagged N-terminal fragment of GRK1 (M1-G183) was monitored using analytical affinity chromatography (pull-down assay) (Zernii et al., 2011). Briefly, 50 µg of the fusion protein was immobilized on Glutathione Sepharose resin in 20 mM Tris (pH 8.0), 150 mM NaCl, 1 mM DTT. Next, 25 µM of NCS-1 was applied to the pellet. This suspension was incubated in a thermostatic shaker (1000 rpm) for 1 h at 4°C in the presence of 1 mM CaCl₂ and 0–100 µM ZnCl₂. After each incubation step non-bound protein was removed by washing the resin with the working buffer containing 0.05% Tween 20. Bound NCS-1 was eluted by SDS-PAGE sample buffer and analyzed by western blotting.

## Precipitation Assay

Precipitation of NCS-1 was monitored in the mixture containing 25 µM protein, 20 mM Tris (pH 8.0), 150 mM NaCl, 1 mM DTT, and 0–500 µM ZnCl₂ with or without addition of 1 mM CaCl₂. NCS-1 was incubated for 30 min at 37°C with mild agitation, then precipitated protein was collected by centrifugation (24000 × g, 15 min) and the pellet was dissolved in SDS-PAGE sample buffer

and analyzed by SDS-PAGE. The ratio of precipitated NCS-1 was estimated by densitometric analysis.

## Transmission Electron Microscopy (TEM)

Four microliters of NCS-1 samples (10 µM) obtained in the presence of 1 mM ZnCl₂ were placed on carbon-coated copper grids (300 mesh) during 1 min. After having been blotted, grids were washed with distilled water, blotted again, negatively stained for 30 s with 2% (wt/vol) uranyl acetate. The grids were then dried and observed with a JEOL 2200FS transmission electron microscope (Tokyo, Japan) operating at 200 kV. Images were recorded using a 4k × 4k slow-scan CCD camera (Gatan, Inc., Pleasanton, United States).

## RESULTS

### Stoichiometry of Zinc Binding to NCS-1

To test the hypothesis of the zinc interaction with NCS-1, we directly assessed the amount of zinc ions that can be bound per protein molecule using a micro-equilibrium dialysis system. To this end, a sample of recombinant myristoylated NCS-1 was prepared by Ca²⁺-dependent hydrophobic and ion exchange chromatographies and subjected, at 25°C, to 20-h dialysis (MWCO of 3.5 kDa) against buffer containing different Zn²⁺ concentrations. Zinc concentrations on both sides of the dialysis membrane were then measured by electrothermal atomization atomic absorption spectroscopy (AAS). The approximation of the resulting experimental data using the Hill equation (Figure 1) revealed half-maximal binding at 4.7 µM [Zn²⁺]ᵣfree. Meanwhile, maximal stoichiometry of the zinc binding reached 1.5. The fractional stoichiometry may have been due to either the manifestation of an intermolecular zinc-binding site (which suggests NCS-1 multimerization) or the inaccessibility to zinc ions for some fraction of the protein molecules. Analytical gel filtration experiments did not reveal NCS-1 multimers in the presence of Zn²⁺ (data not shown). Moreover, under these solution conditions, NCS-1 exhibited an even lower Stokes radius than in the presence of Ca²⁺ or EGTA (which also confirms zinc-binding to NCS-1). Therefore, we supposed that some fraction of the zinc-binding sites of NCS-1 remain shielded from zinc. Since AAS analysis of the NCS-1 sample revealed a calcium-to-protein molar ratio of 0.42, we suggest that at least one of zinc-binding sites of NCS-1 overlapped with its active EF-hands. Therefore, we further decalcified NCS-1 samples by stepwise dialysis against EDTA, deionized water and a reaction buffer, as described earlier for the non-myristoylated protein (Blachford et al., 2009). As evidenced by AAS, this procedure decreased the fraction of residual calcium in NCS-1 sample down to 0.17, which means that the protein remains Ca²⁺-bound only by 5.6% of the saturation. Although the resulting protein sample (apo-NCS-1) bound Ca²⁺ with a stoichiometry of 3 (data not shown), we failed to detect zinc binding to apo-NCS-1 by AAS, as its long-term incubation during equilibrium dialysis in the presence of Zn²⁺ was accompanied by NCS-1 aggregation and accumulation on the dialysis membrane. Given this observation,

<footnote n="1">¹http://www.gelanalyzer.com/</footnote>
<figure><img src="image_4.png" /><figcaption>**FIGURE 1** | Zn²⁺ binding to NCS-1 according to equilibrium dialysis experiments. The NCS-1 (5 µM) sample was pre-equilibrated with (2–150) µM Zn²⁺ by equilibrium dialysis (3.5 kDa MWCO) at 25°C. The concentrations of NCS-1-bound zinc ([Zn²⁺]ᵣbound) and free zinc ([Zn²⁺]ᵣfree) in the resulting solutions were measured by electrothermal atomization AAS, using the absorption bands at 213.9 nm or 307.6 nm. The solid curve approximates the experimental data by Hill equation.</figcaption></figure>

further characterization of the cation-binding properties of NCS-1 was performed by ITC. Yet, the AAS data represent direct evidence of zinc binding to NCS-1, revealing its dependence on calcium binding.

## Thermodynamics of Calcium and Magnesium Binding to NCS-1

We employed ITC to determine Ca²⁺/Mg²⁺-binding parameters of the decalcified myristoylated NCS-1 (apo-NCS-1) sample, given that previous data on these properties were contradictory. Apo-NCS-1 (25 µM) was titrated by CaCl₂ in 50 mM Tris-HCl pH 7.5 buffer in the presence of 1 mM TCEP (Figure 2A, top panel). The use of “one set of sites” or “two sets of sites” models did not allow for a correct fit of the titration curve. Meanwhile, the experimental data were well fitted using the “sequential binding” model assuming three calcium sites (Figure 2A, bottom panel): the respective equilibrium association constants are 4.3 × 10⁶ M⁻¹, 2.0 × 10⁵ M⁻¹, and 3.5 × 10⁶ M⁻¹ (Table 1). It should be noted that calcium binding to the two high-affinity sites was enthalpy-driven (Table 1), while calcium binding to the lower affinity site had an unfavorable enthalpy of 1.4 kcal/mol, indicating significant rearrangement of hydrophobic amino acids upon calcium binding to this site.

The ITC data on magnesium binding to apo-NCS-1 were well described by the “one set of sites” model (Figure 2C), revealing 2.7 Mg²⁺ bound per protein molecule with an equilibrium association constant of 5.2 × 10⁵ M⁻¹ (Table 1). Considering that apo-NCS-1 contained small fraction of residual calcium (0.17, see previous section), one can suppose that actual stoichiometry

of Mg²⁺ binding tends to 3. Indeed, magnesium ions compete with calcium for the same sites as no Mg²⁺ binding was observed for NCS-1 saturated with Ca²⁺ (1 mM). Consistent with this suggestion, in the presence of 5 mM Mg²⁺ NCS-1 exhibited decreased affinity with calcium (Table 1). At the same time, the number of Ca²⁺ bound per NCS-1 molecule in excess of magnesium reached 2.7, indicating that calcium completely replaced the bound magnesium ions.

## Thermodynamics of Zinc Binding to NCS-1

The apo-form of myristoylated NCS-1 (25 µM) in 50 mM Tris-HCl pH 7.5 buffer, in the presence of 1 mM TCEP, was titrated by ZnCl₂ using ITC (Figure 2B). The fitting of the resulting titration curve using the “two sets of sites” model indicated that NCS-1 bound Zn²⁺ in two equal high-affinity sites and one low-affinity site. The corresponding equilibrium association constants were calculated as 9.2 × 10⁶ M⁻¹ and 2.3 × 10⁵ M⁻¹, respectively (Table 1). The saturation of NCS-1 with Ca²⁺ or Mg²⁺ abolished the low-affinity Zn²⁺ binding, while zinc affinity of the other two sites was almost unaffected. Meanwhile, the enthalpy changes (ΔH) accompanying zinc interaction with these forms differed: ΔH was negative for Zn²⁺ binding to apo- and Mg²⁺-loaded NCS-1, and positive for Zn²⁺ binding to Ca²⁺-loaded NCS-1, thereby reflecting significant conformational differences between Mg²⁺- and Ca²⁺-loaded NCS-1 states. Finally, Zn²⁺-loaded NCS-1 was unable to bind magnesium, but coordinated 1.5 Ca²⁺ ions per protein molecule with increased affinity of the first site ($K_a = 5.9 \times 10^7 M^{-1}$).
<figure><img src="image_5.png" /><figcaption>**FIGURE 2** | Thermodynamics of calcium, zinc and magnesium binding to NCS-1. Typical ITC curves (upper panels) and binding isotherms (lower panels) representing titration of NCS-1 (25 µM) by different cations. (A) Binding of Ca²⁺. (B) Binding of Zn²⁺. (C) Binding of Mg²⁺. (D) Binding of Ca²⁺ in the presence of 250 µM Zn²⁺. (E) Binding of Zn²⁺ in the presence of 1 mM Ca²⁺. (F) Binding of Mg²⁺ in the presence of 250 µM Zn²⁺. (G) Binding of Ca²⁺ in the presence of 5 mM Mg²⁺. (H) Binding of Zn²⁺ in the presence of 5 mM Mg²⁺. (I) Binding of Mg²⁺ in the presence of 1 mM Ca²⁺. Best fits are shown as solid curves (see Table 1).</figcaption></figure>
Overall, the different modes of zinc binding to NCS-1, as revealed by our data, indicate that structural and functional consequences of this interaction depend on NCS-1 conformation. Thus, the ITC data argue for the existence of Mg²⁺-bound, Ca²⁺-bound and Zn²⁺-bound conformers of NCS-1 as well as its Zn²⁺(Mg²⁺)-bound, Zn²⁺(Ca²⁺)-bound and Ca²⁺(Zn²⁺)-bound forms, where Zn²⁺ or Ca²⁺ are bound to the background of the excess of Mg²⁺, Ca²⁺ and Zn²⁺, respectively (hereinafter the cation that is taken in excess is indicated in parentheses).

## Conformational Properties of Zinc-Bound NCS-1

The ITC experiments suggested the existence of several distinct states of NCS-1 with two different metal ions bound simultaneously. To explore structural differences between these NCS-1 states, we measured the intrinsic fluorescence spectra of NCS-1 in the presence of various combinations of the metals studied, which enabled examination of the mobility and polarity of the microenvironment of Trp30 and Trp103 residues located in N- and C-terminal domains of the protein.

The fluorescence spectra of 15 µM NCS-1 were measured either under metal-free conditions (1 mM EDTA) or in the presence of 1 mM Mg²⁺, 0.1 mM Ca²⁺ or 0.1 mM Zn²⁺, or their combinations. Apo-NCS-1 exhibited a characteristic tryptophan fluorescence emission spectrum with a maximum position (λmax) at 338 nm (Figures 3A,B). Mg²⁺ binding significantly increased the maximal intensity (Imax) of the fluorescence emission spectrum of NCS-1 without affecting the maximum position. Calcium binding to NCS-1 increased its Imax value and shifted its λmax to 334 nm, indicating movement of the emitting Trp residue(s) to a less polar and/or mobile environment. Zinc binding to apo-protein increased Imax without affecting λmax resembling the effect of magnesium in this respect (Figures 3A,B). Meanwhile, the presence of zinc only moderately affected the fluorescence spectra of Mg²⁺- and Ca²⁺-saturated NCS-1, indicating minor structural rearrangements near to the emitting Trp residue(s) under these experimental conditions.

Apo-NCS-1 (8 µM) represented far-UV circular dichroism (CD) spectra typical for an α-helical fold with characteristic minima at 208 nm and 222 nm (Figure 3C). Binding of all examined cations, including zinc, to NCS-1 was accompanied by a similar increase in its α-helical content and a decrease in the content of β-sheets and unordered regions (Figure 3D and Supplementary Table S1). Importantly, NCS-1 contains only two short antiparallel β-sheets, which connect the Ca²⁺-binding loops of EF1-EF2 and EF3-EF4 pairs of EF-hand motifs (Heidarsson et al., 2012; Pandalaneni et al., 2015). Therefore, the revealed similar changes in the β-structure content of NCS-1 upon binding of Ca²⁺/Mg²⁺ and Zn²⁺ suggest that zinc binds to sites overlapping with the EF-hand loops. It should be noted that single and double ion-bound NCS-1 forms also exhibited certain differences in their secondary structure (Figure 3D and Supplementary Table S1). Thus, Ca²⁺-bound and Zn²⁺(Ca²⁺)-bound forms, as well as Mg²⁺-bound and Zn²⁺(Mg²⁺)-bound forms, were identical in α-helical content, but differed in β-structure content 1.46-fold and 1.27-fold, respectively (Supplementary Table S1). These data suggest that zinc binding to Ca²⁺-saturated or Mg²⁺-saturated NCS1 does not significantly alter its overall secondary structure, but still affects its EF-hands.

To gain further insight into conformational differences between single and double ion-bound forms of NCS-1, we measured concentration dependencies of the ratio of its fluorescence intensities at 350 nm and 330 nm (I350/I330) for apo-NCS-1 or NCS-1, saturated by Ca²⁺, Mg²⁺ or Zn²⁺. In the presence of increasing concentrations of Mg²⁺, Ca²⁺ or Zn²⁺, a gradual reduction in the I350/I330 ratio of NCS-1 was observed, confirming interaction of the protein with these cations (Figure 4A). The decrease in the ratio was most pronounced for Ca²⁺ ions, suggesting that Ca²⁺-loaded NCS-1 is structurally different from its Mg²⁺-bound or Zn²⁺-bound states (Figure 4A). In the case of Mg²⁺-saturated NCS-1, low concentrations of Zn²⁺ decreased the I350/I330 ratio, while, at a more than threefold molar excess of zinc, the ratio increased toward the level of the Zn²⁺-bound form (Figure 4B), probably

<table><caption>**TABLE 1** | Thermodynamic parameters of zinc, calcium and magnesium binding to NCS-1 in 50 mM Tris-HCl buffer (pH 7.5) in the presence of 1 mM TCEP at 25°C, estimated from ITC data (see Figure 2).</caption><thead><tr><th>Ion</th><th>Competitor</th><th>K<sub>A</sub><sup>1</sup>, M<sup>-1</sup></th><th>ΔH<sup>1</sup>, kcal M<sup>-1</sup></th><th>K<sub>A</sub><sup>2</sup>, M<sup>-1</sup></th><th>ΔH<sup>2</sup>, kcal M<sup>-1</sup></th><th>K<sub>A</sub><sup>3</sup>, M<sup>-1</sup></th><th>ΔH<sup>3</sup>, kcal M<sup>-1</sup></th></tr></thead><tbody><tr><td>Ca<sup>2+</sup>**</td><td>-</td><td>4.3 × 10<sup>6</sup></td><td>-10.1</td><td>2.0 × 10<sup>5</sup></td><td>1.4</td><td>3.5 × 10<sup>6</sup></td><td>-17.8</td></tr></tbody></table><table><thead><tr><th>Ion</th><th>Competitor</th><th>N<sup>1</sup></th><th>K<sub>A</sub><sup>1</sup>, M<sup>-1</sup></th><th>ΔH<sup>1</sup>, kcal M<sup>-1</sup></th><th>N<sup>2</sup></th><th>K<sub>A</sub><sup>2</sup>, M<sup>-1</sup></th><th>ΔH<sup>2</sup>, kcal M<sup>-1</sup></th></tr></thead><tbody><tr><td rowspan="2">Ca<sup>2+</sup>***</td><td>Zn<sup>2+</sup></td><td>0.6</td><td>5.9 × 10<sup>7</sup></td><td>-3.9</td><td>0.9</td><td>3.5 × 10<sup>6</sup></td><td>-6.5</td></tr><tr><td>Mg<sup>2+</sup></td><td>1.5</td><td>2.2 × 10<sup>4</sup></td><td>-17.6</td><td>1.2</td><td>4.3 × 10<sup>6</sup></td><td>-5.7</td></tr><tr><td rowspan="3">Zn<sup>2+</sup>***</td><td>-</td><td>0.7</td><td>2.3 × 10<sup>5</sup></td><td>-11.8</td><td>2.0</td><td>9.2 × 10<sup>6</sup></td><td>-7.3</td></tr><tr><td>Ca<sup>2+</sup></td><td></td><td></td><td></td><td>1.7</td><td>2.9 × 10<sup>6</sup></td><td>1.0</td></tr><tr><td>Mg<sup>2+</sup></td><td></td><td></td><td></td><td>1.7</td><td>4.2 × 10<sup>6</sup></td><td>-6.3</td></tr><tr><td rowspan="2">Mg<sup>2+</sup>***</td><td>-</td><td></td><td></td><td></td><td>2.7</td><td>5.2 × 10<sup>5</sup></td><td>-2.1</td></tr><tr><td>Ca<sup>2+</sup></td><td>No binding</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>Zn<sup>2+</sup></td><td>No binding</td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

<footnote>*Data were fitted using "sequential binding" model. **Data were fitted using "two sets of sites" model. The concentrations of competitive ions Zn²⁺, Ca²⁺, and Mg²⁺ were 250 µM, 1 mM, and 5 mM, respectively.</footnote>
<figure><img src="image_6.png" /><figcaption>B</figcaption></figure>

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

<figure><img src="image_8.png" /><figcaption>D</figcaption></figure>

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

**FIGURE 3** | Conformational properties of NCS-1 under saturating concentrations of calcium, magnesium and zinc. (A) Representative tryptophan fluorescence spectra of 14 µM NCS-1 under metal-free conditions (1 mM EDTA) or in the presence of Ca²⁺ (100 µM CaCl₂), Mg²⁺ (1 mM MgCl₂), Zn²⁺ (100 µM ZnCl₂) or their combinations. (B) The histogram demonstrating maximum position λmax (left bar) and maximal intensity Imax (right bar) of fluorescence spectrum for the data shown in panel (A). (C) Representative CD spectra of 8 µM NCS-1 under metal-free conditions or in the presence of Ca²⁺ (100 µM CaCl₂), Mg²⁺ (1 mM MgCl₂), Zn²⁺ (100 µM ZnCl₂) or their combinations. (D) The histogram representing secondary structure contents (in %) for each NCS-1 spectrum shown in panel (C).

reflecting the replacement of magnesium by zinc (see previous section). In contrast, the binding of Zn²⁺ to Ca²⁺-NCS-1 produced a highly moderate increasing effect on the ratio (Figure 4C). Monitoring of the I₃₅₀/I₃₃₀ ratio for Zn²⁺-saturated NCS-1 in the presence of increasing calcium levels revealed signs of the Ca²⁺-bound-like conformation of the protein at lower calcium concentrations than in the case of Ca²⁺ binding to apo-NCS-1 (Figure 4D), which agreed with the increased Ca²⁺ affinity of Zn²⁺-saturated protein (see Table 1). It should be noted that differences in conformational changes induced by zinc binding to apo, Mg²⁺-saturated and Ca²⁺-saturated NCS-1 were the most striking at low zinc levels, when they were likely correlated with the stoichiometry of the metals' binding. Thus, the I₃₅₀/I₃₃₀ ratio for apo-NCS-1 decreased even in the case of a onefold excess of zinc (one Zn²⁺ bound), while the same value for Ca²⁺-NCS-1 exhibited a moderate increase only when Zn²⁺ concentration exceeded the protein concentration by three times (two Zn²⁺ bound). In contrast, the fluorescence of Mg²⁺-NCS-1 remained unchanged until reaching a 2.5-fold excess of zinc, whereas the further elevation of Zn²⁺ concentration resulted in a sequential decrease (presumably one Zn²⁺ and two Mg²⁺ bound to the protein) and an increase (presumably two Zn²⁺ and one Mg²⁺ bound to the protein) in the ratio.

Summing up, the spectral measurements reveal certain structural differences between apo, Mg²⁺-bound, Ca²⁺-bound, Zn²⁺-bound, Zn²⁺(Mg²⁺)-bound, Zn²⁺(Ca²⁺)-bound and Ca²⁺(Zn²⁺)-bound conformers of NCS-1.

## Thermal Stability of NCS-1 in the Presence of Zinc

Thermal unfolding of NCS-1 is accompanied by a red shift in its tryptophan fluorescence spectrum, implying that λmax can be used for monitoring thermal denaturation of the protein (Baksheeva et al., 2015). We compared the thermal unfolding profiles of NCS-1 (15 µM) in the presence of 1 mM EDTA, 1 mM Mg²⁺, 100 µM Ca²⁺, 100 µM Zn²⁺ or their combinations. Inspection of the experimental curves revealed that apo-NCS-1 was relatively unstable with a mid-transition temperature (Tₘ) of 40°C, whereas, in the presence of magnesium, Tₘ increased up to 70°C (Figure 5A). In the presence of Zn²⁺, or Mg²⁺ and Zn²⁺, NCS-1 exhibited similar temperature profiles without a clear transition over the experimental temperature range (Figure 5A). In both cases, the dispersion of λmax values observed at temperatures above 60°C indicated protein aggregation. Ca²⁺-saturated NCS-1 demonstrated blue-shifted emission spectra and a Tₘ value exceeding 80°C (Figure 5B). Meanwhile, zinc binding
<figure><img src="image_10.png" /></figure>

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

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

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

**FIGURE 4** | Calcium, magnesium and zinc dependences of conformational properties of NCS-1. Tryptophan fluorescence intensity at 350/330 nm (*I*<sub>350</sub>/*I*<sub>330</sub>) for apo (A), Mg<sup>2+</sup>-saturated (B), Ca<sup>2+</sup>-saturated (C), or Zn<sup>2+</sup>-saturated (D) NCS-1 (25 µM) in the presence of increasing concentrations of the other cations. Standard deviation of *I*<sub>350</sub>/*I*<sub>330</sub> values did not exceed 0.02.

to calcium-loaded NCS-1 shifted the thermal transition of the protein toward lower temperatures, thereby reflecting structural differences between Ca²⁺-bound and Zn²⁺(Ca²⁺)-bound NCS1.

More information was obtained upon monitoring NCS-1 (25 µM) denaturation by registering temperature dependences of the *I*<sub>350</sub>/*I*<sub>330</sub> ratio at different excesses of the cations. A nanoDSF instrument was used since it allows us to monitor, in parallel, the aggregation of the protein by measuring the LS of the sample at 350 nm upon heating. The binding of any of the three cations to apo-NCS-1 increased the stability of the protein, but with different efficacy and within different concentration ranges (Figure 6A). Indeed, in the presence of the fourfold excess of calcium (100 µM), the *T*<sub>m</sub> of the protein increased to > 80°C, while, in the case of the same Mg²⁺ and Zn²⁺ concentrations, the increase was moderate (48 and 42°C, respectively). Interestingly, the use of higher calcium or magnesium concentrations further improved protein stability without affecting the aggregative state, whereas zinc, at more than a fourfold excess, increased susceptibility of the protein to aggregation as indicated by LS (Figure 6C). The binding of zinc to Mg²⁺-NCS-1 had no effect on its stability until a 2.5-fold excess of Zn²⁺ was used (Figure 6B). At this point, the denaturation temperature increased to 80°C and then began to drop, apparently reflecting the formation of Zn²⁺/2Mg²⁺ NCS-1 intermediate and 2Zn²⁺/Mg²⁺ NCS-1 conformer, respectively (see above). The drop was associated with a reduction in

aggregation temperature (*T*<sub>agg</sub>), indicating increased propensity of the protein to aggregation (Figure 6D). We were technically unable to monitor the impact of low zinc concentrations (one- to fivefold excess) on Ca²⁺-saturated NCS-1 (1 mM Ca²⁺) as *T*<sub>m</sub> of the latter exceeded 90°C, which is beyond the detection limit of the method (Figure 6E). Yet, at higher levels, zinc produced a gradual destabilizing effect on NCS-1 and enhanced its susceptibility to aggregation (Figures 6E,G, see also Supplementary Figure S2A). Finally, the presence of Ca²⁺ inhibited aggregation of Zn²⁺-saturated NCS-1 and increased its thermal stability as soon as the protein bound the first calcium ion (onefold excess of Ca²⁺). However, at high calcium concentrations, *T*<sub>m</sub> did not exceed 78°C, indicating that the resulting NCS-1 conformer represents a Ca²⁺(Zn²⁺)-bound form rather than a Ca²⁺-saturated form of the protein (*T*<sub>m</sub> > 90°C) (Figures 6F,H).

Overall, at low zinc concentrations, corresponding to full saturation of Zn²⁺-binding sites in each of the NCS-1 forms, the binding of the cation slightly destabilizes Ca²⁺-loaded NCS1 and enhances the thermal stability of Ca²⁺-free protein. Meanwhile, upon elevation of Zn²⁺ levels, all NCS-1 forms become gradually destabilized and prone to aggregation. Ca²⁺-NCS-1 is the most resistant to the destabilizing effects of high zinc. Consistently, the binding of calcium to Zn²⁺-saturated NCS-1 improves its structure by forming a relatively stable Ca²⁺(Zn²⁺)-bound conformer.
<figure><img src="image_14.png" /></figure>

## Putative Zinc-Binding Sites in Calcium-Saturated NCS-1

The results of spectroscopic and thermal stability studies revealed the existence of a Zn²⁺(Ca²⁺)-bound form of NCS-1, which structurally differs from the well-recognized Ca²⁺-saturated conformer of the protein. Therefore, we next attempted to predict Zn²⁺-binding site locations in Ca²⁺-NCS-1 *in silico*, based on the available crystal structure of this form of the protein [PDB 5AEQ (Pandalaneni et al., 2015)]. Considering the averaged geometry of zinc coordination in all Zn²⁺-binding proteins presented in PDB, the density of Zn²⁺-binding probability in the NCS-1 structure was built in grid with a step of 0.1 Å (Figure 7A). It was found that areas with the required number of chelating groups for Zn²⁺ was located only in the loops of three functional EF-hand sites, namely EF2 (the highest score), EF3 and EF4. Interestingly, the size of these areas was around 4.5 Å, suggesting that they could simultaneously accommodate calcium and zinc ions. Furthermore, such a configuration would compensate for a negative charge (-2 in EF2, -1 in EF3 and EF4), which remained in EF hand loops upon binding of single Ca²⁺. In order to check this suggestion, we performed QM/MM simulations of molecular dynamics associated with Zn²⁺ binding in each Ca²⁺-occupied EF-hand motif. It was found that, in EF2, the number of coordinators around calcium ions in the presence of zinc decreased from seven to six, but most of the metal-chelating residues of the loop (Asp 73, Asn 75, Asp77, Arg79, Glu 81), as well as a water molecule, remained involved in the binding. In this case, the coordination of zinc was maintained by four oxygen atoms from Asp 73 (two atoms, from α-carbonyl and β-carboxyl groups), Asn 75 and Glu 84 (Figure 7B). In EF3, calcium lost three chelators coordinated by Asp109, Asp 111, Glu120 and a water molecule, whereas zinc possessed less favorable coordination due to three oxygen atoms from Tyr115, Asp 109 and Asp 111 (Figure 7C). As for EF4, both cations bound simultaneously in relatively optimal coordination. Thus, calcium was chelated by Glu168, Asp157, Asn159, Asp161, Lys163 (backbone) and a water molecule, whereas zinc was bound

to Met156 (backbone), Asp157, Glu168 and a water molecule (Figure 7D). Taken together, these data indirectly support our suggestion of zinc coordination in EF2, EF3 and EF4 in apo and Mg²⁺-bound NCS-1, as well as provide a rationale for the prediction of Zn²⁺-binding sites in the second and fourth EF-hands of the Ca²⁺-saturated form of the protein.

## Functional Properties of NCS-1 in the Presence of Zinc

Previous *in vitro* and *in vivo* studies revealed that NCS-1 can regulate a number of targets including membrane-associated proteins. Consistently, an important feature of NCS-1 is its Ca²⁺-induced interaction with cellular membranes via the N-terminal myristoyl group of the protein (Baksheeva et al., 2015). Thus, we next explored whether zinc binding affects the affinity of NCS-1 to photoreceptor membranes. Among the detected forms of protein, we focused on the Zn²⁺(Mg²⁺)-bound and Zn²⁺(Ca²⁺)-bound NCS-1 conformers as they might dominate under physiological conditions. According to the data from the modified equilibrium centrifugation assay, Ca²⁺-saturated NCS-1 (25 µM) at 25°C in 20 mM Tris-HCl pH 8.0 buffer, bound to urea-washed photoreceptor membranes, and the binding decreased approximately twofold in the case of Ca²⁺-free/Mg²⁺-saturated NCS-1. Meanwhile, the presence of 0-100 µM Zn²⁺ did not affect the membrane association of both NCS-1 forms (Supplementary Figure S1).

In order to further address the possible effects of zinc on the functional activity of NCS-1, we monitored the interaction of the protein, with D2R and GRK1 representing its well-established Ca²⁺-dependent targets (Pandalaneni et al., 2015). Firstly, the binding of NCS-1 to the complementary D2R peptide N430-R443 (50 µM) was monitored at 25°C in 20 mM Tris-HCl pH 8.0 buffer using ITC. Without calcium, no interaction between NCS-1 and the peptide was registered, regardless of the presence of zinc (data not shown). Meanwhile, Ca²⁺-loaded NCS-1 bound two moles of D2R peptide with a dissociation constant of 30.12 µM (Figure 8A). Remarkably, in the presence of zinc, Ca²⁺-NCS-1
<figure><img src="image_15.png" /></figure>

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

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

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

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

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

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

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

**FIGURE 6** | Calcium, magnesium and zinc dependences of thermal stability of NCS-1. Mid-transition temperatures of NCS-1 (25 µM) denaturation determined from tryptophan fluorescence at 350/330 nm (*I*<sub>350/330</sub>) for apo (A), Mg<sup>2+</sup>-saturated (B), Ca<sup>2+</sup>-saturated (E), or Zn<sup>2+</sup>-saturated (F) protein in the presence of increasing concentrations of the alternative cations. Mid-transition temperatures of NCS-1 (25 µM) aggregation determined from light scattering at 350 nm for apo (C), Mg<sup>2+</sup>-saturated (D), Ca<sup>2+</sup>-saturated (G), or Zn<sup>2+</sup>-saturated (H) protein in the presence of increasing concentrations of the other cations.
<figure><img src="image_23.png" /><figcaption>**FIGURE 7** | Prediction of zinc-binding sites in Ca²⁺-saturated NCS-1 using molecular modeling. (A) Density of putative Zn²⁺ positions in NCS-1 from geometry search visualized as volumetric data, from low (yellow) to high (black) values. (B–D) Positions of chelators for zinc and calcium in Ca²⁺-binding loops of EF2 (B), EF3 (C) and EF4 (D) according to QM/MM molecular dynamics simulations.</figcaption></figure>

interacted with D2R peptide with the same stoichiometry, but with a 3.5-fold increase in affinity (Figure 8B and Table 2). Similar observations were made upon monitoring the interaction of NCS-1 (25 µM) with N-terminal domain of GRK1 (M1-G183), fused with glutathione-S-transferase (GST-N-GRK1) at 25°C in 20 mM of Tris-HCl pH 8.0 buffer, by means of a pull-down assay. Thus, 1 mM Ca²⁺ GST-N-GRK1 was bound to NCS-1 and the binding was enhanced twofold in the presence of 25 µM of Zn²⁺ (Figure 8C). Interestingly, a further increase in zinc concentration resulted in the gradual destabilization of the NCS-1-GRK1 complex.

NCS-1 is known to interact with D2R and GRK1 via hydrophobic sites, which become available in response to

Ca²⁺ binding (Pandalaneni et al., 2015). Therefore, we next investigated effects of zinc on the accessibility of such sites in Mg²⁺-saturated and Ca²⁺-saturated NCS-1, using fluorescent dye bis-ANS. The interaction of bis-ANS with hydrophobic cavities of a protein is accompanied by increased intensity and a shift in the λmax of the fluorescence of the dye. It was found that, in the presence of zinc, bis-ANS binding to both Mg²⁺-saturated and Ca²⁺-saturated NCS-1 was moderately enhanced, suggesting increased surface hydrophobicity of these forms (Figure 8D). Thus, it is this effect that may partially account for the increased affinity of Zn²⁺(Ca²⁺)-bound NCS-1 to D2R and GRK1.

Taken together, our data demonstrate that, at low physiological concentrations, zinc cannot substitute calcium in relation to NCS-1 activation; rather, it affects the structure and stability of Ca²⁺-saturated protein, thereby improving its normal functionality.

<table><caption>**TABLE 2** | Thermodynamic parameters of binding of D2R peptide to NCS-1 in 20 mM Tris-HCl buffer (pH 8.0), 150 mM NaCl, 5 mM CaCl₂ in the presence or in the absence of 100 µM ZnCl₂.</caption><thead><tr><th></th><th>Ca²⁺</th><th>Ca²⁺+Zn²⁺</th></tr></thead><tbody><tr><td>N</td><td>0.676 ± 0.091</td><td>0.508 ± 0.102</td></tr><tr><td>K<sub>A</sub>, M⁻¹</td><td>(3.320 ± 0.568) × 10⁴</td><td>(11.800 ± 2.550) × 10⁴</td></tr><tr><td>K<sub>D</sub>, M</td><td>30.12 × 10⁻⁶</td><td>8.47 × 10⁻⁶</td></tr><tr><td>ΔH, kcal M⁻¹</td><td>-3.0 ± 0.5</td><td>-3.2 ± 0.8</td></tr><tr><td>ΔS, cal K⁻¹ M⁻¹</td><td>10.6</td><td>12.4</td></tr><tr><td>ΔG, kcal M⁻¹</td><td>-6.2</td><td>-6.9</td></tr></tbody></table>

## Abnormal Behavior of NCS-1 in the Presence of Excessive Zinc Concentrations

Although the estimated intracellular concentration of free zinc is considerably low, it is entirely possible that, under certain pathological conditions, it can abnormally increase. As such, we further analyzed behavior of different forms of NCS-1 in the presence of “pathological” amounts of zinc. According to
<figure><img src="image_24.png" /></figure>

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

**FIGURE 8** | Target-binding properties of NCS-1 in the presence of zinc. (A,B) Typical ITC curves (upper panels) and binding isotherms (lower panels) representing titration of 50 µM D2R peptide with 5–150 µM NCS-1 in the presence of 5 mM Ca²⁺ (A) or 5 mM Ca²⁺ and 100 µM Zn²⁺ (B). (C) Binding of 25 µM NCS-1 to GST-N-GRK1 at 1 mM Ca²⁺ in the presence of 0, 25, 50, or 75 µM Zn²⁺ (i.e., at [Zn²⁺]/[NCS-1] ratio of 0-3), monitored by pull-down assay. (D) Representative fluorescence spectra of bis-ANS (1.2 µM) and NCS-1 (5 µM) complexes formed in the presence of either 1 mM Mg²⁺ or 100 µM Ca²⁺ with or without addition of 100 µM Zn²⁺.

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

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

LS data, the susceptibility of the protein to aggregation, in the presence of high zinc concentrations, decreased in the following order: apo-NCS-1 > Mg²⁺-bound NCS-1 > Ca²⁺-bound NCS-1 (Supplementary Figure S2A). For apo and Mg²⁺-bound NCS-1, the decrease in the temperature of aggregation started when Zn²⁺ concentration exceeded the concentration required for full saturation of the protein by 25 µM. In contrast, Ca²⁺-bound NCS-1 can sustain up to 100 µM free Zn²⁺. At physiological temperatures, the signs of aggregation of apo and Mg²⁺-bound NCS-1 were observed at 300–500 µM free Zn²⁺ (data not shown). Meanwhile, in the presence of calcium, no NCS-1 aggregation was detected under these conditions. Since aggregation includes multimeric associations, which can produce insoluble precipitates of the protein, we also monitored the

Zn²⁺-induced precipitation of Mg²⁺-loaded and Ca²⁺-loaded NCS-1 (25 µM) at 25°C (Supplementary Figure S2B). The formation of NCS-1 precipitates was initiated at 200 and 325 µM of free zinc for Mg²⁺-loaded and Ca²⁺-loaded NCS-1 forms, respectively. To visualize the shape and arrangement of the insoluble NCS-1 conglomerates formed in the presence of zinc, we further examined the respective protein precipitates by means of TEM. As can be seen from TEM data (Supplementary Figures S2C,D), Zn²⁺-bound NCS-1 constitutes fibrilic twisted rope-like structures resembling the aggregates of another neuronal protein, TDP-43, found in the presence of zinc ions (Garnier et al., 2017).

We concluded that, at high concentrations, zinc might bind to NCS-1 non-specifically, thereby deteriorating the structure of the
protein and promoting its aggregation and precipitation, which are most prominent in the absence of calcium.

DISCUSSION

Previous *in vitro* studies reported the existence of three major forms of NCS-1 in terms of metal binding, namely, apo, Mg²⁺-bound and Ca²⁺-bound. Meanwhile, the data concerning the stoichiometry and affinity of calcium binding to the protein are contradictory. Thus, according to flow dialysis, non-myristoylated NCS-1 (nNCS-1) cooperatively binds two calcium ions with nanomolar and micromolar affinities (Cox et al., 1994). Subsequent ITC experiments also suggested the binding of two Ca²⁺, but in a non-cooperative manner and with a dissociation constant of 1.8 µM for both sites (Jeromin et al., 2004). Meanwhile, refinement of the data using NMR studies revealed that nNCS-1 actually coordinates three calcium ions in EF2-EF4 (Chandra et al., 2011; Heidarsson et al., 2012). Myristoylated NCS-1 (mNCS-1) was reported to bind three calcium ions. However, two sets of ITC studies conducted by the same authors report different modes of calcium binding, which likely depend on the preparation of the protein samples and the model applied for fitting of the ITC data (Jeromin et al., 2004; Aravind et al., 2008). Thus, in the first study, the use of a “three sequential binding sites” model revealed binding constants of a micromolar, nanomolar and submicromolar order (Jeromin et al., 2004), whereas, in the second study, the “two sets of sites” model was applied, which allowed for identifying two similar sites with submicromolar affinity and one site with nanomolar affinity (Aravind et al., 2008). The ITC data on Ca²⁺ binding to mNCS-1, obtained in our current study, are generally in agreement with the data reported by Jeromin et al. (2004) including the revealed positive enthalpy of Ca²⁺ binding to the low-affinity site. Thus, we confirmed different calcium affinities of three EF-hands of mNCS-1 and the sequential mode of their filling (Figure 2 and Table 1). Such a mechanism agrees with previous NMR studies, according to which Ca²⁺-binding sites become occupied in the following order EF2→EF3→EF4 (Chandra et al., 2011). Considering the evaluations of the Ca²⁺ affinity of individual EF hands reported by Chandra et al., we can attribute calcium-binding constants KA¹, KA² and KA³, as calculated in the current work (Table 1), to EF3, EF2 and EF4, respectively.

In early magnesium binding experiments, nNCS-1 exhibited the non-cooperative coordination of two Mg²⁺ with a dissociation constant of 12 µM (Cox et al., 1994). Similar findings were reported for mNCS-1, based on ITC, NMR and mutagenesis studies (Aravind et al., 2008). According to our ITC data, the amount of Mg²⁺ bound to mNCS-1 tends toward three (Table 1). We speculate that such stoichiometry is a specific feature of myristoylated protein, where Mg²⁺ binds to EF2-EF4. It should be noted that the actual amount of protein-associated magnesium, which binds with low affinity, might be highly sensitive to the quality of the protein sample (i.e., the content of the nNCS-1 admixture or residual calcium) and may therefore be differently evaluated. Yet, all three studies including ours agree that Mg²⁺ antagonizes Ca²⁺ binding by reducing the

affinity of the respective sites of the protein. These data confirm competition between the ions for the same binding sites with a preference for calcium (Cox et al., 1994; Aravind et al., 2008).

Our brand-new finding is that myristoylated NCS-1 is capable of coordinating up to three zinc ions. The mechanism of zinc binding to the protein and the exact Zn²⁺-binding sites are yet to be determined. For the moment, based on our metal competition analysis, CD studies and molecular modeling, we can hypothesize that zinc binds to functional EF-hands of the protein. Indeed, the ability of EF-hands to coordinate Zn²⁺ was previously reported for another ubiquitous Ca²⁺-binding protein calmodulin by X-ray crystallographic studies (Warren et al., 2007). Based on the analysis of all Zn²⁺-binding proteins presented in PDB, we found that, in NCS-1, the density of chelating groups required for Zn²⁺ binding is located only in the loops of EF2 (the highest score), EF3 and EF4 (Figure 7A). According to our CD measurements, the interaction of zinc with apo-NCS-1 induces a decrease in the content of β-sheets and an increase in α-helical content, exactly as in the case of the binding of Ca²⁺ or Mg²⁺ to EF-hands (Figure 3D and Supplementary Table S1). Since NCS-1 contains only two short antiparallel β-sheets, which connect Ca²⁺-binding loops of EF1-EF2 and EF3-EF4 (Heidarsson et al., 2012; Pandalaneni et al., 2015), one can suggest that zinc binds to EF-hands of the protein. This conclusion is further supported by our ITC data, indicating that Zn²⁺-saturated NCS-1 does not bind magnesium and exhibits reduced stoichiometry of Ca²⁺-binding (Figure 2 and Table 1). Interestingly, Ca²⁺ binding to one of the sites in Zn²⁺-saturated NCS-1 is one order of magnitude higher in affinity than any of the sites in apo-NCS-1 (Table 1). Given the proposed model for sequential filling of the mNCS-1 by calcium in the order EF2→EF3→EF4 (Chandra et al., 2011), we hypothesize that Zn²⁺-bound EF2 may adopt a conformation that facilitates the binding of calcium to the remaining two sites. Consistently, Zn²⁺-bound EF-hands of calmodulin resembled an intermediate state in the chain of conformational transitions induced by Ca²⁺-binding (Warren et al., 2007).

The unique mode of zinc binding to Ca²⁺-saturated NCS-1 is predicted by QM/MM simulations of the associated molecular dynamics, based on the crystal structure of the respective NCS-1 form [PDB 5AEQ (Pandalaneni et al., 2015)]. In the absence of zinc, EF3 possesses the most favourable environment for the coordination of calcium among EF-hands of the protein, which agrees with its maximal affinity for Ca²⁺ (Chandra et al., 2011). At the same time, coordination of both cations in this site seems unlikely due to the absence of the required number of chelating groups (Figure 7C). Therefore, EF3 can bind strictly to one ion with a preference for calcium and the latter can replace zinc from the site but not vice versa. In contrast, EF2 and EF4 can accommodate both Ca²⁺ and Zn²⁺, at least under our *in silico* conditions. In both cases, calcium loses one chelator, in turn becoming coordinated by six oxygen atoms (Figures 7B,D). Yet, such a configuration is common for proteins (Pidcock and Moore, 2001). Furthermore, such a configuration would completely compensate for the high negative
charge in the EF2 (-2 in Ca²⁺-bound NCS-1). Thus, based on these observations, we suggest that Ca²⁺-saturated NCS-1 can accommodate up to two zinc ions, one in EF4 and the other one in EF2.

The proposed binding modes for Zn²⁺ and Ca²⁺ are generally in agreement with our ITC and spectroscopic data. Thus, Ca²⁺-saturated NCS-1 coordinates one or two zinc ions (Table 1), apparently in terms of EF4/EF2 yielding the Zn²⁺(Ca²⁺)-bound protein form, which does not significantly differ from the “open” Ca²⁺-bound conformer in the overall protein fold (Figure 4C) but possesses enhanced thermal stability (Figures 5B, 6E-H). In contrast, Zn²⁺-saturated NCS-1 coordinates two calcium ions (Table 1) yielding a Ca²⁺(Zn²⁺)-bound conformer. In this case, EF2 likely remains occupied by zinc, which could facilitate calcium binding to EF3 (and consequently to EF4) as suggested by the absence of a low-affinity Ca²⁺-binding site and an increased binding constant for the high-affinity site in the ITC data (Table 1). At the same time, calcium replaces zinc from EF3 and could replace or temporary co-reside with zinc in EF4. It cannot be ruled out, however, that one of the EF-hands, being occupied with zinc, might adopt a conformation that is favorable for calcium binding, thereby exhibiting increased Ca²⁺ affinity as seen in our ITC studies. In any case, the resulting Ca²⁺(Zn²⁺)-bound conformer possesses only a small difference with the Zn²⁺(Ca²⁺)-bound from of the protein in the I₃₅₀/I₃₃₀ ratio (Figures 4C,D), but significantly differs from it in thermal stability (> 20°C, Figures 6E,F). It should be emphasized that, despite being highly consistent with the experimental and literature data, the above mechanisms of Zn²⁺/Ca²⁺ binding are mostly speculative and require additional confirmations.

In the aggregate, our *in vitro* studies suggest the existence of Zn²⁺-bound, Zn²⁺(Mg²⁺)-bound, Zn²⁺(Ca²⁺)-bound and Ca²⁺(Zn²⁺)-bound conformers of NCS-1 in addition to previously recognized apo, Mg²⁺-bound and Ca²⁺-bound forms of the protein. It should be mentioned that structural differences between the two latter forms, as observed in this study, are generally in accord with the reported data. Thus, the binding of both Mg²⁺ and Ca²⁺ increases the α-helical content of NCS-1, whereas only Ca²⁺ binding notably increases its surface hydrophobicity, as originally described by Jeromin et al. (2004). In addition, Mg²⁺ binding induced a more pronounced increase in the intensity of intrinsic fluorescence of the protein than Ca²⁺ binding, in agreement with previous observations (Aravind et al., 2008). It has been suggested that the unique mode of Ca²⁺/Mg²⁺ binding and resulting structural alterations govern the target recognition by NCS-1. Indeed, the NCS-1, preloaded with Mg²⁺, binds D2R in response to Ca²⁺ elevation more efficiently when compared to apo-protein, indicating that magnesium can serve as a physiological co-factor with calcium in this interaction (Woll et al., 2011).

Alongside NCS-1, the Ca²⁺/Mg²⁺ interplay was shown to regulate the structure and function of the other NCS proteins belonging to all five classes of the NCS family. Interestingly, the mechanisms of this regulation are quite distinct. Thus, magnesium and calcium bind to different EF-hand motifs of these proteins and the binding differently affects their functional

specificity. In recoverin, Mg²⁺ binds to functional EF2 and EF3, which reduces the Ca²⁺ affinity of the protein (at high magnesium concentration), but only slightly affects its secondary and tertiary structure, does not lead to activation of its myristoyl switch and is not required for its interaction with GRK1 (Ozawa et al., 2000; Ames et al., 2006; Marino et al., 2015). The cooperative sequential binding of calcium to the EF3 and EF2 of recoverin increases its thermal stability and α-helical content, as well as leads to exposure of its myristoyl group and hydrophobic pocket residues, thereby providing the protein with a capability to interact with membranes and GRK1 (Zozulya and Stryer, 1992; Ames et al., 1997, 2006; Permyakov et al., 2000; Zernii et al., 2015). A similar mechanism of Ca²⁺ binding (cooperative binding to EF2 and EF3), structural alterations (increase in α-helical content in the presence of Ca²⁺ but not Mg²⁺) and a Ca²⁺-myristoyl switch were recognized in the case of another NCS protein, VILIP1. However, unlike recoverin, VILIP1 only coordinates magnesium in EF3 with a relatively high affinity (K_D = 20 µM), suggesting the functional significance of this complex. Furthermore, VILIP-1 forms a stable dimer, which is not dependent on Ca²⁺ or Mg²⁺, but seems to be required for proper target recognition (Jheng et al., 2006; Li et al., 2011). In the proteins belonging to another class of the NCS family, GCAPs, Mg²⁺ and Ca²⁺ play a crucial role in tuning their activity toward target enzymes, i.e., retinal guanylate cyclases (GCs). For instance, in GCAP1, Mg²⁺ binds to EF2 with micromolar affinity (EF3 and EF4 exhibit only low affinity with the cation) and the binding stabilizes a tertiary structure of the protein, which otherwise represents a molten globule incapable of regulating GCs (Lim et al., 2009; Dell'Orco et al., 2010). Thus, the presence of magnesium in EF2 is necessary for maintaining a GC-activator state of GCAP1 (Peshenko and Dizhoor, 2004; Lim et al., 2016). Calcium binds to EF2, EF3 and EF4 of GCAP1 in a non-cooperative manner, which drastically increases the thermal stability of the protein without altering its secondary structure and triggering exposure of its myristoyl group (Lim et al., 2009; Marino et al., 2015). Instead, the binding converts GCAP1 into a GC-inhibitor state by inducing local conformational changes via the Ca²⁺-myristoyl tug mechanism (Peshenko et al., 2012; Lim et al., 2016). Finally, a rather different mechanism for Ca²⁺/Mg²⁺-dependent regulation was reported for NCS protein of the KChIP class, i.e., KChIP3, also known as the transcriptional repressor DREAM. In the absence of magnesium, this protein binds Ca²⁺ non-cooperatively in the following sequence EF3→EF4→EF2. Interestingly, the apo-form of KChIP3 coordinates Mg²⁺ with high affinity (K_D = 13 µM) in EF2 (EF3 and EF4 bind Mg²⁺ in the millimolar range), and this bound magnesium cannot be replaced by calcium, suggesting that, under cellular conditions, the protein will exist in either Mg²⁺-bound, or 2Ca²⁺/Mg²⁺-bound forms. Consistently, Mg²⁺-bound KChIP3 exists as a monomer and can specifically recognize target DNA elements, whereas Ca²⁺ binding to EF3 and/or EF4 induces dimerization of the protein and suppresses DNA binding. Similar to GCAP1, apo-KChIP3 represents a molten globule and Ca²⁺/Mg²⁺ binding enhances its stability (Osawa et al., 2005).
Overall, Ca²⁺/Mg²⁺ interplay governs the structural and functional properties of the majority of NCS proteins, although they exhibit different modes of regulation. Meanwhile, the involvement of zinc ions in this regulation so far has only been determined for recoverin. Similar to NCS-1, recoverin binds Zn²⁺, regardless of the presence of calcium, while the binding only slightly affects the secondary structure of the protein and destabilizes its Ca²⁺-saturated form. In recoverin, Zn²⁺ was proposed to be coordinated outside EF-hands since it binds to “inactivated” mutant with E→Q substitutions in the 12th position of the loop of functional EF2 and EF3 (E85Q/E121Q). However, this conclusion does not seem to be strict, as our current calculations indicate that E→Q mutation in such a position does not necessarily prevent the four-chelator coordination of Zn²⁺, which becomes bound by the other chelators in the loop. Consistently, such mutation does not prevent the six-chelator coordination of Mg²⁺ in EF-hands (Cates et al., 1999). Thus, it cannot be excluded that, similar to NCS1, Ca²⁺-loaded recoverin binds Zn²⁺ in one of the functional EF-hands. In this case, the reduced stoichiometry of zinc binding to NCS1 (2 Zn²⁺ per protein), compared to recoverin (1 Zn²⁺ per protein), can be explained by the fact that the latter contains a smaller amount of functional EF-hands: its EF4 is naturally non-functional due to substitutions of the metal coordinating residues in the first and third positions of the EF-hand loop. It should be noted that, unlike NCS-1, recoverin exhibits an increased affinity with photoreceptor membranes in the presence of zinc (Permyakov et al., 2003). Thus, although the coordination of Zn²⁺ may be a common property of NCS proteins, it produces somewhat different effects concerning their function, which are similar to those observed in the case of Ca²⁺ and Mg²⁺. Zinc binding may therefore additionally diversify specific regulation of NCS proteins.

It still remains an open question as to which of the discovered Zn²⁺-bound conformers of NCS-1 (see above) dominate under physiological conditions. In contrast to the well-recognized physiological role of calcium in cell signaling, zinc has long been considered as a solely structural component of proteins. Thus, being bound with picomolar to nanomolar affinities, presumably to sulfur- and nitrogen-containing ligands in tetrahedral coordination, zinc normally serves to maintain the structure and function of enzymes, transcription factors, receptors and signaling proteins (Maret and Li, 2009). According to our data, NCS-1 binds zinc transiently with a much lower affinity and likely to the sites in EF-hands. Assuming that, in neurons, the binding will occur against the background of a constantly high [up to 1-2 mM (Romani and Scarpa, 1992)] magnesium concentration and recurring elevations [up to 1-2 µM (Sabatini et al., 2002)] of calcium concentration, one might suggest physiological relevance only for Zn²⁺(Mg²⁺)-bound and Zn²⁺(Ca²⁺)-bound conformers of NCS-1 in addition to the well-known Mg²⁺-bound and Ca²⁺-bound forms. Nevertheless, even the formation of two additional forms might extend the functional repertoire of the protein.

Our results suggest that the binding of zinc to NCS-1 required micromolar concentration of the free cation. However, it has been generally accepted that, in contrast to magnesium and calcium,

both extracellular and intracellular free zinc concentration is low. Indeed, cytosolic zinc levels are regulated by a complex Zn²⁺-buffering system and the so-called “muffling reactions”, involving buffer proteins of metallothionein (picomolar affinity with zinc) class, as well as transporters, such as ZnTs, ZIPs and DCTs, which shuttle Zn²⁺ outside the cell or into subcellular stores including mitochondria, Golgi apparatus and lysosomes (Cousins et al., 2006; Colvin et al., 2010) [for review, see (Colvin et al., 2000)]. As a result, although the total concentration of zinc in cells reaches 0.2 mM (Colvin et al., 2008), the levels of free zinc in the cytoplasm were estimated as picomolar to low micromolar (Krezel and Maret, 2006). In this case, what are the physiological conditions in which the binding of zinc to NCS-1 can occur? The growing evidence indicates that, under certain conditions, the intracellular zinc levels can transiently increase, while zinc can perform signaling functions by playing complementary signaling roles with calcium (Maret, 2001). This is especially valid for the nervous system, as it is characterized by the highest extracellular zinc concentration and Zn²⁺ is known to be specifically accumulated in neurons (Frederickson et al., 2005). The hallmark of neuronal Zn²⁺ is its neurotransmitter function, along with glutamate in so-called “gluzinergic” neurons of forebrain. In presynaptic terminals, the cation is accumulated in ZnT3-loaded synaptic vesicles and undergoes a Ca²⁺-induced release into a synaptic cleft, where it can modulate various ionotropic and metabotropic receptors. The resulting high zinc concentration in the cleft (raised from 0.5 to 300 µM) can be pumped back to the presynaptic cell by ZnT3, or permeate into the postsynaptic neurons through calcium channels, thereby increasing the local cytosolic level of the cation (Frederickson and Bush, 2001). Furthermore, under certain conditions, zinc can be released from intracellular sources. For instance, in the hippocampal neurons exposure to glutamate-induced Ca²⁺ influx triggers cytosolic acidification and intracellular Zn²⁺ release (Kiedrowski, 2012). The resulting zinc signals could be recognized by specialized Zn²⁺-binding sites in neuronal proteins in order to conduct biochemical stimuli (Maret, 2006; Maret and Li, 2009). The levels of zinc are also high in neuroretina, where the maximal amount of total zinc was found in inner segments and synaptic terminals of photoreceptors cells, suggesting that it may participate in photoreceptor metabolism and neurotransmission (Ugarte et al., 2012; Ugarte and Osborne, 2014). Similar to gluzinergic neurons, photoreceptors have been suggested as releasing zinc together with glutamate in synapses. Indeed, both neurotransmitters were found in synaptic vesicles of the outer plexiform layer, which also contained ZnT3 (Akagi et al., 2001; Ugarte and Osborne, 2014). Importantly, photoreceptors contain considerable amounts of free zinc (or so-called “loosely bound zinc”), while its concentration varies depending on the light conditions, which, in turn, are likely to generate zinc signals (Ugarte and Osborne, 2014). Our data suggest that such signals in CNS and retinal neurons can be detected and transmitted by NCS-1. In ITC studies, the presence of zinc enhances the binding of Ca²⁺-NCS-1 to D2R (Figures 8A,B and Table 2), a process known to suppress the desensitization of the receptor. Therefore, upon receiving a joint Ca²⁺/Zn²⁺ stimulus, NCS-1 can modulate dopamine signaling in a specific
enhanced manner. In addition, Zn²⁺(Ca²⁺)-bound NCS-1
can specifically regulate the desensitization of D2R or other
homologous receptors by GRKs. Indeed, the presence of low zinc
levels improved the binding of Ca²⁺-NCS-1 to GRK1, whereas
the subsequent elevations in zinc concentration produced the
opposite effect (Figure 8C). Finally, Zn²⁺ binding can affect
the well-recognized function of NCS-1 in the Ca²⁺-dependent
regulation of neurotransmission and synaptic plasticity, as the
most pronounced increase in intracellular zinc is expected to be
in synaptic terminals, while NCS-1 is known to be specifically
accumulated in this part of the neurons including photoreceptors
(De Raad et al., 1995; Tsujimoto et al., 2002; Sippy et al., 2003;
Negyessy and Goldman-Rakic, 2005).

On the other side, uncontrolled elevations of zinc in neurons
and the consequential impairment of Zn²⁺ and Ca²⁺ ion
interplay can produce pathological effects. One of the principal
causes of these elevations is thought to be oxidative stress
(Wood and Osborne, 2001; Sensi et al., 2003; Sheline et al.,
2010b). Thus, in CNS neurons, zinc becomes released from
metallothioneins in response to the oxidation or nitrosylation
of their cysteine residues (Bossy-Wetzel et al., 2004). In
photoreceptors, the major source for pathological free zinc is
rhodopsin, which coordinates seven zinc ions per dimer and can
lose them in response to light-induced oxidative stress (Sheline
et al., 2010b). All these events have close ties to neurological
and neuro-ophthalmological disorders. Thus, increased zinc
concentration in postsynaptic neurons was shown to promote
excitotoxic cell death after seizures and mechanical brain
trauma, while zinc chelators were found to be neuroprotective
(Suh et al., 2000). Consistently, exposure of neurons to high
zinc concentrations or induction of intracellular zinc release
promotes their apoptosis (Manev et al., 1997; Bossy-Wetzel et al.,
2004). Elevations of retinal zinc are associated with ischemia,
trophic deprivation or hypoglycemia leading to neuronal death
(Yoo et al., 2004; Suh et al., 2008; Sheline et al., 2010a).
The mechanisms of zinc toxicity generally involve massive
Zn²⁺-induced aggregation of neuronal proteins (Cuajungco
and Faget, 2003). For instance, an increase in extracellular
zinc is a key factor in the aggregation of amyloid plaques
in Alzheimer's disease. Furthermore, overexpression of zinc
transporter proteins was observed in patients in the early
stages of this disorder (Lyubartseva and Lovell, 2012). In the
current study, we have demonstrated that non-specific low-
affinity binding of excessive zinc promotes the aggregation
and precipitation of NCS-1, which are associated with the
formation of fibrilic twisted rope-like structures of the protein
(Supplementary Figure S2). These structures resemble the
Zn²⁺-induced aggregates of another neuronal protein, TDP-
43, the aggregation of which is associated with amyotrophic
lateral sclerosis and frontotemporal lobar degeneration (Garnier
et al., 2017). There are no direct indications concerning the
involvement of NCS-1 in the pathogenesis of neurodegenerative
diseases. Meanwhile, its transcription levels were found to be
altered in Alzheimer's disease, which is known to be associated
with altered zinc homeostasis (Karim et al., 2014). Furthermore,
a large body of evidence supports the neuroprotective role of
NCS-1 (Nakamura et al., 2006; Yip et al., 2010), which might be

suppressed upon the Zn²⁺-induced loss of the protein structure. In addition, it has been proposed that NCS-1 misfolding, together with calcium dysregulation, contributes to neurodegeneration (Heidarsson et al., 2014). These observations are in accord with the fact that the partial proteolytic degradation of NCS-1 and the loss of intracellular calcium signaling induce peripheral neuropathy associated with chemotherapy by paclitaxel (Boeckel and Ehrlich, 2018).

In summary, our study suggests that the complex interplay
between magnesium, calcium and zinc ions results in the
appearance of multiple conformations of NCS-1, thereby
modulating its functional status. It also indicates that the extreme
elevation of zinc levels peculiar to some neurodegenerative and
neuro-ophthalmological disorders may cause the formation of
unstable Zn²⁺-bound conformers of NCS-1 and promote its
aggregation. Further studies are required for unraveling the
molecular mechanism and exact sites of zinc binding to NCS-1
and firmly establishing of physiological and pathological roles of
this phenomenon.

AUTHOR CONTRIBUTIONS

PT, AR, and FD performed ITC, DSF, and TEM studies. VB, AZ, and EZ performed functional assays and ITC studies. AN, MS, and SP performed fluorimetric, CD, contributed to equilibrium dialysis, and AAS experiments. VV and DZ performed expression and purification of the proteins and analytical gel-filtration. MB and AG performed molecular modeling and QM/MM molecular dynamics simulations. PT, SP, and EZ wrote the article.

FUNDING

This work was supported by grants from the Russian Foundation for Basic Research (15-04-07963 and 18-04-01250) to EZ.

ACKNOWLEDGMENTS

The research is carried out using the equipment of the
shared research facilities of HPC computing resources at
Lomonosov Moscow State University supported by the project
RFMEFI62117X0011. ITC and nanoDSF measurements were
performed in Plateau Microcalorimetry Timone.

SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online
at: https://www.frontiersin.org/articles/10.3389/fnmol.2018.
00459/full#supplementary-material

**FIGURE S1** | Binding of NCS-1 to photoreceptor membranes in the presence of zinc. Weight fractions of NCS-1 (25 µM) bound to urea-washed bovine photoreceptor membranes in the presence of 0, 25, 50, 75, or 100 µM Zn²⁺ [i.e., at (Zn²⁺)/(NCS-1) ratio of 0-4] under Ca²⁺-free conditions (–Ca²⁺) or on the background of 1 mM Ca²⁺ (+Ca²⁺) according to the data from equilibrium centrifugation assay.
**FIGURE S2** | Structural properties of NCS-1 in the presence of excessive zinc concentrations. (A) Mid-transition temperatures of NCS-1 (25 µM) aggregation in the presence of 0.1–2000 µM Zn²⁺ on the background of either 1 mM Mg²⁺ or 1 mM Ca²⁺ determined form light scattering at 350 nm. (B) Weight fractions of NCS-1 (25 µM) precipitated in the presence of 0–500 µM Zn²⁺ on the background of either 1 mM Mg²⁺ or 1 mM Ca²⁺. (C,D) Electron microphotographs

of NCS-1 aggregates, obtained in the presence of 5 mM Zn²⁺; scale: 0.5 µm (C), 0.2 µm (D).

**TABLE S1** | Secondary structure content (in %) of NCS-1 estimated from CD data obtained at 25°C (10 mM Hepes-KOH buffer (pH 7.6), 100 mM KCl) under metal-free conditions or in the presence of Ca²⁺ (100 µM CaCl₂), Mg²⁺ (1 mM MgCl₂), Zn²⁺ (100 µM ZnCl₂) or their combinations.

## REFERENCES

Abraham, M. J., Murtola, T., Schulz, R., Päll, S., Smith, J. C., Hess, B., et al. (2015). GROMACS: high performance molecular simulations through multi-level parallelism from laptops to supercomputers. *SoftwareX* 1, 19–25. doi: 10.1016/j.softx.2015.06.001

Akagi, T., Kaneda, M., Ishii, K., and Hashikawa, T. (2001). Differential subcellular localization of zinc in the rat retina. *J. Histochem. Cytochem.* 49, 87–96. doi: 10.1177/002215540104900109

Ames, J. B., Ishima, R., Tanaka, T., Gordon, J. I., Stryer, L., and Ikura, M. (1997). Molecular mechanics of calcium-myristoyl switches. *Nature* 389, 198–202. doi: 10.1038/38310

Ames, J. B., Levay, K., Wingard, J. N., Lusin, J. D., and Slepak, V. Z. (2006). Structural basis for calcium-induced inhibition of rhodopsin kinase by recoverin. *J. Biol. Chem.* 281, 37237–37245. doi: 10.1074/jbc.M606913200

Aravind, P., Chandra, K., Reddy, P. P., Jeromin, A., Chary, K., and Sharma, Y. (2008). Regulatory and structural EF-hand motifs of neuronal calcium sensor-1: Mg2⁺ modulates Ca2⁺ binding, Ca2⁺-induced conformational changes, and equilibrium unfolding transitions. *J. Mol. Biol.* 376, 1100–1115. doi: 10.1016/j.jmb.2007.12.033

Augustine, G. J., Santamaria, F., and Tanaka, K. (2003). Local calcium signaling in neurons. *Neuron* 40, 331–346. doi: 10.1016/S0896-6273(03) 00639-1

Baksheeva, V., Nazipova, A., Zinchenko, D., Serebryakova, M., Senin, I. I., Permyakov, S., et al. (2015). Ca2⁺-Myristoyl Switch in Neuronal Calcium Sensor-1: a Role of C-Terminal Segment. *CNS Neurol. Disord. Drug Targets* 14, 437–451. doi: 10.2174/1871527314666150225143403

Banker, M. J., Clark, T. H., and Williams, J. A. (2003). Development and validation of a 96-well equilibrium dialysis apparatus for measuring plasma protein binding. *J. Pharm. Sci.* 92, 967–974. doi: 10.1002/jps.10332

Barwinska-Sendra, A., and Waldron, K. J. (2017). "Chapter Eight – The role of intermetal competition and mis-metalation in metal toxicity," in *Advances in Microbial Physiology*, ed. R.K. Poole (Cambridge, MA: Academic Press), 315–379.

Blachford, C., Ćelić, A., Petri, E. T., and Ehrlich, B. E. (2009). Discrete proteolysis of neuronal calcium sensor-1 (NCS-1) by µ-calpain disrupts calcium binding. *Cell Calcium* 46, 257–262. doi: 10.1016/j.cca.2009.08.002

Boeckel, G. R., and Ehrlich, B. E. (2018). NCS-1 is a regulator of calcium signaling in health and disease. *Biochim. Biophys. Acta Mol. Cell Res.* [Epub ahead of print], doi: 10.1016/j.bbamcr.2018.05.005

Bossy-Wetzel, E., Talantova, M. V., Lee, W. D., Schölzke, M. N., Harrop, A., Mathews, E., et al. (2004). Crosstalk between nitric oxide and zinc pathways to neuronal cell death involving mitochondrial dysfunction and p38-activated K⁺ channels. *Neuron* 41, 351–365. doi: 10.1016/S0896-6273(04) 00015-7

Burgoyne, R. D. (2007). Neuronal calcium sensor proteins: generating diversity in neuronal Ca2⁺ signalling. *Nat. Rev. Neurosci.* 8, 182–193. doi: 10.1038/nrn2093

Burgoyne, R. D., and Haynes, L. P. (2012). Understanding the physiological roles of the neuronal calcium sensor proteins. *Mol. Brain* 5:2. doi: 10.1186/1756-66 06-5-2

Burstein, E. A., and Emelyanenko, V. I. (1996). Log-normal description of fluorescence spectra of organic fluorophores. *Photochem. Photobiol.* 64, 316– 320. doi: 10.1111/j.1751-1097.1996.tb02464.x

Bussi, G., Donadio, D., and Parrinello, M. (2007). Canonical sampling through velocity rescaling. *J. Chem. Phys.* 126:014101.

Cates, M. S., Berry, M. B., Ho, E. L., Li, Q., Potter, J. D., Phillips, G. N. Jr. (1999). Metal-ion affinity and specificity in EF-hand proteins: coordination geometry and domain plasticity in parvalbumin. *Structure* 7, 1269–1278. doi: 10.1016/ S0969-2126(00)80060-X

Chandra, K., Ramakrishnan, V., Sharma, Y., and Chary, K. (2011). N-terminal myristoylation alters the calcium binding pathways in neuronal calcium sensor- 1. *J. Biol. Inorg. Chem.* 16, 81–95. doi: 10.1007/s00775-010-0705-3

Choi, D. W., and Koh, J. Y. (1998). Zinc and brain injury. *Annu. Rev. Neurosci.* 21, 347–375. doi: 10.1146/annurev.neuro.21.1.347

Colvin, R. A., Bush, A. I., Volitakis, I., Fontaine, C. P., Thomas, D., Kikuchi, K., et al. (2008). Insights into Zn2⁺ homeostasis in neurons from experimental and modeling studies. *Am. J. Physiol. Cell Physiol.* 294, C726–C742. doi: 10.1152/ ajpcell.00541.2007

Colvin, R. A., Davis, N., Nipper, R. W., and Carter, P. A. (2000). Zinc transport in the brain: routes of zinc influx and efflux in neurons. *J. Nutr.* 130, 1484S–1487S. doi: 10.1093/jn/130.5.1484S

Colvin, R. A., Holmes, W. R., Fontaine, C. P., and Maret, W. (2010). Cytosolic zinc buffering and muffling: their role in intracellular zinc homeostasis. *Metallomics* 2, 306–317. doi: 10.1039/b926662c

Cousins, R. J., Liuzzi, J. P., and Lichten, L. A. (2006). Mammalian zinc transport, trafficking, and signals. *J. Biol. Chem.* 281, 24085–24089. doi: 10.1074/jbc. R600011200

Cox, J. A., Durussel, I., Comte, M., Nef, S., Nef, P., Lenz, S. E., et al. (1994). Cation binding and conformational changes in VILIP and NCS-1, two neuron-specific calcium-binding proteins. *J. Biol. Chem.* 269, 32807–32813.

Cuajungco, M. P., and Faget, K. Y. (2003). Zinc takes the center stage: its paradoxical role in Alzheimer's disease. *Brain Res. Rev.* 41, 44–56. doi: 10.1016/ S0165-0173(02)00219-9

Dapprich, S., Komáromi, I., Byun, K. S., Morokuma, K., and Frisch, M. J. (1999). A new ONIOM implementation in Gaussian98. Part I. The calculation of energies, gradients, vibrational frequencies and electric field derivatives1. *J. Mol. Struct. THEOCHEM* 461–462, 1–21. doi: 10.1016/S0166-1280(98)00 475-8

De Raad, S., Comte, M., Nef, P., Lenz, S. E., Gundelfinger, E. D., and Cox, J. A. (1995). Distribution pattern of three neural calcium-binding proteins (NCS-1, VILIP and recoverin) in chicken, bovine and rat retina. *Histochem. J.* 27, 524–535. doi: 10.1007/BF00174325

DeLano, W. L. (2002). *The PyMOL Molecular Graphics System*. Available at: http://www.pymol.org

Dell'Orco, D., Behnen, P., Linse, S., and Koch, K. W. (2010). Calcium binding, structural stability and guanylate cyclase activation in GCAP1 variants associated with human cone dystrophy. *Cell. Mol. Life Sci.* 67, 973–984. doi: 10.1007/s00018-009-0243-8

Dudev, T., and Lim, C. (2003). Principles governing Mg, Ca, and Zn binding and selectivity in proteins. *Chem. Rev.* 103, 773–788. doi: 10.1021/cr020467n

Frederickson, C. J., and Bush, A. I. (2001). Synaptically released zinc: physiological functions and pathological effects. *Biometals* 14, 353–366. doi: 10.1023/A: 1012934207456

Frederickson, C. J., Koh, J. Y., and Bush, A. I. (2005). The neurobiology of zinc in health and disease. *Nat. Rev. Neurosci.* 6, 449–462. doi: 10.1038/nrn1671

Gardner, J. M. M., Powell, C. A., Baker-Henningham, H., Walker, S. P., Cole, T. J., and Grantham-McGregor, S. M. (2005). Zinc supplementation and psychosocial stimulation: effects on the development of undernourished Jamaican children. *Am. J. Clin. Nutr.* 82, 399–405. doi: 10.1093/ajcn/82.2.399

Garnier, C., Devred, F., Byrne, D., Puppo, R., Roman, A. Y., Malesinski, S., et al. (2017). Zinc binding to RNA recognition motif of TDP-43 induces the formation of amyloid-like aggregates. *Sci. Rep.* 7:6812. doi: 10.1038/s41598- 017-07215-7

Gaus, M., Cui, Q., and Elstner, M. (2012). DFTB3: extension of the self-consistent-charge density-functional tight-binding method (SCC-DFTB). *J. Chem. Theory Comput.* 7, 931–948. doi: 10.1021/ct100684s

Grigoriev, I. I., Senin, I. I., Tikhomirova, N. K., Komolov, K. E., Permyakov, S. E., Zernii, E. Y., et al. (2012). Synergetic effect of recoverin and calmodulin on
regulation of rhodopsin kinase. *Front. Mol. Neurosci.* 5:28. doi: 10.3389/fnmol. 2012.00028

Grimme, S., Antony, J., Ehrlich, S., and Krieg, H. (2010). A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. *J. Chem. Phys.* 132:154104. doi: 10.1063/1.3382344

Heidarsson, P. O., Bjerrum-Bohr, I. J., Jensen, G. A., Pongs, O., Finn, B. E., Poulsen, F. M., et al. (2012). The C-terminal tail of human neuronal calcium sensor 1 regulates the conformational stability of the Ca2+-activated state. *J. Mol. Biol.* 417, 51–64. doi: 10.1016/j.jmb.2011.12.049

Heidarsson, P. O., Naqvi, M. M., Otazo, M. R., Mossa, A., Kragelund, B. B., and Cecconi, C. (2014). Direct single-molecule observation of calcium-dependent misfolding in human neuronal calcium sensor-1. *Proc. Natl. Acad. Sci. U.S.A.* 111, 13069–13074. doi: 10.1073/pnas.1401065111

Jeromin, A., Muralidhar, D., Parameswaran, M. N., Roder, J., Fairwell, T., Scarlata, S., et al. (2004). N-terminal myristoylation regulates calcium-induced conformational changes in neuronal calcium sensor-1. *J. Biol. Chem.* 279, 27158–27167. doi: 10.1074/jbc.M312172200

Jheng, F. F., Wang, L., Lee, L., and Chang, L. S. (2006). Functional contribution of Ca2+ and Mg2+ to the intermolecular interaction of visinin-like proteins. *Protein J.* 25, 250–256. doi: 10.1007/s10930-006-9008-5

Karim, S., Mirza, Z., Ansari, S. A., Rasool, M., Iqbal, Z., Sohrab, S. S., et al. (2014). Transcriptomics study of neurodegenerative disease: emphasis on synaptic dysfunction mechanism in Alzheimer's disease. *CNS Neurol. Disord. Drug Targets* 13, 1202–1212. doi: 10.2174/1871527313666140917113446

Kazakov, A. S., Sokolov, A. S., Vologzhannikova, A. A., Permyakova, M. E., Khorn, P. A., Ismailov, R. G., et al. (2017). Interleukin-11 binds specific EF-hand proteins via their conserved structural motifs. *J. Biomol. Struct. Dyn.* 35, 78–91. doi: 10.1080/07391102.2015.1132392

Kiedrowski, L. (2012). Cytosolic acidification and intracellular zinc release in hippocampal neurons. *J. Neurochem.* 121, 438–450. doi: 10.1111/j.1471-4159. 2012.07695.x

Koch, K. W., and Dell'Orco, D. (2015). Protein and signaling networks in vertebrate photoreceptor cells. *Front. Mol. Neurosci.* 8:67. doi: 10.3389/fnmol.2015.00067

Komolov, K. E., Senin, I. I., Kovaleva, N. A., Christoph, M. P., Churumova, V. A., Grigoriev, I. I., et al. (2009). Mechanism of rhodopsin kinase regulation by recoverin. *J. Neurochem.* 110, 72–79. doi: 10.1111/j.1471-4159.2009.06118.x

Krezel, A., and Maret, W. (2006). Zinc-buffering capacity of a eukaryotic cell at physiological pZn. *J. Biol. Inorg. Chem.* 11, 1049–1062. doi: 10.1007/s00775- 006-0150-5

Kubař, T., Welke, K., and Groenhof, G. (2015). New QM/MM implementation of the DFTB3 method in the gromacs package. *J. Comput. Chem.* 36, 1978–1989. doi: 10.1002/jcc.24029

Laitaoja, M., Valjakka, J., and Jänis, J. (2013). Zinc coordination spheres in protein structures. *Inorg. Chem.* 52, 10983–10991. doi: 10.1021/ic401072d

Leung, K. W., Gvritishvili, A., Liu, Y., and Tombran-Tink, J. (2012). ZIP2 and ZIP4 mediate age-related zinc fluxes across the retinal pigment epithelium. *J. Mol. Neurosci.* 46, 122–137. doi: 10.1007/s12031-011-9536-0

Li, C., Pan, W., Braunewell, K. H., and Ames, J. B. (2011). Structural analysis of Mg2+ and Ca2+ binding, myristoylation, and dimerization of the neuronal calcium sensor and visinin-like protein 1 (VILIP-1). *J. Biol. Chem.* 286, 6354– 6366. doi: 10.1074/jbc.M110.173724

Lim, S., Peshenko, I., Dizhoor, A., and Ames, J. B. (2009). Effects of Ca2+, Mg2+, and myristoylation on guanylyl cyclase activating protein 1 structure and stability. *Biochemistry* 48, 850–862. doi: 10.1021/bi801897p

Lim, S., Peshenko, I. V., Olshevskaya, E. V., Dizhoor, A. M., and Ames, J. B. (2016). Structure of guanylyl cyclase activator protein 1 (GCAP1) mutant V77E in a Ca2+-free/Mg2+-bound activator state. *J. Biol. Chem.* 291, 4429–4441. doi: 10.1074/jbc.M115.696161

Lindorff-Larsen, K., Piana, S., Palmo, K., Maragakis, P., Klepeis, J. L., Dror, R. O., et al. (2010). Improved side-chain torsion potentials for the Amber ff99SB protein force field. *Proteins* 78, 1950–1958. doi: 10.1002/prot.22711

Lyubartseva, G., and Lovell, M. A. (2012). A potential role for zinc alterations in the pathogenesis of Alzheimer's disease. *Biofactors* 38, 98–106. doi: 10.1002/biof.199

Manev, H., Kharlamov, E., Uz, T., Mason, R. P., and Cagnoli, C. M. (1997). Characterization of zinc-induced neuronal death in primary cultures of rat cerebellar granule cells. *Exp. Neurol.* 146, 171–178. doi: 10.1006/exnr.1997.6510

Maret, W. (2001). Crosstalk of the group IIa and IIb metals calcium and zinc in cellular signaling. *Proc. Natl. Acad. Sci. U.S.A.* 98, 12325–12327. doi: 10.1073/ pnas.231481398

Maret, W. (2006). Zinc coordination environments in proteins as redox sensors and signal transducers. *Antioxid. Redox Signal.* 8, 1419–1441. doi: 10.1089/ars. 2006.8.1419

Maret, W., and Li, Y. (2009). Coordination dynamics of zinc in proteins. *Chem. Rev.* 109, 4682–4707. doi: 10.1021/cr800556u

Marino, V., Sulmann, S., Koch, K. W., and Dell'orco, D. (2015). Structural effects of Mg(2)(+) on the regulatory states of three neuronal calcium sensors operating in vertebrate phototransduction. *Biochim. Biophys. Acta* 1853, 2055–2065. doi: 10.1016/j.bbamcr.2014.10.026

Moroz, O. V., Wilson, K. S., and Bronstein, I. B. (2011). The role of zinc in the S100 proteins: insights from the X-ray structures. *Amino acids* 41, 761–772. doi: 10.1007/s00726-010-0540-4

Nakamura, T. Y., Jeromin, A., Smith, G., Kurushima, H., Koga, H., Nakabeppu, Y., et al. (2006). Novel role of neuronal Ca2+ sensor-1 as a survival factor up-regulated in injured neurons. *J. Cell Biol.* 172, 1081–1091. doi: 10.1083/jcb. 200508156

Negyessy, L., and Goldman-Rakic, P. S. (2005). Subcellular localization of the dopamine D2 receptor and coexistence with the calcium-binding protein neuronal calcium sensor-1 in the primate prefrontal cortex. *J. Comp. Neurol.* 488, 464–475. doi: 10.1002/cne.20601

Olshevskaya, E. V., Ermilov, A. N., and Dizhoor, A. M. (1999). Dimerization of guanylyl cyclase-activating protein and a mechanism of photoreceptor guanylyl cyclase activation. *J. Biol. Chem.* 274, 25583–25587. doi: 10.1074/jbc.274.36. 25583

Osawa, M., Dace, A., Tong, K. I., Valiveti, A., Ikura, M., and Ames, J. B. (2005). Mg2+ and Ca2+ differentially regulate DNA binding and dimerization of DREAM. *J. Biol. Chem.* 280, 18008–18014. doi: 10.1074/jbc.M500338200

Ozawa, T., Fukuda, M., Nara, M., Nakamura, A., Komine, Y., Kohama, K., et al. (2000). How can Ca2+ selectively activate recoverin in the presence of Mg2+? Surface plasmon resonance and FT-IR spectroscopic studies. *Biochemistry* 39, 14495–14503. doi: 10.1021/bi001930y

Pandalaneni, S., Karuppiah, V., Saleem, M., Haynes, L. P., Burgoyne, R. D., Mayans, O., et al. (2015). Neuronal Calcium Sensor-1 binds the D2 dopamine receptor and G-protein coupled receptor kinase 1 (GRK1) peptides using different modes of interactions. *J. Biol. Chem.* 290, 18744–18756. doi: 10.1074/ jbc.M114.627059

Penland, J. G., Sandstead, H. H., Alcock, N. W., Dayal, H. H., Chen, X. C., Li, J. S., et al. (1997). A preliminary report: effects of zinc and micronutrient repletion on growth and neuropsychological function of urban Chinese children. *J. Am. Coll. Nutr.* 16, 268–272. doi: 10.1080/07315724.1997.10718684

Permyakov, S. E., Cherskaya, A. M., Senin, I. I., Zargarov, A. A., Shulga-Morskoy, S. V., Alekseev, A. M., et al. (2000). Effects of mutations in the calcium-binding sites of recoverin on its calcium affinity: evidence for successive filling of the calcium binding sites. *Protein Eng.* 13, 783–790. doi: 10.1093/protein/13.11.783

Permyakov, S. E., Cherskaya, A. M., Wasserman, L. A., Khokhlova, T. I., Senin, I. I., Zargarov, A. A., et al. (2003). Recoverin is a zinc-binding protein. *J. Proteome Res.* 2, 51–57. doi: 10.1021/pr025553i

Permyakov, S. E., Zernii, E. Y., Knyazeva, E. L., Denesyuk, A. I., Nazipova, A. A., Kolpakova, T. V., et al. (2012). Oxidation mimicking substitution of conservative cysteine in recoverin suppresses its membrane association. *Amino acids* 42, 1435–1442. doi: 10.1007/s00726-011-0843-0

Peshenko, I. V., and Dizhoor, A. M. (2004). Guanylyl cyclase-activating proteins (GCAPs) are Ca2+/Mg2+ sensors: implications for photoreceptor guanylyl cyclase (RetGC) regulation in mammalian photoreceptors. *J. Biol. Chem.* 279, 16903–16906. doi: 10.1074/jbc.C400065200

Peshenko, I. V., Olshevskaya, E. V., Lim, S., Ames, J. B., and Dizhoor, A. M. (2012). Calcium-myristoyl Tug is a new mechanism for intramolecular tuning of calcium sensitivity and target enzyme interaction for guanylyl cyclase-activating protein 1: dynamic connection between N-fatty acyl group and EF-hand controls calcium sensitivity. *J. Biol. Chem.* 287, 13972–13984. doi: 10.1074/jbc.M112.341883

Pidcock, E., and Moore, G. R. (2001). Structural characteristics of protein binding sites for calcium and lanthanide ions. *J. Biol. Inorg. Chem.* 6, 479–489. doi: 10.1007/s007750100214
Romani, A., and Scarpa, A. (1992). Regulation of cell magnesium. *Arch. Biochem. Biophys.* 298, 1–12. doi: 10.1016/0003-9861(92)90086-C

Sabatini, B. L., Oertner, T. G., and Svoboda, K. (2002). The life cycle of Ca(2+)-ions in dendritic spines. *Neuron* 33, 439–452. doi: 10.1016/S0896-6273(02)00573-1

Senin, I. I., Tikhomirova, N., Churumova, V., Grigoriev, I., Kolpakova, T., Zinchenko, D., et al. (2011). Amino acid sequences of two immune-dominant epitopes of recoverin are involved in Ca 2+/recoverin-dependent inhibition of phosphorylation of rhodopsin. *Biochemistry* 76, 332–338.

Sensi, S. L., Ton-That, D., Sullivan, P. G., Jonas, E. A., Gee, K. R., Kaczmarek, L. K., et al. (2003). Modulation of mitochondrial function by endogenous Zn2+ pools. *Proc. Natl. Acad. Sci. U.S.A.* 100, 6157–6162. doi: 10.1073/pnas.1031598100

Sheline, C. T., Cai, A. L., Zhu, J., and Shi, C. (2010a). Serum or target deprivation-induced neuronal death causes oxidative neuronal accumulation of Zn2+ and loss of NAD+. *Eur. J. Neurosci.* 32, 894–904. doi: 10.1111/j.1460-9568.2010.07372.x

Sheline, C. T., Zhou, Y., and Bai, S. (2010b). Light-induced photoreceptor and RPE degeneration involve zinc toxicity and are attenuated by pyruvate, nicotinamide, or cyclic light. *Mol. Vis.* 16, 2639–2652.

Sippy, T., Cruz-Martín, A., Jeromin, A., and Schweizer, F. E. (2003). Acute changes in short-term plasticity at synapses with elevated levels of neuronal calcium sensor-1. *Nat. Neurosci.* 6, 1031–1038. doi: 10.1038/nn1117

Sreerama, N., and Woody, R. W. (2000). Estimation of protein secondary structure from circular dichroism spectra: comparison of CONTIN, SELCON, and CDSSTR methods with an expanded reference set. *Anal. Biochem.* 287, 252–260. doi: 10.1006/abio.2000.4880

Suh, S. W., Chen, J. W., Motamedi, M., Bell, B., Listiak, K., Pons, N. F., et al. (2000). Evidence that synaptically-released zinc contributes to neuronal injury after traumatic brain injury. *Brain Res.* 852, 268–273. doi: 10.1016/S0006-8993(99)02095-8

Suh, S. W., Hamby, A. M., Gum, E. T., Shin, B. S., Won, S. J., Sheline, C. T., et al. (2008). Sequential release of nitric oxide, zinc, and superoxide in hypoglycemic neuronal death. *J. Cereb. Blood Flow Metab.* 28, 1697–1706. doi: 10.1038/jcbfm.2008.61

Tsujimoto, T., Jeromin, A., Saitoh, N., Roder, J. C., and Takahashi, T. (2002). Neuronal calcium sensor 1 and activity-dependent facilitation of P/Q-type calcium currents at presynaptic nerve terminals. *Science* 295, 2276–2279. doi: 10.1126/science.1068278

Tsvetkov, P., Barbier, P., Breuzard, G., Peyrot, V., and Devred, F. (2013). "Microtubule-associated proteins and tubulin interaction by isothermal titration calorimetry," in *Methods in cell biology*, eds J. J. Correia and L. Wilson (Amsterdam: Elsevier), 283–302.

Tsvetkov, P., Devred, F., and Makarov, A. (2010). Thermodynamics of zinc binding to human S100A2. *Mol. Biol.* 44, 832–835. doi: 10.1134/S0026893310050213

Tsvetkov, P. O., Makarov, A. A., Malesinski, S., Peyrot, V., and Devred, F. (2012). New insights into tau-microtubules interaction revealed by isothermal titration calorimetry. *Biochimie* 94, 916–919. doi: 10.1016/j.biochi.2011.09.011

Ugarte, M., Grime, G. W., Lord, G., Geraki, K., Collingwood, J. F., Finnegan, M. E., et al. (2012). Concentration of various trace elements in the rat retina and their distribution in different structures. *Metallomics* 4, 1245–1254. doi: 10.1039/c2mt20157g

Ugarte, M., and Osborne, N. N. (2014). Recent advances in the understanding of the role of zinc in ocular tissues. *Metallomics* 6, 189–200. doi: 10.1039/c3mt00291h

Vladimirov, V. I., Zernii, E. Y., Baksheeva, V. E., Wimberg, H., Kazakov, A. S., Tikhomirova, N. K., et al. (2018). Photoreceptor calcium sensor proteins in detergent-resistant membrane rafts are regulated via binding to caveolin-1. *Cell Calcium* 73, 55–69. doi: 10.1016/j.ceca.2018.04.003

Warren, J. T., Guo, Q., and Tang, W. J. (2007). A 1.3-Å structure of zinc-bound N-terminal domain of calmodulin elucidates potential early ion-binding step. *J. Mol. Biol.* 374, 517–527. doi: 10.1016/j.jmb.2007.09.048

Waters, N. J., Jones, R., Williams, G., and Sohal, B. (2008). Validation of a rapid equilibrium dialysis approach for the measurement of plasma protein binding. *J. Pharm. Sci.* 97, 4586–4595. doi: 10.1002/jps.21317

Weiergräber, O. H., Senin, I. I., Zernii, E. Y., Churumova, V. A., Kovaleva, N. A., Nazipova, A. A., et al. (2006). Tuning of a neuronal calcium sensor. *J. Biol. Chem.* 281, 37594–37602. doi: 10.1074/jbc.M603700200

Wills, N., Ramanujam, V. S., Kalariya, N., Lewis, J., and Van Kuijk, F. (2008). Copper and zinc distribution in the human retina: relationship to cadmium accumulation, age, and gender. *Exp. Eye Res.* 87, 80–88. doi: 10.1016/j.exer.2008.04.013

Woll, M. P., De Cotiis, D. A., Bewley, M. C., Tacelosky, D. M., Levenson, R., and Flanagan, J. M. (2011). Interaction between the D2 dopamine receptor and neuronal calcium sensor-1 analyzed by fluorescence anisotropy. *Biochemistry* 50, 8780–8791. doi: 10.1021/bi200637e

Wood, J. P., and Osborne, N. N. (2001). The influence of zinc on caspase-3 and DNA breakdown in cultured human retinal pigment epithelial cells. *Arch. Ophthalmol.* 119, 81–88.

Yanez, M., Gil-Longo, J., and Campos-Toimil, M. (2012). Calcium binding proteins. *Adv. Exp. Med. Biol.* 740, 461–482. doi: 10.1007/978-94-007-2888-2_19

Yip, P. K., Wong, L. F., Sears, T. A., Yáñez-Muñoz, R. J., and Mcmahon, S. B. (2010). Cortical overexpression of neuronal calcium sensor-1 induces functional plasticity in spinal cord following unilateral pyramidal tract injury in rat. *PLoS Biol.* 8:e1000399. doi: 10.1371/journal.pbio.1000399

Yoo, M. H., Lee, J. Y., Lee, S. E., Koh, J. Y., and Yoon, Y. H. (2004). Protection by pyruvate of rat retinal cells against zinc toxicity in vitro, and pressure-induced ischemia in vivo. *Invest. Ophthalmol. Vis. Sci.* 45, 1523–1530. doi: 10.1167/iovs.03-1315

Zernii, E. Y., Komolov, K. E., Permyakov, S. E., Kolpakova, T., Dell'orco, D., Poetzsch, A., et al. (2011). Involvement of the recoverin C-terminal segment in recognition of the target enzyme rhodopsin kinase. *Biochem. J.* 435, 441–450. doi: 10.1042/BJ20110013

Zernii, E. Y., Nazipova, A. A., Gancharova, O. S., Kazakov, A. S., Serebryakova, M. V., Zinchenko, D. V., et al. (2015). Light-induced disulfide dimerization of recoverin under ex vivo and in vivo conditions. *Free Radic. Biol. Med.* 83, 283–295. doi: 10.1016/j.freeradbiomed.2015.03.001

Zernii, E. Y., Tikhomirova, N., Philippov, P., and Senin, I. I. (2003). Detection of annexin IV in bovine retinal rods. *Biochemistry* 68, 129–134.

Zozulya, S., and Stryer, L. (1992). Calcium-myristoyl protein switch. *Proc. Natl. Acad. Sci. U.S.A.* 89, 11569–11573. doi: 10.1073/pnas.89.23.11569

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2018 Tsvetkov, Roman, Baksheeva, Nazipova, Shevelyova, Vladimirov, Buyanova, Zinchenko, Zamyatnin, Devred, Golovin, Permyakov and Zernii. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.