Halogen–Aromatic π Interactions Modulate Inhibitor Residence Times

Prolonged drug residence times may result in longer‐lasting drug efficacy, improved pharmacodynamic properties, and “kinetic selectivity” over off‐targets with high drug dissociation rates. However, few strategies have been elaborated to rationally modulate drug residence time and thereby to integrate this key property into the drug development process. Herein, we show that the interaction between a halogen moiety on an inhibitor and an aromatic residue in the target protein can significantly increase inhibitor residence time. By using the interaction of the serine/threonine kinase haspin with 5‐iodotubercidin (5‐iTU) derivatives as a model for an archetypal active‐state (type I) kinase–inhibitor binding mode, we demonstrate that inhibitor residence times markedly increase with the size and polarizability of the halogen atom. The halogen–aromatic π interactions in the haspin–inhibitor complexes were characterized by means of kinetic, thermodynamic, and structural measurements along with binding‐energy calculations.     

Ab initio calculation of the interaction energy in the P2 binding pocket of HIV‐1 protease

Interactions at the P2 binding pocket of human immunodeficiency virus type 1 (HIV‐1) protease have been studied using calculated interaction energies for model systems that mimic this binding pocket. Models were built for the P2 pocket of HIV‐1 protease in complex with TMC114, nelfinavir, and amprenavir. A two‐step procedure was applied. In the first step, the size of the model system was confined to ～40 atoms, and the interaction energy was calculated at different computational levels. In the second step, the size of the system was increased to 138 atoms, and the calculations were only performed at the HF/6‐31G** level. The interaction energy of the HIV‐1 protease/TMC114 complex was found to be more favorable than the interaction energies of the other complexes because of the additional hydrogen bond interaction this inhibitor is able to make with the HIV‐1 protease backbone. The results of the calculations are supported by stockholder charges and electrostatic potential maps. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Surfing π Clouds for Noncovalent Interactions: Arenes versus Alkenes

A comparative study of molecular balances by NMR spectroscopy indicates that noncovalent functional‐group interactions with an arene dominate over those with an alkene, and that a π‐facial intramolecular hydrogen bond from a hydroxy group to an arene is favored by approximately 1.2 kJ mol^?1. The strongest interaction observed in this study was with the cyano group. Analysis of the series of groups CH₂CH₃, CH?CH₂, C?CH, and C?N shows a correlation between conformational free‐energy differences and the calculated charge on the C^α atom of these substituents, which is indicative of the electrostatic nature of their π interactions. Changes in the free‐energy differences of conformers show a linear dependence on the solvent hydrogen bond acceptor parameter β.     

Study on the Cation-π Interactions between Ammonium Ion and Aromatic π Systems

The nature and strength of the cation-π interactions between NH4^＋ and toluene, p-cresol, or Me-indole were studied in terms of the topological properties of molecular charge density and binding energy decomposition. The results display that the diversity in the distribution pattern of bond and cage critical points reflects the profound influence of the number and nature of substituent on the electron density of the aromatic rings. On the other hand, the energy decomposition shows that dispersion and repulsive exchange forces play an important role in the organic cation （NH4^＋）-π interaction, although the electrostatic and induction forces dominate the interaction. In addition, it is intriguing that there is an excellent correlation between the electrostatic energy and ellipticity at the bond critical point of the aromatic π systems, which would be helpful to further understand the electrostatic interaction in the cation-π complexes.     

Observation of CH⋅⋅⋅π Interactions between Methyl and Carbonyl Groups in Proteins

Protein structure and function is dependent on myriad noncovalent interactions. Direct detection and characterization of these weak interactions in large biomolecules, such as proteins, is experimentally challenging. Herein, we report the first observation and measurement of long‐range “through‐space” scalar couplings between methyl and backbone carbonyl groups in proteins. These J couplings are indicative of the presence of noncovalent C−H⋅⋅⋅π hydrogen‐bond‐like interactions involving the amide π network. Experimentally detected scalar couplings were corroborated by a natural bond orbital analysis, which revealed the orbital nature of the interaction and the origins of the through‐space J couplings. The experimental observation of this type of CH⋅⋅⋅π interaction adds a new dimension to the study of protein structure, function, and dynamics by NMR spectroscopy.     

