Tailoring approach for exploring electron densities and electrostatic potentials of molecular crystals

Experimental and theoretical studies of electron densities and the corresponding derived entities such as electrostatic potentials have been the primary means of understanding the chemical nature and electronic properties of crystalline substances. Conventional crystal calculation methods such as the embedded cluster models are capable of performing calculations on small and medium-sized molecules, while periodic ab initio methods can treat crystals with up to 200 atoms per unit cell. A linear scaling method, viz. the molecular tailoring approach, has recently been developed for obtaining ab initio quality one-electron properties. In the present study, the molecular tailoring approach is employed to generate electron density, electrostatic potential and interaction density maps with the ibuprofen crystal as a test case. The interaction density and electrostatic potential maps produced in the present work succinctly bring out the actual crystalline environment around a given reference molecule by including the interactions with atoms in its neighborhood. The results obtained from the molecular tailoring approach may thus be expected to enhance our understanding of the environment in the crystalline material with reasonably small computational effort.Contribution to the Jacopo Tomasi Honorary Issue     

Ab Initio Calculation of structures and properties of halogenated general anesthetics: halothane and sevoflurane

To correctly analyze the effects of general anesthetics on their potential targets by large‐scale molecular simulation, the structural parameters and partial atomic charges of the anesthetics are of determinant importance. Geometric optimizations using the Hartree–Fock and the B3LYP density functional theory methods with the large 6‐311+G(2d,p) basis set were performed to determine the structures and charge distributions of two halogenated anesthetics, 2‐bromo‐2‐chloro‐1,1,1‐trifluoroethane (halothane) and fluoromethyl‐2,2,2,‐trifluoro‐1‐(trifluoromethyl) ethyl ether (sevoflurane). The calculated bond lengths and angles are within 3% of the corresponding experimental values reported for the similar molecular groups. Charges are assigned using the Mulliken population analysis and the electrostatic potential (ESP) based on the Merz–Kollman–Singh scheme. The atoms‐in‐molecules (AIM) theory is also used to assign the charges in halothane. The dipole moments calculated with the Mulliken population analysis and ESP for the structures optimized by B3LYP/6‐311+(2d,p) were respectively 1.355 and 1.430 D for halothane and 2.255 and 2.315 D for sevoflurane. These are in excellent agreement with the experimental values of 1.41 and 2.33 D for halothane and sevoflurane, respectively. The calculated structures and partial charge distributions can be readily parameterized for molecular mechanics and molecular dynamics simulations involving these halogenated agents. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 436–444, 2001     

Accurate conformation‐dependent molecular electrostatic potentials for high‐throughput in silico drug discovery

The atom‐centered partial charges‐approximation is commonly used in current molecular modeling tools as a computationally inexpensive alternative to quantum mechanics for modeling electrostatics. Even today, the use of partial charges remains useful despite significant advances in improving the efficiency of ab initio methods. Here, we report on new parameters for the EEM and SFKEEM electronegativity equalization‐based methods for rapidly determining partial charges that will accurately model the electrostatic potential of flexible molecules. The developed parameters cover most pharmaceutically relevant chemistries, and charges obtained using these parameters reproduce the B3LYP/cc‐pVTZ reference electrostatic potential of a set of FDA‐approved drug molecules at best to an average accuracy of 13 ± 4 kJ mol^?1; thus, equipped with these parameters electronegativity equalization‐based methods rival the current best non‐quantum mechanical methods, such as AM1‐BCC, in accuracy, yet incur a lower computational cost. Software implementations of EEM and SFKEEM, including the developed parameters, are included in the conformer‐generation tool BALLOON , available free of charge at http://web.abo.fi/fak/mnf/bkf/research/johnson/software.php . © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

VESPA: A new,fast approach to electrostatic potential-derived atomic charges from semiempirical methods

