Dimitrios Spiliotopoulos

Discovery of CREBBP Bromodomain Inhibitors by High-Throughput Docking and Hit Optimization Guided by Molecular Dynamics.

We have identified two chemotypes of CREBBP bromodomain ligands by fragment-based high-throughput docking. Only 17 molecules from the original library of two-million compounds were tested in vitro. Optimization of the two low-micromolar hits, the 4-acylpyrrole 1 and acylbenzene 9, was driven by molecular dynamics results which suggested improvement of the polar interactions with the Arg1173 side chain at the rim of the binding site. The synthesis of only two derivatives of 1 yielded the 4-acylpyrrole 6 which shows a single-digit micromolar affinity for the CREBBP bromodomain and a ligand efficiency of 0.34 kcal/mol per non-hydrogen atom. Optimization of the acylbenzene hit 9 resulted in a series of derivatives with nanomolar potencies, good ligand efficiency and selectivity (see Unzue, A.; Xu, M.; Dong, J.; Wiedmer, L.; Spiliotopoulos, D.; Caflisch, A.; Nevado, C.Fragment-Based Design of Selective Nanomolar Ligands of the CREBBP Bromodomain. J. Med. Chem. 2015, DOI: 10.1021/acs.jmedchem.5b00172). The in silico predicted binding mode of the acylbenzene derivative 10 was validated by solving the structure of the complex with the CREBBP bromodomain.

Discovery of BRD4 bromodomain inhibitors by fragment-based high-throughput docking.

Bromodomains (BRDs) recognize acetyl-lysine modified histone tails mediating epigenetic processes. BRD4, a protein containing two bromodomains, has emerged as an attractive therapeutic target for several types of cancer as well as inflammatory diseases. Using a fragment-based in silico screening approach, we identified two small molecules that bind to the first bromodomain of BRD4 with low-micromolar affinity and favorable ligand efficiency (0.37 kcal/mol per non-hydrogen atom), selectively over other families of bromodomains. Notably, the hit rate of the fragment-based in silico approach is about 10% as only 24 putative inhibitors, from an initial library of about 9 million molecules, were tested in vitro.

Fragment-Based Design of Selective Nanomolar Ligands of the CREBBP Bromodomain

Novel ligands of the CREBBP bromodomain were identified by fragment-based docking. The in silico discovered hits have been optimized by chemical synthesis into selective nanomolar compounds, thereby preserving the ligand efficiency. The selectivity for the CREBBP bromodomain over other human bromodomain subfamilies has achieved by a benzoate moiety which was predicted by docking to be involved in favorable electrostatic interactions with the Arg1173 side chain, a prediction that could be verified a posteriori by the high-resolution crystal structure of the CREBBP bromodomain in complex with ligand 6 and also by MD simulations (see Xu, M.; Unzue, A.; Dong, J.; Spiliotopoulos, D.; Nevado, C.; Caflisch, A. Discovery of CREBBP bromodomain inhibitors by high-throughput docking and hit optimization guided by molecular dynamics. J. Med. Chem. 2015, DOI: 10.1021/acs.jmedchem.5b00171).

dMM-PBSA: A New HADDOCK Scoring Function for Protein-Peptide Docking

Molecular-docking programs coupled with suitable scoring functions are now established and very useful tools enabling computational chemists to rapidly screen large chemical databases and thereby to identify promising candidate compounds for further experimental processing. In a broader scenario, predicting binding affinity is one of the most critical and challenging components of computer-aided structure-based drug design. The development of a molecular docking scoring function which in principle could combine both features, namely ranking putative poses and predicting complex affinity, would be of paramount importance. Here, we systematically investigated the performance of the MM-PBSA approach, using two different Poisson–Boltzmann solvers (APBS and DelPhi), in the currently rising field of protein-peptide interactions (PPIs), identifying the correct binding conformations of 19 different protein-peptide complexes and predicting their binding free energies. First, we scored the decoy structures from HADDOCK calculation via the MM-PBSA approach in order to assess the capability of retrieving near-native poses in the best-scoring clusters and of evaluating the corresponding free energies of binding. MM-PBSA behaves well in finding the poses corresponding to the lowest binding free energy, however the built-in HADDOCK score shows a better performance. In order to improve the MM-PBSA-based scoring function, we dampened the MM-PBSA solvation and coulombic terms by 0.2, as proposed in the HADDOCK score and LIE approaches. The new dampened MM-PBSA (dMM-PBSA) outperforms the original MM-PBSA and ranks the decoys structures as the HADDOCK score does. Second, we found a good correlation between the dMM-PBSA and HADDOCK scores for the near-native clusters of each system and the experimental binding energies, respectively. Therefore, we propose a new scoring function, dMM-PBSA, to be used together with the built-in HADDOCK score in the context of protein-peptide docking simulations.

