Benchmark of different charges for prediction of the partitioning coefficient through the hydrophilic/lipophilic index

A few different theoretical methods for assigning the partial atomic charges were benchmarked for calculation of the hydrophilic/lipophilic index (HLI). The coefficients were selected to produce the best correlation of the HLI values with the experimental octanol-water partition. Different parameters were checked in calculations of partial charges to get the best performance of the HLI values obtained. Thus, four partitioning schemes (Coulson, Mulliken, Merz-Kollman, Ford-Wang) were benchmarked for calculations of atomic charges with six semiempirical methods (AM1, PM3, RM1, PM6, PM6-D3H4, PM7). Moreover, five distinct types of partial atomic charges (Mulliken, Hirshfeld, Löwdin, CHELPG, NPA), obtained at the Hartree–Fock and DFT levels of theory with three basis sets, were tested for their ability to produce the HLI values with the best correlation to experimental logP coefficients of 50 mono-charged organic anions. In the case of the semiempirical methods, the best correlation between the HLI and logP values (the correlation coefficient r?=?0.9216) was obtained with the AM1 Ford–Wang parametric electrostatic potential charges. The Mulliken and Coulson charges calculated with the PM7 method can be used as an alternative to AM1, with the r values of 0.9107 and 0.8984, respectively. In the case of the DFT, the PBE/def2-TZVP natural population analysis charges produce the best correlation (r?=?0.9220). Nevertheless, in spite of a marginally lower performance (r?=?0.9159), the NPA charges computed at the PBE/def2-SVP level are more robust and can be regarded as the optimum choice for calculating the HLI values.

Open image in new window

Graphical abstract The hydrophilic/lipophilic index (HLI)

     

A Comparative Semiempirical, Ab initio and DFT Study of Decarbonylation of Cyclic Ketones Part 1: Thermal Fragmentation of Cyclopropanone

