Progress in the understanding of drug-receptor interactions, part 2: experimental and theoretical electrostatic moments and interaction energies of an angiotensin II receptor antagonist (C30H30N6(O)3S)

A combined experimental and theoretical charge density study of an angiotensin II receptor antagonist (1) is presented focusing on electrostatic properties such as atomic charges, molecular electric moments up to the fourth rank and energies of the intermolecular interactions, to gain an insight into the physical nature of the drug-receptor interaction. Electrostatic properties were derived from both the experimental electron density (multipole refinement of X-ray data collected at T=17 K) and the ab initio wavefunction (single molecule and fully periodic calculations at the DFT level). The relevance of SO and SN intramolecular interactions on the activity of 1 is highlighted by using both the crystal and gas-phase geometries and their electrostatic nature is documented by means of QTAIM atomic charges. The derived electrostatic properties are consistent with a nearly spherical electron density distribution, characterised by an intermingling of electropositive and -negative zones rather than by a unique electrophilic region opposed to a nucleophilic area. This makes the first molecular moment scarcely significant and ill-determined, whereas the second moment is large, significant and highly reliable. A comparison between experimental and theoretical components of the third electric moment shows a few discrepancies, whereas the agreement for the fourth electric moment is excellent. The most favourable intermolecular bond is show to be an NHN hydrogen bond with an energy of about 50 kJ mol(-1). Key pharmacophoric features responsible for attractive electrostatic interactions include CHX hydrogen bonds. It is shown that methyl and methylene groups, known to be essential for the biological activity of the drug, provide a significant energetic contribution to the total binding energy. Dispersive interactions are important at the thiophene and at both the phenyl fragments. The experimental estimates of the electrostatic contribution to the intermolecular interaction energies of six molecular pairs, obtained by a new model proposed by Spackman, predict the correct relative electrostatic energies with no exceptions.     

Novel properties from experimental charge densities: an application to the zwitterionic neurotransmitter taurine

The charge distribution of taurine (2-aminoethane-sulfonic acid) is revisited by using an orbital-based method that describes the density in a fixed molecular orbital basis with variable orbital occupation numbers. A new neutron data set is also employed to explore whether this improves the deconvolution of thermal motion and charge density. A range of molecular properties that are novel for experimentally determined charge densities are computed, including Weinhold population analysis, Mayer bond orders, and local kinetic energy densities, in addition to charge topological analysis and quantum theory of atoms-in-molecules (QTAIM) integrated properties. The ease with which a distributed multipole analysis can be performed on the fitted density matrix makes it straightforward to compute molecular moments, the lattice energy, and the electrostatic interaction energies of molecules removed from the crystal. Results are compared with high-level (QCISD) gas-phase calculations and band structure calculations employing density functional theory. Finally, the avenues available for extending the range of molecular properties that can be calculated from experimental charge densities still further using this approach are discussed.     

Polarisable multipolar electrostatics from the machine learning method Kriging: an application to alanine

We present a polarisable multipolar interatomic electrostatic potential energy function for force fields and describe its application to the pilot molecule MeNH-Ala-COMe (AlaD). The total electrostatic energy associated with 1, 4 and higher interactions is partitioned into atomic contributions by application of quantum chemical topology (QCT). The exact atom–atom interaction is expressed in terms of atomic multipole moments. The machine learning method Kriging is used to model the dependence of these multipole moments on the conformation of the entire molecule. The resulting models are able to predict the QCT-partitioned multipole moments for arbitrary chemically relevant molecular geometries. The interaction energies between atoms are predicted for these geometries and compared to their true values. The computational expense of the procedure is compared to that of the point charge formalism.     

A divide and conquer real-space approach for all-electron molecular electrostatic potentials and interaction energies

