Interplay of histidine residues of the Alzheimer’s disease Aβ peptide governs its Zn-induced oligomerization

Conformational changes of Aβ peptide result in its transformation from native monomeric state to the toxic soluble dimers, oligomers and insoluble aggregates that are hallmarks of Alzheimer’s disease (AD). Interactions of zinc ions with Aβ are mediated by the N-terminal Aβ1–16 domain and appear to play a key role in AD progression. There is a range of results indicating that these interactions trigger the Aβ plaque formation. We have determined structure and functional characteristics of the metal binding domains derived from several Aβ variants and found that their zinc-induced oligomerization is governed by conformational changes in the minimal zinc binding site 6HDSGYEVHH14. The residue H6 and segment 11EVHH14, which are part of this site are crucial for formation of the two zinc-mediated interaction interfaces in Aβ. These structural determinants can be considered as promising targets for rational design of the AD-modifying drugs aimed at blocking pathological Aβ aggregation.

physiological conditions in the presence of Zn 2+ the metal binding domains of several natural Aβ variants form homo-and hetero-dimeric complexes [35][36][37][38] . Residues 11-14 of the two interacting subunits compose the dimer interface wherein two pairs of E11 and H14 residues coordinate a zinc ion 36,38 .
Along with the intact Aβ isoforms which are heterogenous at their N-termini and/or C-termini, the amyloid plaques involve a variety of chemically modified Aβ variants 39 . The Aβ species extracted from the plaques can initiate pathological aggregation of endogenous Aβ upon intracerebral injections into animal models of AD 3,40,41 . Several post-translational modifications have been discovered to increase the aggregation rate of Aβ 42 . Some chemical modifications and amino acid changes within the metal binding domain of Aβ (e.g. isomerization of D7, phosphorylation of S8, and the H6R English familial mutation associated with early onset AD) facilitate zinc-dependent dimerization and/or oligomerization of the domain 36,38,43 , thus suggesting their potential role in initiating the pathological aggregation process. Indeed, peripheral injections of the synthetic Aβ species bearing isomerized D7 (isoD7-Aβ ) was shown to trigger cerebral amyloidosis in vivo 44 . Little is known about the Aβ metal binding sites in the aggregated state (oligomers or fibrils), but in general, it seems that the binding sites are similar to those in the monomeric peptide 11,45 .
Addition of Zn 2+ ions to all studied peptides leads to substantial line broadening of the NMR signals (Fig.  S7). The nature of the peptide signal broadening primarily is associated with the exchange between multiple conformational states of the complex 46-48 . Potentially the aggregation processes can also cause signal broadening. However, the extent of line broadening for peptide H6R-Аβ 1-16 that never aggregated, and the peptides that aggregate is virtually the same. This means that aggregation does not affect NMR line widths of these peptides under the given conditions. Additionally, an extra set of resonances (Figs 3 and 4, S2 and S3, S8) has been found in NMR spectra of the peptides Аβ 1-16 , H6R-Аβ 1-16 , isoD7-Аβ 1-16 and Аβ [6][7][8][9][10][11][12][13][14][15][16] in the presence of zinc ions. This set has been assigned to the dimeric peptide complexes (see below).  Zinc-induced dimers of Аβ 1-16 , Аβ 6-16 , H6R-Аβ 1-16 and isoD7-Аβ 1-16 have common zinc-mediated dimerization interface. We previously showed by a combination of NMR, isothermal titration calorimetry (ITC) and surface plasmon resonance methods that a stable dimeric form of the peptide H6R-Аβ 1-16 is formed when bound to a zinc ion 36 . Residues E11 and H14 of the two peptide subunits formed  A comparison of the NMR spectra of the peptides Аβ 1-16 , Аβ 6-16 , H6R-Аβ 1-16 and isoD7-Аβ 1-16 in the presence of zinc ions obtained in this study has shown a set of characteristic signals similar to that observed in the work cited above, for all these peptides (Fig. 3, S3, S8). At low concentration of the peptides Аβ 1-16 , H6R-Аβ 1-16 and isoD7-Аβ 1-16 (~0.2 mM) fraction of the dimeric form is nearly identical (Fig. 3). Fraction of the dimeric forms rises with concentration of the peptides as expected for any dimerization process. Such tendency has been shown for H6R-Аβ 1- 16 36 , which has the highest dimer abundance among the studied peptides. The common set of characteristic signals found for the dimers of Аβ 1-16 , Аβ 6-16 , H6R-Аβ 1-16 and isoD7-Аβ [1][2][3][4][5][6][7][8][9][10][11][12][13][14][15][16] indicates that all the dimers have the same zinc-induced dimerization interface, namely, 11 EVHH 14 . Equilibrium between zinc-bound monomers and dimers of Аβ 1-16 isoforms. Equilibrium between free peptide and its zinc-bound monomer complex is fast on the NMR time scale. In contrast, equilibrium between zinc-bound monomers and dimers is slow on the NMR time scale, as evidenced by the two separate sets of NMR signals (Fig. 3, 4, S2,S3,S8). The exchange between monomers and dimers is unambiguously confirmed by the rotating frame nuclear Overhauser effect spectroscopy (ROESY) since cross-peaks originating from the through-space dipole-dipole interaction (NOE) have opposite signs to the cross-peaks derived from chemical exchange. Figure 4 illustrates NOEs between S8 Hα and Hβ protons (negative, blue), and exchange cross peaks between several Hα signals of monomeric and dimeric forms (red, positive). Exchange between zinc-bound monomers and dimers of H6R-Аβ 1-16 and isoD7-Аβ 1-16 is also confirmed by characteristic cross-peaks in the NOESY spectra ( Fig. S8) demonstrating exchange between V12 Hγ 1* signals. Both peptides show identical patterns of NMR signals and differ only by fraction of the dimeric forms.
Rate constants k m → d and k d → m at equilibrium between the monomeric and dimeric complexes of H6R-Аβ 1-16 with Zn 2+ have been measured using the magnetization transfer NMR experiments (see Supporting information pp. S24-S28 for details). It has been found that an effective rate constant k m→d measured at the total peptide concentration of 2.3 mM in the presence of half molar equivalence of ZnCl 2 equals to 8.6 ± 0.6 s −1 and k d→m equals to 31.3 ± 2.3 s −1 .