乙烯、乙炔和苯与氨基锂间的π型锂键相互作用

分别在DFT-B3LYP和MP2/6-311++G**水平上求得乙烯、乙炔和苯与氨基锂锂键复合物势能面上的3个几何构型. 频率分析表明,3个构型均为稳定构型. 计算结果表明,形成锂键复合物后,质子供体N(2)- Li(4)的键长明显增大,且其伸缩振动的频率发生了不同程度的红移. 分别在乙烯…氨基锂、乙炔…氨基锂和苯…氨基锂三种复合物中,经MP2/6-311++G**水平计算的同时含基组重叠误差(BSSE)和零点振动能校正的单体间锂键相互作用能分别为-26.04、-24.86 和 -30.02 kJ·mol-1. 自然键轨道理论（NBO）分析表明单体间的弱相互作用属于π-s型锂键.     

Interplay between Metal⋅⋅⋅π Interactions and Hydrogen Bonds: Some Unusual Synergetic Effects of Coinage Metals and Substituents

The ternary systems of C₂H₄ (C₂H₂ or C₆H₆)‐MCN‐HF (M=Cu, Ag, Au) and the respective binary systems were investigated to study the interplay between metal???π interactions and hydrogen bonds. The metal???π interactions in C₂H₄‐MCN become stronger with the irregular order Ag<Cu<Au, while the hydrogen bonds in MCN‐HF become weaker following the same order. The metal???π interactions are weakened as the H atoms in the π system are replaced with electron‐withdrawing groups and enhanced by electron‐donating groups. Type 1 of these ternary systems, in which MCN acts as Lewis base and acid simultaneously, is more stable than type 2, in which C₂H₄ acts as a double Lewis base. Negative cooperativity is present in type 2 ternary systems with a weakening of the metal???π interactions and the hydrogen bonds. Positive cooperativity is found in type 1 ternary systems with an enhancement of the metal???π interactions and the hydrogen bonds, except for C₂(CN)₄‐AuCN‐HF‐1. The weaker metal???π interaction in C₆H₆‐AuCN has a greater enhancing effect on the hydrogen bond in AuCN‐HF than those in C₂H₄‐AuCN and C₂H₂‐AuCN. These synergetic effects were analyzed with the natural bond orbital and energy decomposition.     

Herringbone Pattern and CH–π Bonding in the Crystal Architecture of Linear Polycyclic Aromatic Hydrocarbons

The herringbone pattern is a pervasive structural motive found in most molecular crystals involving aromatic compounds. A plot of the experimental sublimation enthalpies of members of increasing size of the acene, phenacene and p‐phenyl families versus the number of carbons uncovers a linear relationship between the two magnitudes, suggesting a major role of CH–π bonding. In this work we undertake the task of evaluating the relevance of the edge‐to‐face interaction (or CH–π bond) in the overall reticular energy of the crystal, to quantitatively assess the importance of this structural element. Following a heuristic approach, we considered the series of acenes, phenacenes and p‐phenyls and analyzed the edge‐to‐face interaction between the molecules as they occur in the experimental crystal network. Isolation of the relevant molecular dimers allows to incorporate some of the most sophisticated tools of quantum chemistry and get a reliable picture of the isolated bond. When compared to the experimental sublimation energy, our results are conclusive: this sole interaction is the largest contribution to the lattice energy, and definitively dictates the crystal architecture in all the studied cases. Elusive enough, the edge‐to‐face interaction is mainly dominated by correlation interactions, specifically in the form of dispersion and, to a less extent, of charge‐transfer terms. A suggestive picture of the bond has been obtained by displaying the differences in local electron densities calculated by either correlated or non‐correlated methods.     

