Adsorption and dimerisation of thiol molecules on Au(111) using a Z-matrix approach in density functional theory

The adsorption energetics of methanethiolate and benzenethiolate on Au(111) have been calculated using periodic density functional theory (DFT), based on the SIESTA methodology, with an internal coordinates implementation for geometry input and structure optimisation. Both molecules are covalently bound with interaction energies of 1.85 and 1.43 eV for methanethiolate and benzenethiolate, respectively. The preferred binding site is slightly offset from the bridge site in both cases towards the fcc-hollow. The potential energy surfaces (PES) have depths of 0.36 and 0.22 eV, the hollow sites are local maxima in both cases, and there is no barrier to diffusion of the molecule at the bridge site. The corresponding dimers are weakly bound for methanethiolate and benzenethiolate, with binding energies of 0.38 and 0.16 eV, respectively, and the preferred binding geometry is with the two sulphur atoms close to adjacent atop sites. The barrier to dissociation of the dimer dimethyl disulphide is estimated to lie between 0.3 and 0.35 eV.     

The anticancer drug bleomycin investigated by density functional theory

Activated bleomycin (ABLM) is an oxygenated iron drug complex which embodies the drug's DNA-cleaving activity. This activity is exercised on DNA, if present, but if DNA is absent, the drug itself is inactivated. We have employed quantum density functional theory (DFT)-based methods to investigate (i) the structure of the Fe(II)BLM complex that is first formed in the human body after drug's administration, and (ii) the activation mechanism of the O–O bond present in the ABLM. We have identified the controversial second axial ligand as the endogenous oxygen atom of the carbamoyl group. Our first principles molecular dynamics (MD) simulations indicate a homolytic cleavage as the mechanism of the O–O bond activation in the ABLM complex.     

Probing the coordination properties of glutathione with transition metal ions (Cr2+,Mn2+,Fe2+,Co2+, Ni2+,Cu2+,Zn2+,Cd2+,Hg2+) by density functional theory

Complexes formed by reduced glutathione (GSH) with metal cations (Cr²⁺, Mn²⁺,Fe²⁺,Co²⁺,Ni²⁺,Cu²⁺,Zn²⁺,Cd²⁺,Hg²⁺) were systematically investigated by the density functional theory (DFT). The results showed that the interactions of the metal cations with GSH resulted in nine different stable complexes and many factors had an effect on the binding energy. Generally, for the same period of metal ions, the binding energies ranked in the order of Cu²⁺>Ni²⁺>Co²⁺>Fe²⁺>Cr²⁺>Zn²⁺>Mn²⁺; and for the same group of metal ions, the general trend of binding energies was Zn²⁺>Hg²⁺>Cd²⁺. Moreover, the amounts of charge transferred from S or N to transition metal cations are greater than that of O atoms. For Fe²⁺,Co²⁺,Ni²⁺,Cu²⁺,Zn²⁺,Cd²⁺ and Hg²⁺ complexes, the values of the Wiberg bond indices (WBIs) of M-S (M denotes metal cations) were larger than that of M-N and M-O; for Cr²⁺ complexes, most of the WBIs of M-O in complexes were higher than that of M-S and M-N. Furthermore, the changes in the electron configuration of the metal cations before and after chelate reaction revealed that Cu²⁺, Ni²⁺,Co²⁺ and Hg²⁺ had obvious tendencies to be reduced to Cu⁺,Ni⁺,Co⁺ and Hg⁺ during the coordination process.     

Effect of microsolvation on zwitterionic glycine: an ab initio and density functional theory study

The effect of microsolvation on zwitterionic glycine, considering both (-NH3(+)) as proton donor and (-COO(-)) as proton acceptor at correlated ab initio (MP2) level and density functional methods (B3LYP, PW91, MPW1PW91 and PBE) using 6-311++G** basis set has been reported. DFT methods have been employed so as to compare the performance/quality of different gradient-corrected correlation functionals (PW91, PBE), hybrid functionals (B3LYP, MPW1PW91) and to predict the near quantitative structural and vibrational properties, at reduced computational cost. B3LYP method outperforms among the different DFT methods for the computed hydrogen bond distances and found closer to the value obtained by correlated MP2 level, whereas MPW1PW91 and PBE methods shows very similar values but approximately 0.03 A less, compared to B3LYP method. MP2 calculation and single point CCSD(T)//MP2 calculation have been considered to decompose the interaction energy, including corrections for basis set superposition error (BSSE). Moreover, charge distribution analysis has also been carried out to understand the long raised questions, how and why the two body energies have significant contribution to the total binding energy.     