Structure of zinc-mediated H6R-Аβ 1-16 dimer in solution.
Significant ratio of stable (on NMR time scale) zinc-bound H6R-Аβ 1-16 dimers allowed to determine the dimer structure in solution. The NMR structure has been determined using the set of distance restraints collected from 2D NOESY and ROESY spectra and a set of constraints between zinc ion and residues E11 and H14 (Table 2). In addition to NOEs between the signals of the dimer, NOEs between the signals of zinc-bound monomer have been also included in the list of distance restraints used for structure calculation (Table 2). Due to the equilibrium between zinc-bound monomer and dimer, more intensive signals of monomer ( Fig. 5) contain all the information about interatomic distances within the dimer via the transferred NOE mechanism. Similar transferred NOE mechanism is observed in the equilibrium between free peptide and its complex with a larger protein 49 . Effectiveness of such mechanism is determined by the substantially more effective cross-relaxation in dimer due to its slower molecular tumbling.
QM/MM calculations have been carried out to determine the restraints that describe geometry of the zinc-mediated interface formed by pairs of E11 and H14 from the interacting subunits. An approach similar to that described earlier for the determination of NMR structure of rat Аβ 1-16 complex with zinc ions 47 has been used in the calculations. Tetrahedral geometry of the zinc ion coordination sphere originated from quantum mechanics calculations (Fig. S12B). Such geometry is in good agreement with the principles governing Zn binding in proteins 50 and in agreement with the previously determined structures of Аβ 1-16 zinc-bound complexes 11 Table 2. NMR restrains and structural statistics for the complex of two H6R-Aβ 1-16 peptides with zinc ion.
Single set of zinc-bound dimer chemical shifts indicates that dimer subunits are chemically equivalent. Therefore, each NOE has been assumed to represent either intra-chain correlation or correlation between two adjacent subunits. All NOEs have been treated in the structure calculations as ambiguous distance restraints allowing optimization protocol to find optimal assignments.
A family of 20 NMR structures has been determined on the basis of 181 experimental restraints (see Table 1 for details) using simulated annealing MD protocol in explicit water environment 51 . For most of the residues, the number of NOEs per residue is between 15 and 40 (Fig. S9). Statistics of the final ensemble are given in Table 1 and superposition of the final family of calculated structures is presented in Fig. S10. A representative structure was selected from the ensemble of calculated structures as being the one that is closest to all the other structures and thus gives the lowest sum of pairwise RMSD for the remainder of the structures in the family. RMSD between the family of calculated structures and the representative structure is about 2.3 Å for the backbone heavy atoms. RMSD between the structures in the final family for heavy atoms of the dimerization interface core (residues 8-15) is about 1.0 Å (Table 1). On Ramachandran plot (Fig. S11), 73.5% of the residues in the whole NMR family are in the most favored regions and none in the disallowed regions. Representative structure has been additionally optimized using the QM/MM method in order to refine the geometry of Zn 2+ environment (Figs 6 and S12). In the interface, pairs of residues E11 and H14 of the interacting subunits coordinate a common zinc ion. In the presence of equimolar amounts of zinc ions only H6R-Аβ 1-16 does not precipitate up to 10 mM, whereas the other peptides rapidly precipitate when their concentration reaches threshold values (> 5 mM for Аβ 1-16 , > 0.2 mM for isoD7-Аβ 1-16 , > 0.8 mM for Аβ 6-16 ) ( Table 1). Apparently, precipitation can occur only if an oligomer subunit possesses at least two sites that sterically can be involved in the formation of zinc-mediated interfaces with other peptide subunits. Accordingly, the Аβ 1-16 , Аβ 6-16 and isoD7-Аβ 1-16 peptides must have a second dimerization interface located within the Аβ minimal zinc binding site 6-14 34 . The fact that aggregation of H6R-Аβ 1-16 is not observed at any peptide concentration as shown in Table 1, clearly indicates the involvement of residue H6 in the second zinc-dependent interface. Considering that side chains of residues Asp, Glu и His are typical zinc chelators and E11 и H14 are part of the primary dimerization interface, one can rationally assume that the second dimerization interface should include residues D7 (isoD7) or H13 in addition to H6.
To probe involvement of isoD7 in the second dimerization interface we have studied a capacity of the model peptide isoD7-Аβ 1-10 to form dimers in the presence of zinc ions. We have performed NMR titration experiments using the method of continuous variations 36,38,52 . The results show that stoichiometry of interaction of isoD7-Аβ 1-10 with Zn 2+ is 1:1 (Fig. S13B, Table 1). The coordination centers of zinc ion in the monomeric complex with isoD7-Аβ 1-10 have been identified from changes of the chemical shifts between free and zinc-bound states (Fig. S14, Tables S5, S6) . In addition to imidazole ring of H6, Zn 2+ is coordinated by the side chain carboxyl group of isoD7 and two backbone carbonyl groups of the residues F4 and H6. Chemical shift changes data (Fig. S15) allowed to calculate the binding constant of Zn 2+ to isoD7-Аβ 1-10 , K a = 1.19 ± 0.06·10 3 M −1 (Fig. S13A, Table 1). The data show that the peptide forms a monomeric complex with zinc ion and thus isoD7 is not involved in the second dimerization interface.
Similarly, to probe involvement of residue H13 in the second dimerization interface we have studied the model peptide isoD7-H13R-Аβ 1-16 . The H13R substitution has been designed to minimize changes in the peptide properties relative to isoD7-Аβ 1-16 . Chemical shifts of the peptide in the presence and absence of zinc ions are given in Supplementary material (Tables S3 and S4 ). NMR studies of the interaction of isoD7-H13R-Аβ 1-16 with Zn 2+ (Fig. S16) have shown zinc binding constant K a = 1.35 ± 0.06·10 3 M −1 (Fig. S17A, Table 1) and stoichiometry (1:1) (Fig. S17B, Table 1). Zinc-induced chemical shift changes (Fig. S18) indicate that side chains of the two histidine residues (H6 and H14) are involved in the coordination of Zn 2+ . Chemical shifts of E11 do not change upon interaction of the peptide with Zn 2+ , suggesting that this residue is not involved in its coordination. RMSD of the 1 H and 13 C chemical shifts between the free and zinc-bound peptides (Fig. S18) indicate that zinc ion is coordinated by backbone carbonyl groups of the residues H6 and R5, similarly to that observed earlier in the S8 phophorylated Aβ 1-16 peptide 38 .
Taken together, the described data show that in zinc-induced oligomers of the Аβ 1-16 , Аβ 6-16 and isoD7-Аβ 1-16 peptides each subunit interacts with the adjacent subunits via residues E11 and H14 from the primary dimerization interface, and residues H6 and H13 from the second one.
Our experimental results indicate that isomerization of D7 as well as truncation of the first five residues facilitate formation of the second Zn-dependent interaction site. To understand why a relatively small change of the peptide structure between Аβ 1-16 and isoD7-Аβ 1-16 causes dramatic change in their zinc-induced oligomerization, MD simulations have been performed for homodimer Aβ 1-16 ·Aβ 1-16 and heterodimer Aβ 1-16 ·isoD7-Aβ 1-16 . Conformational behavior of Aβ 1-16 and isoD7-Aβ 1-16 chains in the respective homodimer and heterodimer, has been examined in MD trajectories. In the MD simulations the primary zinc-dependent interface structure has been kept in its initial conformation, and distances between the imidazole rings of residues H6 and H13 in the second interface in both Aβ 1-16 and isoD7-Aβ 1-16 chains have been constrained in order to keep them close enough for zinc ion coordination.
MD simulations have shown that residues D7 and S8 in the Aβ 1-16 chain form a stable bend (RMSD over N, C and Cα atoms of residues 7 and 8 was 0.13 Å). Such bend is stabilized by hydrogen bond between the D7 side-chain carboxyl and S8 side-chain hydroxyl (Fig. 7a). In contrast, in the isoD7-Aβ 1-16 chain residues D7 and S8 do not form a bend, but demonstrate a propensity to form extended conformation (Fig. 7b). Thus, structural changes associated with the isomerization of D7 (elongation of the peptide backbone by one additional CH 2 group and reduction of the side-chain) lead to disruption of the interaction between residues D7 and S8.
MD results show that in both peptides Аβ 1-16 and isoD7-Аβ 1-16 , residues H6 and H13 adopt a conformation supporting interaction with Zn 2+ . However, probability to adopt the most favorable conformation is higher for the isoD7-Аβ 1-16 peptide. On the contrary, such conformational space in Аβ 1-16 peptide is restricted due to the bent region 7-8 (Fig. S19). Thus, isomerization of D7 increases capacity of the residues H6 and H13 to form the second zinc-dependent interface interaction with zinc ion, which facilitates oligomerization of the peptide. Truncation of the first five residues of Аβ 1-16 forms the Аβ 6-16 peptide with increased conformational freedom of the N-terminal residue H6 and similarly facilitates zinc-induced oligomerization.

