A Markov-chain model description of binding funnels to enhance the ranking of docked solutions

Abstract

Within the crowded, seemingly chaotic environment of the cell, proteins are still able to find their binding partners. This is achieved via an ensemble of trajectories, which funnel them towards their functional binding sites, the binding funnel. Here, we characterize funnel-like energy structures on the global energy landscape using time-homogeneous finite state Markov chain models. These models are based on the idea that transitions can occur between structurally similar docking solutions, with transition probabilities determined by their difference in binding energy. Funnel-like energy structures are those containing solutions with very high equilibrium populations. Although these are found surrounding both near-native and false positive binding sites, we show that the removal of nonfunnel-like energy structures, by filtering away solutions with low maximum equilibrium population, can significantly improve the ranking of docked poses.