[Paper Review] Atomic-resolution structures of prion AGAAAAGA amyloid fibrils

Jiapu Zhang|arXiv (Cornell University)|Jun 9, 2011

Prion Diseases and Protein Misfolding14 references3 citations

TL;DR

This study computationally determines atomic-resolution structures of prion AGAAAAGA amyloid fibrils using advanced optimization techniques, including steepest descent, conjugate gradient, simulated annealing, and genetic algorithms, combined with canonical dual optimization and structural bioinformatics. The key contribution is the identification of stable, energetically favorable fibril models for the AGAAAAGA palindrome (residues 113–120), a critical region in prion disease pathogenesis that resists traditional experimental characterization due to its insoluble, noncrystalline nature.

ABSTRACT

To the best of the author's knowledge, there is little structural data available on the AGAAAAGA palindrome in the hydrophobic region (113-120) of prion proteins due to the unstable, noncrystalline and insoluble nature of the amyloid fibril, although many experimental studies have shown that this region has amyloid fibril forming properties and plays an important role in prion diseases. In view of this, the present study is devoted to address this problem from computational approaches such as local optimization steepest descent, conjugate gradient, discrete gradient and Newton methods, global optimization simulated annealing and genetic algorithms, canonical dual optimization theory, and structural bioinformatics. The optimal atomic-resolution structures of prion AGAAAAGA amyloid fibils reported in this Chapter have a value to the scientific community in its drive to find treatments for prion diseases or at least be useful for the goals of medicinal chemistry.

Motivation & Objective

• To address the lack of experimental structural data for the AGAAAAGA palindrome (113–120) in prion proteins due to its insoluble, noncrystalline, and unstable nature.
• To develop and apply novel computational methods to predict stable, atomic-resolution structures of prion amyloid fibrils where traditional X-ray and NMR methods fail.
• To identify the minimal fibril-forming sequence and confirm the role of AGAAAAGA in amyloidogenesis and prion toxicity.
• To provide high-resolution structural models useful for medicinal chemistry and drug design targeting prion diseases.

Proposed method

• Employed local optimization methods including steepest descent (SD), conjugate gradient (CG), discrete gradient, and Newton methods for energy minimization.
• Applied global optimization techniques such as simulated annealing (SA) and genetic algorithms to escape local minima and locate global energy minima.
• Used canonical dual optimization theory to reformulate the molecular structure determination as a global optimization problem.
• Constructed initial fibril models using templates from the Protein Data Bank (PDB), including 2OMP.pdb, 1YJP.pdb, 3FVA.pdb, and 3NHC–3NVH.pdb, via molecular modeling tools (Insight II and Swiss-PdbViewer).
• Refined models using AMBER 10 with hybrid SDCG-SA protocols to achieve structural stability and equilibrium.
• Formulated the molecular structure determination as a distance geometry problem (DGP) and solved it via a global optimization framework minimizing the function $ P_{\epsilon}(X) = \sum_{(i,j)\in S} w_{ij}(||x_i - x_j||^2 - d_{ij}^2)^2 - \epsilon^T X $, with $ d_{ij} $ defined as the sum of van der Waals radii for non-bonded atoms.

Experimental results

Research questions

• RQ1 What is the atomic-resolution structure of the prion AGAAAAGA amyloid fibril segment (residues 113–120), which is critical for prion disease pathogenesis but inaccessible to X-ray and NMR methods?
• RQ2 Which computational optimization strategies can reliably predict stable, low-energy fibril conformations for this highly insoluble and noncrystalline segment?
• RQ3 How do van der Waals and hydrogen bonding interactions contribute to the stability of the predicted fibril models?
• RQ4 Can the minimal fibril-forming sequence within AGAAAAGA be identified and validated through computational energy minimization and structural analysis?
• RQ5 What is the role of the AGAAAAGA palindrome in mediating prion fibril formation and neurotoxicity, based on structural modeling?

Key findings

• The AGAAAAGA segment (113–120) was computationally confirmed as a strong amyloid fibril-forming region, with energy below the threshold of -26, indicating high fibrilogenic potential.
• Twenty distinct atomic-resolution models (MODEL01–MODEL19) were generated, including 12-chain AGAAAA (MODEL01), 10-chain AGAAAAG (MODEL02), and 10-chain GAAAAGA (MODEL03) structures, all stabilized via SDCG-SA refinement.
• Models based on 3NHC, 3NVE, 3NVF, 3NVG, and 3NVH templates yielded stable fibril structures after mutation and optimization, with specific van der Waals contact pairs preserved (e.g., A.ALA3.CB–G.ALA4.CB in 3NHC models).
• The 3NHC, 3NVE/F/G/H models were successfully refined into stable, low-energy configurations, with MODEL06–MODEL19 representing the final optimized fibril structures.
• The study demonstrates that mathematical optimization frameworks can effectively model amyloid fibrils where experimental methods fail, providing high-resolution models for drug targeting.
• The results validate that AGAAAAGA is not only necessary for fibril formation but also capable of inhibiting PrPSc neurotoxicity, supporting its role as a therapeutic target.