Next Article in Journal
Targeting SARS-CoV-2 Variants with Nucleic Acid Therapeutic Nanoparticle Conjugates
Next Article in Special Issue
Pharmacometabolomics by NMR in Oncology: A Systematic Review
Previous Article in Journal
Network Pharmacology of Ginseng (Part II): The Differential Effects of Red Ginseng and Ginsenoside Rg5 in Cancer and Heart Diseases as Determined by Transcriptomics
Previous Article in Special Issue
Semaphorin 3A—Glycosaminoglycans Interaction as Therapeutic Target for Axonal Regeneration
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Zinc Binding to NAP-Type Neuroprotective Peptides: Nuclear Magnetic Resonance Studies and Molecular Modeling

by
Ancuta-Veronica Lupaescu
1,
Cosmin Stefan Mocanu
2,
Gabi Drochioiu
2 and
Catalina-Ionica Ciobanu
3,*
1
Integrated Center for Research, Development and Innovation in Advanced Materials, Nanotechnologies and Distributed Systems for Fabrication and Control (MANSiD), Stefan cel Mare University of Suceava, 720229 Suceava, Romania
2
Faculty of Chemistry, Alexandru Ioan Cuza University, 11 Carol I, 700506 Iasi, Romania
3
CERNESIM Centre, Institute of Interdisciplinary Research, Alexandru Ioan Cuza University of Iasi, 11 Carol I, 700506 Iasi, Romania
*
Author to whom correspondence should be addressed.
Pharmaceuticals 2021, 14(10), 1011; https://doi.org/10.3390/ph14101011
Submission received: 31 July 2021 / Revised: 25 September 2021 / Accepted: 28 September 2021 / Published: 1 October 2021
(This article belongs to the Special Issue High Field NMR and Ultra-High Field NMR in Medicinal Chemistry)

Abstract

:
Aggregation of amyloid-β peptides (Aβ) is a hallmark of Alzheimer’s disease (AD), which is affecting an increasing number of people. Hence, there is an urgent need to develop new pharmaceutical treatments which could be used to prevent the AD symptomatology. Activity-dependent neuroprotective protein (ADNP) was found to be deficient in AD, whereas NAP, an 8-amino-acid peptide (1NAPVSIPQ8) derived from ADNP, was shown to enhance cognitive function. The higher tendency of zinc ion to induce Aβ aggregation and formation of amorphous aggregates is also well-known in the scientific literature. Although zinc binding to Aβ peptides was extensively investigated, there is a shortage of knowledge regarding the relationship between NAP peptide and zinc ions. Therefore, here, we investigated the binding of zinc ions to the native NAP peptide and its analog obtained by replacing the serine residue in the NAP sequence with tyrosine (1NAPVYIPQ8) at various molar ratios and pH values by mass spectrometry (MS) and nuclear magnetic resonancespectroscopy (NMR). Matrix-assisted laser desorption/ionization time-of-flight (MALDI ToF) mass spectrometry confirmed the binding of zinc ions to NAP peptides, while the chemical shift of Asp1, observed in 1H-NMR spectra, provided direct evidence for the coordinating role of zinc in the N-terminal region. In addition, molecular modeling has also contributed largely to our understanding of Zn binding to NAP peptides.

Graphical Abstract

1. Introduction

Amyloid-β peptides (Aβ) are protein fragments obtained following the abnormal cleavage assisted by β- and γ-secretase of β-amyloid precursor protein (APP), responsible for the proteinopathy in Alzheimer’s disease (AD) [1,2]. AD is characterized by extensive atrophy of the brain caused by a series of neuropathological changes, such as neuronal loss, the formation of amyloid plaques, development of neurofibrillary tangles, and synaptic loss [3,4]. The misfolding process associated with Aβ aggregation, which is a cause of AD, is largely influenced by metal ions such as Al, Cu, Fe, and Zn that are found in the proximity of senile plaques [5,6].
Zinc is an essential element involved in brain function. This metal can be found mostly tightly bound to zinc-dependent enzymes and other proteins. However, approximately 10% of zinc is not associated with proteins or amino acid ligands. Under pathological conditions, the “free” or “chelatable” zinc accumulates in neurons in neurotoxic amounts. For instance, during intense synaptic activity, zinc accumulates in the synaptic cleft and enters postsynaptic spines, sharing entry routes with calcium [7].
Zinc ions play a critical role in Aβ aggregation according to recent publications and, therefore, become of particular interest in the study of AD [8,9]. Moreover, zinc is the most prevalent transition metal in the brain and performs several functions in physiological processes in both normal and pathological conditions [10,11,12]. Thus, it was observed that high concentrations of zinc coincide with the regions where amyloid plaques form according to the post-mortem investigation of AD patients [13]. Therefore, the induction of Aβ oligomerization is considered as the main contribution of zinc ions to the pathogenesis of AD [14,15].
Furthermore, different fragments derived from activity-dependent neuroprotective protein (ADNP) have been studied as potential anti-AD drugs [15,16]. For instance, the neuroprotective 1NAPVSIPQ8 peptide (NAP) contains two prolines, which confer β conformation. The potential activity of this peptide consist of the inhibition of Aβ(1–40) and Aβ(25–35) fibrillogenesis, protection against Aβ-induced neurotoxicity, and prevention of AD-associated microtubule disruption [17].
The NAP peptide plays an important role in immune regulation, and it was well recognized as a neuroprotective moiety of ADNP [18,19]. Scanning of peptide activity highlighted the neuroprotective role of octapeptide NAP against the toxicity associated with amyloid-β peptide and zinc overload, which causes the breakdown of microtubules in astrocytes and neurons [20,21]. In addition, recent studies showed that the treatment of astrocytes with ZnCl2, which results in microtubule disassembly and cell death, was protected by the NAP peptide [20].
While the relationship between heavy metal ions, especially zinc ions, and amyloid-β peptides was recently investigated, little is known about the metal ion binding to neuroprotective molecules, such as NAP-derivate peptides. Therefore, here, we investigated the possible in vitro interaction of zinc ions with two NAP analogs.
Our previous research focused on the design and characterization of NAP-like peptides derived from the native one in order to investigate their behavior in the presence or absence of various metal ions [22,23,24]. Thus, several animal reports investigated the evaluation of neuroprotective strategies. For example, oral administration of Tyr-containing dipeptides, such as Ile-Tyr and Ser-Tyr, in a mouse model influenced catecholamine synthesis and transmission in the brain [25]. Another study suggested that Tyr-Pro dipeptide, capable of crossing the blood–cerebrospinal fluid barrier, can improve impaired cognitive deficits in both working and long-term memories in Aβ-injected AD model mice [26]. In a similar study, a Tyr based Aβ electrochemistry analysis was employed to study the effect of metal ions on Aβ aggregation kinetics. Using the direct oxidation signal of Tyr residue, the authors were able to monitor the metal-induced conformational changes on amyloid-beta peptides in vitro. Thus, the impact of Tyr residue oxidation was correlated to the formation of Zn(II)-Aβ (1–16) complex within a wide pH range from 6 to 9 [27].
The present report focuses mainly on the role of hydroxyl groups, alcoholic and phenolic, in the chelation of Zn2+ metal ions. Phenols are known to possess a higher acidic character than alcohols because the negative charge in the phenoxide ion is not localized on the oxygen atom, as it is in an alkoxide ion, but shared by a number of carbon atoms in the benzene ring. Therefore, phenolic hydroxyl of tyrosine is more capable of losing protons and coordinating to the metal ion through oxygen donors than aliphatic hydroxyl of serine [28]. In addition to influencing the binding of metal ions, the presence of an aromatic structure leads to peptide structural rearrangements that may enhance or reduce the peptide neuroprotection activity. Furthermore, according to recent publications, it has been observed that Tyr plays an important role in the complexation and structural modifications of peptides [29]. Oxidative stress is another process involved in the AD pathology [30]. Hence, the Tyr choice was made for this reason as well, since it allowed us to follow the influence on Zn complex formation. In addition, Tyr can be easily oxidized to dihydroxyphenylalanine (DOPA), a precursor used in the treatment of Parkinson’s disease [31]. Thus, following the oxidation process, a new hydroxyl group liable to intervene in metal complexation is introduced. Additionally, when the oxidation process is prolonged, the final reaction product becomes dopaquinone which can easily form metal cluster compounds. Therefore, the complexation process compared to the initial one becomes highly accentuated, with direct consequences in the peptide conformation.
Moreover, it has recently been shown that increasing the hydrophobicity of NAP-derivate peptides enhances the therapeutic effect on microtubule stabilization in AD [32]. Therefore, the addition of a new phenyl group automatically increases the hydrophobicity of the peptide, and therefore, the NAPY (H2N-1NAPVYIPQ8-COOH)-microtubule interaction becomes much more stable due to the hydrophobic system. Additionally, the presence of Tyr residue gives the molecules a much more pronounced rigidity. According to previous studies, replacement of histidine residues from the amyloid beta 1–16 peptide sequence with serine has been shown to increase flexibility. Thus, in aluminum complexation of serine-modified peptides, the aggregation process accelerates with the formation of novel fibrils [33]. Therefore, following this process, the effectiveness of NAP therapeutic treatment decreases drastically.
The molecular mechanisms of neurodegeneration proved to be more complex than a simple association of metal ions with peptides. It also involved the presence of multi-component systems and their time behavior under different pH, ionic strength, and concentration conditions. Consequently, we are searching for a new potential drug for the neurodegenerative disorders due to fibril formation in Aβ plaques in order to address either the prevention or improvement of AD pathophysiology.
In this regard, two synthesized peptides, native NAP fragment and tyrosine-modified peptide NAPY (1NAPVYPQ8), were incubated at different molar ratios with zinc ions. Their behavior was investigated by NMR spectroscopy and mass spectrometry (MS). Furthermore, the effects of pH on the peptide–metal interaction were also highlighted by NMR spectroscopy, suggesting a higher affinity toward metal ions at basic pH. Moreover, the theoretical calculations were used to confirm the experimental approaches.

