Fuzzy fragment selection strategies,basis set dependence and HF–DFT comparisons in the applications of the ADMA method of macromolecular quantum chemistry

The Adjustable Density Matrix Assembler (ADMA) method is examined in this paper. This method approximates (first order) density matrices by taking only those interactions into account which are present between fragments separated by a preset distance parameter (d). The accuracy of this approximation is tested using different basis sets, distance parameters and exchange‐correlation functionals. As an illustration of the applications of the method, the electron density of the hemoglobin molecule is presented in its oxy, deoxy and carbon‐monoxy form. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Toward similarity measures for macromolecular bodies: Medla test calculations for substituted benzene systems

A new method is proposed for the evaluation of numerical similarity measures for large molecules, defined in terms of their electron density (ED) distributions. The technique is based on the Molecular Electron Density Lego Assembler (MEDLA) approach, proposed earlier for the generation of ab initio quality electron densities for proteins and other macromolecules. The reliability of the approach is tested using a family of 13 substituted aromatic systems for which both standard ab initio electron density computations and the MEDLA technique are applicable. These tests also provide additional examples for evaluating the accuracy of the MEDLA technique. Electron densities for a series of 13 substituted benzenes were calculated using the standard ab initio method with STO-3G, 3-21G, and 6-31G** basis sets as well as the MEDLA approach with a 6-31G** database of electron density fragments. For each type of calculation, pairwise similarity measures of these compounds were calculated using a point-by-point numerical comparison of the EDs. From these results, 2D similarity maps were constructed, serving as an aid for quick visual comparisons for the entire molecular family. The MEDLA approach is shown to give virtually equivalent numerical similarity measures and similarity maps as the standard ab initio method using a 6-31G** basis set. By contrast, significant differences are found between the standard ab initio 6-31G** results and the standard ab initio results obtained with smaller STO-3G and 3-21G basis sets. These tests indicate that the MEDLA-based similarity measures faithfully mimic the actual, standard ab initio 6-31G** similarity measures, suggesting the MEDLA method as a reliable technique to assess the shape similarities of proteins and other macromolecules. The speed of the MEDLA computations allows rapid, pairwise comparisons of the actual EDs for a series of molecules, requiring no more computer time than other simplified, less detailed representations of molecular shape. The MEDLA method also reduces the need to store large volumes of numerical density data on disk, as these densities can be quickly recomputed when needed. For these reasons, the proposed MEDLA similarity analysis technique is likely to become a useful tool in computational drug design. © 1995 John Wiley & Sons, Inc.     

Holographic electron density shape theorem and its role in drug design and toxicological risk assessment

Each complete, boundaryless molecular electron density is fully determined by any nonzero volume piece of the electron density cloud. This inherent feature of molecules, called the "holographic" property of molecular electron densities, provides a strong foundation for the local, quantum chemical shape analysis of various functional groups, pharmacophores, and other local molecular moieties. A proof is presented for the relevant molecular shape theorem, the "holographic electron density shape theorem", and the role of this theorem in quantum chemical, quantitative shape-activity relations (QShAR) is discussed. The quantum chemical methods of molecular shape analysis can be extended to ab initio quality electron densities of macromolecules, such as proteins, as well as to local molecular moieties, such as functional groups or pharmacophores, based on the transferability and additivity of local, fuzzy density fragments and the associated local density matrixes within the framework of the ADMA (Adjustable Density Matrix Assembler) approach. In addition to new results on chemical bonding and the development of macromolecular force methods, the new methodologies are also applicable to QShAR studies in computer-aided drug discovery and in toxicological risk assessment.     

An alternative to the “Star Path” enhancement of the ADMA linear scaling method for protein modeling

