CYP 2D6 binding affinity predictions using multiple ligand and protein conformations

Citations of this article
Mendeley users who have this article in their library.


Because of the large flexibility and malleability of Cytochrome P450 enzymes (CYPs), in silico prediction of CYP binding affinities to drugs and other xenobiotic compounds is a true challenge. In the current work, we use an iterative linear interaction energy (LIE) approach to compute CYP binding affinities from molecular dynamics (MD) simulation. In order to improve sampling of conformational space, we combine results from simulations starting with different relevant protein-ligand geometries. For calculated binding free energies of a set of thiourea compounds binding to the flexible CYP 2D6 isoform, improved correlation with experiment was obtained by combining results of MD runs starting from distinct protein conformations and ligand-binding orientations. This accuracy was obtained from relatively short MD simulations, which makes our approach computationally attractive for automated calculations of ligand-binding affinities to flexible proteins such as CYPs.




Perić-Hassler, L., Stjernschantz, E., Oostenbrink, C., & Geerke, D. P. (2013). CYP 2D6 binding affinity predictions using multiple ligand and protein conformations. International Journal of Molecular Sciences, 14(12), 24514–24530.

Register to see more suggestions

Mendeley helps you to discover research relevant for your work.

Already have an account?

Save time finding and organizing research with Mendeley

Sign up for free