2. Results

2.1. Mass Spectrometric Studies

In this study, we examined the interaction of zinc ions (Zn2+) with the native NAP peptide and its tyrosine analog NAPY obtained by substituting the fifth amino acid (serine) with tyrosine, in order to investigate the structural changes regarding Zn binding. Thus, the structure of the obtained compound was confirmed by Matrix-assisted laser desorption/ionization time-of-flight (MALDI ToF) mass spectrometry. The detailed characterization of synthesized peptides was previously described in our research [34]. A matrix-assisted laser desorption/ionization time-of-flight mass spectrometer was used to examine the affinity of NAP-like peptides for Zn2+ ions (MALDI-ToF MS) [35]. Consequently, Figure 1 and Figure 2 show the NAP-derived peptides positive-ion MALDI-ToF mass spectra before and after adding metal ions (Zn2+). The theoretical data from the ChemCalc software was well correlated with the m/z values found experimentally by the MALDI ToF mass spectrometry (Table 1).
Figure 1 shows the MS spectra obtained for the NAP peptide in the presence and absence of zinc ions. The addition of zinc ions in the sample generated an intense signal at m/z 887.6 characteristic of Zn2+-peptide complexes [M+Zn-H]+. The intense molecular peak observed at m/z 809.6 in both spectra was assigned to the deamidated peptide [M-16+H]+ ion. The structure rearrangement that leads to peptide deamidation was favored by the C-terminus glutamine residue. This process takes place during peptide ionization and thus can be considered as a post-source decay mechanism. The most intense signal observed at m/z 847.6, in the spectrum of native octapeptide generated in the absence of metal ions, was assigned to the sodium adducts ions [M+Na]+. Other signals observed in the spectra at m/z 863.6, m/z 869.6, and m/z 885.6 were attributed to the adduct with potassium [M+K]+ and two ions of sodium [M+2Na-H]+ as well as to [M+K+Na-H]+ ions.
In the case of the tyrosine-modified peptide, the MS spectra obtained (Figure 2) showed an intense signal at m/z 885.7, characteristic of the molecular [M-16+H]+ ion. As a result, the presence of Zn2+ did not enhance peptide deamination, which was shown in both spectra of peptides with and without metal ions. In addition to the signals corresponding to the deaminated ion, the NAPY spectra clearly showed typical signals of the molecular ion [M+H]+ at m/z 901.7. The signals obtained by the addition of +22 Da and +38 Da were attributed to the adducts of sodium, such as [M+Na]+, [M+2Na-H]+ and potassium, [M+K]+, and [M+K+Na-H]+. For the sample incubated with zinc ions, the obtained spectrum showed a distinctive signal at m/z 963.6, characteristic of the [M+Zn-H]+ ions generated by the Zn2+-peptide complex.

2.2. NMR Study

The synthesized peptide NAP (H2N-1NAPVSIPQ8-COOH) and NAPY (H2N-1NAPVYIPQ8-COOH) were characterized by NMR spectroscopy in order to perform atom-specific chemical shift assignments and structural information regarding specific metal binding sites. The NMR technique is the method of choice for studying the interaction of metal ions with peptides [36,37,38]. The perturbations generated in the proton NMR spectra, upon Zn2+ binding to NAP-like peptides, provide insights about the interaction established between metal ions and various amino acid residues. Furthermore, a comparison between the displayed chemical shifts and the assigned values identified in various peptide studies offers relevant information regarding the ligands involved in coordinating zinc ions.
In this regard, one-dimensional proton spectra of the native NAP peptide and tyrosine-modified NAPY peptide dissolved in deuterium oxide were recorded. All the signals measured were assigned to the corresponding protons according to the observed chemical shift, homonuclear correlation (1H–1H correlation), and heteronuclear direct correlation (1H–13C correlation). Additionally, the 13C spectrum presented a suitable signal dispersion, which allowed the assignment of carbon atoms that are directly bound to the protons. Although some signals were affected by an exchange process with the deuterium atoms of the solvent and are missing in the spectra, such as amide, amines, hydroxyl, and carboxylic protons, the NMR experiments confirmed the amino acid sequence. The observed protons are defined in the following structure (Figure 3).
The 1H-NMR spectrum of NAP peptide (Figure 4) showed signals in the aliphatic region between 0.2 and 4.7 ppm. The most deshielded signals belonged to the eight alpha protons Hα defined in the peptide structure (amidic protons are not visible in the spectrum). The characteristic signals of Hβ protons of serine (2H) and Hδ pyrrolidine protons of proline (4H) were observed in the range of 3.6–3.9 ppm. Due to the spatial arrangement of the atoms in the molecule, the Hβ protons of Asn generated two signals at 2.98 ppm (1H) and 2.85 ppm (1H). Other Hβ protons corresponding to the Val (1H) and Ile (1H) along with Hγ protons (4H) from the pyrrolidine ring generated signals between 2.1 ppm and 1.8 ppm. In the case of glutamine, the Hβ protons were observed at 2.20 and 1.98 ppm (2H), while Hγ protons generated signals identified at 2.41 ppm (2H). The intense signal observed as a doublet at 1.36 ppm was easily attributed to the methyl protons (3H) of alanine, and the signals found at 1.47 ppm and 1.15 ppm corresponded to the Hγ isoleucine protons (2H). Additionally, the most shielded protons were the methyl group of valine (Hγ) and isoleucine (Hγ, Hδ) that generated signals in the range of 1.5–0.8 ppm. Furthermore, the intense signal observed at 2.06 ppm was due to acetic acid traces from the sample reagent used in the peptide elution step.
In the case of NAPY, the proton NMR spectrum (Figure 5) contains structural similarities compared to the previously investigated one (Figure 4). However, small differences caused by the replacement of the fifth amino acid were observed in the 1H-NMR spectra. For example, signals corresponding to the protons attached to the phenolic nucleus of tyrosine were observed in the aromatic region at 7.07 ppm (2H) and 6.78 ppm (2H). Additionally, the signal given by the tyrosine alpha proton was identified at 4.59 ppm (1H), while the Hβ proton (2H) overlapped the signals generated by the beta protons of asparagine in the region of 3.03–2.77 ppm. Another notable change was the rearrangement of signals generated by gamma protons found in the structure of the Ile and Val amino acids. Thus, the observed aspects highlight the structural changes suffered by the peptide following the replacement of serine with tyrosine and the influence that the new amino acid has on the peptide structural rearrangement.
After characterization of pure peptides, the metal ion binding capacity to NAP octapeptides was investigated by NMR titration. Thus, solutions of peptide and Zn(II) ions at a molar ratio equal to 1:0.1, 1:0.5, 1:1, 1:1.5, 1:2, and 1:10 were prepared by gradually adding 5 µL of metal ion stock solutions (60, 480, 2400, or 9600 mM in H2O) to the peptide solution (10 mM in D2O). The mixture series of metal ions and peptides were transferred into standard 5 mm NMR tubes and all experiments recorded at 25 °C.
The figures which show the 1D proton NMR spectrum of NAP (Figure 6) and NAPY (Figure 7) peptides titrated with increasing amounts of Zn2+ at pH 5 are presented below. The addition of Zn2+ ions did not cause a perturbation of 1H-NMR resonances at the natural pH of the solution upon adding ZnSO4 to peptides. Thus, no chemical shifts were observed upon metal increasing, suggesting the presence of only intermediate exchange regime.
A similar pattern was observed in Figure 7. According to these results, the 1H-NMR spectrum incorporating Zn2+ titrations of tyrosine modified NAPY peptide did not cause any perturbation to the system.
On the other hand, increasing the pH value of solution caused perturbation of NMR resonances (Figure 8 and Figure 9). Thus, at pH 9.5 and a 1:10 molar ratio of peptide:Zn2+ titration, new exchange cross peaks and shifts appeared in 1H- NMR spectra. The pH was adjusted by adding crystals of potassium carbonate K2CO3 and was measured before recording each spectrum. Both peptides presented a significant chemical shift deviation for the side-chain peak of asparagine residues (Hβ protons). Thus, upon addition of zinc ions, Hα protons from Asn and Gln residues generated a noticeable chemical shift.
All changes caused by affected protons appear as upfield signals in the spectrum. Additionally, both peptides provided a direct indication regarding the involvement of the N-terminus and C-terminus in zinc binding. Therefore, the nitrogen atoms of the amino group oxygen from the carbonyl group can be envisaged as possible coordination sites for metal ions. The chemical shift changes are presented in Table 2 for NAP and Table 3 for the NAPY-modified peptide.
In the case of the NAP peptide, in addition to the mentioned modifications, the Hα protons from alanine generated chemical shift changes at pH 9.5. Thus, considering the position of Ala in the peptide structure (near the N-terminal), it can be concluded that the changes are not purely caused by conformational alterations, but by a direct (de)shielding effect of the positively charged zinc ion found nearby. Moreover, there were no changes in the chemical displacement of protons in the range of 0.8–1.6 ppm (Hβ in Ala, Hγ in Ile and Val, Hδ in Ile) and 3.6–4.7 ppm.
An evident change of the peptides NAP and NAPY conformation in solution was visibly recorded in the 1H-NMR spectra at pH 7.5. Chemical shifts appear in protons from asparagine and glutamine residues. In experiments conducted with zinc ions and peptides at pH 7.5, new modifications were observed in the alpha protons from the asparagine amino acid.
At pH 9.5, the tyrosine beta resonances from the NAPY peptide were unchanged upon zinc addition, while the aromatic protons displayed light chemical shift changes (Figure 9) upon shielding effect. Thus, Zn2+ slightly affects the aromatic protons from Tyr5, appearing to lower the chemical shift. These observations highlight the important role played by Y5 hydroxyl in binding zinc ions. Other chemical shift changes observed in the 1H-NMR spectra of the NAPY peptide were similar to that suffered by the NAP peptide.