Virtual screen to NMR (VS2NMR): Discovery of fragment hits for the CBP bromodomain.

Overexpression of the CREB-binding protein (CBP), a bromodomain-containing transcription coactivator involved in a variety of cellular processes, has been observed in several types of cancer with a correlation to aggressiveness. We have screened a library of nearly 1500 fragments by high-throughput docking into the CBP bromodomain followed by binding energy evaluation using a force field with electrostatic solvation. Twenty of the 39 fragments selected by virtual screening are positive in one or more ligand-observed nuclear magnetic resonance (NMR) experiments. Four crystal structures of the CBP bromodomain in complex with in silico screening hits validate the pose predicted by docking. Thus, the success ratio of the high-throughput docking procedure is 50% or 10% if one considers the validation by ligand-observed NMR spectroscopy or X-ray crystallography, respectively. Compounds 1 and 3 show favorable ligand efficiency in two different in vitro binding assays. The structure of the CBP bromodomain in the complex with the brominated pyrrole 1 suggests fragment growing by Suzuki coupling.

Structural Analysis of the Binding of Type I, I1/2, and II Inhibitors to Eph Tyrosine Kinases.

We have solved the crystal structures of the EphA3 tyrosine kinase in complex with nine small-molecule inhibitors, which represent five different chemotypes and three main binding modes, i.e., types I and I1/2 (DFG in) and type II (DFG out). The three structures with type I1/2 inhibitors show that the higher affinity with respect to type I is due to an additional polar group (hydroxyl or pyrazole ring of indazole) which is fully buried and is involved in the same hydrogen bonds as the (urea or amide) linker of the type II inhibitors. Overall, the type I and type II binding modes belong to the lock-and-key and induced fit mechanism, respectively. In the type II binding, the scaffold in contact with the hinge region influences the position of the Phe765 side chain of the DFG motif and the orientation of the Gly-rich loop. The binding mode of Birb796 in the EphA3 kinase does not involve any hydrogen bond with the hinge region, which is different from the Birb796/p38 MAP kinase complex. Our structural analysis emphasizes the importance of accounting for structural plasticity of the ATP binding site in the design of type II inhibitors of tyrosine kinases.

Fragment-based in silico screening of bromodomain ligands.

We review the results of fragment-based high-throughput docking to the N-terminal bromodomain of BRD4 and the CREBBP bromodomain. In both docking campaigns the ALTA (anchor-based library tailoring) procedure was used to reduce the size of the initial library by selecting for flexible docking only the molecules that contain a fragment with favorable predicted binding energy. Ranking by a force field-based energy with solvation has resulted in small-molecule hits with low-micromolar affinity and favorable ligand efficiency. Importantly, the binding modes predicted by docking have been validated by X-ray crystallography. One of the hits for the CREBBP bromodomain has been optimized by medicinal chemistry into a series of potent and selective ligands.

Chemical Space Expansion of Bromodomain Ligands Guided by in Silico Virtual Couplings (AutoCouple)

Expanding the chemical space and simultaneously ensuring synthetic accessibility is of upmost importance, not only for the discovery of effective binders for novel protein classes but, more importantly, for the development of compounds against hard-to-drug proteins. Here, we present AutoCouple, a de novo approach to computational ligand design focused on the diversity-oriented generation of chemical entities via virtual couplings. In a benchmark application, chemically diverse compounds with low-nanomolar potency for the CBP bromodomain and high selectivity against the BRD4(1) bromodomain were achieved by the synthesis of about 50 derivatives of the original fragment. The binding mode was confirmed by X-ray crystallography, target engagement in cells was demonstrated, and antiproliferative activity was showcased in three cancer cell lines. These results reveal AutoCouple as a useful in silico coupling method to expand the chemical space in hit optimization campaigns resulting in potent, selective, and cell permeable bromodomain ligands.