A computational scheme to perform accurate numerical calculations of electrostatic potentials and interaction energies for molecular systems has been developed and implemented. Molecular electron and energy densities are divided into overlapping atom-centered atomic contributions and a three-dimensional molecular remainder. The steep nuclear cusps are included in the atom-centered functions making the three-dimensional remainder smooth enough to be accurately represented with a tractable amount of grid points. The one-dimensional radial functions of the atom-centered contributions as well as the three-dimensional remainder are expanded using finite element functions. The electrostatic potential is calculated by integrating the Coulomb potential for each separate density contribution, using our tensorial finite element method for the three-dimensional remainder. We also provide algorithms to compute accurate electron-electron and electron-nuclear interactions numerically using the proposed partitioning. The methods have been tested on all-electron densities of 18 reasonable large molecules containing elements up to Zn. The accuracy of the calculated Coulomb interaction energies is in the range of 10(-3) to 10(-6) E(h) when using an equidistant grid with a step length of 0.05 a(0).     

A nonempirical anisotropic atom-atom model potential for chlorobenzene crystals

A nearly nonempirical, transferable model potential is developed for the chlorobenzene molecules (C6ClnH6-n, n = 1 to 6) with anisotropy in the atom-atom form of both electrostatic and repulsion interactions. The potential is largely derived from the charge densities of the molecules, using a distributed multipole electrostatic model and a transferable dispersion model derived from the molecular polarizabilities. A nonempirical transferable repulsion model is obtained by analyzing the overlap of the charge densities in dimers as a function of orientation and separation and then calibrating this anisotropic atom-atom model against a limited number of intermolecular perturbation theory calculations of the short-range energies. The resulting model potential is a significant improvement over empirical model potentials in reproducing the twelve chlorobenzene crystal structures. Further validation calculations of the lattice energies and rigid-body k = 0 phonon frequencies provide satisfactory agreement with experiment, with the discrepancies being primarily due to approximations in the theoretical methods rather than the model intermolecular potential. The potential is able to give a good account of the three polymorphs of p-dichlorobenzene in a detailed crystal structure prediction study. Thus, by introducing repulsion anisotropy into a transferable potential scheme, it is possible to produce a set of potentials for the chlorobenzenes that can account for their crystal properties in an unprecedentedly realistic fashion.     

Ab initio based polarizable force field generation and application to liquid silica and magnesia

We extend the program potfit, which generates effective atomic interaction potentials from ab initio data, to electrostatic interactions and induced dipoles. The potential parametrization algorithm uses the Wolf direct, pairwise summation method with spherical truncation. The polarizability of oxygen atoms is modeled with the Tangney-Scandolo interatomic force field approach. Due to the Wolf summation, the computational effort in simulation scales linearly in the number of particles, despite the presence of electrostatic interactions. Thus, this model allows to perform large-scale molecular dynamics simulations of metal oxides with realistic potentials. Details of the implementation are given, and the generation of potentials for SiO(2) and MgO is demonstrated. The approach is validated by simulations of microstructural, thermodynamic, and vibrational properties of liquid silica and magnesia.     

How dependent are molecular and atomic properties on the electronic structure method? Comparison of Hartree‐Fock,DFT, and MP2 on a biologically relevant set of molecules

This article compares molecular properties and atomic properties defined by the quantum theory of atoms in molecules (QTAIM) obtained from three underlying levels of theory: MP2(full), density functional theory (DFT) (B3LYP), and Hartree‐Fock (H‐F). The same basis set (6‐311++G(d,p)) has been used throughout the study. The calculations and comparisons were applied to a set of 30 small molecules representing common fragments of biological molecules. The molecular properties investigated are the energies and the electrostatic moments (up to and including the quadrupoles), and the atomic properties include electron populations (and atomic charge), atomic dipolar and quadrupolar polarizations, atomic volumes, and corrected and raw atomic energies. The Cartesian distance between dipole vectors and the Frobenius distance between the quadrupole tensors calculated at the three levels of theory provide a measure of their correlation (or lack thereof). With the exception of energies (atomic and molecular), it is found that both DFT and H‐F are in excellent agreement with MP2, especially with regards to the electrostatic mutipoles up to the quadrupoles, but DFT and MP2 agree better in almost all studied properties (with the exception of molecular geometries). QTAIM properties whether obtained from H‐F, DFT(B3LYP), or MP2 calculations when used in the construction of empirical correlations with experiment such as quantitative structure‐activity‐(or property)‐relationships (QSAR/QSPR) are equivalent (because the properties calculated at the three levels are very highly correlated among themselves with r² typically >0.95, and therefore preserving trends). These results suggest that the massive volume of results that were published in the older literature at the H‐F level is valid especially when used to study trends or in QSAR or QSPR studies, and, as long as our test set of molecules is representative, there is no pressing need to re‐evaluate them at other levels of theory except when inadequate basis sets were used by today's standards. Extensive tabulation of molecular and atomic properties at the three theoretical levels is available in the Supporting Information, including optimized geometries, molecular energies, virial ratios, molecular electrostatic moments up to and including hexadecapoles, atomic populations, atomic volumes, atomic electrostatic moments up to and including the quadrupoles, and atomic energies. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Nonempirical Atom-Atom Potentials for Main Components of Intermolecular Interaction Energy

Atom-atom potentials representing separate contributions to the nonempirical interaction energy have been derived in the SCF decomposition scheme corrected for basis set superposition error by the counterpoise method. The nontransferable long-range electrostatic multipole and classical induction terms have been evaluated directly from cumulative atomic multipole expansions, whereas the short-range exchange, charge-transfer, and electrostatic penetration contributions have been represented by simplified potentials of the form (β + δR^?1) exp(?δR) fitted to the corresponding ab initio results for 336 dimer configurations formed by HF, H₂O, NH₃, CH₄, CO, and CO₂. The dominant anisotropic character of electrostatic multipole atom-atom potentials and much more isotropic nature of the potentials representing short-range terms is illustrated in the Appendix for head-on interactions in CO ‥ OC and HF ‥ FH dimers.     

Electrostatic field-adapted molecular fractionation with conjugated caps for energy calculations of charged biomolecules

An electrostatic field-adapted molecular fractionation with conjugated caps (EFA-MFCC) approach is implemented for treating macromolecules with several charge centers. The molecular fragmentation is performed in an "electrostatic field," which is described by putting point charges on charge centers, directly affecting the Hamiltonians of both fragments and conjugated caps. So the present method does not need truncation during the calculation of electrostatic interactions. Our test calculations on a series of charged model systems and biological macromolecules using the HF and B3LYP methods have demonstrated that this approach is capable of describing the electronic structure with accuracy comparable to other fragment-based methods. The EFA-MFCC approach is an alternative way for predicting the total energies of charged macromolecules with acyclic, loop, and intersectional loop structures and interaction energies between two molecules.     

On the interplay between CH...O and OH...O interactions in determining crystal packing and molecular conformation: an experimental and theoretical charge density study of the fungal secondary metabolite austdiol (C12H12O5)

The total experimental electron density rho(r), its Laplacian inverted delta(2)rho(r), the molecular dipole moment, the electrostatic potential phi(r), and the intermolecular interaction energies have been obtained from an extensive set of single-crystal X-ray diffracted intensities, collected at T = 70(1) K, for the fungal metabolite austdiol (1). The experimental results have been compared with theoretical densities from DFT calculations on the isolated molecule and with fully periodic calculations. The crystal structure of (1) consists of zigzag ribbons extended along one cell axis and formed by molecules connected by both OH...O and CH...O interactions, while in a perpendicular direction, adjacent molecules are linked by short CH...O intermolecular contacts. An extensive, quantitative study of all the intra- and intermolecular H...O interactions, based not only on geometrical criteria, but also on the topological analysis of rho(r), as well as on the evaluation of the pertinent energetics, allowed us (i) to assess the mutual role of OH...O and CH...O interactions in determining molecular conformation and crystal packing; (ii) to identify those CH...O contacts which are true hydrogen bonds (HBs); (iii) to determine the relative hydrogen bond strengths. An experimental, quantitative evidence is given that CH...O HBs are very similar to the conventional OH...O HBs, albeit generally weaker. The comparison between experimental and theoretical electric dipole moments indicates that a noticeable charge rearrangement occurs upon crystallization and shows the effects of the mutual cooperation of HBs in the crystal. The total intermolecular interaction energies and the electrostatic energy contribution obtained through different theoretical methods are reported and compared with the experimental results. It is found that the new approach proposed by Spackman, based on the use of the promolecular charge density to approximate the penetration contribution to intermolecular electrostatic energies, predicts the correct relative electrostatic interaction energies in most of the cases.     