2.3. Molecular Modeling

Oftentimes, the atomic-level details of the structures of pharmaceutically relevant receptors are not available. Three-dimensional alignment of potential ligands can be utilized to derive structural requirements for biological activity in such instances. Another strategy is pharmacophore elucidation, in which several ligands are aligned, and a small collection of essential molecular features required for biological activity is derived from the alignment [39]. Additionally, another method is to search a database for ligands with one or more conformations which superpose well with a query molecule using 3D molecular alignment [40]. Methodologies based on 3D alignment for discovering biologically active ligands employ the qualitative premise that if two ligands have similar biological activity and bind in similar modes, then their bound conformations match well, and inferences about the nature of the receptor may be drawn.
In order to develop and confirm the interaction of zinc ions with the sequence of the NAP and NAPY peptide, we resorted to a series of theoretical approaches. First, we applied Flexible Alignment using Molecular Operating Environment 2016.02 (MOE) in order to find the coordination pattern of Zn2+ for the NAP and NAPY peptide. According to the submitted simulation, the zinc ion binds to the NAP through Asn1 or Gln8 at the C-terminal or N-terminal of the peptide (Figure 10A,B). These bindings were also noticed in the NMR spectra, confirming the experimental procedure. However, these representations suggest that the interaction between Zn2+ and the NAP sequence is not well defined, as in the case of the NAPY-Zn2+ complex (Figure 10C).
Furthermore, according to Figure 10C, the zinc ions bind strongly with the NAPY peptide structure through a cluster formed by Asn1, Ala2, Pro3, Val4, and Tyr5. Moreover, besides the oxygen atom, Tyr5 seems to contribute even with its backbone through the π-cation interaction, acting as a stabilizer to the Zn2+-NAPY complex. Hence, in this case, the binding potential between the zinc ion and the tyrosine analogue becomes stronger.
In addition, the NAP/NAPY-Zn2+ complex was subjected to molecular dynamics (MD) in order to investigate the stability and the binding interaction between zinc ions and the corresponding amino acid under thermic conditions with the same software. The purpose of this approach was to study the influence of body temperature (37 °C) regarding the distances between Zn2+ and the atom involved in the coordination process. Hence, a series of parameters described in the Section 4.2.3 was applied using three different methods regarding MD algorithms: NPA (Nosé-Poincaré-Andersen), NHA (Nosé-Hoover-Andersen), and BER (Berendsen equations) (Figure S1, Supplementary Materials).
The coordination binding, as well as the structure of the peptide under 25 °C (298.15 K) for complex stabilization and 37 °C (310.15 K) for the body temperature effect, is represented in Figure S2 in the case of NAP-Zn2+-Asn1, in Figure S3 for NAP-Zn2+-Gln8, and in Figure S4 for the case of NAPY and zinc ion binding (all in Supplementary Materials).
According to the MD runnings, there are few modifications regarding the binding pattern of NAP and NAPY peptides under a thermic condition at 25 degrees (room temperature) and 37 degrees (body temperature). In the case of the NAPY peptide, under a thermic condition, it was noticed that the Tyr plays an important role in metal complexation. According to Figure S3, the oxygen atom of Tyr influences the stability of Zn2+ complexation. At 37 degrees, it seems that the coordination binding is stronger than at 25 °C. This suggests that in the human body, the NAPY-Zn2+ complex becomes more stable. Moreover, these data also suggest that the peptide has a strong stability which did not influence the binding cluster with zinc ions.
Consequently, significant changes were also analyzed regarding binding sites and the distance between zinc ions and the corresponding atom from the amino acid involved in the binding process (Table S1, Supplementary Materials).
According to these data, relevant modifications were observed regarding the distance between zinc ions and the amino acid from the NAP/NAPY sequence. As can be noticed, under a thermic condition compared with the static momentum, the bond between Zn2+ ions and the atom involved in the coordination process becomes stronger at 37 degrees (310.15 K). This suggests that the temperature plays an important role in the metal ion binding in our case. Moreover, it also can be observed that in the coordination process, the oxygen in the sp2 state strengthens the bond more than in the sp3 state. Therefore, according to these data, the more oxygen the peptide has in the sp2 state, the stronger the bond between the metal and the peptide becomes regarding the temperature.
In addition, it was noticed that in the case of NAP-Zn2+-Asn1, there was a slight dissociation of the zinc ion when it was applied to a higher temperature. Furthermore, in the NAP-Zn2+Gln8 complex, there were no major differences in the distance between the zinc ion and the corresponding atom when the temperature increased from 25 to 37 degrees—just a minor drop of binding energy. This suggests that in this situation, at body temperature, the bonds in the coordination process become stronger, even if the distances remain almost the same. In the case of tyrosine-modified peptides, the bond between Osp2(Asn1) and zinc increases while the distance that implies alanine and proline becomes shorter, and the bond energies are constant.
Thus, according to these data, at 37 degrees (310.15 K), there is a higher probability to form a NAP-Zn2+Gln8 complex than NAP-Zn2+-Asn1. Moreover, the addition of tyrosine increases the stability of the zinc ions complex though alanine and proline residues, as well as the oxygen and the π-cations of this amino acid (stability is confirmed in Figure 10 and Figure S4).
Hence, this approach could have an important implication in drug design with direct consequences in peptide stability and blood–brain barrier permeability. Furthermore, these observations were well correlated with NMR investigations, suggesting that the data obtained in our study were accurate according to the experimental and theoretical approaches. Therefore, our experimental and theoretical design helped us to reveal the underlying structural modification under thermic and metal conditions.

3. Discussion

