Exploration of 3D activity cliffs on the basis of compound binding modes and comparison of 2D and 3D cliffs

Activity cliffs are formed by pairs or groups of structurally similar compounds having large differences in potency and are focal points of structure-activity relationship (SAR) analysis. The choice of molecular representations is a critically important aspect of activity cliffs analysis. Thus far, activity cliffs have predominantly been defined on the basis of molecular graph or fingerprint representations. Herein we introduce 3D activity cliffs derived from comparisons of experimentally determined compound binding modes. The analysis of 3D activity cliffs is generally applicable to target proteins for which structures of multiple ligand complexes are available. For two popular targets, β-secretase 1 (BACE1) and factor Xa (FXa), public domain X-ray structures with bound inhibitors were collected. Crystallographic binding modes of inhibitors were systematically compared using a 3D similarity method taking conformational, positional, and atomic property differences into account. In addition, standard 2D similarity relationships were also determined. SAR information associated with individual compounds substantially changed when either bioactive conformations or 2D molecular graphs were used for similarity evaluation. 3D activity cliffs were identified for BACE1 and FXa inhibitor sets and systematically compared to 2D cliffs. It was found that less than 40% of 3D activity cliffs were conserved when 2D similarity was applied. The limited conservation of 3D and 2D cliffs provides further evidence for the strong molecule representation dependence of activity cliffs. Moreover, 3D cliffs represent a new class of activity cliffs that convey SAR information in ways that differ from graph-based similarity measures. In cases where sufficient structural information is available, the comparison of 3D and 2D cliffs is expected to aid in SAR analysis and mapping of critical binding determinants.