q-bio.QM q-bio.BM
Weighted persistent homology for osmolyte molecular aggregation and hydrogen bonding network analysis. It has long been observed that trimethylamin N oxide (TMAO) and urea demonstrate dramatically different properties in a protein folding process. Even with the enormous theoretical and experimental research work of the two osmolytes, various aspects of their underlying mechanisms still remain largely elusive. In this paper, we propose to use the weighted persistent homology to systematically study the osmolytes molecular aggregation and their hydrogen bonding network from a local topological perspective. We consider two weighted models, i.e., localized persistent homology (LPH) and interactive persistent homology (IPH). From the localized persistent homology models, we have found that TMAO and urea have very different local topology. TMAO shows local network structures. With the concentration increase, the circle elements in these networks show a clear increase in their total numbers and a decrease in their relative sizes. In contrast, urea shows two types of local topological patterns, i.e., local clusters around 6 AA and a few global circle elements at around 12 AA. From the interactive persistent homology models, it has been found that our persistent radial distribution function (PRDF) from the global scale IPH has same physical properties as the traditional radial distribution function (RDF). Moreover, PRDFs from the local scale IPH can also be generated and used to characterize the local interaction information. Other than the clear difference of the first peak value of PRDFs at filtration size 4 AA, TMAO and urea also shows very different behaviors at the second peak region from filtration size 5 AA to 10 AA.