Structural and electronic analysis of the octarepeat region of prion protein with four Cu2+ by polarizable MD and QM/MM simulations

Prions have been linked to neurodegenerative diseases that affect various species of mammals including humans. The prion protein, located mainly in neurons, is believed to play the role of metal ion transporter. High levels of copper ions have been related to structural changes. A 32-residue region of the N-terminal domain, known as octarepeat, can bind up to four copper ions. Different coordination modes have been observed and are strongly dependent on Cu2+ concentration. Many theoretical studies carried out so far have focused on studying the coordination modes of a single copper ion. In this work we investigate the octarepeat region coordinated with four copper ions. Molecular dynamics (MD) and hybrid quantum mechanics/molecular mechanics (QM/MM) simulations using the polarizable AMOEBA force field have been carried out. The polarizable MD simulations starting from a fully extended conformation indicate that the tetra-Cu2+/octarepeat complex forms a globular structure. The globular form is stabilized by interactions between Cu2+ and tryptophan residues resulting in some coordination sites observed to be in close proximity, in agreement with experimental results. Subsequent QM/MM simulations on several snapshots suggests the system is in a high-spin quintet state, with all Cu2+ bearing one single electron, and all unpaired electrons are ferromagnetically coupled. NMR simulations on selected structures provides insights on the chemical shifts of the first shell ligands around the metals with respect to inter-metal distances.