Cl-π Interactions in Protein-Ligand Complexes

Experimental and computational evidence for a stabilising C–Cl(lone-pair)⋯π(chelate-ring) interaction

Stabilising C–Cl(lone-pair)⋯π(chelate ring) interactions are described.

Theoretical investigation of AhR binding property with relevant structural requirements for AhR-mediated toxicity of polybrominated diphenyl ethers

Polybrominated diphenyl ethers (PBDEs) are more frequently suspected with the induction of toxicity via signal transduction pathway of cytosolic aryl hydrocarbon receptor (AhR), the initial binding to which is assumed to be an essential prerequisite during the ligand-dependent activation. However, the AhR binding property and associated toxicity of PBDEs is yet to be clearly known for lacking insights into the structural requirements at molecular level. To understand the AhR binding property of PBDEs, the ligand binding domain (LBD) of AhR was simulatively developed on homologous protein after basic validation of geometrical rationality and the binding interaction profile was visually described using molecular docking approach. For AhR binding, the offset or edge-on π-π stackings with aromatic motifs including Phe289, Phe345 and His285 were shown to be structurally required whereas the electrostatic attraction validated for AhR binding to dioxins might be less effective for 2,2',3,4,4'-pentabromodiphenyl ether (BDE-85). Besides the demands of less steric hindrance from alanines and weak formulation of hydrogen bonds, the dispersion force through large contact and polarization of S-π electrons seemed to be impactful when BDE-85 were closer to Cys327, Met334 or Met342. With theoretical computation of AhR binding energies, the more significant correlativity with bioassays was derived especially for the lowly/moderately brominated congeners, and could be used to predict the AhR binding affinity on certain degree. The informative results would thus not only help well understand the molecular basis of AhR-mediated toxicity but give an approach for accelerative evaluation of AhR binding and toxicity of PBDEs.

Synthesis, Structural Characterization, and Crystal Packing of the Elusive Pentachlorinated Azafullerene C59 NCl5

We report on the synthesis and the structure elucidation of the elusive azafullerene pentachloride C₅₉ NCl₅ , which was obtained by high temperature halogenation of (C₅₉ N)₂ . The exceptionally strong host-guest interaction of the title compound in the solid is discussed.

Halogen and Hydrogen Bonding Benzothiophene Diol Derivatives: A Study Using ab initio Calculations and X-Ray Crystal Structure Measurements

The aim of this study is to describe and compare the supramolecular interactions, in the solid state, of chloro-, bromo-, and iodobenzothiophene diols. The compounds were obtained through organo-catalyzed reactions starting from 3-substituted halobenzothiophene carbaldehydes. Energies of the noncovalent interactions were obtained by density functional theory calculations. Bond distances and angles were found to be in accordance with those determined by X-ray structure analysis. anti-Bromobenzothiophene derivatives showed strong halogen⋅⋅⋅π interactions between bromine and the heterocyclic phenyl ring, corresponding to an energy of 7.5 kcal mol(-1). syn-Bromo and syn-iodo derivatives appeared to be isostructural, showing X⋅⋅⋅O (carbonyl) interactions, π stacking, and formation of extended hydrogen bonding networks. In contrast, the chloro derivatives displayed no halogen bonding interactions.

High-resolution crystal structures of factor XIa coagulation factor in complex with nonbasic high-affinity synthetic inhibitors

Factor XI (FXI) is a key enzyme in the coagulation pathway and an attractive target for the development of anticoagulant drugs. A small number of high-resolution crystal structures of FXIa in complex with small synthetic inhibitors have been published to date. All of these ligands have a basic P1 group and bind exclusively in the nonprime side of the active site of FXIa. Here, two structures of FXIa in complex with nonbasic inhibitors that occupy both the prime and nonprime sides of the active site are presented. These new structures could be valuable in the design and optimization of new FXIa synthethic inhibitors.

Compound Design by Fragment-Linking

The linking together of two fragment compounds that bind to distinct protein sub-sites can lead to a superadditivity of binding affinities, in which the binding free energy of the linked fragments exceeds the simple sum of the binding energies of individual fragments (linking coefficient E<1). However, a review of the literature shows that such events are relatively rare and, in the majority of the cases, linking coefficients are far from optimal being much greater than 1. It is critical to design a linker that does not disturb the original binding poses of each fragment in order to achieve successful linking. However, such an ideal linker is often difficult to design and even more difficult to actually synthesize. We suggest that the chance of achieving successful fragment linking can be significantly improved by choosing a fragment pair that consists of one fragment that binds by strong H-bonds (or non-classical equivalents) and a second fragment that is more tolerant of changes in binding mode (hydrophobic or vdW binders). We also propose that the fragment molecular orbital (FMO) calculations can be used to analyse the nature of the binding interactions of the fragment hits for the selection of fragments for evolution, merging and linking in order to optimize the chance of success.

Prediction of cyclin-dependent kinase 2 inhibitor potency using the fragment molecular orbital method

BACKGROUND

The reliable and robust estimation of ligand binding affinity continues to be a challenge in drug design. Many current methods rely on molecular mechanics (MM) calculations which do not fully explain complex molecular interactions. Full quantum mechanical (QM) computation of the electronic state of protein-ligand complexes has recently become possible by the latest advances in the development of linear-scaling QM methods such as the ab initio fragment molecular orbital (FMO) method. This approximate molecular orbital method is sufficiently fast that it can be incorporated into the development cycle during structure-based drug design for the reliable estimation of ligand binding affinity. Additionally, the FMO method can be combined with approximations for entropy and solvation to make it applicable for binding affinity prediction for a broad range of target and chemotypes.

RESULTS

We applied this method to examine the binding affinity for a series of published cyclin-dependent kinase 2 (CDK2) inhibitors. We calculated the binding affinity for 28 CDK2 inhibitors using the ab initio FMO method based on a number of X-ray crystal structures. The sum of the pair interaction energies (PIE) was calculated and used to explain the gas-phase enthalpic contribution to binding. The correlation of the ligand potencies to the protein-ligand interaction energies gained from FMO was examined and was seen to give a good correlation which outperformed three MM force field based scoring functions used to appoximate the free energy of binding. Although the FMO calculation allows for the enthalpic component of binding interactions to be understood at the quantum level, as it is an in vacuo single point calculation, the entropic component and solvation terms are neglected. For this reason a more accurate and predictive estimate for binding free energy was desired. Therefore, additional terms used to describe the protein-ligand interactions were then calculated to improve the correlation of the FMO derived values to experimental free energies of binding. These terms were used to account for the polar and non-polar solvation of the molecule estimated by the Poisson-Boltzmann equation and the solvent accessible surface area (SASA), respectively, as well as a correction term for ligand entropy. A quantitative structure-activity relationship (QSAR) model obtained by Partial Least Squares projection to latent structures (PLS) analysis of the ligand potencies and the calculated terms showed a strong correlation (r2 = 0.939, q2 = 0.896) for the 14 molecule test set which had a Pearson rank order correlation of 0.97. A training set of a further 14 molecules was well predicted (r2 = 0.842), and could be used to obtain meaningful estimations of the binding free energy.

CONCLUSIONS

Our results show that binding energies calculated with the FMO method correlate well with published data. Analysis of the terms used to derive the FMO energies adds greater understanding to the binding interactions than can be gained by MM methods. Combining this information with additional terms and creating a scaled model to describe the data results in more accurate predictions of ligand potencies than the absolute values obtained by FMO alone.