Ab initio prediction of optical rotation: comparison of density functional theory and Hartree-Fock methods for three 2,7,8-trioxabicyclo[3.2.1]octanes

We report ab initio calculations of the frequency-dependent electric dipole-magnetic dipole polarizabilities, beta(nu), at the sodium D line frequency and, thence, of the specific rotations, [alpha](D), of 2,7,8-trioxabicyclo[3.2.1]octane, 1, and its 1-methyl derivative, 2, using the Density Functional Theory (DFT) and Hartree-Fock/Self-Consistent Field (HF/SCF) methodologies. Gauge-invariant (including) atomic orbitals (GIAOs) are used to ensure origin-independent [alpha](D) values. Using large basis sets which include diffuse functions DFT [alpha](D) values are in good agreement with experimental values (175.8 degrees and 139.2 degrees for (1S,5R)-1 and -2, respectively); errors are in the range 25-35 degrees. HF/SCF [alpha](D) values, in contrast, are much less accurate; errors are in the range 75-95 degrees. The use of small basis sets which do not include diffuse functions substantially lowers the accuracy of predicted [alpha](D) values, as does the use of the static limit approximation: beta(nu) approximately beta(o). The use of magnetic-field-independent atomic orbitals, FIAOs, instead of GIAOs, leads to origin-dependent, and therefore nonphysical, [alpha](D) values. We also report DFT calculations of [alpha](D) for the 1-phenyl derivative of 1, 3. DFT calculations find two stable conformations, differing in the orientation of the phenyl group, of very similar energy, and separated by low barriers. Values of [alpha](D) predicted using two different algorithms for averaging over phenyl group orientations are in good agreement with experiment. In principle, the absolute configuration (AC) of a chiral molecule can be assigned by comparison of the optical rotation predicted ab initio to the experimental value. Our results demonstrate the critical importance of the choice of ab initio methodology in obtaining reliable optical rotations and, hence, ACs, and show that, at the present time, DFT constitutes the method of choice.     

Interaction energies in non-covalently bound intermolecular complexes derived using the subsystem formulation of density functional theory

Interaction energies for a representative sample of 39 intermolecular complexes are calculated using two computational approaches based on the subsystem formulation of density functional theory introduced by Cortona (Phys. Rev. B 44:8454, 1991), adopted for studies of intermolecular complexes (Wesolowski and Weber in Chem. Phys. Lett. 248:71, 1996). The energy components (exchange-correlation and non-additive kinetic) expressed as explicit density functionals are approximated by means of gradient-free- (local density approximation) of gradient-dependent- (generalized gradient approximation) approximations. The sample of the considered intermolecular complexes was used previously by Zhao and Truhlar to compare the interaction energies derived using various methods based on the Kohn-Sham equations with high-level quantum chemistry results considered as the reference. It stretches from rare gas dimers up to strong hydrogen bonds. Our results indicate that the subsystem-based methods provide an interesting alternative to that based on the Kohn-Sham equations. Local density approximation, which is the simplest approximation for the relevant density functionals and which does not rely on any empirical data, leads to a computational approach comparing favorably with more than twenty methods based on the Kohn-Sham equations including the ones, which use extensively empirical parameterizations. For various types of non-bonding interactions, the strengths and weaknesses of gradient-free and gradient-dependent approximations to exchange-correlation and non-additive kinetic energy density functionals are discussed in detail.     

Conformational preferences of 1,4,7-trithiacyclononane: a molecular mechanics and density functional theory study

Conformational preferences of 1,4,7-trithiacyclononane were studied using a highly efficient sampling technique based on local nonstochastic deformations and the MM2(91) force field. The results show that conformers that the molecule adopts in the crystal state were found to be low-energy conformers (LECs) within 5 kcal mol(-1) of the global minimum. A conformation with C1 symmetry was the global minimum and the C3 and C2 conformations were calculated to be 0.03 and 1.78 kcal mol(-1) higher in energy, respectively. The structures were further minimized using Density Functional Theory (DFT) calculations with two different functionals. The C2 and the C1 conformations were found to be LECs with the C3 conformation more than 4.0 kcal mol(-1) above the global minimum. The relative energies and structural ordering obtained using the BP86 functional are in agreement with the previously reported relative energies calculated using second-order Moller-Plesset (MP2) ab initio calculations. With the energy ordering being dependent on the molecular mechanics force field used, the approach of MM-->DFT (searching exhaustively the available conformational space at the MM level followed by generating the energy ordering through DFT calculations) appears to be appropriate for thiacrown ethers.     

Adsorption of CO on Cu (110) and (100) surfaces using COSMO-based DFT