With the aim of improving the performance of macromolecular quantum chemistry conformation analysis and reaction path following methods, the Adjustable Density Matrix Assembler (ADMA) method has already been combined with some faster although less accurate density matrix extrapolation methods, such as the Löwdin‐Inverse‐Löwdin (LIL) extrapolation along a potential energy surface, and a strategically arranged back‐and‐forth switching between these methods has been proven to be advantageous. Here, an alternative approach is proposed and investigated, based on several actual test calculations, where the “inexpensive” LIL density matrix extrapolation steps are replaced by only somewhat more expensive, but still ADMA‐based calculations, where in the “rough‐search stage,” only interactions of shorter distances within the macromolecule are considered. It is shown that this approach is viable, as an alternative to the “Star Path” method including both ADMA and LIL steps. © 2017 Wiley Periodicals, Inc.     

Integration of Graph Theory and Quantum Chemistry for Structure-Activity Relationships

Abstract

The objective of this article is to outline both graph-theoretically based and quantum chemically based structural indices of potential use in quantitative structure activity correlations. We consider graph-theoretical indices such as the connectivity index, topological index, Wiener index and molecular ID indices. Several structural and geometry-dependent indices can be derived from semiempirical and ab initio quantum calculations based on the charge densities, overlap matrices, frontier orbitals, molecular hardness, free valence, density matrices, quantum spectral difference indices, quantum spectral indices and bond matrices. Finally, the use of electrostatic potentials and charge densities for the prediction of reactive sites will be discussed.     

Theoretical studies of molecular conformation. Derivation of an additive procedure for the computation of intramolecular interaction energies. Comparison withab initio SCF computations

An additive procedure (SIBFA) is developed for the rapid computation of conformational energy variations in very large molecules. The macromolecule is built out of constitutive molecular fragments and the intramolecular energy is computed as a sum of interaction energies between the fragments. The electrostatic and the polarization components are calculated using multicenter multipole expansions of theab initio SCF electron density of the fragments. The repulsion component is obtained as a sum of bond and lone pair interactions.Tests of the procedure on a series of model compounds containing ether oxygens and pyridine-like nitrogens are reported and compared with the results of correspondingab initio SCF calculations. The resulting methodology is compatible with the simultaneous computation of intermolecular interactions.     

Search for a density-based alternative quantum mechanics of many-electron systems

There are three reasons for seeking an alternative density-based quantum mechanics of many-electron systems, incorporating both interpretive and basic quantum mechanical aspects: (i) failure of popularad hoc chemical concepts underab initio scrutiny; (ii) failure ofab initio calculations to provide simple concepts; and (iii) highly attractive concepts and pictures generated by the electron density in three-dimensional space. At present the three interlinked pillars for such a density mechanics (in contrast to wave mechanics) are: (a) density functional theory; (b) quantum fluid dynamics; and (c) property densities in three-dimensional space. This article describes several studies dealing with these aspects. Although a density mechanics may well be an impossible ideal to realize, the search for it is indeed rejuvenating the whole of quantum chemistry.     

Bond fukui indices: Comparison of frozen molecular orbital and finite differences through mulliken populations

Bond Fukui functions and matrices are introduced for ab initio levels of theory using a Mulliken atoms in molecules model. It is shown how these indices may be obtained from first‐order density matrix derivatives without need for going to second‐order density matrices as in a previous work. The importance of taking into account the nonorthogonality of the basis in ab initio calculations is shown, contrasting the present results with previous work based on Hückel theory. It is shown how the extension of Fukui functions to Fukui matrices allows getting more insight into the nature of bond Fukui functions. All presently introduced indices respect the necessary normalization conditions and include the classical single atom condensed Fukui functions. © 2013 Wiley Periodicals, Inc.     

Accurate density functional theory description of binding constants and NMR chemical shifts of weakly interacting complexes of C60 with corannulene‐based molecular bowls