Calculations using different quantum mechanical methods including semiempirical (MNDO,AM1 and PM3), ab initio (RHF and MP2 calculations using the 6-311G and 6-311++G^** basis sets), and density functional theory (LSDA, BP, MIXBP and B3LYP, i.e., B3LYP/6-311+G**//B3LYP/6-31G*) have been performed on the thermal fragmentation of cyclopropanone to ethylene and carbon monoxide. All RHF calculations predict a concerted single step mechanism for this conversion. The estimated activation energies vary from 34.4 to 54.6 kcal·mol^-1, mainly localized around 37±2 kcal·mol^-1, depending on the method. Whereas the calculated RHF reaction energies also varied from 14.5 to -33.3 kcal·mol^-1, the B3LYP/6-311+G**//B3LYP/6-31G* method predicts the experimental value (-17.7 kcal·mol^-1) within experimental uncertainties. Remarkably, semiempirical AM1 and PM3 methods and simple DFT calculations, LSDA, predict comparable results to the more advanced methods. UHF ab initio calculations predict the same single step mechanism, whereas a multistep biradical mechanism with an unrealistically low activation energy is favored by the semiempirical methods. Structures of the activated complex of the single step mechanism, estimated by different methods, are very similar and consistent with a nonlinear cheletropic [_2s + _2a] reaction, as predicted by the orbital symmetry rules and earlier EHT calculations.Electronic Supplementary Material available.     

A search for quantum mechanical methods suitable for the conformational analysis of models related to some inverse agonists of the benzodiazepine receptor

The target of this work consists in a search, among several methods of quantum mechanical calculations, for the most suitable one (both with regard to accuracy and speed) for doing conformational analysis on models related to molecules belonging to a group of indol-3-yl glyoxylyl derivatives, many of which behave as ligands of the benzodiazepine receptor (BzR). Among the wide variety of ligands of the BzR, -carboline derivatives are included, many of which show a pharmacological profile similar to the one exhibited by compounds of interest to us. A structure analogy between our compounds and -carboline derivatives could be hypothesized, depending upon the preferred arrangement of the investigated molecules with respect to the dihedral angles of the indol glyoxylyl moiety. Therefore a potential energy surface scan has been carried out on selected models either by the quantum mechanical semiempirical methods MNDO, AM1, PM3 (by using the programs included in themopac package) or by a quantum mechanical SCFab initio method (by using thegaussian-80 program), both in order to find support for the choice of the more suitable semiempirical method, and in order to verify a structural analogy between our compounds and the -carboline ring. The results showed that AM1 and PM3, unlike MNDO, seem to be suitable semiempirical methods for doing potential energy surface scans on this kind of molecule. The preferred arrangement showed by the models with respect to the indol glyoxylyl moiety appeared to support the hypothesis about the structure analogy between the investigated ligands and -carboline derivatives.     

PM3 and AM1 Study on β-N-acetyl-muramic Acid and 3 Murein Related Derivatives

The energetically favoured conformations of -N-acetyl-Muramic acid, its C6-O-acetylated form, the methylamide and the methyl-glycoside have been investigated using the semiempirical PM3 and AM1 methods. All these compounds are either components or fragmentary structures of the murein network. The atomic coordinates of the starting set of the -N-acetyl-Muramic acid molecule have been obtained by a PM3 minimization of one saccharide molecule cut out from the murein single strand model proposed by Barnickel at al. [1]. The sidegroups of the derivatives have been introduced by a molecular editor. These conformations served as starting points in conformational space for a grid search by scanning all sidechain torsional angles for non-hydrogen atoms with exception of the N-acetyl group which was held in cisoid position (i.e. N2-H bond is parallel to C1-H and C3-H bond) and only minimized. The PM3 method with an additional amide correction potential and the AM1 method were used. The torsional angle distributions of the lactyl sidechain (free acid and methylamide), the C6-O-acetylated sidechain and the C1-methoxy sidechain have been investigated, showing distinct energetically favoured torsional angle regions. The results are compared to earlier studies on -N-acetyl-Muramic acid by J.S. Yadav et al. [2,3] who were using the MNDO and PCILO methods and by P.N.S. Yadav et al. [4] who were using the empirical MM2 force-field.Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1007/s0089460020190     

Hydrogen bonding interactions in E- or Z-2-phenyl-3-( X'-pyridyl)propenoic acid ( X=2, 3 or 4) assemblies--a molecular modeling study

The agglomeration properties of stereoisomeric 2-phenyl-3-(X-pyridyl)propenoic acids (X=2, 3 or 4) were studied by the PM3 semiempirical quantum chemical method. Calculations revealed that dimers kept together by the intermolecular hydrogen bonding interactions of the carboxylic groups could be built from both stereoisomers irrespective of the position of the nitrogen heteroatom. The dimers of the Z-isomers could also be built through (aromatic)C–H...N hydrogen bonds between the dimer units. The longest agglomerate was the pentamer of the dimers when the nitrogen was in the 2 position. Longer hydrogen-bonded agglomerates than dimers could only be constructed from the E-isomer with the nitrogen in position 4. Here, the trimer of the dimers proved to be the longest hydrogen-bonded entity and similarly to the Z-isomers, the dimer units are kept together by (aromatic)C–H...N intermolecular hydrogen bonds.Figure The helical structure of the pentamer of dimers for the Z-2-phenyl-3-(2-pyridyl)propenoic acid molecule, when in the initial geometry (aromatic)C(4)–H...N close contact was enforced.     

Normal coordinate analysis and vibrational spectra of 9-β-D-arabinofuranosyladenine hydrochloride (ara-A.HCl)

The vibrational spectra of a synthetic purine nucleoside with known antiviral activity, 9--D-arabino-furanosyladenine hydrochloride (ara-A.HCl) are reported. The Fourier transform infrared (FT-IR) and Fourier transform Raman (FT-Raman) spectra were recorded in the 4000-30 cm^–1 spectral region. The harmonic frequencies and potential energy distributions (PED) of the vibrational modes of ara-A.HCI were calculated by two different methods: a classical molecular mechanics method and a semiempirical molecular orbital (MO) method, PM3. The results of both computational methods, based on the Wilson GF method, are compared with observed spectra, and an assignment of the vibrational modes of ara-A.HCl is proposed on the basis of the potential energy distributions (PED). It is found that the wavenumbers can be calculated with remarkable accuracy (1% deviation in most cases), with the classical mechanics method, by transferring a sufficiently large set of available harmonic force constants, thus permitting a reliable assignment. The semiempirical MO method, PM3, is found to be useful for the assignment of experimental frequencies although it is less accurate (10% deviation). IR intensities calculated by this method did not coincide with the experimental values. Certain out-of-plane vibrations in the base, not reported in previous studies, have been observed. The performance of both methods was related to the crystallographic and ab initio data available. Previous normal coordinate calculations for the adenine base and the nucleoside 5-dGMP are compared with our results and discussed, in relation to the crystal structure of Ara-A.HCl. Correspondence to: A. Hernanz     

The origin and recycling of sedimented biogenic debris in a subalphine eutrophic lake (Lake Bled, Slovenia)

The areal distribution of organic C contents, ¹³C values, total N and P and biogenic Si contents in surficial sediments were used to study the distribution, origin and diagenetic transformations of sedimented biogenic debris in the eutrophic subalpine Lake Bled (Slovenia), which for most of the yearhas an anoxic hypolimnion. The influence of an allochthonous input, restricted to the western basin, was clearly traced by higher organic C and total N and P contents, higher ¹³C values, and higher sedimentation rate in comparison to the eastern basin. The low ¹³C values of sedimentary organic matter in the major part of the lake, lower than the ¹³C values of different types of organic matter, suggest that this sedimentary organic matter is most probably the product of a microbial community and not a residue of primary production.The temporal variation of benthic diffusive fluxes of NH₄, Si and PO₄, derived from modelling the pore water profiles, was related to sedimentation of phytoplanktonic blooms, while the PO₄ fluxes were also dependent on changing redox conditions at the sediment-water interface in the period of the winter-spring overtum. The removal of PO₄ in pore waters is probably due to the adsorption of phosphate and precipitation of apatite and vivianite. The budget of C, N and P at the sediment-water interface revealed a high recycling efficiency (>70%), also confirmed by the rather uniform (or only slightly decreasing) vertical profiles of organic C, total N and P in sediment cores and C/N and C/P ratios. The percentage of biogenic Si recycling is low (<10%), suggesting its removal in sediments.     

Validation of semiempirical methods for modeling of corrinoid systems

Several semiempirical methods (MNDO-d, PM3tm, PM3-d, PM5, PM6, and AM1-d) have been tested against experimental data and density functional theory (DFT) results in search for the best methods that can be used for quantum-mechanical-molecular mechanics (QM/MM) modeling of corrinoid systems of vitamin B(12) co-factor. It has been found that the PM6 parametrization in its present form gives results closest to hybrid DFT calculations that are most widely used thus far. In comparison with pure DFT and experimental data the best agreement is obtained for PM3tm parametrization, while PM6 yields slightly worse results. AM1-d yields bad geometry of the corrin moiety. The worst performance was observed for MNDO-d, which has severe problem with position and orientation of the alpha-ligands.     

用分子轨道法探讨油田废水中取代苯系物的毒性效应

采用半经验的分子轨道AM1方法中的MOPAC软件包计算了55种在油田废水中易于检出的取代苯系物的分子轨道能（EHOMO,ELUMO,ENHOMO,ENLUMO）、分子生成熟（△H^f0)和偶极矩（μ）等量化参数,结合一阶价分子连接性指数（^1X^V)和正辛醇／水分配系数(logP)与实验所得发光菌的半数活性浓度（EC50）成功建立了多参数定量－结构活性关系模式,在分类建模的基础上,又获得了仅包含1^X^V和EHOMO两个参数的55种取代苯系物的定量结构－活性相关模式,探讨了不同取代基的毒性作用机制,结果表明,量化参数与物化参数结合能够很好地预测油田废水中具有不同取代基团的取代苯系物的生物活性,预测模式中量子化学参数的出现有利于深入探讨油田废水中有机污染物的毒性效应。     

