Boltzmann Machine Learning and Regularization Methods for Inferring Evolutionary Fields and Couplings From a Multiple Sequence Alignment
| Autor: | Sanzo Miyazawa |
|---|---|
| Rok vydání: | 2022 |
| Předmět: |
FOS: Computer and information sciences
Computer Science - Machine Learning 0206 medical engineering Monte Carlo method Boltzmann machine FOS: Physical sciences Inverse Machine Learning (stat.ML) 02 engineering and technology Regularization (mathematics) Machine Learning (cs.LG) Machine Learning Statistics - Machine Learning Genetics Amino Acid Sequence Statistical physics Quantitative Biology - Populations and Evolution Condensed Matter - Statistical Mechanics Mathematics Quantitative Biology::Biomolecules Multiple sequence alignment Statistical Mechanics (cond-mat.stat-mech) Applied Mathematics Populations and Evolution (q-bio.PE) Proteins Biomolecules (q-bio.BM) Function (mathematics) Boltzmann distribution Quantitative Biology - Biomolecules FOS: Biological sciences Pairwise comparison Sequence Alignment 020602 bioinformatics Biotechnology |
| Zdroj: | IEEE/ACM Transactions on Computational Biology and Bioinformatics. 19:328-342 |
| ISSN: | 2374-0043 1545-5963 |
| Popis: | The inverse Potts problem to infer a Boltzmann distribution for homologous protein sequences from their single-site and pairwise amino acid frequencies recently attracts a great deal of attention in the studies of protein structure and evolution. We study regularization and learning methods and how to tune regularization parameters to correctly infer interactions in Boltzmann machine learning. Using $L_2$ regularization for fields, group $L_1$ for couplings is shown to be very effective for sparse couplings in comparison with $L_2$ and $L_1$. Two regularization parameters are tuned to yield equal values for both the sample and ensemble averages of evolutionary energy. Both averages smoothly change and converge, but their learning profiles are very different between learning methods. The Adam method is modified to make stepsize proportional to the gradient for sparse couplings and to use a soft-thresholding function for group $L_1$. It is shown by first inferring interactions from protein sequences and then from Monte Carlo samples that the fields and couplings can be well recovered, but that recovering the pairwise correlations in the resolution of a total energy is harder for the natural proteins than for the protein-like sequences. Selective temperature for folding/structural constrains in protein evolution is also estimated. In this version (arXiv:1909.05006v3), the values of selective temperature for protein PF00153, $T_s$ in Table 5 and in the section 2.8, and folding free energy for PF00595 are corrected and shown in red. The version 2 (arXiv:1909.05006v2) was published in the IEEE/ACM Transactions on Computational Biology and Bioinformatics. The program is available from https://gitlab.com/sanzo.miyazawa/BM/ |
| Databáze: | OpenAIRE |
| Externí odkaz: |
|