Density functional calculations on “catch and release” complexes of C₆₀ with corannulene derived molecular bowls show that computationally obtained ¹H nuclear magnetic resonance (NMR) chemical shifts can be used as a reliable predictor of binding constants. A wide range of functionals was benchmarked against accurate ab initio calculations to ensure a credible representation of the weak forces that dominate the interactions in these systems. The most reliable density functional theory (DFT) results were then calibrated using experimentally observed NMR data. Careful analysis and comparison of a wide range of commonly used density functionals shows that the explicit inclusion of dispersion corrections is currently the only reliable way to accurately describe the systems investigated in our study. Moreover, we are able to show that the B97‐D and ωB97X‐D functionals are not only able to reproduce ab initio benchmark calculations, but they do so accurately with a moderately sized basis sets and without the problems of numerical integration we encountered with other functionals in this study. © 2013 Wiley Periodicals, Inc.     

Dipole moment derivatives and integrated intensities for the vibrational transitions of N2 … HF

The integrated intensity of vibrational transitions depends on the magnitude of the derivatives of the dipole with respect to nuclear motion. These derivatives are usually obtained by time-consuming ab initio calculations. In this paper we apply a long-range model, based on distributed schemes for describing the charge densities and polarizabilities of molecules, to the prediction of dipole derivatives and infrared intensities for the N₂ … HF complex. The results are found to agree qualitatively with full ab initio self-consistent field calculations. © 1996 by John Wiley & Sons, Inc.     

Quantum chemical and physicochemical influences on structure–activity relations and drug design

Quantum chemical results will be presented on drugs, carcinogens, teratogens, and endogenous biomolecules using our new nonempirical ab initio MODPOT /VRDDO method, which incorporates as options to our ab initio LCAO -MO -SCF /CI programs ab initio effective core model potentials (MODPOT ) permitting one to calculate only the valence electrons explicitly yet accurately and an integral prescreening technique (VRDDO , variable retention of diatomic differential overlap) especially effective for spatially extended molecules. For molecules of the size of those of interest the MODPOT /VRDDO calculations run an order-of-magnitude faster than with our own fast ab initio programs and still retain accuracy to the third decimal place for the valence electron properties. We have also just implemented a new efficient MERGE technique which allows us to reuse integrals from a common skeletal fragment and only to have to recalculate those for a new atom or group or a change in its position. Examples will be presented of the use of this technique on a carcinogenic polycyclic aromatic hydrocarbon and its metabolites. The pK_a's, oil-water partition, and drug distribution coefficients as a sensitive function of pH have been measured for a number of drugs as well as for relevant endogenous biomolecules. The pH dependence of the lipophilicities of such molecules has profound implication on appropriate use of such data in QSAR studies.     

Electron Density Calculations for Crystals with a NaCl Lattice

The ab initio local DFT method combined with the ab initio pseudopotential technique is used to calculate valence electron densities in series of crystal substances MA (M = Li, Na, K, Rb, Ag, Mg, Ca; A = F, Cl, Br, O, S) with a NaCl lattice. Systematic variations of valence electron density depending on the atomic number of anion and cation have been found. According to electron density distribution, the compounds are classified into three groups: a) oxides and fluorides; b) sulfides, chlorides, and bromides; c) noble metal halides. In oxides and fluorides, the maximum of valence electron density is in the middle of the M–M bond. In sulfides, chlorides, and bromides, the minimum density is in the middle of the M–M bond, with two symmetric maxima or two shoulders (depending on the atomic number of the cation) lying away from the center of the bond. In noble metal halides, the maximum of valence density is on the metal due to the presence of metal d-states, and the density map is rotated through 90° relative to the map of alkali and alkali earth metals, so that the Hal–Hal bond becomes an analog of the M–M bond.     

Molecular symmetry. IV. The coupled perturbed Hartree–Fock method

Symmetry methods employed in the ab initio polyatomic program HONDO are extended to the coupled perturbed Hartree–Fock (CPHF) formalism, a key step in the analytical computation of energy first derivatives for configuration interaction (CI) wavefunctions, and energy second derivatives for Hartree–Fock (HF) wavefunctions. One possible computational strategy is to construct Fock-like matrices for each nuclear coordinate in which the one- and two-electron integrals of the usual Fock matrix are replaced by the integral first derivatives. “Skeleton” matrices are constructed from the unique blocks of electron-repulsion integral derivatives. The correct matrices are generated by applying a symmetrization operator. The analysis is valid for many wavefunctions, including closed- or open-shell spin-restricted and spin-unrestricted HF wavefunctions. To illustrate the method, we compare the computer time required for setting up the coupled perturbed HF equations for eclipsed ethane using D_3h symmetry point group and various subgroups of D_3h. Computational times are roughly inversely proportional to the order of the point group.     