Short-range potential functions in computer simulations of water and aqueous solutions

A review and comparative analyses of methods for restricting the range of molecular interactions within the concept of atom-atom potentials are presented. Emphasis is placed on the problem of calculating the electrostatic energy in models with periodic boundary conditions. Numerous calculations of the thermodynamic and structural characteristics of water using parallel Monte Carlo computations have shown that the use of functional forms simulating the electric potentials of “screened charges” provides very good results.     

Accuracy of distributed multipoles and polarizabilities: comparison between the LoProp and MpProp models

Localized multipole moments up to the fifth moment as well as localized dipole polarizabilities are calculated with the MpProp and the newly developed LoProp methods for a total of 20 molecules, predominantly derived from amino acids. A comparison of electrostatic potentials calculated from the multipole expansion obtained by the two methods with ab initio results shows that both methods reproduce the electrostatic interaction with an elementary charge with a mean absolute error of approximately 1.5 kJ/mol at contact distance and less than 0.1 kJ/mol at distances 2 A further out when terms up to the octupole moments are included. The polarizabilities are tested with homogenous electric fields and are found to have similar accuracy. The MpProp method gives better multipole moments unless diffuse basis sets are used, whereas LoProp gives better polarizabilities.     

Influence of metal ion on chelate–aryl stacking interactions

CCSD(T)/CBS and DFT methods are employed to study the stacking interactions of acetylacetonate‐type (acac‐type) chelates of nickel, palladium, and platinum with benzene. The strongest chelate–aryl stacking interactions are formed by nickel and palladium chelate, with interaction energies of −5.75 kcal mol⁻¹ and −5.73 kcal mol⁻¹, while the interaction of platinum chelate is weaker, with interaction energy of −5.36 kcal mol⁻¹. These interaction energies are significantly stronger than stacking of two benzenes, −2.73 kcal mol⁻¹. The strongest nickel and palladium chelate–aryl interactions are with benzene center above the metal area, while the strongest platinum chelate–aryl interaction is with the benzene center above the C2 atom of the acac‐type chelate ring. These preferences arise from very different electrostatic potentials above the metal ions, ranging from very positive above nickel to slightly negative above platinum. While the differences in electrostatic potentials above metal atoms cause different geometries with the most stable interaction among the three metals, the dispersion (correlation energy) component is the largest contribution to the total interaction energy for all three metals.     

Atom-atom partitioning of total (super)molecular energy: the hidden terms of classical force fields

Classical force fields describe the interaction between atoms that are bonded or nonbonded via simple potential energy expressions. Their parameters are often determined by fitting to ab initio energies and electrostatic potentials. A direct quantum chemical guide to constructing a force field would be the atom-atom partitioning of the energy of molecules and van der Waals complexes relevant to the force field. The authors used the theory of quantum chemical topology to partition the energy of five systems [H2, CO, H2O, (H2O)2, and (HF)2] in terms of kinetic, Coulomb, and exchange intra-atomic and interatomic contributions. The authors monitored the variation of these contributions with changing bond length or angle. Current force fields focus only on interatomic interaction energies and assume that these purely potential energy terms are the only ones that govern structure and dynamics in atomistic simulations. Here the authors highlight the importance of self-energy terms (kinetic and intra-atomic Coulomb and exchange).     

Anion⋅⋅⋅Anion Interactions Involving σ-Holes of Perrhenate,Pertechnetate and Permanganate Anions

In this communication experimental and theoretical results are reported affording strong evidence that interactions between electron rich atoms and the metal of tetroxide anions of group 7 elements are a new case of attractive and σ-hole interactions. Single crystal X-ray analyses, molecular electrostatic potentials, quantum theory of atoms-in-molecules, and noncovalent interaction plot analyses show that in crystalline permanganate and perrhenate salts the metal in Mn/ReO₄⁻ anion can act as electron acceptors, the oxygen of another Mn/ReO₄⁻ anion can act as the donor and supramolecular anionic dimers or polymers are formed. The name matere bond (MaB) is proposed to categorize these noncovalent interactions and to differentiate them from the classical metal-ligand coordination bond.