Skip to main content
  • Poster presentation
  • Open access
  • Published:

IterTunnel; a method for predicting and evaluating ligand EgressTunnels in proteins with buried active sites

The computational prediction of ligand entry and egress paths in proteins has become an emerging topic in computational biology due to the potential for estimating kinetic properties of drug binding. These properties are related to important pharmacological quantities such as the kon and koff rate of drugs [1, 2].We have investigated the influence of protein flexibility on tunnel prediction using geometric methods by comparing tunnels identified in static structures with those found in structural ensembles of three CYP isozymes. We found drastic differences between tunnels predicted in the crystal structures as opposed to those predicted in the ensembles [3]. Furthermore, we found significant differences between tunnels identified in the apo versus the holo protein ensembles [3].

While geometric prediction provides a good starting point for tunnel prediction, in order to estimate kinetic properties, more detailed investigations of the ligand binding process are required. We have developed a tunnel prediction methodology, IterTunnel, which predicts tunnels in proteins and estimates the free energy of ligand unbinding using a combination of geometric tunnel prediction with steered molecular dynamics and umbrella sampling [4]. Applying this new method to cytochrome P450 2B6 (CYP2B6), we demonstrate that the ligand itself plays an important role in reshaping tunnels as it traverses through a protein. This process results in the exposure of new tunnels and the closure of pre-existing tunnels as the ligand migrates from the active site. We found that many of the tunnels that are exposed due to ligand-induced conformational changes are amongst the most energetically favorable tunnels for ligand egress in CYP2B6 [4].


  1. Colizzi F, Perozzo R, Scapozza L, Recanatini M, Cavalli A: Molecule Pulling Simulations Can Discern Active from Inactive Enzyme Inhibitors. J Am Chem Soc. 2010, 132: 7361-7371. 10.1021/ja100259r.

    Article  CAS  Google Scholar 

  2. Jorgensen WL: Drug discovery: Pulled from a protein's embrace. Nature. 2010, 466: 42-43. 10.1038/466042a.

    Article  CAS  Google Scholar 

  3. Kingsley LJ, Lill MA: The Influence of Protein Flexibility on Geometric Tunnel Prediction in Cytochrome P450 Enzymes. Submitted

  4. Kingsley LJ, Lill MA: Including Protein Flexibility Induced by Ligand Migration into Protein Tunnel Prediction. In preparation

Download references

Author information

Authors and Affiliations


Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits use, duplication, adaptation, distribution, and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Reprints and permissions

About this article

Cite this article

Kingsley, L.J., Lill, M.A. IterTunnel; a method for predicting and evaluating ligand EgressTunnels in proteins with buried active sites. J Cheminform 6 (Suppl 1), P62 (2014).

Download citation

  • Published:

  • DOI: