Modelling of interactions between volatile anaesthetics (halothane,enflurane) and aromatic compounds,ab initio study

For many years halothane and enflurane have been used clinically as volatile anaesthetics, however, their mechanism of action is still not fully understood. Recently, it has been suggested that they can act by a direct bonding to neuroreceptors containing the aromatic groups. In this work, the halothane?benzene and enflurane?benzene complexes were studied by the ab initio MP2 and CCSD(T) methods. All possible structures of the complexes were calculated by means of the counterpoise CP-corrected gradient optimization technique. It has been found that among these species, the C–H?π hydrogen bonded complexes are the most stable. The CCSD(T)/CBS calculated stabilization energies for halothane and enflurane complexes are: −10.56 and −9.72 kcal mol⁻¹, respectively. The interaction energy is mainly dominated by the dispersion attraction. In the case of enflurane, the C–H bond shows a very small contraction (by −0.0008 Å) upon complexation. This change is accompanied by the blue-shift (20 cm⁻¹) of the C–H stretching frequency and an increase of the infrared intensity of the corresponding mode by 7 km mol⁻¹. Similar results were obtained for the halothane complex: a small contraction of the C–H bond; an increase of the C–H stretching frequency by 11 cm⁻¹ (blue-shift); and an increase of the infrared intensity by 37 km mol⁻¹. In order to explain the nature of these effects, the halothane and enflurane molecules were studied in the electric field generated by benzene atoms, and Natural Bond Orbital (NBO) analyses were performed. The molecular dipole moments of these molecules were calculated with respect to the C–H bond changes. The positive dipole moment derivative obtained for halothane is in agreement with the literature data, while, in the case of enflurane, an unusual effect is observed, the blue-shift of the C–H stretching frequency is accompanied by the positive dipole moment derivative for one C–H bond and the negative for the other C–H bond. The mechanisms responsible for contraction and strengthening of the C–H bonds are discussed.     

Computational study of structures and proton transfer in hydrogen‐bonded ammonia complexes using semiempirical valence‐bond approach

In this study, the seGVB method was implemented for the N H bonding system, specifically for hydrogen‐bonded ammonia complexes, and the model well reproduces the MP2 geometries and energetics. A comparison between the ammonia dimer and water dimer is given from the viewpoint of valance‐bond structures in terms of the calculated bond energies and pair–pair interactions. The linear hydrogen bond is found to be stronger than the bent bonds in both cases, with the difference in energy between the linear and cyclic structures being comparable in both cases although the NH bonds are generally weaker. The energy decomposition clearly demonstrates that the changes in electronic energy are quite different in the two cases due to the presence of an additional lone pair on the water molecule, and it is this effect which leads to the net stabilization of the cyclic structure for the ammonia dimer. Proton‐transfer profiles for hydrogen‐bonded ammonia complexes [NH₂ H NH₂]⁻ and [NH₃ H NH₃]⁺ were calculated. The barrier for proton transfer in [NH₃ H NH₃]⁺ is larger than that in [NH₂ H NH₂]⁻, but smaller than that in the protonated water dimer. The different bonding structures substantially affect the barrier to proton transfer, even though they are isoelectronic systems. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 73: 357–367, 1999     

A comparative molecular field analysis (CoMFA) study using semiempirical, density functional, ab initio methods and pharmacophore derivation using DISCOtech on sigma 1 ligands