Discussion
In human tissues, Aβ is a group of peptides, heterogenous at their N-termini and/or C-termini and produced by the β -and γ -secretase-dependent cleavage of amyloid precursor protein 53 . Structural characterization of the full lenght Aβ monomer in solution is difficult due to the tendency of the peptide to aggregate. It was found that the soluble Aβ 1-42 peptide is an intrinsically disordered polypeptide in aqueous solution, having β -strand secondary structure within the segments 2-7, 16-23, 28-32, and 34-36 54 . Solid-state NMR spectroscopy and electron microscopy allowed to identify diverse morphologies and structures in Aβ fibrils (reviewed in 55,56 ). The metal-binding domain 1-16 was found to be flexible, situated outside of the hydrophobic core formed by the residues 17-42 20,57 . It is worth noting that the absence of this domain in the so called p3 peptide, which comprises the Aβ region 17-42, prevents the formation of oligomers 57,58 . Another fact that three amino acid substitutions that discern Aβ of rats and mice from all other mammals and protect rodents against AD are situated within the metal binding domain, indicating its involvement in initialization of the pathological Aβ aggregation. Several AD-causing mutations are also located in the metal binding domain of Aβ 59,60 . A role of Aβ 1-16 in the AD pathology is that it is required for the formation of zinc-induced oligomers. Antibodies targeting the N-terminal region of Aβ are able to block the formation of Aβ 40 amyloids 61 , thus suggesting that the N-terminal domain contributes to the formation and/or the maturation of Aβ fibrils.
Our results indicate that conformation of the 10-15 region of the metal binding domain in absence of zinc ions is identical in all studied variants of Aβ [1][2][3][4][5][6][7][8][9][10][11][12][13][14][15][16] , and that it is likely to be pre-formed for an effective zinc ion trapping. The region 11-14 was shown to play a principle role in Zn 2+ binding in the monomeric complex 17,34 . Such interaction includes two consecutive steps, (i) primary zinc ion binding by side chains of the residues E11, H13, H14 and a water molecule, (ii) substitution of water by the residue H6 due to movement of the fragments 1-8 and 11-14 towards each other. It was suggested that a zinc ion binds to monomeric Aβ and then the zinc-peptide complex undergoes conformational changes that lead to the formation in the formation of zinc-induced Aβ oligomers 11,21,62,63 . The role of 11 EVHH 14 fragment in the formation of Aβ dimers was reported previously 11,[36][37][38]64 . In the current study we have established the three-dimensional structure of a zinc-linked dimer of H6R-Аβ 1-16 , with the dimerization interface formed by the side chains of residues E11 and H14. Residues 10-15 also contribute to structure stabilization by hydrophobic interactions between Y10 and V12, and by forming hydrogen bonds between the HN protons of H14 and Gln15 and the backbone carbonyl groups of E11 and V12 respectively (Fig. 6, Fig. S12). Low field shifts of these two amide protons are in agreement with this conclusion. Characteristic high field shift of the resonance of Hγ 1*of V12 originates from its proximity to the aromatic ring of Y10 in the hydrophobic core and additionally validates the structure. Thus, binding of a zinc ion to the 11 EVHH 14 fragment leads to formation of the peptide dimer where one zinc ion is coordinated by the residues E11 and H14 from the two interacting peptide chains.
Miller et al. 62 discuss two possible scenarios of the zinc-induced amyloidogenesis. In the first one, metal ions bind amyloid monomers and induce their assembly into oligomers via interactions of the zinc-bound complexes. Alternative mechanism presumes that metal ions bind to pre-formed apo-oligomers of Aβ . Lim et al. 13 showed that zinc binding to Aβ initiates a cascade of conformational changes leading to intermolecular interactions of Aβ via residues 12-21. After an initial structure rearrangement caused by zinc binding, the C-terminal residues readjust their conformation to support effective intermolecular interaction 18,21,47,[57][58][59]61,65 . Similar structural changes of the metal-free Aβ 1-40 peptides were also observed in the presence of the preformed oligomers 13 , suggesting that such conformational transitions may constitute a general molecular mechanism of the Aβ amyloidogenesis.
We have shown here for Aβ 1-16 and isoD7-Aβ 1-16 that formation of the primary zinc-dependent dimerization interface by residues E11 and H14 initiates conformational rearrangement leading to formation of the second dimerization interface by residues H6 and H13 (Fig. 7). Emergence of this second interface is a key event enabling zinc-induced oligomerization of the metal binding domain (Fig. 8).
If a zinc-induced dimer is formed by native Aβ 1-16 via primary interface, it restricts conformational mobility of the peptide. We hypothesize that this can result in disruption of the reciprocal aligning of residues 6 and 13 necessary to form the second interface, leading to lower fraction of peptides with both interfaces (Fig. S19). This explains moderate susceptibility of the native peptide Aβ 1-16 to oligomerization (Aβ 1-16 precipitates at concentrations > 5 mM, Table 1). In comparison with Aβ 1-16 , isoD7-Aβ 1-16 undergoes oligomerization at substantially higher rate (isoD7-Aβ 1-16 precipitates at concentrations > 0.3 mM). Isomerization of D7 considerably increases conformational space of the peptide backbone and facilitates favorable aligning of the side chains of residues H6 and H13 that form the second dimerization interface. Another way to increase conformational freedom of the H6 side chain is to remove the first five residues of Aβ , including residues E1 and D3 capable of electrostatically interacting with H6 34 . Indeed, peptide Aβ 6-16 precipitates at ~0.8 mM (Table 1). Notably, enzymatic removal of the first five Aβ residues is performed by ACE 66 as we showed earlier, which can potentially link ACE to Alzheimer's disease.
The results obtained in this study allow to propose the molecular mechanism of zinc-dependent oligomerization of Aβ metal binding domain (Fig. 8). Oligomerization starts with formation of zinc-peptide monomeric complex (Fig. 8b), and subsequent transformation of the complex to a dimer (Fig. 8c), where zinc ion is coordinated by the side chains of residues E11 and H14 from the interacting peptide molecules. After this, conformational rearrangement in the segments 6-14 of each subunit leads to formation of the second zinc-dependent dimerization interface composed of residues H6 and H13. The dimer becomes a seed of subsequent zinc-dependent oligomerization leading to formation of higher order soluble oligomers that are transformed into insoluble aggregates (Fig. 8d). In contrast with the earlier concepts of polymorphism of Aβ within zinc-dependent oligomers, our data suggest that the Aβ metal binding domain has a distinct three-dimensional structure, allowing the domain to simultaneously interact with two other Aβ molecules. This could trigger a "chain reaction" of zinc-induced Aβ oligomerization. This mechanism is in line with our recent in vivo studies showing that synthetic peptides

