Elucidation of the coordination of metal ions to Aβ is essential to understand their role in its aggregation and to rationally design new chelators with potential therapeutic applications in Alzheimer disease. Because of that, in the last 10 years several studies have focused their attention in determining the coordination properties of Cu(2+) interacting with Aβ. However, more important than characterizing the first coordination sphere of the metal is the determination of the whole Cu(2+)-Aβ structure. In this study, we combine homology modeling (HM) techniques with quantum mechanics based approaches (QM) to determine plausible three-dimensional models for Cu(2+)-Aβ(1-16) with three histidines in their coordination sphere. We considered both ε and δ coordination of histidines 6, 13, and 14 as well as the coordination of different possible candidates containing oxygen as fourth ligand (Asp1, Glu3, Asp7, Glu11, and CO(Ala2)). Among the 32 models that enclose COO(-), the lowest energy structures correspond to [O(E3),N(δ)(H6),N(ε)(H13),N(ε)(H14)] (1), [O(E3),N(δ)(H6),N(δ)(H13),N(δ)(H14)] (2), and [O(D7),N(ε)(H6),N(δ)(H13),N(δ)(H14)] (3). The most stable model containing CO(Ala2) as fourth ligand in the Cu(2+) coordination sphere is [O(c)(A2),N(ε)(H6),N(δ)(H13),N(ε)(H14)] (4). An estimation of the relative stability between Glu3 (1) and CO(Ala2) (4) coordinated complexes seems to indicate that the preference for the latter coordination may be due to solvent effects. The present results also show the relationship between the peptidic and metallic moieties in defining the overall geometry of the complex and illustrate that the final stability of the complexes results from a balance between the metal coordination site and amyloid folding upon complexation.

Oxidative stress induced by redox-active metal cations such as Cu(2+) is a key event in the development of Alzheimer's disease. A detailed knowledge of the structure of Cu(2+)-Aβ complex is thus important to get a better understanding of this critical process. In the present study, we use a computational approach that combines homology modeling with quantum-mechanics-based methods to determine plausible 3D structures of Cu(2+)-Aβ(1-16) complexes that enclose the different metal coordination spheres proposed experimentally at different pH values. With these models in hand, we determine their standard reduction potential (SRP) with the aim of getting new insights into the relation between the structure of these complexes and their redox behavior. Results show that in all cases copper reduction induces CObackbone decoordination, which, for distorted square planar structures in the oxidized state (Ia_δδ, IIa_εδε, IIa_εεε, and IIc_ε), leads to tricoordinated species. For the pentacoordinated structural candidate Ib_δε with Glu11 at the apical position, the reduction leads to a distorted tetrahedral structure. The present results highlight the importance of the nature of the ligands on the SRP. The computed values (with respect to the standard hydrogen electrode) for complexes enclosing negatively charged ligands in the coordination sphere (from -0.81 to -0.12 V) are significantly lower than those computed for models involving neutral ligands (from 0.19 to 0.28 V). Major geometry changes induced by reduction, on both the metal site and the peptide configuration, are discussed as well as their possible influence in the formation of reactive oxygen species.

The imidazole ring is part of the lateral chain of histidine. One of the main features of this amino acid is the ability to coordinate copper, especially Cu(2+), because of the intermediate base nature of its imidazole ring, which has a great biological relevance. Proteins such as cytochrome c oxidase, a crucial enzyme in the respiratory chain, and β-amyloid peptide, implicated in the pathology of Alzheimer's disease, are examples of proteins containing histidines in their coordination sphere. Several studies indicate that the presence of this metal ion produces a decrease in the pK(a) of the imidazole ring of histidine. However, there are no reports of systematic studies of pK(a) variation in these types of metal cation complexes. In this work we use density functional theory to study the dependence of imidazole pK(a) with the number of imidazole rings in Cu(2+) coordination environments. The pK(a) of isolated imidazole (ImH), and the pK(a) of imidazole in Cu(2+)(ImH)(m)(H(2)O)(4-m) (m=1-3) complexes have been studied using two different functionals, B3LYP and MPWB1K, which have different percentage of exact exchange, and the highly-correlated CCSD(T) method. Results show that imidazole pK(a) decreases between 2 and 7 units depending on the method employed and the number of imidazole rings coordinating the metal cation. Taking into account that the pK(a) of imidazole is 14, this decrease could be relevant in biological processes.

Iron is one of the most abundant metals found in senile plaques of post mortem patients with Alzheimer's disease. However, the interaction mode between iron ions and β-amyloid peptide as well as their precise affinity is unknown. In this study we apply ab initio computational methodology to calculate binding energies of Fe(2+/3+) with the His13-His14 sequence of Aβ, as well as other important ligands such as His6 and Tyr10. Calculations were carried out at the "MP2/6-311+G(2df,2p)"//B3LYP/6-31+G(d) level of theory and solvent effects included by the IEFPCM procedure. Several reaction paths for the binding of imidazole, phenol, and the His13-His14 fragment (modeled by N-(2-(1H-imidazol-4-yl)ethyl)-3-(1H-imidazol-4-yl)propanamide) were sequentially explored. The results show that the most stable complexes containing His13-His14 and phenolate of Tyr10 are the pentacoordinated [Fe(2+)(O-HisHis)(PhO(-))(H(2)O)](+) and [Fe(3+)(N-HisHis)(PhO(-))(H(2)O)](+) compounds and that simultaneous coordination of tyrosine and His13-His14 to Fe(2+/3+) is thermodynamically favorable in water at physiological pH. Computed Raman spectra confirm the conclusion obtained by Miura et al. ( Biochemistry 2000 , 39 , 7024 ) that tyrosine is coordinated to Fe(3+) but do not exclude coordination of imidazoles. Finally, calculations of standard reduction potentials indicate that phenol coordination reduces the redox activity of the iron/Aβ complexes.

We investigate, by means of density-functional theory, the binding of dioxygen to Cu(I)-amyloid β (Aβ), one of the first steps in the oxidation of ascorbate by dioxygen. Cu, Aβ, ascorbate and dioxygen are all present in the synapse during neurodegeneration, when the above species can trigger an irreversible oxidative stress inducing the eventual death of neurons. The binding of dioxygen to Cu(I) is possible and its role in dioxygen activation of Cu ligands and of residues in the first coordination sphere is described in atomic detail. Dioxygen is activated when a micro-environment suitable for a square-planar Cu(2+) coordination is present and a negatively charged group like Asp 1 carboxylate takes part in the Cu coordination anti to O2.