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.     

Pyrene: hydrogenation, hydrogen evolution, and π-band model

We present a theoretical investigation of the hydrogenation of pyrene and of the subsequent molecular hydrogen evolution. Using density functional theory (DFT) at the GGA-PBE level, the chemical binding of atomic hydrogen to pyrene is found to be exothermic by up to 1.6 eV with a strong site dependence. The edge C atoms are found most reactive. The barrier for the formation of the hydrogen-pyrene bond is small, down to 0.06 eV. A second hydrogen binds barrierless at many sites. The most stable structure of dihydrogenpyrene is more stable by 0.64 eV than pyrene plus a molecular hydrogen molecule and a large barrier of 3.7 eV for the molecular hydrogen evolution is found. Using a simple tight-binding model we demonstrate that the projected density of π-states can be used to predict the most stable binding sites for hydrogen atoms and the model is used to investigate the most favorable binding sites on more hydrogenated pyrene molecules and on coronene. Some of the DFT calculations were complemented with hybrid-DFT (PBE0) showing a general agreement between the DFT and hybrid-DFT results.     

A study of the phase behavior of a simple model of chiral molecules and enantiomeric mixtures

We present a study of the solid-fluid and solid-solid phase equilibrium for molecular models representative of chiral molecules and enantiomeric mixtures. The models consist of four hard sphere interaction sites of different diameters in a tetrahedral arrangement with the fifth hard sphere interaction site at the center of the tetrahedron. The volumetric properties and free energies of the pure enantiomers and binary mixtures were calculated in both fluid and solid phases using isobaric Monte Carlo simulations. The models exhibit essentially ideal solution behavior in the fluid phase with little chiral discrimination. In the solid phase the effects of chirality are much greater. Solid-fluid phase behavior involving the pure enantiomer solids and also racemic compounds was calculated. The calculations indicate that, depending on the relative sizes of the hard sphere interaction sites, packing effects alone can be sufficient to stabilize a racemic compound with respect to the pure enantiomer solids.     

Direct correlation functions of binary mixtures of hard Gaussian overlap molecules

We study the direct correlation function (DCF) of a classical fluid mixture of nonspherical molecules. The components of the mixture are two types of hard ellipsoidal molecules with different elongations, interacting through the hard Gaussian overlap (HGO) model. Two different approaches are used to calculate the DCFs of this fluid, and the results are compared. Here, the Pynn approximation [J. Chem. Phys. 60, 4579 (1974)] is extended to calculate the DCF of the binary mixtures of HGO molecules, then we use a formalism based on the weighted density functional theory introduced by Chamoux and Perera [J. Chem. Phys. 104, 1493 (1996)]. These results are fairly in agreement with each other. The pressure of this system is also calculated using the Fourier zero components of the DCF. The results are in agreement with the Monte Carlo molecular simulation.     

A solvation‐free‐energy functional: A reference‐modified density functional formulation

The three‐dimensional reference interaction site model (3D‐RISM) theory, which is one of the most applicable integral equation theories for molecular liquids, overestimates the absolute values of solvation‐free‐energy (SFE) for large solute molecules in water. To improve the free‐energy density functional for the SFE of solute molecules, we propose a reference‐modified density functional theory (RMDFT) that is a general theoretical approach to construct the free‐energy density functional systematically. In the RMDFT formulation, hard‐sphere (HS) fluids are introduced as the reference system instead of an ideal polyatomic molecular gas, which has been regarded as the appropriate reference system of the interaction‐site‐model density functional theory for polyatomic molecular fluids. We show that using RMDFT with a reference HS system can significantly improve the absolute values of the SFE for a set of neutral amino acid side‐chain analogues as well as for 504 small organic molecules. © 2015 Wiley Periodicals, Inc.     

Adlayer structures of aza- and/or oxo-bridged calix[2]arene[2]triazines on Au(111) investigated by scanning tunneling microscopy (STM)

The adlayers formed by a series of aza- and/or oxo-bridged calix[2]arene[2]triazines on Au(111) surfaces were investigated by scanning tunneling microscopy (STM) and density functional theory (DFT) calculations. 1,3-Alternate configurations of these molecules are preserved on gold surfaces as in their three-dimensional crystals. STM images show that the cavity sizes of these molecules are finely tuned by substituting the bridging nitrogen atom with oxygen atoms, which change the strengths and densities of the intermolecular hydrogen bonds. Hydrogen bond interaction influences the molecular orientation and conformation in the adlayers, and it plays a key role in the formation of these two-dimensional supramolecular architectures. Coadsorption of calix[2]arene[2]triazine with 1,3,5-tris(5-carboxyamyloxy)benzene (TCAB) intervenes with the intermolecular hydrogen bond formations among the calix[2]arene[2]triazine molecules and consequently causes a conformational transition of the calixarene molecules from rhombic to square. These results demonstrate the role of hydrogen bonds in molecular assembly formations.     