An improved semiempirical method for computing electrostatic potential-derived atomic charges is described. It includes a very fast algorithm for the generation of the grid points around the molecule and the calculation of the electrostatic potential at these points. The dependency of the atomic point charges obtained on the number of grid points used in the fitting procedure is examined. For “buried” atoms a high density grid is necessary. It is possible to obtain 6–31G*-quality atom-centered point charges, even for phosphorus compounds, using AM1 or PM3. This approach can therefore be recommended for general use in QSAR or molecular mechanics for any organic and bioorganic system up to about 200 atoms. © 1997 by John Wiley & Sons, Inc. J Comput Chem 18: 744–756, 1997     

On the representation of electrostatic fields around ab initio charge distributions

Summary We compare two methods (Mulliken charges and a distributed multipole analysis, DMA) of representing an ab initio charge distribution for calculating the electrostatic field and potential outside the molecule, using pyrimidine and the RNA base uracil as examples. This is done using a 3-D graphical display of the electrostatic fields, which, when used with real-time rotation, zooming and clipping, has many advantages for qualitatively assessing the electrostatic interactions of a molecule. The errors involved in using Mulliken point charges may be of similar magnitude to the total electrostatic field in regions which are important in recognition processes. The DMA representation automatically includes the anisotropic electrostatic effects of non-spherical features in the charge distribution of each atom, and yet the displayed electrostatic fields around the atoms which have lone-pair density do not show marked anisotropy.     

Comparison of methods for point‐charge representation of electrostatic fields

The calculation of the electrostatic potential resulting from an infinite or extended array of charges in the interior of a region of interest is a frequent task in computational chemistry. In case of a periodic potential this can, for example, be done by Ewald summation or by multipole methods. An important alternative are those methods where arrays of auxiliary point charges are optimized with respect to charge and/or position to reproduce the original electrostatic potential. In the literature different variations are reported. We compare the performance of some of these with respect to their ability to reproduce the original potential and the computational effort required. Between (1) surface charges determined by the conductor‐boundary condition, (2) optimized surface charges, and (3) surface charges floating on the surface we find that (2) offers good quality with small computational costs involved. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

Strategic approaches to drug design. I. An integrated software framework for molecular modelling

Summary An integrated molecular graphics and computational chemistry framework is described which has been designed primarily to handle small molecules of up to 300 atoms. The system provides a means of integrating software from any source into a single framework. It is split into two functional subsystems. The first subsystem, called COSMIC. runs on low-cost, serial-linked colour graphics terminals and allows the user to prepare and examine structural data and to submit them for extensive computational chemistry. Links also allow access to databases, other modelling systems and user-written modules. Much of the output from COSMIC cannot be examined with low level graphics. A second subsystem, called ASTRAL, has been developed for the high-resolution Evans & Sutherland PS300 colour graphics terminal and is designed to manipulate complex display structures. The COSMIC minimisers, geometry investigators, molecular orbital displays, electrostatic isopotential generators and various interfaces and utilities are described.     

A practical procedure for the determination of electrostatic charges of large molecules

Summary A practical procedure for the precise determination of electrostatic charges, which are evaluated by fitting the rigorous quantum mechanical molecular electrostatic potential to a monopole-monopole expression, is presented. The proposal of this procedure arises from the study of the minimum requirements necessary to obtain reliable electrostatic charges. Such a study is focused on: (i) the dependence of the electrostatic charges on the set of points where the quantum mechanical and the monopole-monopole molecular electrostatic potentials are fitted; thus, both the influence of the number of points and their distribution in layers located out of the van der Waals radii of the atoms are examined, and (ii) the reliability of the use of fractional models for the evaluation of electrostatic charges of large molecules. Results point out that the optimum number of points is defined by a density of points ranging from 0.45 to 0.60 points/Ų when four layers (separated by 0.2 Å) are considered. Nevertheless, the use of only two layers (separated by 0.4 Å) for large molecules is recommended, thus enabling one to obtain reliable charges at a reduced computational cost. Moreover, results justify the use of fractional models for the determination of electrostatic charges of extremely large molecules, even when aromatic structures exist.     

Fully polarizable QM/MM calculations: An application to the nonbonded iodine–oxygen interaction in dimethyl‐2‐iodobenzoylphosphonate