Conclusions
We have demonstrated that in the Aβ metal binding domains of intact, H6R mutant and the isoD7 Aβ isoforms in the absence and presence of zinc ions, the dominant backbone conformation of the fragment 10-15 is identical. This fragment forms primary zinc-mediated dimerization interface of all studied metal binding domains. Solution structure of zinc-mediated H6R-Аβ 1-16 dimer has been determined, providing insight into the mechanism of formation of the dimerization interface. Zinc-induced oligomerization of synthetic peptides representing the 1-16 metal binding domains of natural Аβ variants has been shown to follow the same molecular mechanism: (i) peptide dimer is formed through the primary zinc-mediated interface 11 EVHH 14 ; (ii) residues H6 and H13 are realigned creating the second zinc-dependent interface in each subunit; (iii) the dimer becomes a seed of subsequent zinc-dependent oligomerization of Aβ 1-16 . Our results indicate that the extent of conformational freedom of residue H6 determines the propensity of Aβ 1-16 isoforms to undergo zinc-induced oligomerization. Targeting of Aβ zinc-mediated interfaces may provide a therapeutic route for AD treatment.

Dynamic light scattering. Dynamic light scattering (DLS) measurements were carried out on a Zetasizer
Nano ZS apparatus (Malvern Instruments Ltd., UK) in accordance with the manufacturer instruction. The 120-μ L aliquots of peptide solutions were placed into a BRAND UV microcuvetter (BRAND GMBH, Germany) and used for the measurements. Measurements of peptides in the presence of Zn 2+ were carried out within 10 minutes after addition of two-fold molar excess of ZnCl 2 to the peptide solutions. The instrument is equipped with a He-Ne laser source (λ = 632.8 nm) and operates in the back-scatting mode, measuring particle size in the range between 0.6 nm and 10 μ m. Particle size distribution was estimated using a CONTIN data analysis utility with spherical approximation of the particles, available as a part of the instrument software.
Turbidity measurements. Measurements were performed on a NanoDrop 1000 spectrophotometer (Thermo Scientific, USA). The optical density of peptide solutions was measured at 350 nm and 405 nm, using 2 μ L aliquots of the peptide solutions. Measurements of peptides in the presence of Zn 2+ were carried out after 30-40 minutes following addition of ZnCl 2 to the peptide solutions.s.

