Exploring the conformation,charge density distribution and the electrostatic properties of galanthamine molecule in the active site of AChE using DFT and AIM theory

Quantum chemical calculations and charge density analysis were carried out to understand the geometry, charge density distribution and the electrostatic properties of isolated galanthamine molecule (form I) and for the same lifted out from the active site (form II) of AChE. The optimized geometry of isolated galanthamine was obtained from a hybrid density functional theory (B3LYP/6‐311G**) calculation. A docking analysis on galanthamine with AChE was performed, and the lowest docked energy structure was selected from the active site of AChE for the further study. A single point energy quantum chemical calculation (B3LYP/6‐311G**) was carried out for the lowest energy structure, which was lifted from the galanthamine–AChE complex from molecular docking analysis. The structural comparison between (I) and (II) helps to understand the conformational modification of the galanthamine molecule in the active site. When the molecule present in the active site, the molecular geometry is seen to be significantly altered, specifically, large changes were observed in the outer core of the molecule while the inner core geometry is intact. The bond topological and electrostatic properties of (I) and (II) were calculated. The dipole moment of the galanthamine molecule also increases from 2.09 to 2.67 D in the process. A large negative electrostatic potential region is found at the vicinity of oxygen and nitrogen atoms of the molecule, which predominantly involve strong hydrophobic and electrostatic interactions with the amino acid residues TRP84, PHE330, GLY118, TYR70, and SER122 present in the active site of AChE. © 2013 Wiley Periodicals, Inc.     

pH‐responsive ampholytic terpolymers of acrylamide,sodium 3‐acrylamido‐3‐methylbutanoate,and (3‐acrylamidopropyl)trimethylammonium chloride. II. Solution properties

The solution properties of low‐charge‐density ampholytic terpolymers of acrylamide, sodium 3‐acrylamido‐3‐methylbutanoate, and (3‐acrylamidopropyl)trimethylammonium chloride were studied as functions of the solution pH, ionic strength, and polymer concentration. Terpolymers with low charge densities, large charge asymmetries, or both exhibited excellent solubility in deionized (DI) water, and higher charge density terpolymers were readily dispersible in DI water; however, the higher charge density terpolymer solutions separated into polymer‐rich and polymer‐poor phases upon standing over time. Charge‐balanced terpolymers exhibited antipolyelectrolyte behavior at pH values greater than or equal to the ambient pH (6.5 ± 0.2); the same terpolymers behaved increasingly as cationic polyelectrolytes with decreasing solution pH because of the protonation of the 3‐acrylamido‐3‐methylbutanoate (AMB) repeat units. Unbalanced terpolymers generally exhibited polyelectrolyte behavior, although the effects of intramolecular electrostatic attractions (i.e., polyampholyte effects) on the hydrodynamic volume of the unbalanced terpolymer coils were evident at certain values of the solution pH and salt concentration. The dilute‐solution behavior of the terpolymers correlated well with the behavior predicted by several polyampholyte solution theories. In the semidilute regime, solution viscosities increased with increasing terpolymer charge density, and this indicated a significant enhancement of the solution viscosity by intermolecular electrostatic associations. Upon the addition of NaCl, semidilute‐solution viscosities tended to decrease because of the disruption of the intermolecular electrostatic associations. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 3252–3270, 2004     

Destructive Effect of an α‐Particle on TEX

By computational studies it was found that interaction of an α‐particle with TEX (an explosive molecule having a cage structure) disturbs the molecular skeleton dramatically. As the α‐particle exerts its positive electrostatic field, the C–C bond between the dioxalane rings present in TEX breaks down forming the dication of TEX. The calculation were done within the constraints of the density functional theory at the level of B3LYP using various bases sets, mainly 6‐31++G(d,p).     

(E)‐Methyl 2‐anilinomethylene‐3‐oxobutanoate: X‐ray and density functional theory studies

Molecules of the title compound, C₁₂H₁₃NO₃, are not planar and are stabilized by electrostatic interactions, as the dipole moment of the molecule is 3.76 D. They are also stabilized by intramolecular hydrogen bonds of N...O and C...O types, and by a complicated network of weak intermolecular hydrogen bonds of the C...O type. This paper also reports the theoretical investigation of the hydrogen bonding and electronic structure of the title compound using natural bond orbital (NBO) analysis.     

Intermolecular complexes of HArF and HKrF with CO