Density functional theory (DFT) combined with the conductor-like solvent model (COSMO) can provide valuable atomistic level insights into CO adsorption on Cu surface interactions in liquid paraffin. The objective of this research was to investigate the solvent effect of liquid paraffin. It was found that both structural parameters and relative energies are very sensitive to the COSMO solvent model. Solvent effects can improve the stability of CO adsorption on Cu (110) and (100) surfaces and the extent of CO activation.     

Molecular simulation of the preferential adsorption of CH4 and CO2 in middle-rank coal

The adsorption of the CO₂/CH₄ mixture in coal affects the CO₂-enhanced coalbed methane recovery project. To gain a better understanding of CH₄ and CO₂ interaction with middle-rank coal, we developed a molecular concept with support for the sorption of CH₄ and CO₂ on Ximing-8 coal (XM-8) (1.8% vitrinite reflectance). A XM-8 coal model was built by using molecular dynamic (MD) simulations. The molecular simulations were established by the Grand Canonical Monte Carlo and MD methods to study the effects of the temperature, pressure, and species bulk mole fraction on the pure component adsorption isotherms, isosteric heat and adsorption selectivity. It turns out that the CO₂ selectivity decreases as the pressure and its own bulk mole fraction increases, but it increases as temperature increases, and the selectivity values are not always greater than 1. The interactions between the small molecules and XM-8 were determined by using density functional theory. It was found that the interactions between the CO₂ and XM-8 surface is greater, particularly for the heteroatoms than CH₄. The adsorption selectivity and interaction were simultaneously used to reveal that the advantageously substituted range is high temperature, low pressure and a high content of heteroatoms.     

Molecular insight into adsorption affinities of Carmustine drug on boron and nitrogen doped functionalized single-walled carbon nanotubes using density functional theory including dispersion correction calculations and molecular dynamics simulation

Abstract

We report a quantum mechanics calculation and molecular dynamics simulation study of Carmustine drug (BNU) adsorption on the surface of nitrogen (N) and boron (B) doped-functionalized single-walled carbon nanotubes. The stability of the optimized complexes is determined on the basis of relative adsorption energy (ΔE_ads). The ΔE_ads results claim that drug molecule tends to adsorb on the nitrogen and boron doped functionalized tubes with the energy values in the range of ?61.177 to ?95.806?kJ/mol. Based on the obtained results, it is observed that N-doping compared with B-doping has improved more effectively drug absorption on the surface of functionalized nanotube. The results of Atoms in Molecule calculations indicate that drug adsorbs molecularly via hydrogen bonds interactions on the surface doped-functionalized carbon nanotubes. Moreover, molecular dynamics simulation is performed to investigate the dynamics behavior of the drug molecules on the nitrogen-doped functionalized carbon nanotube (f-NNT) and functionalized carbon nanotube (f-CNT). The higher average calculated electrostatic and van der Waals energies as well as higher number of intermolecular hydrogen bonds in BNU-f-NNT in comparison with BNU-f-CNT model suggest the more effectual interaction between drug molecules and nitrogen-doped functionalized carbon nanotube.

Communicated by Ramaswamy H. Sarma     

Density functional theory for the selective adsorption of small molecules on a surface modified with polymer brushes

A density functional method based on weighted density approximation is extended to study the selective adsorption of small molecules on a surface modified with end-grafted square-well chains. The excess part of the Helmholtz free energy functional is divided into two components: the hard sphere repulsion and the square-well attraction. The equation of state for hard sphere chain fluids developed by Liu et al. is used to calculate the repulsive part of the excess Helmholtz free energy functional, and the equation of state for square-well chain fluid with variable range developed by Li et al. is employed to calculate the attractive part. With this theoretical model, we examine the physical properties of the grafted polymer and the selective adsorption of small molecules on the modified surface.     

A density functional theory study on the role of His-107 in arylamine N-acetyltransferase 2 acetylation

Arylamine N-acetyltransferases (NATs, EC 2.3.1.5) catalyze an acetyl group transfer from acetyl coenzyme A (AcCoA) to primary arylamines, and are responsible for the biotransformation and metabolism of drugs, carcinogens, etc. Structure analysis revealed that His-107 was likely the residue accountable for mediating acetyl transfer. We have examined the full catalytic mechanism of this system by means of DFT method. The results indicate that if the acetyl group directly transferred from the donor, p-nitrophenyl acetate, to the acceptor, cysteine, the high activation energy will be a great hindrance. These energies have dropped a little in a range of 20-25 kJ/mol when His-107 is assisting the transfer process. However, when protonated His-107 is mediating the reaction, the activation energies have dropped about 70-85 kJ/mol. Our calculations strongly support an enzymatic acetylation mechanism that experiences a thiolate-imidazolium pair, which have verified the presumption from experiments.     