The antioxidant octapeptide NAP acts as a therapeutic agent at femtomolar concentrations, enhancing neuronal development, learning, and memory [41]. This peptide has been applied in different clinical studies, showing cognitive and functional protection in research trials involving patients with AD-specific cognitive impairment [42]. Moreover, it has been shown to inhibit Aβ aggregation, while the neuroprotective mechanism exerted by NAP assumes the folding modulation of toxic proteins in the extracellular environment [18].
The etiology of this pathology is related to the formation of fibrils and protein aggregates associated with metal binding, which leads to conformational changes in the morphology of the entire brain [43].
Although the homeostasis of metal ions, such as iron, copper, zinc, and calcium, in the brain is crucial for maintaining normal physiological functions, the metabolism alteration of these ions in the brain is closely linked to the onset and progression of AD [44].
Various studies showed that zinc is involved in several AD-associated processes, such as neurotransmission, nonamyloidogenic APP processing, enzymatic degradation of the Aβ peptide, the Aβ–DNA interaction, and inducing Aβ oligomerization [12,14]. Furthermore, the zinc ion concentration in amyloid plaques was demonstrated to be high, reaching 1 mM [45]. Thus, chelating agents capable of linking together metal ions to form complex ring-like structures are of particular interest for reducing metal toxicity.
The zinc binding site to the amyloid peptide is located at the disordered N-termini reach in aspartate (Asp1, Asp7), glutamate (Glu3, Glu11), and histidine residues (His6, His13, and His14) [46]. NMR studies performed in aqueous solution at pH 6.5 and 7.4 using the 1–16 amyloid fragment showed that Zn2+ was bound to the three histidine residues and to the two carboxylic groups of Glu11 [47]. Another NMR study, performed on Aβ 1–28, reveled a penta-coordination of Zn2+ by His6, Glu11, His13, His14, and Asp1 at the N-terminus [48].
At low pH, where the majority of the amino groups are protonated, the affinity towards the metal ions was weak, and the stability of the complex was plain. As the pH rises, the affinity and stability of the peptide–metal complexes increase. This suggests that at low pH, almost no metal ions are complexed. However, once the pH increases, the peptide–metal affinity reaches a maximal value. This process indicates the possibility of a recovery process by readjusting the pH [49].
There are five potentially active amino acid residues in human Aβ peptides: three His at positions 6, 13, and 14; one Met at position 35; and one Tyr at position 10 [50]. The results obtained after increasing the pH value showed a higher affinity toward metal ions. Moreover, similar studies on zinc chelation agents showed that water-soluble macroligands had the best affinity at basic pH (pH 9), while acidic pH promotes competition between H+ and divalent metal ions [49].
In addition, the structural changes of the Zn-NAP/NAPY complex could also be attributed to the solvation of the system. Ben-Naim highlights this process by establishing that the effect of a pair of functional groups (hydroxyl or amino), either intermolecular or intramolecular hydrogen bonds, depends on the total standard free energy of the peptide folding [51]. Moreover, a similar study reached the same conclusions. The relevance of the solute–solvent coupling, as well as the differences in behavior of hydrophilic and hydrophobic moieties, depends on the amino acid residues [52].
Additionally, these interferences could become significant in our case because they might interfere with the aggregation process in AD. A model has been used to investigate peptide aggregation at hydrophobic–hydrophilic interfaces and compare it to aggregation in a clear water solvent [53]. It was also relevant in a study which investigated the size and fractal character of the monomers, oligomers, and fibrils of the various amyloid-beta peptides using the small angle neutron scattering method (SANS) [54]. Using molecular dynamics simulation, new insight regarding the influence of water molecules on the folding landscape of amyloid peptides was obtained. Thus, in the presence of water, hydrophobic regions were shown to accelerate fibril construction while hydrophilic sequences reduced the growth of fibrils by water-stabilizing metastable intermediates [55]. The relation between hydrogen bonding and the hydrophobic sequence leading to protein folding is predominantly determined by the backbone intermolecular hydrogen bonding interactions [56]. Hydrogen bonds have significant roles in protein assemblies and in the organization of biomolecular protein–ligand complexes. All amino acids can generate hydrogen bonds and serve as both donors and acceptors. Additionally, the hydroxyl group of tyrosine can undergo hydrogen bonding with the backbone amide of nearby tyrosine residue, leading to structure stabilization [57]. Hydrogen bonding events can vary at different pH conditions. For example, collagen-like peptides exhibit a reduced intramolecular hydrogen bonding behavior at basic condition and a minimal event of H-bond between amino acid residues and water molecules at acid conditions [58]. Moreover, it was proved that in different protein aggregations, water performs an important role [55]. Therefore, the Zn-NAP/NAPY complexations and the obtained results could be attributed to the solvent, due to the hydrophobic and structure stabilization.

4. Materials and Methods

4.1. Peptide Synthesis

The octapeptides NAP (NH2-1NAPVSIPQ8-COOH) and NAPY (NH2-1NAPVYIPQ8-COOH) were synthesized using the solid phase Fmoc/tBu solid phase synthesis (SPPS). The peptide synthesis was carried out manually in a linear fashion on a Fmoc-Gln(Trt)-Wang Resin and was described in our previous work [28]. Summarily, the manual peptides synthesis was performed in a fritted plastic reactor attached to a vacuum pump in order to remove the washing solutions. This synthesis was performed in dimethylformamide (DMF) solvent. Protecting groups such as Fmoc or tert-butyl, used to protect the amino acids at the N-terminus or the side chain, were removed with 20% piperidine in DMF. Final cleavage of the peptides from the resin was performed with a TFA: TIS: H2O solution, in a ratio of 95:2.5:2.5 (v/v/v). The precipitation of peptides was performed using diethyl ether (−20 °C). The resulting fractions, eluted with 5% acetic acid solution in ultrapure water, were lyophilized. The ZnSO4·7H20 salt used as metal ion source was procured from Sigma Aldrich (Germany). All solutions were prepared using deionized water (18.2 MΩ·cm) from a Milli-Q system (Millipore, Bedford, MA, USA).

4.2. Instruments

4.2.1. Mass Spectrometry

The matrix-assisted laser desorption ionization time-of-flight (MALDI-ToF) MS analysis was performed on a Bruker Ultraflex MALDI ToF/ToF mass spectrometer operated in positive reflectron mode and equipped with a pulsed nitrogen UV laser (λmax 337 nm). For MS studies, metal complexes of NAP-like peptides obtained by mixing the corresponding octapeptides with metal solution at physiological pH (pH 7.4) at a 1:10 peptide-to-metal ratio were used. MilliQ-grade water was used to make the stock solutions of peptides (8 mM) and metal ions (80 mM). The samples were loaded onto a 384-spot target plate of the MALDI-ToF device using the dried-droplet method. On the target, the sample and matrix solution (α-cyano-4-hydroxycinnamic acid dissolved in a solution comprising 2:1 ACN: 0.1% TFA in MilliQ) were mixed and allowed to dry. After co-crystallization, the metal plate was introduced into the mass spectrometer and subjected to short laser pulses. The acquired signals were expressed as atomic mass units (Da). Finally, the obtained spectra were processed using Bruker’s FlexAnalysis 3.4 software.

4.2.2. NMR Spectroscopy

Proton NMR experiments were recorded using a Bruker Avance III, 500 MHz frequency spectrometer, equipped with a 5 mm Pabbo detection probe and operating at 500 MHz for 1H nucleus and 125 MHz for 13C nucleus. Chemical shifts are reported in δ units (ppm), relative to the internal standard (TSP, δ 0.0 ppm). The experiments were recorded without introduction of a pre-saturation pulse for solvent signal suppression. All spectra were recorded at 298 K with 12 mg/mL peptide dissolved in D2O with 5 mM TSP 3-(trimethylsilyl)propionic-2,2,3,3-d4 acid, sodium salt in D2O added, as the internal chemical shift calibration standard. Due to the dissolution of the peptide in deuterium oxide, in the 1H–NMR spectra, some signals are affected by a process of exchange with deuterium in the solvent and are missing in the spectra; for example, the exchange of amidic protons, amino protons from asparagine and glutamine residues, and hydroxyl hydrogens from serine and tyrosine residues. Based on the peak integral measured in the quantitative 1D 1H-NMR spectrum, the number of protons in the peptide was verified. NMR data were processed and analyzed using BrukerTopSpin3.2. software.

4.2.3. Data Analysis

