As the main sequence responsible for metal ion coordination in the amyloid beta (Aβ) peptide, Aβ_1–16 plays a key role in the understanding of the aggregation of Aβ induced by Cu²⁺ ions. There is no consensus on the nature of the coordination sphere of the Cu²⁺–Aβ complex so far due to the amorphous conformation of the Aβ_1–16 peptide itself and the pH dependence of Cu²⁺–Aβ coordination. The simulation reported here reveals that human Aβ_1–16 monomer has a U-shape morphology, which is preserved at any pH. This morphology accommodates Cu²⁺ ions with several binding sites and is also the basis for establishing a center-distance statistical method (CDSM). Based on this CDSM, specific histidine residues for a Cu²⁺-coordinated sphere are identified and the corresponding accurate pH range is established, indicating that the CDSM can be used as a reference to predict the potential coordination sites of metal ions in other amorphous peptides. By contrast, mouse Aβ_1–16 monomer has a more open and random morphology than human Aβ_1–16 due to the differences of three sequence positions. These mutations not only reduce the number of binding sites required by a stable Cu²⁺-binding sphere but also diminish the capacity to generate salt bridges compared to the human peptide. These observations offer insights into the roles of three residues that differ in the mouse Aβ_1–16 and perhaps into the reasons mice seldom develop Alzheimer's disease.

In this study, structural and mechanical properties of a series of models of Aβ₄₂ (one- and two-fold) and Aβ₄₀ (two- and three-fold) fibrils have been computed by using all-atom molecular dynamics simulations. Based on calculations of the twist angle (θ) and periodicity (v=360d/θ), oligomers formed by 20, 11, and 13 monomers were found to be the smallest realistic models of three-fold Aβ₄₀, one-fold Aβ₄₂, and two-fold Aβ₄₂ fibrils, respectively. Our results predict that the Aβ₄₀ fibrils initially exist in two staggered conformations [STAG(+2) and STAG(+1)] and then undergo a [STAG(+2)STAG(+1)] transformation in a size-dependent manner. The length of the loop region consisting of the residues 23–29 shrinks with the elongation of both Aβ₄₀ and Aβ₄₂ fibrils. A comparison of the computed potential energy suggests that a two-fold Aβ₄₀ aggregate is more stable than its three-fold counterpart, and that Aβ₄₂ oligomers can exist only in one-fold conformation for aggregates of more than 11 monomers in length. The computed Young′s modulus and yield strengths of 50 GPa and 0.95 GPa, respectively, show that these aggregates possess excellent material properties.

Stability orderings of 150 stable complexes formed by metal ions (Na⁺, K⁺, Ca²⁺, Mg²⁺, and Zn²⁺) and 13 stable thymine tautomers in both solvent and gas phases are obtained, and the optimal binding site for a metal ion in a specific thymine tautomer is identified. Results indicate that the complex with the canonical thymine tautomer (T1) is more stable than those with the rare ones, and the monodentate complex M–T1o4(o2) are their ground-state form in the solvent phase. The ground-state thymine complexes bound by Ca²⁺, Mg²⁺, or Zn²⁺ become bidentate M–T3o4lo2,n3, which is derived from a rare thymine tautomer T3o4l, whereas those bound by Na⁺ and K⁺ are still monodentate complexes M–T1o4(o2), however, in the gas phase. The differences in stability are discussed in detail from the binding strength of metal ions, relative energy of the corresponding thymine tautomers, and solution effect. Copyright © 2011 John Wiley & Sons, Ltd.