No abstract

No abstract

Binding affinity prediction of potential drugs to target and off-target proteins is an essential asset in drug development. These predictions require the calculation of binding free energies. In such calculations, it is a major challenge to properly account for both the dynamic nature of the protein and the possible variety of ligand-binding orientations, while keeping computational costs tractable. Recently, an iterative Linear Interaction Energy (LIE) approach was introduced, in which results from multiple simulations of a protein-ligand complex are combined into a single binding free energy using a Boltzmann weighting-based scheme. This method was shown to reach experimental accuracy for flexible proteins while retaining the computational efficiency of the general LIE approach. Here, we show that the iterative LIE approach can be used to predict binding affinities in an automated way. A workflow was designed using preselected protein conformations, automated ligand docking and clustering, and a (semi-)automated molecular dynamics simulation setup. We show that using this workflow, binding affinities of aryloxypropanolamines to the malleable Cytochrome P450 2D6 enzyme can be predicted without a priori knowledge of dominant protein-ligand conformations. In addition, we provide an outlook for an approach to assess the quality of the LIE predictions, based on simulation outcomes only.

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.

No abstract