Free energy calculations using flexible-constrained, hard-constrained and non-constrained molecular dynamics simulations.

A comparison of different treatments of bond-stretching interactions in molecular dynamics simulation is presented. Relative free energies from simulations using rigid bonds maintained with the SHAKE algorithm, using partially rigid bonds maintained with a recently introduced flexible constraints algorithm, and using fully flexible bonds are compared in a multi-configurational thermodynamic integration calculation of changing liquid water into liquid methanol. The formula for the free energy change due to a changing flexible constraint in a flexible constraint simulation is derived. To allow for a more direct comparison between these three methods, three different pairs of models for water and methanol were used: a flexible model (simulated without constraints and with flexible constraints), a rigid model (simulated with standard hard constraints), and an alternative flexible model (simulated with flexible constraints and standard hard constraints) in which the ideal or constrained bond lengths correspond to the average bond lengths obtained from a short simulation of the unconstrained flexible model. The particular treatment of the bonds induces differences of up to 2 % in the liquid densities, whereas (excess) free energy differences of up to 5.7 (4.3) kJ mol(-1) are observed. These values are smaller than the differences observed between the three different pairs of methanol/water models: up to 5 % in density and up to 8.5 kJ mol(-1) in (excess) free energy.     

Chemical reactivity and selectivity using Fukui functions: basis set and population scheme dependence in the framework of B3LYP theory

The use of Fukui functions for the site selectivity of the formaldehyde molecule for nucleophilic, electrophilic and radical attacks has been made with special emphasis to the dependence of Fukui values on the basis sets as well as population schemes in the framework of B3LYP theory. Out of the five population schemes selected viz., Mulliken population analysis, natural population analysis, CHELP, CHELPG and atoms in molecules (AIM), it is found that the CHELPG and AIM schemes predict precise reactive site with less dependency on the basis sets. Charges derived from Hirshfeld partitioning, calculated using the BLYP/dnd method (implemented in the DMOL³ package), provide non-negative Fukui values for all the molecular systems considered in this study. Supporting results have been obtained for acetaldehyde and acetone molecules at the 6-31+G** basis set level. These results support the fact that high Fukui values correspond to soft–soft interaction sites. On the other hand, the correlation of the low Fukui value to the hard–hard interaction site merits further investigation. Received: 10 November 2001 / Accepted: 6 March 2002 / Published online: 13 June 2002     

Carbon chains and the (5,5) single-walled nanotube: structure and energetics versus length

Reliable thermochemistry is computed for infinite stretches of pure-carbon materials including acetylenic and cumulenic carbon chains, graphene sheet, and single-walled carbon nanotubes (SWCNTs) by connection to the properties of finite size molecules that grow into the infinitely long systems. Using ab initio G3 theory, the infinite cumulenic chain (:C[double bond]C[double bond]C[double bond]C:) is found to be 1.9+/-0.4 kcal/mol per carbon less stable in free energy at room temperature than the acetylenic chain (.C[triple bond]C-C[triple bond]C.) which is 24.0 kcal/mol less stable than graphite. The difference between carbon-carbon triple, double, and single bond lengths (1.257, 1.279, and 1.333 A, respectively) in infinite chains is evident but much less than with small hydrocarbon molecules. These results are used to evaluate the efficacy of similar calculations with the less rigorous PM3 semiempirical method on the (5,5) SWCNT, which is too large to be studied with high-level ab initio methods. The equilibrium electronic energy change for C(g)-->C[infinite (5,5) SWCNT] is -166.7 kcal/mol, while the corresponding free energy change at room temperature is -153.3 kcal/mol (6.7 kcal/mol less stable than graphite). A threefold alternation (6.866, 6.866, and 6.823 A) in the ring diameter of the equilibrium structure of infinitely long (5,5) SWCNT is apparent, although the stability of this structure over the constant diameter structure is small compared to the zero point energy of the nanotube. In general, different (n,m) SWCNTs have different infinite tube energetics, as well as very different energetic trends that vary significantly with length, diameter, and capping.     

Hydrogen bond energies and cooperativity in substituted calix[n]arenes (n = 4, 5)

Hydrogen-bonded interactions in para-substituted calix[n]arenes (CX[n]) (n = 4, 5) and their thia analogues are analyzed using the recently proposed molecular tailoring approach. The cooperative contribution toward the hydrogen-bonding network within the CX[5] host is shown to be nearly 5 times larger than that in its thia analogue. Hydrogen bond strengths in the O-H···O network are enhanced on substitution of an electron-donating group. The cooperativity contributions are reflected in the electron density at the bond critical point in the quantum theory of atoms in molecules.     

Silicon radicals in silicon oxynitride: a theoretical ESR study