Accurate modeling of the intramolecular electrostatic energy of proteins

The (?, ψ) energy surface of blocked alanine (N-acetyl–N′-methyl alanineamide) was calculated at the Hartree-Fock (HF)/6-31G* level using ab initio molecular orbital theory. A collection of six electrostatic models was constructed, and the term electrostatic model was used to refer to (1) a set of atomic charge densities, each unable to deform with conformation; and (2) a rule for estimating the electrostatic interaction energy between a pair of atomic charge densities. In addition to two partial charge and three multipole electrostatic models, this collection includes one extremely detailed model, which we refer to as nonspherical CPK. For each of these six electrostatic models, parameters—in the form of partial charges, atomic multipoles, or generalized atomic densities—were calculated from the HF/6-31G* wave functions whose energies define the ab initio energy surface. This calculation of parameters was complicated by a problem that was found to originate from the locking in of a set of atomic charge densities, each of which contains a small polarization-induced deformation from its idealized unpolarized state. It was observed that the collective contribution of these small polarization-induced deformations to electrostatic energy differences between conformations can become large relative to ab initio energy differences between conformations. For each of the six electrostatic models, this contribution was reduced by an averaging of atomic charge densities (or electrostatic energy surfaces) over a large collection of conformations. The ab initio energy surface was used as a target with respect to which relative accuracies were determined for the six electrostatic models. A collection of 42 more complete molecular mechanics models was created by combining each of our six electrostatic models with a collection of seven models of repulsion + dispersion + intrinsic torsional energy, chosen to provide a representative sample of functional forms and parameter sets. A measure of distance was defined between model and ab initio energy surfaces; and distances were calculated for each of our 42 molecular mechanics models. For most of our 12 standard molecular mechanics models, the average error between model and ab initio energy surfaces is greater than 1.5 kcal/mol. This error is decreased by (1) careful treatment of the nonspherical nature of atomic charge densities, and (2) accurate representation of electrostatic interaction energies of types 1—2 and 1—3. This result suggests an electrostatic origin for at least part of the error between standard model and ab initio energy surfaces. Given the range of functional forms that is used by the current generation of protein potential functions, these errors cannot be corrected by compensating for errors in other energy components. © 1995 by John Wiley & Sons, Inc.     

Multiobjective evolutionary algorithm with many tables for purely ab initio protein structure prediction

This article focuses on the development of an approach for ab initio protein structure prediction (PSP) without using any earlier knowledge from similar protein structures, as fragment‐based statistics or inference of secondary structures. Such an approach is called purely ab initio prediction. The article shows that well‐designed multiobjective evolutionary algorithms can predict relevant protein structures in a purely ab initio way. One challenge for purely ab initio PSP is the prediction of structures with β‐sheets. To work with such proteins, this research has also developed procedures to efficiently estimate hydrogen bond and solvation contribution energies. Considering van der Waals, electrostatic, hydrogen bond, and solvation contribution energies, the PSP is a problem with four energetic terms to be minimized. Each interaction energy term can be considered an objective of an optimization method. Combinatorial problems with four objectives have been considered too complex for the available multiobjective optimization (MOO) methods. The proposed approach, called “Multiobjective evolutionary algorithms with many tables” (MEAMT), can efficiently deal with four objectives through the combination thereof, performing a more adequate sampling of the objective space. Therefore, this method can better map the promising regions in this space, predicting structures in a purely ab initio way. In other words, MEAMT is an efficient optimization method for MOO, which explores simultaneously the search space as well as the objective space. MEAMT can predict structures with one or two domains with RMSDs comparable to values obtained by recently developed ab initio methods (GAPF_CG, I‐PAES, and Quark) that use different levels of earlier knowledge. © 2013 Wiley Periodicals, Inc.     