The Comparative Molecular Field Analysis (CoMFA) was developed to investigate a three-dimensional quantitative structure activity relationship (3D-QSAR) model of ligands for the sigma 1 receptor. The starting geometry of sigma-1 receptor ligands was obtained from the Tripos force field minimizations and conformations were decided from DISCOtech using the SYBYL 6.8. program. The structures of 48 molecules were fully optimized at the ab initio HF/3-21G* and semiempirical AM1 calculations using GAUSSIAN 98. The electrostatic charges were calculated using several methods such as semiempirical AM1, density functional B3LYP/3-21G*, and ab initio HF/3-21G*, MP2/3-21G* calculations within GAUSSIAN 98. Using the optimized geometries, the CoMFA results derived from the HF/3-21G method were better than those from AM1. The best CoMFA was obtained from HF/3-21G* optimized geometry and charges (R2 = 0.977). Using the optimized geometries, the CoMFA results derived from the HF/3-21G methods were better than those from AM1 calculations. The training set of 43 molecules gave higher R2 (0.989-0.977) from HF/3-21G* optimized geometries than R2 (0.966-0.911) values from AM1 optimized geometries. The test set of five molecules also suggested that HF/3-21G* optimized geometries produced good CoMFA models to predict bioactivity of sigma 1 receptor ligands but AM1 optimized geometries failed to predict reasonable bioactivity of sigma 1 receptor ligands using different calculations for atomic charges.     

Quantum chemical study of the bond orders in the ruthenium,diruthenium and dirhodium nitrosyl complexes

Density functional calculations with the B3LYP functional were carried out for the [Ru(NO)Cl₅]²⁻, [Ru(NO)(NH₃)₅]³⁺, [Ru(NO)(CN)₅]²⁻, [Ru(NO)(CN)₅]³⁻, [Ru(NO)(hedta)]^q (hedta = N-(hydroxyethyl)ethylenediaminetriacetate triple-charged anion; q = 0, −1, −2), Rh₂(O₂CR)₄, Rh₂(O₂CR)₄(NO)₂, Ru₂(O₂CR)₄, Ru₂(O₂CR)₄(NO)₂, Ru₂(dpf)₄, and Ru₂(dpf)₄(NO)₂ (dpf = N,N′-diphenylformamidinate ion; R = H, CH₃, CF₃) complexes. The electronic structure was analyzed in terms of Mayer and Wiberg bond order indices. The technique of bond order indices decomposition into σ-, π-, and δ-contributions was proposed.     

Theoretical study of the strong intramolecular hydrogen bond and metal-ligand interactions in group 10 (Ni, Pd, Pt) bis(dimethylglyoximato) complexes

The structural and bonding characteristics of the bis(dimethylglyoximato) complexes of group 10 transition metals ([M(dmg)₂], where M = Ni, Pd and Pt) were investigated by means of quantum chemical computations. The equilibrium geometries, energetic and bonding properties were computed using the B3P86 exchange-correlation density functional in conjunction with a 6-311+(+)G^∗∗ basis set. The computations revealed that the strong O⁻?H-O hydrogen bond exists only in the presence of the metal cations. The free (dmg)₂²⁻ ligand has significantly different geometry in which the O⁻?H-O interaction is replaced by N?O-H bonds. The characteristics of the metal-ligand interactions were determined by natural bond orbital analysis.     

Interactions of the “piano‐stool” [ruthenium(II) (η6‐arene)(en)CL]+ complexes with water and nucleobases; ab initio and DFT study

Piano stool ruthenium complexes of the composition [Ru(II)(η⁶‐arene)(en)Cl]^+/2+ (en = ethylenediamine) represent an emerging class of cisplatin‐analogue anticancer drug candidates. In this study, we use computational quantum chemistry to characterize the structure, stability and reactivity of these compounds. All these structures were optimized at DFT(B3LYP)/6‐31G(d) level and their single point properties were determined by the MP2/6‐31++G(2df,2pd) method. Thermodynamic parameters and rate constants were determined for the aquation process, as a replacement of the initial chloro ligand by water and subsequent exchange reaction of aqua ligand by nucleobases. The computations were carried out at several levels of DFT and ab initio theories (B3LYP, MP2 and CCSD) utilizing a range of bases sets (from 6‐31G(d) to aug‐cc‐pVQZ). Excellent agreement with experimental results for aquation process was obtained at the CCSD level and reasonable match was achieved also with the B3LYP/6‐31++G(2df,2pd) method. This level was used also for nucleobase‐water exchange reaction where a smaller rate constant for guanine exchange was found in comparison with adenine. Although adenine follows a simple replacement mechanism, guanine complex passes by a two‐step mechanism. At first, Ru‐O6(G) adduct is formed, which is transformed through a chelate TS2 to the Ru‐N7(G) final complex. In case of guanine, the exchange reaction is more favorable thermodynamically (releasing in total by about 8 kcal/mol) but according to our results, the rate constant for guanine substitution is slightly smaller than the analogous constant in adenine case when reaction course from local minimum is considered. © 2008 Wiley Periodicals, Inc. J Comput Chem, 2009     

