The identification of intermediate states for folding and aggregation is important from a fundamental standpoint and for the design of novel therapeutic strategies targeted at conformational disorders. Protein human β2-microglobulin (HB2m) is classically associated with dialysis-related amyloidosis, but the single point mutant D76N was recently identified as the causative agent of a hereditary systemic amyloidosis affecting visceral organs. Here, we use D76N as a model system to explore the early stage of the aggregation mechanism of HB2m by means of an integrative approach framed on molecular simulations. Discrete molecular dynamics simulations of a structured-based model predict the existence of two intermediate states populating the folding landscape. The intermediate I1 features an unstructured C-terminus, while I2 , which is exclusively populated by the mutant, exhibits two unstructured termini. Docking simulations indicate that I2 is the key species for aggregation at acidic and physiological pH contributing to rationalize the higher amyloidogenic potential of D76N relative to the wild-type protein and the ΔN6 variant. The analysis carried out here recapitulates the importance of the DE-loop in HB2m self-association at a neutral pH and predicts a leading role of the C-terminus and the adjacent G-strand in the dimerization process under acidic conditions. The identification of aggregation hot-spots is in line with experimental results that support the importance of Phe56, Asp59, Trp60, Phe62, Tyr63, and Tyr66 in HB2m amyloidogenesis. We further predict the involvement of new residues such as Lys94 and Trp95 in the aggregation process.
Keywords: aggregation; amyloids; biophysics; computational biology; computer modeling; conformational change; conformational disorders; protein folding.
© 2017 Wiley Periodicals, Inc.