Determination of the absolute configurations of chiral organometallic complexes via density functional theory calculations of their vibrational circular dichroism spectra: The chiral chromium tricarbonyl complex of N-pivaloyl-tetrahydroquinoline

The racemate of the chiral tricarbonyl-η⁶-arene-chromium(0) complex, tricarbonyl-η⁶-N-pivaloyl-tetrahydroquinoline-chromium(0), 1, has been synthesized and resolved using chromatography on a (R,R)-Whelk-O1 column. The Absolute Configuration (AC) of 1 has been determined using vibrational circular dichroism (VCD) spectroscopy. The VCD spectrum of 1 has been predicted using the Stephens equation for vibrational rotational strengths, implemented using density functional theory (DFT) in the gaussian program. Using the B3PW91 functional and the 6-311++G (2d,2p) basis set, the predicted VCD spectrum of S-1 is in excellent agreement with the experimental VCD spectrum of (+)-1, leading unambiguously to the AC S-(+). It is concluded that VCD is a useful technique for determining the ACs of chiral organometallic complexes, given the use of optimum functionals and basis sets.     

Deformability of Fe(II)CO and Fe(III)CN groups in heme protein models: nonlocal density functional theory calculations

 High-quality nonlocal density functional calculations have been carried out on the deformability of Fe(II)CO and Fe(III)CN units in model compounds of heme proteins. The results confirm our previous finding with a local functional that the Fe(II)CO unit is significantly deformable with respect to tilting and bending. This deformability stems in large part from a large, negative interaction constant between the FeC tilt and FeCO bend coordinates. The Fe(III)CN unit is also significantly deformable, but in this case the deformability results from a very small Fe(III)CN bend force constant and the ability of the cyano nitrogen to act as a hydrogen bond acceptor. The prediction that the energetic penalty associated with deforming the Fe(II)CO unit is relatively modest indicates that such deformations are unlikely to be the dominant contributor to myoglobin's discrimation against CO in favor of O₂. Received, accepted: 16 June 1997     

The influence of nicotine on pioglitazone encapsulation into carbon nanotube: the investigation of molecular dynamic and density functional theory

In this work, molecular dynamics simulations of the insertion of pioglitazone into the nanotube with chirality (10, 10) at 400 K and 1 bar in the presence and absence of nicotine molecules and in different drug concentrations have been studied. The main aim is consideration of the effect of nicotine in the drug encapsulation process. The results indicate that encapsulation of pioglitazone could be attributed to the water flow via van der Waals and hydrophilic interactions. Because of the existence of the partial π–π interactions between aromatic rings of pioglitazone and the conjugated aromatic rings of nanotube, pioglitazone molecule can enter inside the nanotube. Some physical properties such as hydrogen bonding, number of contacts, also, the diffusion coefficient of the pioglitazone and water molecules, and variation of the center of mass have been calculated during the simulation. Furthermore, computing the electronic structure has also been done on model systems for quantitative determination of the adsorption energy (E_ads). The B3LYP/6-31G* level calculations on four different configurations of pioglitazone/carbon nanotube (CNT) and nicotine/CNT show that the interaction of drug with the inside of the nanotube is stronger than the other forms.     

Syntheses, characterization, density functional theory calculations, and activity of tridentate SNS zinc pincer complexes

A series of tridentate SNS ligand precursors were metallated with ZnCl₂ to give new tridentate SNS pincer zinc complexes. The zinc complexes serve as models for the zinc active site in liver alcohol dehydrogenase (LADH) and were characterized with single crystal X-ray diffraction, ¹H, ¹³C, and HSQC NMR spectroscopies and electrospray mass spectrometry. The bond lengths and bond angles of the zinc complexes correlate well to those in horse LADH. The zinc complexes feature SNS donor atoms and pseudotetrahedral geometry about the zinc center, as is seen for liver alcohol dehydrogenase. The SNS ligand precursors were characterized with ¹H, ¹³C, and HSQC NMR spectroscopies and cyclic voltammetry, and were found to be redox active. Gaussian calculations were performed and agree quite well with the experimentally observed oxidation potential for the pincer ligand. The zinc complexes were screened for the reduction of electron poor aldehydes in the presence of a hydrogen donor, 1-benzyl-1,4-dihydronicotinamide (BNAH). The zinc complexes enhance the reduction of electron poor aldehydes. Density functional theory calculations were performed to better understand why the geometry about the zinc center is pseudo-tetrahedral rather than pseudo-square planar, which is seen for most pincer complexes. For the SNS tridentate pincer complexes, the data indicate that the pseudo-tetrahedral geometry was 43.8 kcal/mol more stable than the pseudo-square planar geometry. Density functional theory calculations were also performed on zinc complexes with monodentate ligands and the data indicate that the pseudo-tetrahedral geometry was 30.6 kcal/mol more stable than pseudo-square planar geometry. Overall, the relative stabilities of the pseudo-tetrahedral and pseudo-square planar systems are the same for this coordination environment whether the ligand set is a single tridentate SNS system or is broken into three separate units. The preference of a d¹⁰ Zn center to attain a tetrahedral local environment trumps any stabilization gained by removal of constraints within the ligand set.     