A mass spectrometry and DFT study of pyrithione complexes with transition metals in the gas phase

2‐Mercaptopyridine N ‐oxide (pyrithione, PTOH) along with several transition metal ions forms coordination compounds displaying notable biological activities. Gas‐phase complexes formed between pyrithione and manganese (II), cobalt (II), nickel (II), copper (II), and zinc (II) were investigated by infusion in the electrospray source of a quadrupole‐time of flight mass spectrometer. Remarkably, positive ion mode spectra displayed the singly charged metal adduct ion [C₁₀H₈MN₂O₂S₂]²⁺ ([M(PTO)₂]^+• or [M(DPTO)]^+•), where DPTO is dipyrithione, 2,2′‐dithiobis(pyridine N ‐oxide), among the most abundant peaks, implying a change in the oxidation state of whether the metal ion or the ligands. In addition, doubly charged ions were recognized as metal adduct ions containing DPTO ligands, [M(DPTO)_n]²⁺. Generation of [M(PTO)₂]^+• / [M(DPTO)]^+• could be traced by CID of [M(DPTO)₂]²⁺, by observation of the sequential losses of a charged (PTO⁺) and a radical (PTO^•) deprotonated pyrithione ligand. The fragmentation pathways of [M(PTO)₂]^+• / [M(DPTO)]^+• were compared among the different metal ions, and some common features were noticed. Density functional theory (DFT) calculations were employed to study the structures of the observed adduct ions, and especially, to decide in the adduct ion [M(PTO)₂]^+• / [M(DPTO)]^+• whether the ligands are 2 deprotonated pyrithiones or a single dipyrithione as well as the oxidation state of the metal ion in the complex. Characterization of gas‐phase pyrithione metal ion complexes becomes important, especially taking into account the presence of a redox‐active ligand in the complexes, because redox state changes that produce new species can have a marked effect on the overall toxicological/biological response elicited by the metal system.     

A theoretical analysis of the explanation of the significant differences in antiferromagnetic interactions between homologous μ-alkoxo and acetate bridged dicopper(II) complexes: ab initio and semi-empirical calculations

A magnetostructural classification of dimmers, containing the Cu (μ-alkoxo) Cu core, based on data obtained from X-ray diffraction analysis reported in the literature has been performed. In these complexes, the local geometry around the copper ions is generally a square planar and each copper ion is surrounded by one N atom and three O atoms. The influence of the overlap interactions between the bridging ligands and the metal (Cu) d orbitals on the super-exchange coupling constant has been studied by means of ab initio Restricted Hatree–Fock molecular orbital calculations. The interaction between the magnetic d orbitals and highest occupied molecular orbitals of the acetate oxygens has been investigated in homologous μ-acetato-bridged dicopper(II) complexes which have significantly different −2J values (the energy separation between the spin-triplet and spin-singlet states). In order to determine the nature of the fronter orbitals, Extended Hückel molecular Orbital calculations are also reported. Ab initio restricted Hartree–Fock calculations have shown that the acetato bridge and the alkoxide bridge contribute to the magnetic interaction countercomplementarily to reduce antiferromagnetic interaction.     

Hydrogen bond interactions in sulfamerazine: DFT study of the O-17, N-14, and H-2 electric field gradient tensors