Comparison of NDDO and quasi-ab initio approaches to compute semiempirical molecular electrostatic potentials

The suitability of the two most widely used strategies to compute semiempirical MEPs is examined. For this purpose, MEP minima, electrostatic charges, and dipoles for a large number of molecules were computed at the AM1, MNDO, and PM3 levels using both the NDDO strategy developed by Ferenczy, Reynolds, and Richards and our own quasi-ab initio method. Results demonstrate that the quasi-ab initio is preferred over the NDDO method for the computation of MEP minima. It is also found that the best set of semiempirical charges and dipoles are obtained using either the AM1 NDDO or the MNDO quasi-ab initio methods. In these two cases, the quality of the results is fully comparable with 6-31G* values. © 1994 by John Wiley & Sons, Inc.     

A density functional treatment of organolithium compounds: Comparison to ab initio,semiempirical, and experimental results

Density functional calculations on several classes of organolithium compounds are described. The compounds studied include lithium bonds to carbon, oxygen, and nitrogen and are representative of most types of organolithium compounds that have appeared in the recent literature. The computational results are compared to those using MNDO, which has been shown to have some serious deficiencies in compounds involving carbon–lithium bonds, and to PM3 results, which offer some improvement over MNDO for many organolithium compounds. Most of the density functional calculations with a large basis set are in good agreement with available ab initio and experimental data. Calculated carbon–lithium bond lengths were slightly shorter than those calculated by other ab initio methods and were substantially longer than those calculated by MNDO, which is known to underestimate carbon–lithium bond lengths severely. Dimerization energies of methyllithium, calculated by DMol, were also in good agreement with those of other ab initio calculations. Lithium–nitrogen bonds in lithium amides were calculated to be slightly shorter by DMol than by MNDO, although the two methods were in qualitative agreement for this type of compound. © 1995 by John Wiley & Sons, Inc.     

A comparative study of the structure and bonding of HOO,HOS, HSO,and HSS radicals by CNDO/2 and INDO methods

Equilibrium geometries, force constants, barriers to linearity, charge distributions, dipole moments, and electron spin density of HOO, HOS, HSO, and HSS radicals are calculated by CNDO/2 and INDO methods using respectively the original and some recently introduced scheme of parametrization. Three sets of calculations, namely, CNDO/2(sp), CNDO/2(spd), and INDO, are performed, and the results are compared with the ab initio and experimental values, wherever available. A good agreement is obtained for geometry in the case of CNDO/2 (sp) and INDO calculations. The performance of CNDO/2 (spd) calculations in this regard is quite unreliable. The stretching force constants are considerably overestimated by all the methods, while the bending force constants are in reasonable agreement with the ab initio values. With respect to dipole moments, the CNDO/2 values are in better agreement with the ab initio results than the INDO values. In all the cases, the dipole moment vector directions are in complete disagreement with the ab initio predictions.     

Ab initio quasirelativistic calculations on electronic transitions in ICl by the multireference many‐body perturbation theory

A quasirelativistic perturbative method of ab initio calculations on ground and excited molecular electronic states and transition properties within the relativistic effective core potential approximation is presented and discussed. The method is based on the construction of a state‐selective many‐electron effective Hamiltonian in the model space spanned by an appropriate set of Slater determinants by means of the second‐order many‐body multireference perturbation theory. The neglect of effective spin–orbit interactions outside of the model space allows the exploitation of relatively high nonrelativistic symmetry during the evaluation of perturbative corrections and therefore dramatic reduction of the cost of computations without any contraction of the model‐space functions. One‐electron transition properties are evaluated via the perturbative construction of spin‐free transition density matrices. Illustrative calculations on the X0⁺ ? A1, B0⁺, and (ii)1 transitions in the ICl molecule are reported. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem, 2002