Stable linear weakly bound hydrogen-bonded complexes of HArF and HKrF with the CO molecule have been predicted by ab initio computations at the MP2/6- 311+ +G(2d,2p) level of theory. The complexes, having stabilities in the order, FArH...CO>FKrH...CO>FArH...OC>FKrH...OC are compared. They exhibit unusual blueshifts of the Ar-H (Kr-H) stretching frequency, as well as contraction of the Ar-H (Kr-H) bond, and these effects are rationalized mainly in terms of shifts in the electron density of HArF (HKrF) on complexation, caused mainly by a combination of the intermolecular electrostatic interaction, electron-electron (Pauli) and nuclear-nuclear repulsion and charge density transferred from the CO molecule to the rare-gas-containing molecule.     

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.     

Towards a force field based on density fitting

Total intermolecular interaction energies are determined with a first version of the Gaussian electrostatic model (GEM-0), a force field based on a density fitting approach using s-type Gaussian functions. The total interaction energy is computed in the spirit of the sum of interacting fragment ab initio (SIBFA) force field by separately evaluating each one of its components: electrostatic (Coulomb), exchange repulsion, polarization, and charge transfer intermolecular interaction energies, in order to reproduce reference constrained space orbital variation (CSOV) energy decomposition calculations at the B3LYP/aug-cc-pVTZ level. The use of an auxiliary basis set restricted to spherical Gaussian functions facilitates the rotation of the fitted densities of rigid fragments and enables a fast and accurate density fitting evaluation of Coulomb and exchange-repulsion energy, the latter using the overlap model introduced by Wheatley and Price [Mol. Phys. 69, 50718 (1990)]. The SIBFA energy scheme for polarization and charge transfer has been implemented using the electric fields and electrostatic potentials generated by the fitted densities. GEM-0 has been tested on ten stationary points of the water dimer potential energy surface and on three water clusters (n = 16,20,64). The results show very good agreement with density functional theory calculations, reproducing the individual CSOV energy contributions for a given interaction as well as the B3LYP total interaction energies with errors below kBT at room temperature. Preliminary results for Coulomb and exchange-repulsion energies of metal cation complexes and coupled cluster singles doubles electron densities are discussed.     

Biotransformation,spectroscopic investigation,crystal structure and electrostatic properties of 3,7α‐dihydroxyestra‐1,3,5(10)‐trien‐17‐one monohydrate studied using transferred electron‐density parameters