The compound dimethyl‐2‐iodobenzoylphosphonate is unusual in that it forms well‐ordered crystals that clearly show short iodine‐oxygen interactions in which both the iodine and the oxygen are in their normal oxidation states. These interactions were studied using a new hybrid quantum mechanical–molecular mechanical approach that employs a polarizable molecular mechanics component. The electric field at the molecular mechanics atoms was calculated from a distributed multipole expansion of the wave function; this induced dipoles on the molecular mechanics atoms. The electrostatic potential in a spherical shell around the induced dipoles was reproduced through induced charges on the atomic center and those bonded to it using an analytical (rather than numerical) procedure. The new atomic charges (induced charges plus permanent charges) were then able to interact with the quantum mechanical entity and polarize the wave function. The procedure was iterated to convergence. The calculations show that the iodine atom becomes more positive in the crystal environment (modeled by a chain of three molecules of dimethyl‐2‐iodobenzoylphosphonate). Thus, while the cooperative effects of the crystal environment may not be the only feature stabilizing this unusual interaction, they do play a significant role in reducing the otherwise unfavourable iodine–oxygen monopole–monopole interaction. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 478–482, 2000     

Toward a new approach for determination of solute's charge distribution to analyze interatomic electrostatic interactions in quantum mechanical/molecular mechanical simulations

We present an alternative approach to determine "density-dependent property"-derived charges for molecules in the condensed phase. In the case of a solution, it is essential to take into consideration the electron polarization of molecules in the active site of this system. The solute and solvent molecules in this site have to be described by a quantum mechanical technique and the others are allowed to be treated by a molecular mechanical method (QM/MM scheme). For calculations based on this scheme, using the forces and interaction energy as density-dependent property our charges from interaction energy and forces (CHIEF) approach can provide the atom-centered charges on the solute atoms. These charges reproduce well the electrostatic potentials around the solvent molecules and present properly the picture of the electron density of the QM subsystem in the solution system. Thus, the CHIEF charges can be considered as the atomic charges under the conditions of the QM/MM simulation, and then enable one to analyze electrostatic interactions between atoms in the QM and MM regions. This approach would give a view of the QM nuclei and electrons different from the conventional methods.     

Near ab initio quality molecular electrostatic potential maps using hybridization displacement charges at PM3 level and effects of geometrical changes in amino groups

Charge distributions, dipole moments, and molecular electrostatic potentials (MEP) around several molecules consisting of carbon, nitrogen, oxygen, fluorine, sulfur, and chlorine atoms were studied using the PM3 semiempirical method and the results compared with those obtained using ab initio calculations at the RHF/6‐31G** level. Thus it is shown that relative MEP values near different atoms can be obtained using hybridization displacement charges (HDC) obtained by employing PM3 density matrices that usually agree quite satisfactorily with the ab initio ones. Further, positions of ab initio MEP minima are correctly located and the corresponding relative MEP values usually correctly predicted using the PM3(HDC) charges distributed continuously in three dimensions according to the forms of squares of valence s atomic orbitals. The necessary parameters for HDC calculations using the PM3 method were optimized. It is shown how within the frameworks of both PM3 and AM1 methods the π electrons or lone pairs associated with amino group nitrogen atoms and ring atoms can be satisfactorily treated in different situations. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 82: 299–312, 2001     

A numerically stable restrained electrostatic potential charge fitting method

Inspired by the idea of charge decomposition in calculation of the dipole preserving and polarization consistent charges (Zhang et al., J. Comput. Chem. 2011, 32, 2127), we have proposed a numerically stable restrained electrostatic potential (ESP)‐based charge fitting method for protein. The atomic charge is composed of two parts. The dominant part is fixed to a predefined value (e.g., AMBER charge), and the residual part is to be determined by restrained fitting to residual ESP on grid points around the molecule. Nonuniform weighting factors as a function of the dominant charge are assigned to the atoms. Because the residual part is several folds to several orders smaller than the dominant part, the impact of ill‐conditioning is alleviated. This charge fitting method can be used in quantum mechanical/molecular mechanical (QM/MM) simulations and similar studies, where QM calculated electronic properties are frequently mapped to partial atomic charges. © 2012 Wiley Periodicals, Inc.     

