CYP 2D6 Binding Affinity Predictions Using Multiple Ligand and Protein Conformations

International Journal of Molecular Sciences

Lovorka Perić-Hassler¹  Eva Stjernschantz¹  Chris Oostenbrink² 

[1] AIMMS Division of Molecular Toxicology, Department of Chemistry and Pharmaceutical Sciences, Faculty of Sciences, VU University Amsterdam, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands; E-Mail:;Institute of Molecular Modeling and Simulation, University of Natural Resources and Life Sciences, Muthgasse 18, 1190 Vienna, Austria; E-Mail:

关键词: binding free energies;    LIE method;    CYP 2D6;    thiourea ligands;   

DOI  :  10.3390/ijms141224514

来源: mdpi

【 摘 要 】

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 ofMDruns 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.

【 授权许可】

CC BY   
© 2013 by the authors; licensee MDPI, Basel, Switzerland