Valence state parameters of all transition metal atoms in metalloproteins—development of ABEEMσπ fluctuating charge force field

To promote accuracy of the atom‐bond electronegativity equalization method (ABEEMσπ) fluctuating charge polarizable force fields, and extend it to include all transition metal atoms, a new parameter, the reference charge is set up in the expression of the total energy potential function. We select over 700 model molecules most of which model metalloprotein molecules that come from Protein Data Bank. We set reference charges for different apparent valence states of transition metals and calibrate the parameters of reference charges, valence state electronegativities, and valence state hardnesses for ABEEMσπ through linear regression and least square method. These parameters can be used to calculate charge distributions of metalloproteins containing transition metal atoms (Sc‐Zn, Y‐Cd, and Lu‐Hg). Compared the results of ABEEMσπ charge distributions with those obtained by ab initio method, the quite good linear correlations of the two kinds of charge distributions are shown. The reason why the STO‐3G basis set in Mulliken population analysis for the parameter calibration is specially explained in detail. Furthermore, ABEEMσπ method can also quickly and quite accurately calculate dipole moments of molecules. Molecular dynamics optimizations of five metalloproteins as the examples show that their structures obtained by ABEEMσπ fluctuating charge polarizable force field are very close to the structures optimized by the ab initio MP2/6–311G method. This means that the ABEEMσπ/MM can now be applied to molecular dynamics simulations of systems that contain metalloproteins with good accuracy. © 2014 Wiley Periodicals, Inc.     

Ab initio through‐space/bond interaction analysis of the long C–C bonds in Bi(anthracene‐9,10‐dimethylene) photoisomers

The bi(anthracene‐9,10‐dimethylene) photoisomer has remarkably long C–C single bonds. To examine the lengthening of the C–C bond, we propose a novel procedure for quantitatively analyzing orbital interactions in a molecule at the level of the ab initio molecular orbital method. In this procedure, we can cut off the specific through‐space/bond interactions in a molecule by artificially increasing the absolute magnitude of the exponents in a Gaussian function. Then, the spatial orbital interactions were perfectly cut off, and, each term that makes up the total energy, that is, the nuclear–electron attractions, the electron–electron repulsions, and the nuclear–nuclear repulsions cancel each other. Several model molecules of the photoisomer were analyzed by this procedure. It was found that the orbital interaction between the p orbital on the benzene ring and the σ* orbital on the C–C bond in question, σ→σ* electron transfer through π orbital, weakens the C–C bond efficiently when these orbitals were located in the “periplanar” conformation. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Mixed implicit/explicit solvation models in quantum mechanical calculations of binding enthalpy for protein–ligand complexes

An approach to quantum mechanical investigation of interactions in protein–ligand complexes has been developed that treats the solvation effect in a mixed scheme combining implicit and explicit solvent models. In this approach, the first solvation shell of the solvent around the solute is modeled with a limited number of hydrogen bonded explicit solvent molecules. The influence of the remaining bulk solvent is treated as a surrounding continuum in the conductor‐like screening model (COSMO). The enthalpy term of the binding free energy for the protein–ligand complexes was calculated using the semiempirical PM3 method implemented in the MOPAC package, applied to a trimmed model of the protein–ligand complex constructed with special rules. The dependence of the accuracy of binding enthalpy calculations on size of the trimmed model and number of optimized parameters was evaluated. Testing of the approach was performed for 12 complexes of different ligands with trypsin, thrombin, and ribonuclease with experimentally known binding enthalpies. The root‐mean‐square deviation (RMSD) of the calculated binding enthalpies from experimental data was found as ～1 kcal/mol over a large range. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

The Cation–π Interaction in Small‐Molecule Catalysis