Exploring aromatic chemical space with NEAT: novel and electronically equivalent aromatic template

In this paper, we describe a lead transformation tool, NEAT (Novel and Electronically equivalent Aromatic Template), which can help identify novel aromatic rings that are estimated to have similar electrostatic potentials, dipoles, and hydrogen bonding capabilities to a query template; hence, they may offer similar bioactivity profiles. In this work, we built a comprehensive heteroaryl database, and precalculated high-level quantum mechanical (QM) properties, including electrostatic potential charges, hydrogen bonding ability, dipole moments, chemical reactivity, and othe properties. NEAT bioisosteric similarities are based on the electrostatic potential surface calculated by Brood, using the precalculated QM ESP charges and other QM properties. Compared with existing commercial lead transformation software, (1) NEAT is the only one that covers the comprehensive heteroaryl chemical space, and (2) NEAT offers a better characterization of novel aryl cores by using high-evel QM properties that are relevant to molecular interactions. NEAT provides unique value to medicinal chemists quickly exploring the largely uncharted aromatic chemical space, and one successful example of its application is discussed herein.     

Protein‐specific force field derived from the fragment molecular orbital method can improve protein–ligand binding interactions

Accurate computational estimate of the protein–ligand binding affinity is of central importance in rational drug design. To improve accuracy of the molecular mechanics (MM) force field (FF) for protein–ligand simulations, we use a protein‐specific FF derived by the fragment molecular orbital (FMO) method and by the restrained electrostatic potential (RESP) method. Applying this FMO‐RESP method to two proteins, dodecin, and lysozyme, we found that protein‐specific partial charges tend to differ more significantly from the standard AMBER charges for isolated charged atoms. We did not see the dependence of partial charges on the secondary structure. Computing the binding affinities of dodecin with five ligands by MM PBSA protocol with the FMO‐RESP charge set as well as with the standard AMBER charges, we found that the former gives better correlation with experimental affinities than the latter. While, for lysozyme with five ligands, both charge sets gave similar and relatively accurate estimates of binding affinities. © 2013 Wiley Periodicals, Inc.     

Quantum chemical study of electronic and structural properties of retinal and some aromatic analogs

The electronic and structural properties of retinal and four analogs were studied using semiempirical, ab initio Hartree-Fock, and density functional theory methods with the aim to evaluate the effects caused by some structural modifications in the ring bound to the polyenic chain and compared with the all-E-trans-retinal molecule. Therefore, some properties such as bond lengths, bond angles, atomic charges derived from electrostatic potential charges from electrostatic potential using grid based method as well as frontier orbitals of the polyenic chain were analyzed. Furthermore, the transition energies of the molecules were also calculated using the Zerner's intermediate neglect of differential overlap-spectroscopic, time-dependent Hartree-Fock, and time-dependent density functional theory methods. The results indicate that in spite of the structural modifications of retinal derivatives in comparison with all-E-trans-retinal, their properties seem similar. Thus, these molecules may behave similarly to all-E-trans-retinal and possibly be attempted in the search of novel molecular devices.     

Building molecular charge distributions from fragments: Application to HIV-1 protease inhibitors

Interaction energies are a function of the molecular charge distribution. In previous work, we found that the set of atomic partial charges giving the best agreement with experimental vacuum dipole moments were from density functional theory calculations using an extended basis set. Extension of such computations to larger molecules requires an atomic partial charge calculation beyond present computational resources. A solution to this problem is the calculation of atomic partial charges for segments of the molecule and reassociation of such fragments to yield partial charges for the entire molecule. Various partitions and reassociation methods for five molecules relevant to HIV-1 protease inhibitors are examined. A useful method of reassociation is introduced in which atomic partial charges for a large molecule are computed by fitting to the combined electrostatic potential calculated from the fragment partial charges. As expected, the best sites for partitions are shown to be carbon—carbon rather than carbon—nitrogen bonds. © 1997 by John Wiley & Sons, Inc.     

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