This is an archive of papers published by the staff and faculty of Fox Chase Cancer Center. For questions about content, please contact Talbot Research Library
Last updated on
Haldane A , Levy RM
Mi3-GPU: MCMC-based Inverse Ising Inference on GPUs for protein covariation analysis
Comput Phys Commun. 2021 Mar;260
PMID: 33716309 PMCID: PMC7944406 URL: https://www.ncbi.nlm.nih.gov/pubmed/33716309
AbstractInverse Ising inference is a method for inferring the coupling parameters of a Potts/Ising model based on observed site-covariation, which has found important applications in protein physics for detecting interactions between residues in protein families. We introduce Mi3-GPU ("mee-three", for MCMC Inverse Ising Inference) software for solving the inverse Ising problem for protein-sequence datasets with few analytic approximations, by parallel Markov-Chain Monte-Carlo sampling on GPUs. We also provide tools for analysis and preparation of protein-family Multiple Sequence Alignments (MSAs) to account for finite-sampling issues, which are a major source of error or bias in inverse Ising inference. Our method is "generative" in the sense that the inferred model can be used to generate synthetic MSAs whose mutational statistics (marginals) can be verified to match the dataset MSA statistics up to the limits imposed by the effects of finite sampling. Our GPU implementation enables the construction of models which reproduce the covariation patterns of the observed MSA with a precision that is not possible with more approximate methods. The main components of our method are a GPU-optimized algorithm to greatly accelerate MCMC sampling, combined with a multi-step Quasi-Newton parameter-update scheme using a "Zwanzig reweighting" technique. We demonstrate the ability of this software to produce generative models on typical protein family datasets for sequence lengths L ~ 300 with 21 residue types with tens of millions of inferred parameters in short running times.
NotesHaldane, Allan Levy, Ronald M P50 GM103368/GM/NIGMS NIH HHS/United States R35 GM132090/GM/NIGMS NIH HHS/United States U54 GM103368/GM/NIGMS NIH HHS/United States Journal Article Comput Phys Commun. 2021 Mar;260:107312. doi: 10.1016/j.cpc.2020.107312. Epub 2020 Apr 17.