Catalysis by small molecules (≤1000 Da, 10^?9 m) that are capable of binding and activating substrates through attractive, noncovalent interactions has emerged as an important approach in organic and organometallic chemistry. While the canonical noncovalent interactions, including hydrogen bonding, ion pairing, and π stacking, have become mainstays of catalyst design, the cation–π interaction has been comparatively underutilized in this context since its discovery in the 1980s. However, like a hydrogen bond, the cation–π interaction exhibits a typical binding affinity of several kcal mol^?1 with substantial directionality. These properties render it attractive as a design element for the development of small‐molecule catalysts, and in recent years, the catalysis community has begun to take advantage of these features, drawing inspiration from pioneering research in molecular recognition and structural biology. This Review surveys the burgeoning application of the cation–π interaction in catalysis.     

Theoretical studies of uracil–(H2O)n (n = 1–7) clusters by ab initio and ABEEMσπ/MM fluctuating charge model

Uracil–(H₂O)_n (n = 1–7) clusters were systemically investigated by ab initio methods and the newly constructed ABEEMσπ/MM fluctuating charge model. Water molecules have been gradually placed in an average plane containing uracil. The geometries of 38 uracil–water complexes were obtained using B3LYP/6-311++G** level optimizations, and the energies were determined at the MP2/6-311++G** level with BSSE corrections. The ABEEMσπ/MM potential model gives reasonable properties of these clusters when comparing with the present ab initio data. For interaction energies, the root mean square deviation is 0.96 kcal/mol, and the linear coefficient reaches 0.997. Furthermore, the ABEEMσπ charges changed when H₂O interacted with the uracil molecule, especially at the sites where the hydrogen bond form. These results show that the ABEEMσπ/MM model is fine giving the overall characteristic hydration properties of uracil–water systems in good agreement with the high-level ab initio calculations.     

New parameterization approaches of the LIE method to improve free energy calculations of PlmII‐Inhibitors complexes

The standard parameterization of the Linear Interaction Energy (LIE) method has been applied with quite good results to reproduce the experimental absolute binding free energies for several protein–ligand systems. However, we found that this parameterization failed to reproduce the experimental binding free energy of Plasmepsin II (PlmII) in complexes with inhibitors belonging to four dissimilar scaffolds. To overcome this fact, we developed three approaches of LIE, which combine systematic approaches to predict the inhibitor‐specific values of α, β, and γ parameters, to gauge their ability to calculate the absolute binding free energies for these PlmII‐Inhibitor complexes. Specifically: (i) we modified the linear relationship between the weighted nonpolar desolvation ratio (WNDR) and the α parameter, by introducing two models of the β parameter determined by the free energy perturbation (FEP) method in the absence of the constant term γ, and (ii) we developed a new parameterization model to investigate the linear correlation between WNDR and the correction term γ. Using these parameterizations, we were able to reproduce the experimental binding free energy from these systems with mean absolute errors lower than 1.5 kcal/mol. © 2010 Wiley Periodicals, Inc. J Comput Chem 2010     

Combining the polarizable Drude force field with a continuum electrostatic Poisson–Boltzmann implicit solvation model

In this work, we have combined the polarizable force field based on the classical Drude oscillator with a continuum Poisson–Boltzmann/solvent‐accessible surface area (PB/SASA) model. In practice, the positions of the Drude particles experiencing the solvent reaction field arising from the fixed charges and induced polarization of the solute must be optimized in a self‐consistent manner. Here, we parameterized the model to reproduce experimental solvation free energies of a set of small molecules. The model reproduces well‐experimental solvation free energies of 70 molecules, yielding a root mean square difference of 0.8 kcal/mol versus 2.5 kcal/mol for the CHARMM36 additive force field. The polarization work associated with the solute transfer from the gas‐phase to the polar solvent, a term neglected in the framework of additive force fields, was found to make a large contribution to the total solvation free energy, comparable to the polar solute–solvent solvation contribution. The Drude PB/SASA also reproduces well the electronic polarization from the explicit solvent simulations of a small protein, BPTI. Model validation was based on comparisons with the experimental relative binding free energies of 371 single alanine mutations. With the Drude PB/SASA model the root mean square deviation between the predicted and experimental relative binding free energies is 3.35 kcal/mol, lower than 5.11 kcal/mol computed with the CHARMM36 additive force field. Overall, the results indicate that the main limitation of the Drude PB/SASA model is the inability of the SASA term to accurately capture non‐polar solvation effects. © 2018 Wiley Periodicals, Inc.     

