A Binary Matrix Method to Enumerate, Hierarchically Order and Structurally Classify Peptide Aggregation
Protein aggregation is a common and complex phenomenon in biological processes, yet a robust analysis of this aggregation process remains elusive. The commonly used methods such as centre-of-mass to centre-of-mass (COM-COM) distance, the radius of gyration (Rg), hydrogen bonding (HB) and solvent accessible surface area (SASA) do not quantify the aggregation accurately. Herein, a new and robust method that uses an aggregation matrix (AM) approach to investigate peptide aggregation in a MD simulation trajectory is presented. A nxn two-dimensional aggregation matrix (AM) is created by using the inter-peptide CA-CA cut-off distances which are binarily encoded (0 or 1). These aggregation matrices are analysed to enumerate, hierarchically order and structurally classify the aggregates. Moreover, the comparison between the present AM method and the conventional Rg, COM-COM and HB methods shows that the conventional methods grossly underestimate the aggregation propensity. Additionally, the conventional methods do not address the hierarchy and structural ordering of the aggregates, which the present AM method does. Finally, the present AM method utilises only nxn two-dimensional matrices to analyse aggregates consisting of several peptide units. To the best of our knowledge, this is a maiden approach to enumerate, hierarchically order and structurally classify peptide aggregation.