Quantum mechanical calculations are performed on a series of silicon radical defects. These are the upward arrow Si[triple bond]O(3-x)Nx, upward arrow Si[triple bond]N(3-x)Si(x), and upward arrow Si[triple bond]Si(3-x)Ox defects, where x takes on values between 0 and 3. The defects under study constitute a central silicon radical, upward arrow Si, with differing first-nearest-neighbor substitution, as may be found at a Si/SiOxNy interface. These first-nearest neighbor atoms are connected to the silicon radical via three single covalent bonds, denoted as " [triple bond] ". A hybrid defect, upward arrow Si[triple bond]ONSi, is also included. Calculations are performed on gas-phase-like cluster models, as well as more-constrained hybrid quantum and molecular mechanical (QM/MM) models. The isotropic hyperfine coupling constants of these defects are calculated via density functional theory (DFT). Trends in these calculated hyperfines are consistent between the different models utilized. Analysis of the electronic structure and geometries of defects correlate well with trends in the electronegativity of the first-nearest-neighbor atoms. Changes in radical hybridization, induced by changes in the first-nearest-neighbor composition, are the primary factor that affects the calculated hyperfines. Furthermore, comparisons to experimental results are encouraging. Agreement is found between experiments on amorphous to crystalline materials.     

Atomic forces for geometry‐dependent point multipole and Gaussian multipole models

In standard treatments of atomic multipole models, interaction energies, total molecular forces, and total molecular torques are given for multipolar interactions between rigid molecules. However, if the molecules are assumed to be flexible, two additional multipolar atomic forces arise because of (1) the transfer of torque between neighboring atoms and (2) the dependence of multipole moment on internal geometry (bond lengths, bond angles, etc.) for geometry‐dependent multipole models. In this study, atomic force expressions for geometry‐dependent multipoles are presented for use in simulations of flexible molecules. The atomic forces are derived by first proposing a new general expression for Wigner function derivatives . The force equations can be applied to electrostatic models based on atomic point multipoles or Gaussian multipole charge density. Hydrogen‐bonded dimers are used to test the intermolecular electrostatic energies and atomic forces calculated by geometry‐dependent multipoles fit to the ab initio electrostatic potential. The electrostatic energies and forces are compared with their reference ab initio values. It is shown that both static and geometry‐dependent multipole models are able to reproduce total molecular forces and torques with respect to ab initio, whereas geometry‐dependent multipoles are needed to reproduce ab initio atomic forces. The expressions for atomic force can be used in simulations of flexible molecules with atomic multipoles. In addition, the results presented in this work should lead to further development of next generation force fields composed of geometry‐dependent multipole models. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Theoretical Study of the Regioselectivity of the Interaction of 3-Methyl-4-pyrimidone and 1-Methyl-2-pyrimidone with Lewis Acids

A density functional theory (DFT) study is performed to determine the stability of the complexes formed between either the N or O site of 3-methyl-4-pyrimidone and 1-methyl-2-pyrimidone molecules and different ligands. The studied ligands are boron and alkali Lewis acids, namely, B(CH(3))(3), HB(CH(3))(2), H(2)B(CH(3)), BH(3), H(2)BF, HBF(2), BF(3), Li(+), Na(+), and K(+). The acids are divided into two groups according to their hardness. The reactivity predictions, according to the molecular electrostatic potential (MEP) map and the natural bond orbital (NBO) analysis, are in agreement with the calculated relative stabilities. Our findings reveal a strong regioselectivity with borane and its derivatives preferring the nitrogen site in both pyrimidone isomers, while a preference for oxygen is observed for the alkali acids in the 3-methyl-4-pyrimidone molecule. The complexation of 1-methyl-2-pyrimidone with these hard alkali acids does not show any discrimination between the two sites due to the presence of a continuous delocalized density region between the nitrogen and the oxygen atoms. The preference of boron Lewis acids toward the N site is due to the stronger B-N bond as compared to the B-O bond. The influence of fluorine or methyl substitution on the boron atom is discussed through natural orbital analysis (NBO) concentrating on the overlap of the boron empty p-orbital with the F lone pairs and methyl hyperconjugation, respectively. The electrophilicity of the boron acids gives a good overall picture of the interaction capabilities with the Lewis base.     

Exact density functional theory for ideal polymer fluids with nearest neighbor bonding constraints

We present a new density functional theory of ideal polymer fluids, assuming nearest-neighbor bonding constraints. The free energy functional is expressed in terms of end site densities of chain segments and thus has a simpler mathematical structure than previously used expressions using multipoint distributions. This work is based on a formalism proposed by Tripathi and Chapman [Phys. Rev. Lett. 94, 087801 (2005)]. Those authors obtain an approximate free energy functional for ideal polymers in terms of monomer site densities. Calculations on both repulsive and attractive surfaces show that their theory is reasonably accurate in some cases, but does differ significantly from the exact result for longer polymers with attractive surfaces. We suggest that segment end site densities, rather than monomer site densities, are the preferred choice of "site functions" for expressing the free energy functional of polymer fluids. We illustrate the application of our theory to derive an expression for the free energy of an ideal fluid of infinitely long polymers.     