Hydrogen bond (HB) interactions are studied in the real crystalline structure of sulfamerazine by density functional theory (DFT) calculations of the electric field gradient (EFG) tensors at the sites of O-17, N-14, and H-2 nuclei. One-molecule (single) and four-molecule (cluster) models of sulfamerazine are created by available crystal coordinates and the EFG tensors are calculated in both models to indicate the influence of HB interactions on the tensors. Directly relate to the experiments, the calculated EFG tensors are converted to the experimentally measurable nuclear quadrupole resonance (NQR) parameters, quadrupole coupling constant (qcc) and asymmetry parameter (η_Q). The evaluated NQR parameters reveal that due to contribution of the target molecule to N–HN and N–HO types of HB interactions, the EFG tensors at the sites of various nuclei are influenced from single model to the target molecule in cluster. Additionally, O2, N4, and H2 nuclei of the target molecule are significantly influenced by HB interactions, consequently, they have the major contributions to HB interactions in cluster model of sulfamerazine. The calculations are performed employing B3LYP method and 6-311++G** basis set using GAUSSIAN 98 suite of program.     

DFT:B3LYP ab initio molecular dynamics study of the Zundel and Eigen proton complexes, H5O2+ and H9O4+, in the triplet state in gas phase and solution

DFT:B3LYP ab initio molecular dynamics (MD) approach is used to elucidate the properties of the Zundel and Eigen, H5O2+ and H9O4+, proton complexes in the triplet state. The simulation considers the complexes in the gas phase (isolated complexes) and inside the clusters composed of 32, 64, and 128 water molecules, mimicking the behavior of aqueous solutions. MD simulations reveal three distinct periods. For the complex in solutions, the periods are smoothed out. The H5O2+ and H9O4+ complexes in the triplet state undergo structural rearrangements, which eventually result in hydrogen elimination. For the H5O2+, the hydrogen is eliminated from the center of the water cluster, whereas for the H9O4+ it is removed from a near-surface water molecule. The rate of hydrogen elimination decreases with increasing the number of water molecules surrounding the complex.     

A comparative study by AM1, PM3 and ab initio HF/3-21G methods on the structure and reactivity of monosubstituted carbanion 1,2,4-triazolium ylides

A full conformational analysis of six 1,2,4-monosubstituted carbanion 1,2,4-triazolium ylides 4 a–f was performed using AM1, PM3 and HF/3-21G methods. The C-type conformers were found as the most stable structures by these different methods. This study also includes a qualitative estimation of the chemical behavior of triazolium ylides 4 a–f as nucleophilic agents on the level of ylide carbon atoms. The ab initio 3-21G method seems to be the most suitable in the characterization of these molecular systems.     

A comparative study of proton conduction between two new Cd(II) and Co(II) complexes and in vitro antibacterial study of the Cd(II) complex

Two new complexes based on 4,4′-[1,3-phenylenebis(oxy)]diphthalic acid (H₄L) ligands were synthesized, namely, [Cd₂(L)(1,10-phen)₂(H₂O)]_n( 1 ) and [Co₂(L)(1,10-phen)₂(H₂O)]_n( 2 ), in which 2D structures transform into 3D supramolecular structures by C H···π interaction. The proton conductivity of complexes 1 and 2 at low temperature is close (σ₁ = 3.12 × 10⁻⁸ S cm⁻¹ and σ₂ = 3.81 × 10⁻⁸ S cm⁻¹ at 30°C), but these two complexes show different conduction mechanisms. The Vehicular mechanism in 1 is caused by the O···H/H···O contact in 1 , which is stronger than 2 , and the Grotthuss mechanism in 2 is caused by the N···H/H···N contact in 2 , which is stronger than 1 . At the same time, complex 1 showed excellent antibacterial properties in vitro, mainly reflected in that five kinds of bacteria (Escherichia coli, Bacillus amyloliquefaciens, Pantoea agglomerans, Pseudomonas putida, and Pectobacterium carotovora) could play an obvious inhibitory effect in the concentration range of 20 μg·ml⁻¹.     

(N,N) vs. (N,S) chelation of palladium in asymmetric allylic substitution using bis(thiazoline) ligands: A theoretical and experimental study