Isothermal titration calorimetry (ITC).
The thermodynamic parameters of zinc binding to Aβ 6-16 were measured using a MicroCal iTC200 System (GE Healthcare Life Sciences, USA) as described previously 36,43,47 . Experiments were carried out at 25 °C in 50 mM Tris buffer, pH 7.3. 2 μ L aliquots of the ZnCl 2 solution were injected into the 0.2 mL cell containing the peptide solution to obtain a complete binding isotherm. The titration curves were fitted using MicroCal Origin software. Association constant (K a ), binding stoichiometry and enthalpy were determined by a non-linear regression fitting procedure (Fig. S20).

NMR titration experiments.
To identify the amino acid residues that coordinate zinc ion in the peptides isoD7-Аβ 1-10 and isoD7-H13R-Аβ 1-16 , and for determining association constants of zinc ions NMR titration Scientific RepoRts | 6:21734 | DOI: 10.1038/srep21734 experiments were carried out. Peptides at concentration of 1.0-1.5 mM at pH 6.8-7.0 were titrated with a solution of ZnCl 2 in a buffer of identical composition at the same pH value. 1D spectra were recorded for each titration point. Figures S15 and S16 show changes of the chemical shifts with increasing molar content of [Zn 2+ ] from 0.05 to 10.0 relatively to the peptide concentration. The volume of solution at the final point increased from 600 to 800 μ L. Chemical shift changes of the representative signals in titration experiments were used for calculation of the K a values (Figs S13A and S17A). The linear nature of Δ δ values Δ δ presented in Scatchard coordinates 52 confirms equilibrium between the free and zinc-bound peptide forms.
Method of continuous variations 52 was used to determine stoichiometry of zinc binding to the peptides isoD7-Аβ 1-10 and isoD7-H13R-Аβ 1-16 . This involved preparation of a series of samples containing both the peptide and ZnCl 2 in varying proportions of the components, but in a fixed total concentration (from 1.0 to 1.7 mM). Changes of the chemical shifts induced by the interaction of the peptide with zinc ions were analyzed. The plots of the product Δ δ ·[P] 0 (Δ δ -change of the chemical shift; [P] 0 -total peptide concentration in the sample) versus the mole fraction of [Zn 2+ ] show maximum at the fraction value, which corresponds to the stoichiometry (Figs S13B and S17B).