Cathepsin B and D activity in stimulated peritoneal macrophages

Summary Highly sensitive and specific synthetic substrates were used to quantitate cathepsin B and D activity in peritoneal macrophages in response to stimulation in vivo with mineral oil and thioglycollate. After intraperitoneal instillation of mineral oil the activity of cathepsin B increased significantly (to 15 300 units/mg protein versus 7 340 in saline controls), reaching values approaching those found in alveolar macrophages (18 400 units/mg protein). Significantly greater stimulation of enzyme activity was obtained after intraperitoneal instillation of thioglycollate (23 600 units/mg protein). Cathepsin D activity also increased significantly after both mineral oil and thioglycollate. However, the increase was moderate (from 806 to about 1 200 units/mg protein), remaining still more than six times lower-than in alveolar macrophages. The data are the first to demonstrate that cathepsin B activity can be stimulated in vivo in peritoneal macrophages by instillation of agents that induce acute inflammation. They also point to a differential control of expression of cathepsin B and D activity in both peritoneal and alveolar macrophages in spite of the common lysosomal origin of the two enzymes.Abbreviations Cbz -N-benxyloxycarbonyl - 2NA 2-naphthylamine - EDTA ethylenediamine tetraacetate - DMSO dimethylsulfoxide - PBS phosphate-buffered saline - PM peritoneal macrophage - AM alveolar macrophage     

A new method for the determination of microbial activity and critical logP in the presence of organic solvents

A new method to determine microbial activity and critical logP of an organism in the presence of organic solvents has been developed which involves direct contact with a solvent, and a measurement of the developing colony size. This technique has been used to estimate the critical logP of Alcaligenes xylosoxidans Y234, and, although the critical logP for this organism is 3.5, solvents with logP values of up to 4.5 can still reduce microbial activity by up to 55% of the uninhibited amount.