New thiazoline-containing ligands including non-symmetric bis(thiazolines) and oxazoline-thiazolines were synthesized and then compared to C₂-symmetric bis(thiazolines) in the palladium-catalyzed allylic substitution. The experimental results obtained in this study support the hypothesis of a competition between the (N,N) and the (N,S) palladium chelation, when sterically hindered bis(thiazolines) are used as ligands. A quantum chemical study performed on the Pd-complexes derived from three selected ligands, two C₂-symmetric bis(thiazolines) and one oxazoline-thiazoline, also supports this hypothesis.     

Analysis of pair interaction in a bipolar mesogen 4‐(4‐hydroxylbutyloxy)‐4′‐cyano‐biphenyl: A comparative study based on semiempirical and DFT methods

Structure of 4‐(4‐hydroxylbutyloxy)‐4′‐cyano‐biphenyl (H4CBP) molecule has been optimized using density functional B3LYP with 6‐31G (d) basis set taking crystallographic geometry as input. Using the optimized geometry, electronic structure of the H4CBP molecule has been evaluated on the basis of semiempirical methods and DFT calculations. Intermolecular interaction energy between a pair of H4CBP molecules has been evaluated by using Rayleigh–Schrodinger perturbation theory modified with multicentered multipole expansion method for the electrostatic part while dispersion and repulsion terms have been calculated using Kitaigorodskii formula. The results obtained through semiempirical and DFT calculations have been compared for various interacting conditions, viz.: (a) stacking, (b) in‐plane, and (c) terminal interactions. A comparative analysis of the results has been carried out with a view to examine suitability of different methods to study molecular aggregations in moderately large organic systems. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Synthesis,structure and in vitro antiproliferative activities of oxamido‐bridged dicopper(II) complexes: A comparative study of experimental evidence and molecular docking of DNA/protein binding

Two μ‐oxamido‐bridged dicopper(II) complexes, namely [Cu₂(hmpoxd)(H₂O)(phen)](ClO₄) ( 1 ) and [Cu₂(papo)(H₂O)(phen)](ClO₄)·2H₂O ( 2 ), where H₃hmpoxd and H₃papo represent N‐(2‐hydroxy‐5‐methylphenyl)‐N′‐[3‐(dimethylamino)propyl]oxamide and N‐(2‐hydroxylphenyl)‐N′‐(3‐aminopropyl)oxamide, respectively, and phen represents 1,10‐phenanthroline, were synthesized. Single‐crystal X‐ray crystallography and other methods revealed that the two copper(II) ions in complex 1 are bridged by the cis‐hmpoxd^3? with Cu···Cu separation of 5.1896(7) Å, in which the inner (Cu1) and outer (Cu2) copper(II) atoms are located in square‐planar and square‐pyramidal geometries, respectively. To evaluate the effects of bridging ligand hydrophobicity on DNA/protein binding and potential anticancer activities, comparative studies of the reactivity towards herring sperm DNA and protein bovine serum albumin (BSA) as well as cytotoxicity of complex 1 with our previously reported complex 2 were conducted theoretically and experimentally. The results indicate that the two complexes can interact interactively with DNA, and bind to BSA via the binding sites Trp213 for 1 and Trp134 for 2 . Interestingly, the in vitro anticancer activities and DNA/protein binding affinities consistently follow the order of 1 > 2 .     

Theoretical investigation of C–HM interactions in organometallic complexes: A natural bond orbital (NBO) study

The nature of C–HM agostic interactions in model metal complexes [M²⁺(CH₂CH₃)(PH₃)_nCl] (where M = Sc, Ti, V, Ti, Cr, Mn, Fe, Co, Ni, Cu, Zn; n = 1, 2, 3, 4) was studied with the natural bond orbital analysis (NBO) approach using density functional theory (DFT) optimized geometries at the B3LYP/6-31G(d,p) level of theory. The effect of nature of metal, coordination number, oxidation state and ligand field effects on the agostic interaction is examined. A set of 20 crystal structures of organometallic complexes taken from the Cambridge Structural Database (CSD) was studied computationally employing AIM theory and NBO analysis, and the applicability of these methods was critically accessed in demarcating the two types of interaction.