Minimization-Aware Recursive K (MARK ): A Novel, Provable Algorithm that Accelerates Ensemble-Based Protein Design and Provably Approximates the Energy Landscape

Published

Conference Paper

© 2019, Springer Nature Switzerland AG. Protein design algorithms that model continuous sidechain flexibility and conformational ensembles better approximate the in vitro and in vivo behavior of proteins. The previous state of the art, iMinDEE- A∗ - K∗, computes provable ε -approximations to partition functions of protein states (e.g., bound vs. unbound) by computing provable, admissible pairwise-minimized energy lower bounds on protein conformations and using the A∗ enumeration algorithm to return a gap-free list of lowest-energy conformations. iMinDEE-A ∗ - K∗ runs in time sublinear in the number of conformations, but can be trapped in loosely-bounded, low-energy conformational wells containing many conformations with highly similar energies. That is, iMinDEE- A∗ - K∗ is unable to exploit the correlation between protein conformation and energy: similar conformations often have similar energy. We introduce two new concepts that exploit this correlation: Minimization-Aware Enumeration and Recursive K∗. We combine these two insights into a novel algorithm, Minimization-Aware Recursive K∗ (MARK∗ ), that tightens bounds not on single conformations, but instead on distinct regions of the conformation space. We compare the performance of iMinDEE- A∗ - K∗ vs. MARK∗ by running the BBK∗ algorithm, which provably returns sequences in order of decreasing K∗ score, using either iMinDEE- A∗ - K∗ or MARK∗ to approximate partition functions. We show on 200 design problems that MARK∗ not only enumerates and minimizes vastly fewer conformations than the previous state of the art, but also runs up to two orders of magnitude faster. Finally, we show that MARK∗ not only efficiently approximates the partition function, but also provably approximates the energy landscape. To our knowledge, MARK∗ is the first algorithm to do so. We use MARK∗ to analyze the change in energy landscape of the bound and unbound states of the HIV-1 capsid protein C-terminal domain in complex with camelid V H H, and measure the change in conformational entropy induced by binding. Thus, MARK∗ both accelerates existing designs and offers new capabilities not possible with previous algorithms.

Full Text

Cited Authors

  • Jou, JD; Holt, GT; Lowegard, AU; Donald, BR

Published Date

  • January 1, 2019

Published In

Volume / Issue

  • 11467 LNBI /

Start / End Page

  • 101 - 119

Electronic International Standard Serial Number (EISSN)

  • 1611-3349

International Standard Serial Number (ISSN)

  • 0302-9743

International Standard Book Number 13 (ISBN-13)

  • 9783030170820

Digital Object Identifier (DOI)

  • 10.1007/978-3-030-17083-7_7

Citation Source

  • Scopus