Magnetization transfer experiments.
In order to measure exchange rates between monomeric and dimeric states of the Zn·H6R-Аβ 1-16 complex, magnetization transfer NMR technique was used. Magnetization transfer experiments involved selective excitation of the signal at 0.9 ppm, which belongs to the methyl group Hγ 1* of V12 in the monomeric form of zinc-peptide complex, with subsequent detection of the intensity of the signal at 0.2 ppm, which belongs to the same group in the dimeric form. Series of experiments with varying delays between the inverting and reading pulses were carried out. All the experimental details and data analysis are provided in the supplementary material (pages S24-S28). NMR restrains. NOE distance restraints used in structure calculation of the dimeric form of Zn-H6R-Аβ 1-16 complex were obtained from 2D NOESY spectra acquired at 274 K in H 2 O or at 278 K in D 2 O. NMR signal assignment of Zn···H6R-Аβ 1-16 complex was reported earlier 36 . Representative fragment of the NOESY spectrum with some principal assignments is shown on Fig. 5. Signal intensities in NOESY spectra were calibrated and converted into the distance restrains using the fixed distance intraresidue NOEs as the reference. Distance and torsion angle restraints representing the coordination site of zinc ion were obtained using the quantum mechanical calculations (see below).

QM/MM and restrained molecular dynamics.
Structures of the complex of H6R-Аβ 1-16 with zinc ions were determined using the GROMACS 3.3.1 software 70 , AMBER 03 force field 71 and optimized protocol of the simulated annealing MD calculations 47,51 with the set of distance restraints listed in Table 2. QM/MM Car-Parrinello simulation approach 72,73 was applied to optimize the zinc binding site and to obtain restrains that describe geometry of Zn 2+ coordination center. All the parameters and protocols of the MD and QM/MM calculations were described earlier in detail 47 . Convergence of the calculations was determined using root mean square deviation (RMSD) of the coordinates of the C, Cα , and N atoms of protein backbone, calculated for the whole family of structures. Quality of the calculated structures was defined on the basis of the percent of hits of the main dihedral angles ϕ and ψ in the most favorable and prohibited areas of Ramachandran map using Procheck_ NMR 74 . Structures were visualized and analyzed using the InsightII or Discovery Studio software (Accelrys Inc.). MD simulations. Molecular modeling has been performed for homodimer Aβ 1-16 ·Aβ 1-16 and heterodimer Aβ 1-16 ·isoD7-Aβ 1-16 using the representative NMR conformer of the H6R-Aβ 1-16 dimer as initial structure. Homology modeling was performed using the Chimera software 75 . Partial charges for the non-standard residue isoD7 were assigned using the AM1-BCC method 76 . Initial models were refined using the restrained molecular dynamic simulation with GROMACS 4.6.5 software 70 and Amber99SB-ILDN force field 77 . Peptide dimers were placed in the cubic cell with a minimum distance (0.8 nm) between a protein and the box wall and soaked with TIP3P water molecules 78 . Total charge of the solution has been neutralized with sodium ions. Energy was minimized using the steepest descent algorithm and the system was further equilibrated during 100 ps of constant volume (NVT) molecular dynamic simulation followed by 100 ps of constant pressure (NPT) molecular dynamics. 10 ns MD simulations were carried out using the NPT ensemble to relax protein chains. 6 Å restraint has been applied to the distance between H6 Nε 2 and H13 Nε 2 atoms. 20 ns NPT MD simulations were performed to follow approaching of imidazole rings of H6 and H13. Position of the side chains of residues E11 and H14 together with the coordinated zinc ion were restrained during all simulation steps. Calculations were performed at 300 K, pressure 1 bar, with a 2 fs integration step using Berendsen barostat and velocity rescale method for thermostat. Particle-mesh Ewald method 79 has been implemented to treat long-range electrostatic interactions, and LINCS algorithm to control the lengths of covalent bonds 80 .