The biologically transformed product of estradiol valerate, namely 3,7α‐dihydroxyestra‐1,3,5(10)‐trien‐17‐one monohydrate, C₁₈H₂₂O₃·H₂O, has been investigated using UV–Vis, IR, ¹H and ¹³C NMR spectroscopic techniques, as well as by mass spectrometric analysis. Its crystal structure was determined using single‐crystal X‐ray diffraction based on data collected at 100 K. The structure was refined using the independent atom model (IAM) and the transferred electron‐density parameters from the ELMAM2 database. The structure is stabilized by a network of hydrogen bonds and van der Waals interactions. The topology of the hydrogen bonds has been analyzed by the Bader theory of `Atoms in Molecules' framework. The molecular electrostatic potential for the transferred multipolar atom model reveals an asymmetric character of the charge distribution across the molecule due to a substantial charge delocalization within the molecule. The molecular dipole moment was also calculated, which shows that the molecule has a strongly polar character.     

Experimental electron density study of a complex between copper(II) and the antibacterial quinolone family member ciprofloxacin

An experimental charge density study of a 1 : 1 complex of Cu-cfx (cfx = ciprofloxacin), 1 [Cu(cfx)(H(2)O)(3)]SO4.2H(2)O, has been performed using single-crystal X-ray diffraction data collected at 100 K using conventional Mo Kalpha radiation. Metal-ligand (ML) bonds and hydrogen bonds (HBs) have been analysed using topological analysis of the electron density with the atoms in molecules (AIM) approach. The copper atom binds to two oxygen atoms in one end of the zwitterionic form of the cfx molecule, in addition to forming bonds with three water molecules, forming a square pyramidal coordination geometry. AIM decomposition of the experimental electron density establishes that the copper atom binds more strongly to the cfx molecule than to the water molecules, suggesting that the latter can be detached leaving behind a reactive, water-free Cu-cfx complex available for interaction with e.g. a macromolecular site. AIM analysis of the extensive hydrogen bond pattern reveals that the positively charged N-end of the zwitterionic cfx forms a relatively strong N-H-O hydrogen bond implying that this region of cfx may play an important role in the docking process in the active site. Visualisation and statistics of selected density derived properties on the molecular surface of the isolated cfx molecule vs its metal complexed counterpart points out regions of potential reactivity. The effect of the fluorine atom is to expand the negative region of the electrostatic potential, while the nitrogen end is heavily electropositive and willingly donates to--for molecular docking purposes--relatively strong hydrogen bonding. The Cu atom is highlighted as a potentially highly reactive site which is likely to interact strongly with any given negative ligand.     

Theoretical Study on Hydrogen Bonding of Mono‐ and Dihydrated Complexes of 7‐(3′‐Pyridyl)indole in Excited States

The intermolecular hydrogen bonds of mono‐ and dihydrated complexes of 7‐(3′‐Pyridyl)indole (7‐3′PI) have been investigated using the time‐dependent density functional theory (TD‐DFT) method. The electrostatic potential analysis of monomer 7‐3′PI and 7‐(3′‐Pyridyl)indole‐water (7‐3′PI‐W) indicates that an intermolecular hydrogen bond between two waters can be formed for 7‐(3′‐Pyridyl)indole‐2water (7‐3′PI‐2W) complex. The calculated bond lengths of the intermolecular hydrogen bonds of 7‐3′PI‐W and 7‐3′PI‐2W in the S₁ state (the first excited singlet state) are all shortened compared to the ground state. By the analysis of bond length, charge population and infrared spectra, it is demonstrated that the intermolecular hydrogen bonds of 7‐3′PI‐W and 7‐3′PI‐2W are all strengthened upon electronic excitation to the S1 state. Moreover, the fluorescence of 7‐3′PI‐W and 7‐3′PI‐2W are all red‐shifted to larger wavelength compared to monomer 7‐3′PI. The red‐shift of fluorescence peak of 7‐3′PI‐W and 7‐3′PI‐2W should be attributed to the change of hydrogen bond interaction before and after photoexcitation. Therefore, it can be concluded that the intermolecular hydrogen bonding strengthening in the excited S₁ state induces the fluorescence weakening of 7‐3′PI.     

Electrostatic Energies in the 1,4‐Dichlorobenzene Polymorph Crystals: the Role of Charge Density Overlap Effects in Crystal Packing Analysis

Electrostatic and polarization energies for the three known polymorphic crystal structures of 1,4‐dichlorobenzene, as well as for one particularly stable virtual crystal structure generated by computer search, were calculated by a new accurate numerical integration method over static molecular charge densities obtained from high level ab initio molecular‐orbital calculations. Results are compared with those from standard empirical atom‐atom force fields. The new electrostatic energies, which include charge density overlap (penetration) effects, are seen to be much larger than and sometimes of opposite sign to those derived from point‐charge models. None of the four polymorphs is substantially more stable than the others on electrostatic‐energy grounds. Molecule‐molecule electrostatic energies have been calculated for the more important molecular pairs in each of the four structures; trends are found to be very different from those indicated by point‐charge energies or by total energies estimated with a parametric atom‐atom force field. Conclusions based exclusively on analysis of intermolecular atom contacts and point‐charge electrostatics may need to be modified in the light of the new kind of calculation.     

Spectroscopic, electronic structure and natural bond orbital analysis of o-fluoronitrobenzene and p-fluoronitrobenzene: a comparative study

Experimental FTIR, FT-Raman and FT-NMR spectroscopic studies of o-fluoronitrobenzene and p-fluoronitrobenzene have been carried out. A detailed quantum chemical calculations have been performed using DFT/B3LYP method with 6-311++G** and 6-31G** basis sets. Complete vibrational analyses of the compounds were performed. The temperature dependence of thermodynamic properties has been analysed. The atomic charges, electronic exchange interaction and charge delocalisation of the molecule have been performed by natural bond orbital (NBO) analysis. Molecular electrostatic surface potential (MESP), total electron density distribution and frontier molecular orbitals (FMOs) are constructed at B3LYP/6-311++G** level to understand the electronic properties. The charge density distribution and site of chemical reactivity of the molecules have been obtained by mapping electron density isosurface with electrostatic potential surfaces (ESP). The electronic properties, HOMO and LUMO energies were measured by time-dependent TD-DFT approach. (1)H and (13)C NMR spectra were recorded and (1)H and (13)C nuclear magnetic resonance chemical shifts of the molecule were calculated. The (1)H and (13)C nuclear magnetic resonance (NMR) chemical shifts of the molecules in chloroform solvent and in gas phase were calculated by using the Gauge-Independent Atomic Orbital (GIAO) method and are found to be in good agreement with experimental values. The theoretical parameters obtained at B3LYP levels have been compared with the experimental values.     

Topological analysis of the electronic charge density in the ethene protonation reaction catalyzed by acidic zeolite

In the present work, the distribution of the electronic charge density in the ethene protonation reaction by a zeolite acid site is studied within the framework of the density functional theory and the atoms in molecules (AIM) theory. The key electronic effects such as topological distribution of the charge density involved in the reaction are presented and discussed. The results are obtained at B3LYP/6-31G(**) level theory. Attention is focused on topological parameters such as electron density, its Laplacian, kinetic energy density, potential energy density, and electronic energy density at the bond critical points (BCP) in all bonds involved in the interaction zone, in the reactants, pi-complex, transition state, and alkoxy product. In addition, the topological atomic properties are determined on the selected atoms in the course of the reaction (average electron population, N(Omega), atomic net charge, q(Omega), atomic energy, E(Omega), atomic volume, v(Omega), and first moment of the atomic charge distribution, M(Omega)) and their changes are analyzed exhaustively. The topological study clearly shows that the ethene interaction with the acid site of the zeolite cluster, T5-OH, in the ethene adsorbed, is dominated by a strong O-H...pi interaction with some degree of covalence. AIM analysis based on DFT calculation for the transition state (TS) shows that the hydrogen atom from the acid site in the zeolitic fragment is connected to the carbon atom by a covalent bond with some contribution of electrostatic interaction and to the oxygen atom by closed shell interaction with some contribution of covalent character. The C-O bond formed in the alkoxy product can be defined as a weaker shared interaction. Our results show that in the transition state, the dominant interactions are partially electrostatic and partially covalent in nature, in which the covalent contribution increases as the concentration and accumulation of the charge density along the bond path between the nuclei linked increases.     

How the oxazole fragment influences the conformation of the tetraoxazocane ring in a cyclohexanespiro‐3′‐(1,2,4,5,7‐tetraoxazocane): single‐crystal X‐ray and theoretical study

Single crystals of (2S,5R)‐2‐isopropyl‐5‐methyl‐7‐(5‐methylisoxazol‐3‐yl)cyclohexanespiro‐3′‐(1,2,4,5,7‐tetraoxazocane), C₁₆H₂₆N₂O₅, have been studied via X‐ray diffraction. The tetraoxazocane ring adopts a boat–chair conformation in the crystalline state, which is due to intramolecular interactions. Conformational analysis of the tetraoxazocane fragment performed at the B3LYP/6‐31G(d,2p) level of theory showed that there are three minima on the potential energy surface, one of which corresponds to the conformation realized in the solid state, but not to a global minimum. Analysis of the geometry and the topological parameters of the electron density at the (3,?1) bond critical points (BCPs), and the charge transfer in the tetraoxazocane ring indicated that there are stereoelectronic effects in the O—C—O and N—C—O fragments. There is a two‐cross hyperconjugation in the N—C—O fragment between the lone electron pair of the N atom (lpN) and the antibonding orbital of a C—O bond (σ*_C—O) and vice versa between lpO and σ*_C—N. The oxazole substituent has a considerable effect on the geometry and the topological parameters of the electron density at the (3,?1) BCPs of the tetraoxazocane ring. The crystal structure is stabilized via intermolecular C—H…N and C—H…O hydrogen bonds, which is unambiguously confirmed with PIXEL calculations, a quantum theory of atoms in molecules (QTAIM) topological analysis of the electron density at the (3,?1) BCPs and a Hirshfeld analysis of the electrostatic potential. The molecules form zigzag chains in the crystal due to intermolecular C—H…N interactions being electrostatic in origin. The molecules are further stacked due to C—H…O hydrogen bonds. The dispersion component in the total stabilization energy of the crystal lattice is 68.09%.     

Host Perturbation in a β‐Hydroquinone Clathrate Studied by Combined X‐ray/Neutron Charge‐Density Analysis: Implications for Molecular Inclusion in Supramolecular Entities

X‐ray/neutron (X/N) diffraction data measured at very low temperature (15 K) in conjunction with ab initio theoretical calculations were used to model the crystal charge density (CD) of the host–guest complex of hydroquinone (HQ) and acetonitrile. Due to pseudosymmetry, information about the ordering of the acetonitrile molecules within the HQ cavities is present only in almost extinct, very weak diffraction data, which cannot be measured with sufficient accuracy even by using the brightest X‐ray and neutron sources available, and the CD model of the guest molecule was ultimately based on theoretical calculations. On the other hand, the CD of the HQ host structure is well determined by the experimental data. The neutron diffraction data provide hydrogen anisotropic thermal parameters and positions, which are important to obtain a reliable CD for this light‐atom‐only crystal. Atomic displacement parameters obtained independently from the X‐ray and neutron diffraction data show excellent agreement with a |ΔU| value of 0.00058 Ų indicating outstanding data quality. The CD and especially the derived electrostatic properties clearly reveal increased polarization of the HQ molecules in the host–guest complex compared with the HQ molecules in the empty HQ apohost crystal structure. It was found that the origin of the increased polarization is inclusion of the acetonitrile molecule, whereas the change in geometry of the HQ host structure following inclusion of the guest has very little effect on the electrostatic potential. The fact that guest inclusion has a profound effect on the electrostatic potential suggests that nonpolarizable force fields may be unsuitable for molecular dynamics simulations of host–guest interaction (e.g., in protein–drug complexes), at least for polar molecules.     

The investigation of proton transfer and fluorescence‐sensing mechanisms of [2‐(2‐hydroxy‐phenyl)‐1H‐benzoimidazol‐5‐yl]‐phenyl‐methanone

Given the tremendous potential of fluorescence sensors in recent years, in this present work, we theoretically explore a novel fluorescence chemosensor [2‐(2‐Hydroxy‐phenyl)‐1H‐benzoimidazol‐5‐yl]‐phenyl‐methanone (HBPM) about its excited state behaviors and probe‐response mechanism. Using density functional theory (DFT) and time‐dependent density functional theory (TDDFT) methods, we explore the S₀‐state and S₁‐state hydrogen bond dynamical behaviors and confirm that the strengthening intramolecular hydrogen bond in the S₁ state may promote the excited state intramolecular proton transfer (ESIPT) reaction. In view of the photoexcitation process, we find that the charge redistribution around the hydroxyl moiety plays important roles in providing driving force for ESIPT. And the constructed potential energy curves further verify that the ESIPT process of HBPM should be ultrafast. That is the reason why the normal HBPM fluorescence cannot be detected in previous experiment. Furthermore, with the addition of fluoride anions, the exothermal deprotonation process occurs spontaneously along with the intermolecular hydrogen bond O–H?F. It reveals the uniqueness of detecting fluoride anions using HBPM molecules. As a whole, the fluoride anions inhibit the initial ESIPT process of HBPM, which results in different fluorescence behaviors. This work presents the clear ESIPT process and fluoride anion‐sensing mechanism of a novel HBPM chemosensor.     

Two crystal forms of mesogenic bis(4′‐cyanobiphenyl‐4‐yl) butanedioate

The title compound, C₃₀H₂₀N₂O₄, exhibits a nematic phase in the wide temperature range between 498.5 and 538.6 K, in spite of the short linker moiety. Two crystal forms have been found. In both forms, the molecule is centrosymmetric. Form I has a planar biphenyl group, while form II has a twisted biphenyl group with a twist angle of 34.75 (6)°. The packing modes are also different. In form I the long molecular axes are tilted with respect to each other at about 30°, while in form II the long molecular axes have an almost parallel arrangement.     

5‐(3,4‐Dimethoxybenzyl)‐7‐isopropyl‐1,3,5‐triazepane‐2,6‐dione acetonitrile solvate refined using a multipolar atom model

The crystal structure of the title compound, C₁₆H₂₃N₃O₄·CH₃CN, was refined using a multipolar atom model transferred from an experimental electron‐density database. The refinement showed some improvement in crystallographic statistical indices compared with the independent atom model. The triazepane ring adopts a twist‐boat conformation. In the crystal structure, the molecule forms intermolecular contacts with 14 different neighbours. There are two N—H...O and one C—H...O intermolecular hydrogen bond.