Synthesis, crystal structures and density functional theory study of dimeric copper(II) complexes of a tridentate Schiff-base ligand

Oxalate- or 4,4′-bipyridine-bridged dimeric copper(II) complexes, [Cu₂L₂(μ-ox)] (1) and [Cu₂L₂(μ-bipy)](BF₄)₂ (2) [where ox = oxalate, bipy = 4,4′-bipyridine, HL = N-(1H-pyrrol-2-ylmethylene)-2-pyridineethanamine, L⁻ = HL−H⁺], have been synthesised and characterised by elemental analysis, IR, UV-Vis and single crystal X-ray diffraction. Crystal structure determinations carried out on 1 and 2 reveal that 1 is an oxalate-bridged centrosymmetrical square pyramidal dimeric copper(II) complex while 2 is a 4,4′-bipyridine-bridged non-centrosymmetric square planar dinuclear copper(II) complex. Comparison of the optimised geometries with the corresponding crystal structures suggests that the B3LYP/LANL2DZ level can reproduce the structures of 1 and 2 on the whole. The electronic spectra of 1 and 2 predicted by B3LYP/LANL2DZ method show some blue shifts compared with their experimental data. Thermal analysis carried out on 1 shows that there is only one exothermal peak at about 260 °C and the residue is presumably Cu₂O₄N₆.     

青蒿素及其类似物抗疟构效关系的DFT研究

为了从原子水平上揭示青蒿素及其类似物的结构与抗疟活性之间的关系,运用密度泛函理论DFT方法,在B3LYP/6-31G*水平上对青蒿素及其类似物二氢青蒿素、蒿甲醚和青蒿琥酯的结构和性质进行了理论计算。从分子的平衡构型、Wiberg键级、溶剂化能、偶极矩和静电势等方面分析了青蒿素及其类似物的抗疟构效关系。结果表明,青蒿素及其类似物结构中七元环上的过氧桥键、醚氧键以及六元环上的内酯结构是其抗疟作用的关键活性位,过氧桥键处负的静电势越多,青蒿素与血红素的相互作用越强,分子的抗疟活性越强。理论预测四个药物分子的抗疟活性顺序为:青蒿素<二氢青蒿素<蒿甲醚<青蒿琥酯,与实验活性结果一致。     

On the oxidation of alkyl and aryl sulfides by [(Me3TACN)MnO(OH)2]: A density functional study

Density functional theory suggests that the formal 2-electron oxidation of sulfides, RR′S, to sulfoxides, by the model Mn^VO catalyst, [(TACN)Mn^V O(OH)₂]⁺, proceeds in two quite distinct 1-electron steps. Transfer of the first electron is barrierless and generates a sulfur radical cation, antiferromagnetically coupled to a Mn^IV centre via a covalent μ-oxo bridge. The second electron-transfer step is accompanied by migration of the oxygen atom to the sulfur centre, and is rate-determining. The absence of a barrier in the first step, where a sulfur radical is formed, means that the presence of electron-donating or withdrawing substituents on the sulfide has only a minor impact on the rate of reaction.     

A density functional theory-based analysis of the structural,topological and electronic properties of gemcitabine drug adsorption on the pyrrolidine functionalized single-walled carbon nanotube

The stability of gemcitabine anticancer drug on the functionalized (8,0) zigzag carbon nanotube as a drug delivery vehicle is studied within the formalism of the density functional theory calculations to understand the role of the pyrrolidine functional group in binding the adsorbed molecule to the drug delivery system as well as improving water solubility. The binding energies, natural bond orbital calculations, and the quantum theory of atoms in molecules results are obtained to provide more evidences related to the intermolecular interaction between gemcitabine drug and the functionalized nanotube. The negative binding energy corresponds to favorable binding of gemcitabine drug to the functionalized nanotube and presence of the active sites available for hydrogen bond formation facilitates better drug binding to the nanotube sidewall. The results presented in this article indicate that pyrrolidine functionalized carbon nanotube seems to be a novel material for drug delivery applications.

Communicated by Heidar Moradnia