The experimental charge-density approach in the evaluation of intermolecular interactions. Application of a new module of the XD programming package to several solids including a pentapeptide

On the calculation of the electrostatic potential, electric field and electric field gradient from the aspherical pseudoatom model. II. Evaluation of the properties in an infinite crystal

The previously reported exact potential and multipole moment (EP/MM) method for fast and precise evaluation of the intermolecular electrostatic interaction energies in molecular crystals using the pseudoatom representation of the electron density [Nguyen, Macchi & Volkov (2020), Acta Cryst. A76, 630-651] has been extended to the calculation of the electrostatic potential (ESP), electric field (EF) and electric field gradient (EFG) in an infinite crystal. The presented approach combines an efficient Ewald-type summation (ES) of atomic multipoles up to the hexadecapolar level in direct and reciprocal spaces with corrections for (i) the net polarization of the sample (the `surface term') due to a net dipole moment of the crystallographic unit cell (if present) and (ii) the short-range electron-density penetration effects. The rederived and reported closed-form expressions for all terms in the ES algorithm have been augmented by the expressions for the surface term available in the literature [Stenhammar, Trulsson & Linse (2011), J. Chem. Phys. 134, 224104] and the exact potential expressions reported in a previous study [Volkov, King, Coppens & Farrugia (2006), Acta Cryst. A62, 400-408]. The resulting algorithm, coded using Fortran in the XDPROP module of the software package XD, was tested on several small molecular crystal systems (formamide, benzene, L-dopa, paracetamol, amino acids etc.) and compared with a series of EP/MM-based direct-space summations (DS) performed within a certain number of unit cells generated along both the positive and negative crystallographic directions. The EP/MM-based ES technique allows for a noticeably more precise determination of the EF and EFG and significantly better precision of the evaluated ESP when compared with the DS calculations, even when the latter include contributions from an array of symmetry-equivalent atoms generated within four additional unit cells along each crystallographic direction. In terms of computational performance, the ES/EP/MM method is significantly faster than the DS calculations performed within the extended unit-cell limits but trails the DS calculations within the reduced summation ranges. Nonetheless, the described EP/MM-based ES algorithm is superior to the direct-space summations as it does not require the user to monitor continuously the convergence of the evaluated properties as a function of the summation limits and offers a better precision-performance balance.

Fast analytical evaluation of intermolecular electrostatic interaction energies using the pseudoatom representation of the electron density. III. Application to crystal structures via the Ewald and direct summation methods

The previously reported exact potential and multipole moment (EP/MM) method for fast and accurate evaluation of the intermolecular electrostatic interaction energies using the pseudoatom representation of the electron density [Volkov, Koritsanszky & Coppens (2004). Chem. Phys. Lett. 391, 170-175; Nguyen, Kisiel & Volkov (2018). Acta Cryst. A74, 524-536; Nguyen & Volkov (2019). Acta Cryst. A75, 448-464] is extended to the calculation of electrostatic interaction energies in molecular crystals using two newly developed implementations: (i) the Ewald summation (ES), which includes interactions up to the hexadecapolar level and the EP correction to account for short-range electron-density penetration effects, and (ii) the enhanced EP/MM-based direct summation (DS), which at sufficiently large intermolecular separations replaces the atomic multipole moment approximation to the electrostatic energy with that based on the molecular multipole moments. As in the previous study [Nguyen, Kisiel & Volkov (2018). Acta Cryst. A74, 524-536], the EP electron repulsion integral is evaluated analytically using the Löwdin α-function approach. The resulting techniques, incorporated in the XDPROP module of the software package XD2016, have been tested on several small-molecule crystal systems (benzene, L-dopa, paracetamol, amino acids etc.) and the crystal structure of a 181-atom decapeptide molecule (Z = 4) using electron densities constructed via the University at Buffalo Aspherical Pseudoatom Databank [Volkov, Li, Koritsanszky & Coppens (2004). J. Phys. Chem. A, 108, 4283-4300]. Using a 2015 2.8 GHz Intel Xeon E3-1505M v5 computer processor, a 64-bit implementation of the Löwdin α-function and one of the higher optimization levels in the GNU Fortran compiler, the ES method evaluates the electrostatic interaction energy with a numerical precision of at least 10^-5 kJ mol^-1 in under 6 s for any of the tested small-molecule crystal structures, and in 48.5 s for the decapeptide structure. The DS approach is competitive in terms of precision and speed with the ES technique only for crystal structures of small molecules that do not carry a large molecular dipole moment. The electron-density penetration effects, correctly accounted for by the two described methods, contribute 28-64% to the total electrostatic interaction energy in the examined systems, and thus cannot be neglected.

Fast analytical evaluation of intermolecular electrostatic interaction energies using the pseudoatom representation of the electron density. I. The Löwdin α-function method

The previously reported [Volkov et al. (2004). Chem. Phys. Lett. 391, 170-175] exact potential and multipole moment (EP/MM) method for evaluation of intermolecular electrostatic interaction energies using the nuclei-centered pseudoatom representation of electron densities is significantly improved in terms of both speed and accuracy by replacing the numerical quadrature integration of the exact potential with a fully analytical technique. The resulting approach, incorporated in the XDPROP module of the software package XD, has been tested on several molecular systems ranging in size from water-water to dodecapeptide-dodecapeptide dimers using electron densities constructed via the University at Buffalo Aspherical Atom Databank. The improved hybrid method provides electrostatic interaction energies within the uncertainty of ≤0.2 kJ mol^-1 for all benchmark systems. The running time for a dimer of a sizable, 225-atom dodecapeptide is under 4 s on a 2012 central processing unit (2.8 GHz AMD Opteron 6348) and under 3 s on a relatively modern processor (2.8 GHz Intel Xeon E3-1505M v5).

Towards the use of experimental electron densities to estimate reliable lattice energies

Lattice energies derived from experimental charge densities are critically assessed, with a view to encouraging further research of this nature.

Energy frameworks: insights into interaction anisotropy and the mechanical properties of molecular crystals

We present an approach to understanding crystal packing via 'energy frameworks', that combines efficient calculation of accurate intermolecular interaction energies with a novel graphical representation of their magnitude. In this manner intriguing questions, such as why some crystals bend with an applied force while others break, and why one polymorph of a drug exhibits exceptional tabletability compared to others, can be addressed in terms of the anisotropy of the topology of pairwise intermolecular interaction energies. This approach is applied to a sample of organic molecular crystals with known bending, shearing and brittle behaviour, to illustrate its use in rationalising their mechanical behaviour at a molecular level.

Contributions of charge-density research to medicinal chemistry

This article reviews efforts in accurate experimental charge-density studies with relevance to medicinal chemistry. Initially, classical charge-density studies that measure electron density distribution via least-squares refinement of aspherical-atom population parameters are summarized. Next, interaction density is discussed as an idealized situation resembling drug-receptor interactions. Scattering-factor databases play an increasing role in charge-density research, and they can be applied both to small-molecule and macromolecular structures in refinement and analysis; software development facilitates their use. Therefore combining both of these complementary branches of X-ray crystallography is recommended, and examples are given where such a combination already proved useful. On the side of the experiment, new pixel detectors are allowing rapid measurements, thereby enabling both high-throughput small-molecule studies and macromolecular structure determination to higher resolutions. Currently, the most ambitious studies compute intermolecular interaction energies of drug-receptor complexes, and it is recommended that future studies benefit from recent method developments. Selected new developments in theoretical charge-density studies are discussed with emphasis on its symbiotic relation to crystallography.

C—H...X(X= O, N or π) interactions in benzyl carbamate

The crystal packing and interaction energy of benzyl carbamate, C8H9NO2, have been analysed in detail by thePIXELmethod. Benzyl carbamate forms layers of hydrogen-bonded molecules, with the layers connected by weaker C—H...π interactions. According to thePIXELanalysis, combinations of C—H...X(X= O, N or π) interactions are comparable in energy with hydrogen bonding. These interactions are necessary for explaining the geometry and the assembly of the layers.

Physicochemical properties of zwitterionic L- and DL-alanine crystals from their experimental and theoretical charge densities

The total experimental electron density distributions rho(r) of zwitterionic L- and DL-alanine crystals, as derived from extensive sets of X-ray diffracted intensities collected at 23 and 19 K, are compared to gain an insight into the different physical properties of the two related chiral compounds in the solid state and to explore the extent of the rho(r) transferability. Relevant parameters that characterize the two crystal forms are obtained, showing differences and similarities in terms of (i) geometric descriptors, (ii) topological indexes, (iii) molecular electrostatic potential Phi(r) distributions, (iv) atomic volumes and charges, (v) molecular electric moments, and (vi) electrostatic interaction energies. To assess the relative stability of the racemate with respect to the pure enantiomer, the crystal lattice energies, as obtained through DFT fully periodic calculations, are also discussed and compared with the experimental sublimation enthalpies after correction for the proton-transfer energies. In-crystal group charges, evaluated with the quantum theory of atoms in molecules, are found to be transferable between the racemic and the pure enantiomer, at variance with group volumes. Similarly, molecular first and third moments are not strictly transferable and indicate that for the zwitterionic alanine molecule the molecular charge distribution in the DL-crystal is more polarized in the c direction by about 10%. By contrast, quantitative agreement is observed for second and fourth moments. Significant differences arise from (1) the crystal packing of the dipole vectors, which are aligned in an antiparallel fashion in the L-crystal, to be compared with a parallel alignment in the racemate, due the polar space group Pna21 of the latter, (2) the strongly attractive electrostatic energy of a homochiral pair in the L-crystal, which is opposed to the corresponding heterochiral pair in the DL-crystal form. The difference between these Ees values amounts to 135-150 kJ mol(-1). Despite this, the two crystal forms are predicted as equally thermodynamically favored by the theoretical P-B3LYP estimates of the crystal lattice energies. Finally, the necessity of an upgrading of the dispersion and exchange-repulsion terms currently adopted within the experimental charge density approach to intermolecular interactions is recognized and discussed.

Calculation of electrostatic interaction energies in molecular dimers from atomic multipole moments obtained by different methods of electron density partitioning

Accurate and fast evaluation of electrostatic interactions in molecular systems is still one of the most challenging tasks in the rapidly advancing field of macromolecular chemistry, including molecular recognition, protein modeling and drug design. One of the most convenient and accurate approaches is based on a Buckingham-type approximation that uses the multipole moment expansion of molecular/atomic charge distributions. In the mid-1980s it was shown that the pseudoatom model commonly used in experimental X-ray charge density studies can be easily combined with the Buckingham-type approach for calculation of electrostatic interactions, plus atom-atom potentials for evaluation of the total interaction energies in molecular systems. While many such studies have been reported, little attention has been paid to the accuracy of evaluation of the purely electrostatic interactions as errors may be absorbed in the semiempirical atom-atom potentials that have to be used to account for exchange repulsion and dispersion forces. This study is aimed at the evaluation of the accuracy of the calculation of electrostatic interaction energies with the Buckingham approach. To eliminate experimental uncertainties, the atomic moments are based on theoretical single-molecule electron densities calculated at various levels of theory. The electrostatic interaction energies for a total of 11 dimers of alpha-glycine, N-acetylglycine and L-(+)-lactic acid structures calculated according to Buckingham with pseudoatom, stockholder and atoms-in-molecules moments are compared with those evaluated with the Morokuma-Ziegler energy decomposition scheme. For alpha-glycine a comparison with direct "pixel-by-pixel" integration method, recently developed Gavezzotti, is also made. It is found that the theoretical pseudoatom moments combined with the Buckingham model do predict the correct relative electrostatic interactions energies, although the absolute interaction energies are underestimated in some cases. The good agreement between electrostatic interaction energies computed with Morokuma-Ziegler partitioning, Gavezzotti's method, and the Buckingham approach with atoms-in-molecules moments demonstrates that reliable and accurate evaluation of electrostatic interactions in molecular systems of considerable complexity is now feasible.

Application of charge density methods to a protein model compound: calculation of Coulombic intermolecular interaction energies from the experimental charge density

A combined experimental and theoretical charge density study of the pentapeptide Boc-Gln-d-Iva-Hyp-Ala-Phol (Boc, butoxycarbonyl; Gln, glutamine; Iva, isovaline; Hyp, hydroxyproline; Ala, ethylalanine; Phol, phenylalaninol) is described. The experimental analysis, based on synchrotron x-ray data collected at 20 K, is combined with ab initio theoretical calculations. The topologies of the experimental and theoretical densities are analyzed in terms of the atoms in molecules quantum theory. Topological parameters, including atomic charges and higher moments integrated over the atomic basins, have been evaluated with the program topxd and are used to calculate the electrostatic interactions between the molecules in the crystal. The interaction energies obtained after adding dispersive and repulsive van der Waals contributions agree quite well with those based on M-B3LYP/6-31G** dimer calculations for two of the three dimers in the crystal, whereas for the third a larger stabilization is obtained than predicted by the calculation. The agreement with theory is significantly better than that obtained with multipole moments derived directly from the aspherical atom refinement. The convergence of the interaction as a function of addition of successively higher moments up to and including hexadecapoles (l = 4) is found to be within 2-3 kJ/mol. Although shortcomings of both the theoretical and experimental procedures are pointed out, the agreement obtained supports the potential of the experimental method for the evaluation of interactions in larger biologically relevant molecules.