Structure of inhomogeneous polymer solutions: a density functional approach

The structure of polymer solutions confined between surfaces is studied using a density functional theory where the polymer molecules have been modeled as a pearl necklace of freely jointed hard spheres and the solvent as hard spheres. The present theory uses the concept of universality of the free energy density functional to obtain the first-order direct correlation function of the nonuniform system from that of the corresponding uniform system, calculated through the Verlet-modified type bridge function. The uniform bulk fluid direct correlation function required as input has been calculated from the reference interaction site model integral equation theory using the Percus-Yevick closure relation. The calculated results on the density profiles of the polymer as well as the solvent are shown to compare well with computer simulation results.     

Molecular model with quantum mechanical bonding information

The molecular structure can be defined quantum mechanically thanks to the theory of atoms in molecules. Here, we report a new molecular model that reflects quantum mechanical properties of the chemical bonds. This graphical representation of molecules is based on the topology of the electron density at the critical points. The eigenvalues of the Hessian are used for depicting the critical points three-dimensionally. The bond path linking two atoms has a thickness that is proportional to the electron density at the bond critical point. The nuclei are represented according to the experimentally determined atomic radii. The resulting molecular structures are similar to the traditional ball and stick ones, with the difference that in this model each object included in the plot provides topological information about the atoms and bonding interactions. As a result, the character and intensity of any given interatomic interaction can be identified by visual inspection, including the noncovalent ones. Because similar bonding interactions have similar plots, this tool permits the visualization of chemical bond transferability, revealing the presence of functional groups in large molecules.     

Characteristics of multiple N,O bonds

The electron distribution for N,O bonds in a wide range of molecules is analyzed using Bader's atoms in molecules theory and measures of bond order. It is shown that electron density derived parameters correlate very well with bond length and bond orders. An unexpectedly large number of bond orders are near to or greater than 2, indicating that nitrogen is effectively pentavalent in a many simple molecules. This conclusion is substantiated by examination of the Laplacian of the electron density. Such hypervalence indicates that nitrogen does not always obey the classical octet rule.     

The elastic constants of nematics

Abstract

It is pointed out that some often cited molecular statistical calculations of the elastic constants of nematics are based on a postulated instead of a derived expression for the distortion free energy density. In particular attention is paid to the contradictory results of Priest and Poniewierski and Stecki for hard rod systems. The appropriate way to calculate the elastic constants of hard rod models from first principles is discussed briefly.     

Reaction mechanism of hydrogen cyanide catalyzed by gas‐phase titanium

To explore the details of the reaction mechanism of Ti atom with HCN, the reactive site and reactivity have been predicted first, the potential energy surfaces have been systematically studied at different theoretical levels. Four different reaction pathways and product distribution are discussed by means of the activation strain model and Curtin–Hammett principle. In addition, the structures, bonding properties and the frontier molecular orbital interaction diagrams of main stationary points were analyzed by atoms in molecules and natural bond orbital. The results show that for this system, there are four reaction pathways, in which path b (HCN＋Ti→IM1→TS1→IM2→T2b→IM4) is the most favorable pathway.     

Theoretical investigation of hyperthermal reactions at the gas-liquid interface: O (3P) and squalane

Hyperthermal collisions (5 eV) of ground-state atomic oxygen [O ((3)P)] with a liquid-saturated hydrocarbon, squalane (C(30)H(62)), have been studied using QM/MM hybrid "on-the-fly" direct dynamics. The surface structure of the liquid squalane is obtained from a classical molecular dynamics simulation using the OPLS-AA force field. The MSINDO semiempirical Hamiltonian is combined with OPLS-AA for the QM/MM calculations. In order to achieve a more consistent and efficient simulation of the collisions, we implemented a dynamic partitioning of the QM and MM atoms in which atoms are assigned to QM or MM regions based on their proximity to "seed" (open-shell) atoms that determine where bond making/breaking can occur. In addition, the number of seed atoms is allowed to increase or decrease as time evolves so that multiple reactive events can be described. The results show that H abstraction is the most important process for all incident angles, with H elimination, double H abstraction, and C-C bond cleavage also being important. A number of properties of these reactive channels, as well as inelastic nonreactive scattering, are investigated, including angular and translational energy distributions, the effect of incident collision angle, variation with depth of the reactive event within the liquid, with the reaction site on the hydrocarbon, and the effect of dynamics before and after reaction (direct reaction versus trapping reaction-desorption).