In this study, we investigated the binding pattern changes of NAP and NAPY peptides under the influence of Zn2+ ions using different theoretical approaches. For this purpose, the Molecular Operating Environment 2016.02 modeling software from Chemical Computing Group (MOE) was used [59]. First, the peptides (1NAPVSIPQ8—NAP and 1NAPVYIPQ8—NAPY) were built and subjected to energy minimization in order to calculate atomic coordinates which are local minima of a potential energy function. The Amber10:EHT forcefield was used with van der Waals and electrostatic forces considered. Furthermore, tethers, distances, angles, and dihedrals were restrained. All atoms were calculated with partial charges and cutoff of (8,10). The Hamiltonian was PM3-RHF with al spin state. The gradient was set at 0.05 RMS kcal/mol/A2, and regarding quantum mechanics, the basis set was AM1, and the multiplicity singlet. These restraints were removed upon termination of the energy minimization [60,61,62].
Subsequently, the NAP/NAPY-Zn2+ system was subjected to Flexible Alignment, used in general for flexibly aligning small molecules. This algorithm takes a collection of small molecules with 3D coordinates as input and generates a set of alignments. Our purpose was to set a coordination binding to the assigned peptide. Each alignment obtained a score that measured the quality of the alignment in terms of both overlap and internal strain of molecular features. The probability density functions used in our case were Gaussians. The iteration limit was set at 1000 and a failure limit at 50. The energy cutoff was set at 15. Moreover, forcefield charges were calculated prior to search and the stochastic conformations [63,64,65].
The ligand interactions were used to diagrammatically depict the active site of the NAP/NAPY-Zn2+ complex. These interactions were obtained according to the method described in the literature [66]. In our case, the hydrogen bonds, metal ligation, and nonbonded residues were taken into account. Regarding this protocol, a series of different steps were followed: the chosen ligand was set as the central focus of the diagram, the ligand’s coordinates were determined, the residues were placed to form strong hydrogen bonds with the ligand, the residues that are close to the ligand with weaker interactions were placed, and proximity contour and solvent accessibility were calculated [66,67,68].
Consequently, the NAP/NAPY-Zn2+ system was further subjected to molecular dynamics (MD) in order to analyze the influence of thermic condition in metal–peptide complexation. Three independent MD runnings were provided using NPA (Nosé-Poincaré-Andersen), NHA (Nosé-Hoover-Andersen), and BER (Berendsen equations) algorithms. For each MD procedure, a run time of 1000 ps (picoseconds, or 1 ns) was set followed by an equilibrium stage of 100 ps at 298.15 K and 310.15 K—a total run time of 2.2 ns for each algorithm. The same forcefield was used as in the minimization of the system energy, Amber10:EHT with a cutoff of (8,10). The constrains were set at light bonds, and tethers, distances, angles, and dihedrals were restrained. The checkpoint of the MD run was set at 250 ps, with a time step of 0.002 ps [69,70,71].

5. Conclusions

In recent years, molecular studies have become a useful tool in medicinal chemistry, furnishing relevant information concerning binding sites and binding energy. According to our study, the data obtained showed that zinc ions bind not only to amyloid-β peptides but also to the neuroprotective molecules, such as the native NAP peptide and its analog obtained by replacing the serine residue from the NAP sequence with tyrosine (1NAPVYIPQ8). Furthermore, it was observed that at different molar ratios and pH values, the interaction between Zn and our investigated peptides changing and the N-terminal sequence of peptides was involved in ion binding. Moreover, the inquiry of mass spectrometry and NMR spectroscopy confirmed the binding of the zinc ion to NAP peptides, while the chemical shift of Asp1, observed in 1H-NMR spectra, provided direct evidence for the coordinating role of zinc in the N-terminal region. In addition, the molecular modeling showed a perspective regarding their interactions, confirming the experimental data obtained during mass spectrometry and NMR investigations.
Hence, our study offers a unique insight into metal ion interactions with the NAP and NAPY peptide by analyzing the influence of Tyr swapping with Ser, more specifically the influence of the benzene nucleus. Our approaches demonstrated that the addition of a phenyl group increases the binding potential of Zn ions to the structure of peptides. Therefore, the investigation of such short NAP peptide fragments and analogues provides improvements in the structure and biocompatibility of peptide-based drugs enhancing the therapeutic pathways of neuropathologies and can also help to assess chemical phenomena in peptide–metal ion complexes to understand the modifications and structural effects of peptides with therapeutic activity. All these bring important information for understanding AD. Thus, further studies aiming to obtain new modified peptides starting from the active site of the NAP peptide, which is represented by the SIP center, are needed in order to better understand the reactivity of the interactions between metal ions and peptides.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ph14101011/s1. Figure S1 —The structures at 310.15 K following molecular dynamics using 3 algorithms of NAP/NAPY-Zn2+ complexes; Figure S2—Molecular dynamics of NAP-Zn2+-Asn1 complex; Figure S3—Molecular dynamics of NAP-Zn2+-Gln8 complex; Figure S4—Molecular dynamics of NAPY-Zn2+; Table S1—The molecular parameters of NAP/NAPY-Zn2+ complex according to the MD simulation.

Author Contributions

Conceptualization, writing—original draft preparation, C.-I.C. and A.-V.L.; in silico studies, C.S.M.; writing—review and editing, G.D. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by a grant of the Ministry of Research, Innovation, and Digitization, CNCS/CCCDI-UEFISCDI, project number: CNFIS-FDI-2021-0501.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in article and Supplementary Materials.

Acknowledgments

