No abstract

No abstract

PHD fingers represent one of the largest families of epigenetic readers capable of decoding post-translationally modified or unmodified histone H3 tails. Because of their direct involvement in human pathologies they are increasingly considered as a potential therapeutic target. Several PHD/histone-peptide structures have been determined, however relatively little information is available on their dynamics. Studies aiming to characterize the dynamic and energetic determinants driving histone peptide recognition by epigenetic readers would strongly benefit from computational studies. Herein we focus on the dynamic and energetic characterization of the PHD finger subclass specialized in the recognition of histone H3 peptides unmodified in position K4 (H3K4me0). As a case study we focused on the first PHD finger of autoimmune regulator protein (AIRE-PHD1) in complex with H3K4me0. PCA analysis of the covariance matrix of free AIRE-PHD1 highlights the presence of a “flapping” movement, which is blocked in an open conformation upon binding to H3K4me0. Moreover, binding free energy calculations obtained through Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) methodology are in good qualitative agreement with experiments and allow dissection of the energetic terms associated with native and alanine mutants of AIRE-PHD1/H3K4me0 complexes. MM/PBSA calculations have also been applied to the energetic analysis of other PHD fingers recognizing H3K4me0. In this case we observe excellent correlation between computed and experimental binding free energies. Overall calculations show that H3K4me0 recognition by PHD fingers relies on compensation of the electrostatic and polar solvation energy terms and is stabilized by non-polar interactions.

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.

No abstract