Agonist–PPARγ interactions: Molecular modeling study with docking approach

Docking simulation of 18 agonists with the ligand binding pocket (LBP) of PPARγ has been performed. The binding conformations and binding affinities of these agonists were obtained by use of the flexible docking protocol FlexX. Test compound calculations indicated that FlexX can reproduce the binding conformation of the crystal structure (root mean square deviation = 1.43 Å); moreover, the predicted binding affinities correlate well with the activities of these agonists. The interaction model and pharmacophore of PPARγ agonists were derived and the difference in biologic activities of these agonists can be well explained. The PPARγ agonists must have both polar head and the hydrophobic tail, which form hydrogen bonds and hydrophobic contacts with hydrophilic and hydrophobic regions of the LBP of PPARγ, respectively. In addition, a suitable linker is also necessary. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem 93: 405–410, 2003     

Noncovalent Lone Pair⋅⋅⋅(No‐π!)‐Heteroarene Interactions: The Janus‐Faced Hydroxy Group

A comparative study using NMR spectroscopy and designed top‐pan molecular balances demonstrates that the noncovalent interaction of a hydroxy group with π‐deficient pyrazine and quinoxaline units involves a lone pair–heteroarene interaction which is much stronger and solvent independent when measured relative to the classical π‐facial hydrogen bond to a benzene ring. Alkyl fluorides also prefer the heteroarene rings over the benzene ring. The attractive interaction between a quinoxaline and a terminal alkyne is also stronger than the intramolecular hydrogen bond to an arene.     

Strong N−X⋅⋅⋅O−N Halogen Bonds: A Comprehensive Study on N‐Halosaccharin Pyridine N‐Oxide Complexes

A study of the strong N?X???^?O?N⁺ (X=I, Br) halogen bonding interactions reports 2×27 donor×acceptor complexes of N‐halosaccharins and pyridine N‐oxides (PyNO). DFT calculations were used to investigate the X???O halogen bond (XB) interaction energies in 54 complexes. A simplified computationally fast electrostatic model was developed for predicting the X???O XBs. The XB interaction energies vary from ?47.5 to ?120.3 kJ mol^?1; the strongest N?I???^?O?N⁺ XBs approaching those of 3‐center‐4‐electron [N?I?N]⁺ halogen‐bonded systems (ca. 160 kJ mol^?1). ¹H NMR association constants (K_XB) determined in CDCl₃ and [D₆]acetone vary from 2.0×10⁰ to >10⁸ m ^?1 and correlate well with the calculated donor×acceptor complexation enthalpies found between ?38.4 and ?77.5 kJ mol^?1. In X‐ray crystal structures, the N‐iodosaccharin‐PyNO complexes manifest short interaction ratios (R_XB) between 0.65–0.67 for the N?I???^?O?N⁺ halogen bond.     

σ–π Continuum in Indole–Palladium(II) Complexes

The intrinsic features of (hetero‐arene)–metal interactions have been elusive mainly because the systematic structure analysis of non‐anchored hetero‐arene–metal complexes has been hampered by their labile nature. We report successful isolation and systematic structure analysis of a series of non‐anchored indole–palladium(II) complexes. It was revealed that there is a σ–π continuum for the indole–metal interaction, while it has been thought that the dominant coordination mode of indole to a metal center is the Wheland‐intermediate‐type σ‐mode in light of the seemingly strong electron‐donating ability of indole. Several factors which affect the σ‐ or π‐character of indole–metal interactions are discussed.