The authors are thankful to the CERNESIM Center, within the Institute for Interdisciplinary Research from Alexandru Ioan Cuza University of Iasi, for the infrastructure used in recording the NMR experiments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jamerlan, A.; An, S.S.A. The influence of Aβ-dependent and independent pathways on TDP-43 proteinopathy in Alzheimer’s disease: A possible connection to LATE-NC. Neurobiol. Aging 2020, 95, 161–167. [Google Scholar] [CrossRef]
  2. Xing, S.; Lu, Y. P3-199: Comprehensive understanding of tau exon 10 alternative splicing regulation. Alzheimer’s Dement. 2018, 14, 1142–1143. [Google Scholar] [CrossRef]
  3. Isik, A.T. Late onset Alzheimer’s disease in older people. Clin. Interv. Aging 2010, 5, 307–311. [Google Scholar] [CrossRef] [Green Version]
  4. Liu, P.P.; Xie, Y.; Meng, X.Y.; Kang, J.S. History and progress of hypotheses and clinical trials for Alzheimer’s disease. Signal Transd. Targ. Ther. 2019, 4, 1–22. [Google Scholar] [CrossRef]
  5. Viles, J.H. Metal ions and amyloid fiber formation in neurodegenerative diseases. Copper, zinc and iron in Alzheimer’s, Parkinson’s and prion diseases. Coord. Chem. Rev. 2012, 256, 2271–2284. [Google Scholar] [CrossRef]
  6. Ben-Shushan, S.; Miller, Y. Neuropeptides: Roles and activities as metal chelators in neurodegenerative diseases. J. Phys. Chem. B 2021, 125, 2796–2811. [Google Scholar] [CrossRef]
  7. Morris, D.R.; Levenson, C.W. Ion channels and zinc: Mechanisms of neurotoxicity and neurodegeneration. J. Toxicol. 2012, 2012, 785647. [Google Scholar] [CrossRef]
  8. Datki, Z.; Galik-Olah, Z.; Janosi-Mozes, E.; Szegedi, V.; Kalman, J.; Hunya, Á.G.; Fulop, L.; Tamano, H.; Takeda, A.; Adlard, P.A.; et al. Alzheimer risk factors age and female sex induce cortical Aβ aggregation by raising extracellular zinc. Mol. Psychiatr. 2020, 25, 2728–2741. [Google Scholar] [CrossRef] [PubMed]
  9. Ayton, S.; Lei, P.; Bush, A.I. Metallostasis in Alzheimer’s disease. Free Radic. Biol. Med. 2013, 62, 76–89. [Google Scholar] [CrossRef] [PubMed]
  10. Chasapis, C.T.; Ntoupa, P.S.A.; Spiliopoulou, C.A.; Stefanidou, M.E. Recent aspects of the effects of zinc on human health. Arch. Toxicol. 2020, 94, 1443–1460. [Google Scholar] [CrossRef]
  11. Olechnowicz, J.; Tinkov, A.; Skalny, A.; Suliburska, J. Zinc status is associated with inflammation, oxidative stress, lipid, and glucose metabolism. J. Physiol. Sci. 2018, 68, 19–31. [Google Scholar] [CrossRef] [Green Version]
  12. Kulikova, A.A.; Makarov, A.A.; Kozin, S.A. Roles of zinc ions and structural polymorphism of β-amyloid in the development of Alzheimer’s disease. Mol. Biol. 2015, 49, 217–230. [Google Scholar] [CrossRef]
  13. Eisele, Y.S.; Bolmont, T.; Heikenwalder, M.; Langer, F.; Jacobson, L.H.; Yan, Z.X.; Roth, K.; Aguzzi, A.; Staufenbiel, M.; Walker, L.C.; et al. Induction of cerebral β-amyloidosis: Intracerebral versus systemic Aβ inoculation. Proc. Natl. Acad. Sci. USA 2009, 106, 12926–12931. [Google Scholar] [CrossRef] [Green Version]
  14. Xie, Z.; Wu, H.; Zhao, J. Multifunctional roles of zinc in Alzheimer’s disease. Neurotoxicology 2020, 80, 112–123. [Google Scholar] [CrossRef]
  15. Isaev, N.K.; Stelmashook, E.V.; Genrikhs, E.E. Role of zinc and copper ions in the pathogenetic mechanisms of traumatic brain injury and Alzheimer’s disease. Rev. Neurosci. 2020, 31, 233–243. [Google Scholar] [CrossRef] [PubMed]
  16. Gozes, I. Activity-dependent neuroprotective protein (ADNP)/NAP (CP201): Autism, schizophrenia, and Alzheimer’s disease. In Neuroprotection in Autism, Schizophrenia and Alzheimer’s Disease; Academic Press: Cambridge, MA, USA, 2019; pp. 3–20. [Google Scholar] [CrossRef]
  17. Cacabelos, R.; Cacabelos, P.; Torrellas, C.; Tellado, I.; Carril, J.C. Pharmacogenomics of Alzheimer’s disease: Novel therapeutic strategies for drug development. In Pharmacogenomics in Drug Discovery and Development; Yan, Q., Ed.; Springer: New York, NY, USA, 2014; Volume 1175, pp. 323–556. [Google Scholar] [CrossRef]
  18. Heimesaat, M.M.; Mousavi, S.; Kløve, S.; Genger, C.; Weschka, D.; Giladi, E. Immune-modulatory properties of the octapeptide nap in campylobacter jejuni infected mice suffering from acute enterocolitis. Microorganisms 2020, 8, 802. [Google Scholar] [CrossRef]
  19. Matsuoka, Y.; Gray, A.J.; Hirata-Fukae, C.; Minami, S.S.; Waterhouse, E.G.; Mattson, M.P.; LaFerla, F.M.; Gozes, I.; Aisen, P.S. Intranasal NAP administration reduces accumulation of amyloid peptide and tau hyperphosphorylation in a transgenic mouse model of Alzheimer’s disease at early pathological stage. J. Mol. Neurosci. 2007, 31, 165–170. [Google Scholar] [CrossRef] [PubMed]
  20. Divinski, I.; Mittelman, L.; Gozes, I. A femtomolar acting octapeptide interacts with tubulin and protects astrocytes against zinc intoxication. J. Biol. Chem. 2004, 279, 28531–28538. [Google Scholar] [CrossRef] [Green Version]
  21. Kitai, S.T.; Shepard, P.D.; Callaway, J.C.; Scroggs, R. Afferent modulation of dopamine neuron firing patterns. Curr. Opin. Neurobiol. 1999, 9, 690–697. [Google Scholar] [CrossRef]
  22. Ciobanu, C.I.; Stefanescu, R.; Niculaua, M.; Teslaru, T.; Gradinaru, R.; Drochioiu, G. Mass spectrometric evidence for iron binding to the neuroprotective peptide NAP and its Cys5 mutant. Eur. J. Mass Spectrom. 2016, 22, 97–104. [Google Scholar] [CrossRef]
  23. Lupaescu, A.V.; Jureschi, M.; Ciobanu, C.I.; Ion, L.; Zbancioc, G.; Petre, B.A.; Drochioiu, G. FTIR and MS evidence for heavy metal binding to anti-amyloidal NAP-like peptides. Int. J. Pept. Res. Ther. 2019, 25, 303–309. [Google Scholar] [CrossRef]
  24. Lupaescu, A.V.; Sandu, I.; Petre, B.A.; Ion, L.; Ciobanu, C.I.; Drochioiu, G. Nap neuroprotective peptide and its analogs: Simultaneously copper and iron binding and reduction. Rev. Chim. 2019, 70, 1784–1790. [Google Scholar] [CrossRef]
  25. Moriyasu, K.; Ichinose, T.; Nakahata, A.; Tanaka, M.; Matsui, T.; Furuya, S. The dipeptides Ile-Tyr and Ser-Tyr exert distinct effects on catecholamine metabolism in the mouse brainstem. Int. J. Pep. 2016, 2016, 6020786. [Google Scholar] [CrossRef] [Green Version]
  26. Tanaka, M.; Kiyohara, H.; Yoshino, A.; Nakano, A.; Takata, F.; Dohgu, S.; Kataoka, Y.; Matsui, T. Brain-transportable soy dipeptide, Tyr-Pro, attenuates amyloid β peptide 25-35-induced memory impairment in mice. NPJ Sci. Food 2020, 4, 1–4. [Google Scholar] [CrossRef]
  27. Suprun, E.V.; Zaryanov, N.V.; Radko, S.P.; Kulikova, A.A.; Kozin, S.A.; Makarov, A.A.; Archakov, A.I.; Shumyantseva, V.V. Tyrosine based electrochemical analysis of amyloid-β fragment (1-16) binding to metal (II) ions. Electrochim. Acta 2015, 179, 93–99. [Google Scholar] [CrossRef]
  28. Yamauchi, O.; Odani, A.; Takani, M. Metal-amino acid chemistry. Weak interactions and related functions of side chain groups. J. Chem. Soc. Dalton Trans. 2002, 18, 3411–3421. [Google Scholar] [CrossRef]
  29. Ohata, J.; Miller, M.K.; Mountain, C.M.; Vohidov, F.; Ball, Z.T. A three-component organometallic tyrosine bioconjugation. Angew. Chem. Int. Ed. 2018, 57, 2827–2830. [Google Scholar] [CrossRef] [PubMed]
  30. Cassidy, L.; Fernandez, F.; Johnson, J.B.; Naiker, M.; Owoola, A.G.; Broszczak, D.A. Oxidative stress in alzheimer’s disease: A review on emergent natural polyphenolic therapeutics. Complement. Ther. Med. 2020, 49, 102294. [Google Scholar] [CrossRef]
  31. Mocanu, C.S.; Bocec, A.S.; Gradinaru, V.R.; Anton-Paduraru, D.T. A biochemical method for tyrosine determination in phenylketonuria using a colorimetric enzymatic approach. Rev. Chim. 2020, 71, 285–294. [Google Scholar] [CrossRef]
  32. Mocanu, C.S.; Drochioiu, G. The interaction of possible anti-AD ASA-NAP peptide conjugate with tubulin: A theoretical and experimental insight. Int. J. Pept. Res. Ther. 2021, 1–17. [Google Scholar] [CrossRef]
  33. Mocanu, C.S.; Jureschi, M.; Drochioiu, G. Aluminium binding to modified amyloid-β peptides: Implications for alzheimer’s disease. Molecules 2020, 25, 4536. [Google Scholar] [CrossRef]
  34. Lupaescu, A.V.; Ciobanu, C.I.; Humelnicu, I.; Petre, B.A.; Murariu, M.; Drochioiu, G. Design and synthesis of new anti-amyloid NAP-based/like peptides. Rev. Roum. Chim. 2019, 64, 535–546. [Google Scholar] [CrossRef]
  35. Ciobanu, C.I.; Lupaescu, A.V.; Drochioiu, G. Study of zinc binding to neuroprotective peptides. In Proceedings of the 19th International Multidisciplinary Scientific GeoConference (SGEM 2019), Albena, Bulgaria, 28 June–7 July 2019; pp. 905–912. [Google Scholar] [CrossRef]
  36. Ekblad, B.; Kristiansen, P.E. NMR structures and mutational analysis of the two peptides constituting the bacteriocin plantaricin S. Sci. Rep. 2019, 9, 2333. [Google Scholar] [CrossRef] [Green Version]
  37. Danielsson, J.; Pierattelli, R.; Banci, L.; Gräslund, A. High-resolution NMR studies of the zinc-binding site of the Alzheimer’s amyloid beta-peptide. FEBS J. 2007, 274, 46–59. [Google Scholar] [CrossRef]
  38. Syme, C.D.; Viles, J.H. Solution 1H NMR investigation of Zn2+ and Cd2+ binding to amyloid-beta peptide (Abeta) of Alzheimer’s disease. Biochim. Biophys. Acta 2006, 1764, 246–256. [Google Scholar] [CrossRef] [PubMed]
  39. Barranco, E.; Martín, N.; Segura, J.L.; Seoane, C.; De la Cruz, P.; Langa, F.; Gonzalez, A.; Pingarrón, J.M. Syntheses, electrochemistry and molecular modeling of N,N’-dicyanoquinonediimine (DCNQI) derivatives of substituted 1,4-anthracenediones: Precursors for organic metals. Tetrahedron 1993, 49, 4881–4892. [Google Scholar] [CrossRef]
  40. Auchus, R.J.; Miller, W.L. Molecular modeling of human P450c17 (17α-Hydroxylase/ 17,20-Lyase): Insights into reaction mechanisms and effects of mutations. J. Mol. Endocrinol. 1999, 13, 1169–1182. [Google Scholar] [CrossRef] [PubMed]
  41. Leker, R.R.; Teichner, A.; Grigoriadis, N.; Ovadia, H.; Brenneman, D.E.; Fridkin, M.; Giladi, E.; Romano, J.; Gozes, I. NAP, a femtomolar-acting peptide, protects the brain against ischemic injury by reducing apoptotic death. Stroke 2002, 33, 1085–1092. [Google Scholar] [CrossRef] [Green Version]
  42. Sragovich, S.; Ziv, Y.; Vaisvaser, S.; Shomron, N.; Hendler, T.; Gozes, I. The autism-mutated ADNP plays a key role in stress response. Transl. Psychiatry 2019, 9, 1–12. [Google Scholar] [CrossRef]
  43. Roychaudhuri, R.; Yang, M.; Hoshi, M.M.; Teplow, D.B. Amyloid β-protein assembly and Alzheimer disease. J. Biol. Chem. 2009, 284, 4749–4753. [Google Scholar] [CrossRef] [Green Version]
  44. Wang, L.; Yin, Y.L.; Liu, X.Z.; Shen, P.; Zheng, Y.G.; Lan, X.R.; Lu, C.B.; Wang, J.Z. Current understanding of metal ions in the pathogenesis of Alzheimer’s disease. Transl. Neurodegener. 2020, 9, 1–13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Sensi, S.L.; Granzotto, A.; Siotto, M.; Squitti, R. Copper and zinc dysregulation in alzheimer’s disease. Trends Pharmacol. Sci. 2018, 39, 1049–1063. [Google Scholar] [CrossRef] [PubMed]
  46. Al-Shammari, N.; Savva, L.; Kennedy-Britten, O.; Platts, J.A. Forcefield evaluation and accelerated molecular dynamics simulation of Zn(II) binding to N-terminus of amyloid-β. Comput. Biol. Chem. 2021, 93, 107540. [Google Scholar] [CrossRef]
  47. Zirah, S.; Kozin, S.A.; Mazur, A.K.; Blond, A.; Cheminant, M.; Ségalas-Milazzo, I.; Debey, P.; Rebuffat, S. Structural changes of region 1-16 of the Alzheimer disease amyloid beta-peptide upon zinc binding and in vitro aging. J. Biol. Chem. 2006, 281, 2151–2161. [Google Scholar] [CrossRef] [Green Version]
  48. Gaggelli, E.; Janicka-Klos, A.; Jankowska, E.; Kozlowski, H.; Migliorini, C.; Molteni, E.; Valensin, D.; Valensin, G.; Wieczerzak, E. NMR studies of the Zn2+ interactions with rat and human beta-amyloid (1–28) peptides in water-micelle environment. J. Phys. Chem. B 2008, 112, 100–109. [Google Scholar] [CrossRef]
  49. Trivunac, K.; Stevanovic, S. Effects of operating parameters on efficiency of cadmium and zinc removal by the complexation–filtration process. Desalination 2006, 198, 282–287. [Google Scholar] [CrossRef]
  50. Suprun, E.V.; Shumyantseva, V.V.; Archakov, A.I. Protein electrochemistry: Application in medicine. A review. Electrochim. Acta 2014, 140, 72–82. [Google Scholar] [CrossRef]
  51. Ben-Naim, A. The role of hydrogen bonds in protein folding and protein association. J. Phys. Chem. 1991, 95, 1437–1444. [Google Scholar] [CrossRef]
  52. Mallamace, D.; Fazio, E.; Mallamace, F.; Corsaro, C. The role of hydrogen bonding in the folding/unfolding process of hydrated lysozyme: A review of recent NMR and FTIR results. Int. J. Mol. Sci. 2018, 19, 3825. [Google Scholar] [CrossRef] [Green Version]
  53. Ganesan, S.J.; Matysiak, S. Interplay between the hydrophobic effect and dipole interactions in peptide aggregation at interfaces. Phys. Chem. Chem. Phys. 2016, 18, 2449–2458. [Google Scholar] [CrossRef]
  54. Festa, G.; Mallamace, F.; Sancesario, G.M.; Corsaro, C.; Mallamace, D.; Fazio, E. Aggregation states of Aβ1–40, Aβ1–42 and Aβp3–42 amyloid beta peptides: A SANS study. Int. J. Mol. Sci. 2019, 20, 4126. [Google Scholar] [CrossRef] [Green Version]
  55. Thirumalai, D.; Reddy, G.; Straub, J.E. Role of water in protein aggregation and amyloid polymorphism. ACC. Chem. Res. 2012, 45, 83–92. [Google Scholar] [CrossRef] [PubMed]
  56. Fitzpatrick, A.W.; Knowles, T.P.; Waudby, C.A.; Vendruscolo, M.; Dobson, C.M. Inversion of the balance between hydrophobic and hydrogen bonding interactions in protein folding and aggregation. PLoS Comput. Biol. 2011, 10, 1002169. [Google Scholar] [CrossRef] [Green Version]
  57. Lee, J.; Ju, M.; Hyun Cho, O.; Kim, Y.; Nam, K. Tyrosine-rich peptides as a platform for assembly and material synthesis. Adv. Sci. 2019, 6, 1801255. [Google Scholar] [CrossRef] [Green Version]
  58. Pu, S.-Z.; Zhang, W.-H.; Shi, B. Effect of pH on structure and stability of collagen-like peptide: Insight from molecular dynamics simulation. J. Theor. Comput. Chem. 2011, 10, 245–259. [Google Scholar] [CrossRef]
  59. Molecular Operating Environment (MOE); 2013.08; Chemical Computing Group ULC: Montreal, QC, Canada, 2018.
  60. Bomble, Y.J. Amsterdam Density Functional 2005 Scientific Computing and Modelling NV, Vrije Universiteit, Theoretical Chemistry, De Boelelaan 1083, 1081 HV Amsterdam, the Netherlands. J. Am. Chem. Soc. 2006, 128, 3103. [Google Scholar] [CrossRef]
  61. Schmidt, M.W.; Baldridge, K.K.; Boatz, J.A.; Elbert, S.T.; Gordon, M.S.; Jensen, J.H.; Koseki, S.; Matsunaga, N.; Nguyen, K.A.; Su, S.; et al. General atomic and molecular electronic structure system. J. Comp. Chem. 1993, 14, 1347–1363. [Google Scholar] [CrossRef]
  62. Gill, P.; Murray, W.; Wright, M. Practical Optimization; Academic Press: London, UK, 1981. [Google Scholar]
  63. Chan, S.L.; Labute, P. Training a scoring function for the alignment of small molecules. J. Chem. Inf. Model. 2010, 50, 1724–1735. [Google Scholar] [CrossRef] [PubMed]
  64. Kearsley, S.K.; Smith, G.M. An alternative method for the alignment of molecular structures: Maximizing electrostatic and steric overlap. Tetrahedron Comput. Methodol. 1990, 3, 615–633. [Google Scholar] [CrossRef]
  65. Labute, P.; Williams, C.; Feher, M.; Sourial, E.; Schmidt, J.M. Flexible alignment of small molecules. J. Med. Chem. 2001, 44, 1483–1490. [Google Scholar] [CrossRef]
  66. Clark, A.M.; Labute, P. 2D depiction of protein−ligand complexes. J. Chem. Inf. Model. 2007, 47, 1933–1944. [Google Scholar] [CrossRef]
  67. Wallace, A.C.; Laskowski, R.A.; Thornton, J.M. LIGPLOT: A program to generate schematic diagrams of protein-ligand interactions. Protein Eng. Des. Sel. 1995, 8, 127–134. [Google Scholar] [CrossRef] [PubMed]
  68. Clark, A.M.; Labute, P.; Santavy, M. 2D structure depiction. J. Chem. Inf. Model. 2006, 46, 1107–1123. [Google Scholar] [CrossRef] [PubMed]
  69. Berendsen, H.J.; Postma, J.V.; van Gunsteren, W.F.; DiNola, A.R.H.J.; Haak, J.R. Molecular dynamics with coupling to an external bath. J. Chem. Phys. 1984, 81, 3684–3690. [Google Scholar] [CrossRef] [Green Version]
  70. Bond, S.D.; Leimkuhler, B.J.; Laird, B.B. The Nosé–Poincaré method for constant temperature molecular dynamics. J. Comput. Phys. 1999, 151, 114–134. [Google Scholar] [CrossRef] [Green Version]
  71. Phillips, J.C.; Braun, R.; Wang, W.; Gumbart, J.; Tajkhorshid, E.; Villa, E.; Chipot, C.; Skeel, R.D.; Kale, L.; Schulten, K. Scalable molecular dynamics with NAMD. J. Comput. Chem. 2005, 26, 1781–1802. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. MALDI ToF MS spectra of the NAP peptide in the presence and absence of zinc ions at a 1:10 peptide to metal ratio.
Figure 1. MALDI ToF MS spectra of the NAP peptide in the presence and absence of zinc ions at a 1:10 peptide to metal ratio.
Pharmaceuticals 14 01011 g001
Figure 2. MALDI ToF MS spectra of the tyrosine-modified NAPY peptide in the presence and absence of zinc ions at a 1:10 peptide to metal ratio.
Figure 2. MALDI ToF MS spectra of the tyrosine-modified NAPY peptide in the presence and absence of zinc ions at a 1:10 peptide to metal ratio.
Pharmaceuticals 14 01011 g002
Figure 3. Two-dimensional structure of native NAP peptide (Asn-Ala-Pro-Val-Ser-Ile-Pro-Gln) with the marked alpha (α), beta (β), gamma (γ), and delta (δ) protons.
Figure 3. Two-dimensional structure of native NAP peptide (Asn-Ala-Pro-Val-Ser-Ile-Pro-Gln) with the marked alpha (α), beta (β), gamma (γ), and delta (δ) protons.
Pharmaceuticals 14 01011 g003
Figure 4. 1H-NMR spectrum of NAP peptide dissolved in D2O solvent.
Figure 4. 1H-NMR spectrum of NAP peptide dissolved in D2O solvent.
Pharmaceuticals 14 01011 g004
Figure 5. 1H-NMR spectrum of NAPY peptide dissolved in D2O solvent.
Figure 5. 1H-NMR spectrum of NAPY peptide dissolved in D2O solvent.
Pharmaceuticals 14 01011 g005
Figure 6. 1H-NMR spectrum in the aliphatic region of NAP peptide titrated with increasing amounts of Zn2+ at pH = 5.
Figure 6. 1H-NMR spectrum in the aliphatic region of NAP peptide titrated with increasing amounts of Zn2+ at pH = 5.
Pharmaceuticals 14 01011 g006
Figure 7. 1H-NMR spectrum in the aliphatic region of NAPY peptide titrated with increasing amounts of Zn2+ at pH = 5.
Figure 7. 1H-NMR spectrum in the aliphatic region of NAPY peptide titrated with increasing amounts of Zn2+ at pH = 5.
Pharmaceuticals 14 01011 g007
Figure 8. 1H-NMR spectra of NAP peptide titrated with Zn2+ ions at pH 7.5 and pH 9.5, at a 1:10 molar ratio.
Figure 8. 1H-NMR spectra of NAP peptide titrated with Zn2+ ions at pH 7.5 and pH 9.5, at a 1:10 molar ratio.
Pharmaceuticals 14 01011 g008
Figure 9. 1H-NMR spectra of NAPY peptide titrated with of Zn2+ ions at pH 7.5 and pH 9.5, at a 1:10 molar ratio.
Figure 9. 1H-NMR spectra of NAPY peptide titrated with of Zn2+ ions at pH 7.5 and pH 9.5, at a 1:10 molar ratio.
Pharmaceuticals 14 01011 g009
Figure 10. Flexible alignment of Zn2+ interaction, with NAP (A,B) and NAPY (C), performed with Molecular Operating Environment software.
Figure 10. Flexible alignment of Zn2+ interaction, with NAP (A,B) and NAPY (C), performed with Molecular Operating Environment software.
Pharmaceuticals 14 01011 g010
Table 1. The primary structure and molecular weight of the NAP/NAPY peptides in the presence of zinc ions (Zn2+), as experimentally determined with a MALDI-ToF instrument or calculated with ChemCalc online software.
Table 1. The primary structure and molecular weight of the NAP/NAPY peptides in the presence of zinc ions (Zn2+), as experimentally determined with a MALDI-ToF instrument or calculated with ChemCalc online software.
PeptideMolecular IonTheoretical (m/z)Experimental (m/z)
NAP
(H2N–NAPVSIPQ-COOH)
[M-16+H]+809.4809.6
[M+H]+825.4825.6
[M+Na]+847.4847.6
[M+K]+863.4863.6
[M+2Na-H]+869.4869.6
[M+Na+K-H]+885.4885.6
[M+Zn-H]+887.4887.6
NAPY
(H2N–NAPVYIPQ-COOH)
[M-16+H]+885.5885.7
[M+H]+901.5901.7
[M+Na]+923.5923.7
[M+K]+939.4939.6
[M+2Na-H]+945.4945.6
[M+Na+K-H]+961.4961.6
[M+Zn-H]+963.4963.6
Table 2. Assignment of signal modifications in the 1H-NMR spectra of NAP peptide at different values of pH and in the presence of zinc ions.
Table 2. Assignment of signal modifications in the 1H-NMR spectra of NAP peptide at different values of pH and in the presence of zinc ions.
Chemical Shift (ppm)NAP pH 5NAP pH 7.5NAP-Zn pH 7.5NAP-Zn pH 9.5
1Hα Asn 4.344.104.043.77
1Hβ Asn 2.992.872.822.65
2.852.722.692.51
1Hα Gln 4.364.174.174.17
1Hβ Gln 2.202.132.132.13
1Hγ Gln 2.412.362.362.36
Table 3. Assignment of signal modifications from 1H-NMR spectrum of NAP peptide at different value of pH and in the presence of zinc ions.
Table 3. Assignment of signal modifications from 1H-NMR spectrum of NAP peptide at different value of pH and in the presence of zinc ions.
Chemical Shift (ppm)NAPY pH 5NAPY pH 7.5NAPY-Zn pH 7.5NAPY-Zn pH 9.5
1Hα Asn4.333.393.873.75
1Hβ Asn2.992.742.732.65
2.852.642.592.51
1Hα Gln4.344.154.154.15
1Hβ Gln2.222.132.132.13
1Hγ Gln2.412.342.342.34
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lupaescu, A.-V.; Mocanu, C.S.; Drochioiu, G.; Ciobanu, C.-I. Zinc Binding to NAP-Type Neuroprotective Peptides: Nuclear Magnetic Resonance Studies and Molecular Modeling. Pharmaceuticals 2021, 14, 1011. https://doi.org/10.3390/ph14101011

AMA Style

Lupaescu A-V, Mocanu CS, Drochioiu G, Ciobanu C-I. Zinc Binding to NAP-Type Neuroprotective Peptides: Nuclear Magnetic Resonance Studies and Molecular Modeling. Pharmaceuticals. 2021; 14(10):1011. https://doi.org/10.3390/ph14101011

Chicago/Turabian Style

Lupaescu, Ancuta-Veronica, Cosmin Stefan Mocanu, Gabi Drochioiu, and Catalina-Ionica Ciobanu. 2021. "Zinc Binding to NAP-Type Neuroprotective Peptides: Nuclear Magnetic Resonance Studies and Molecular Modeling" Pharmaceuticals 14, no. 10: 1011. https://doi.org/10.3390/ph14101011

APA Style

Lupaescu, A. -V., Mocanu, C. S., Drochioiu, G., & Ciobanu, C. -I. (2021). Zinc Binding to NAP-Type Neuroprotective Peptides: Nuclear Magnetic Resonance Studies and Molecular Modeling. Pharmaceuticals, 14(10), 1011. https://doi.org/10.3390/ph14101011

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop