Density functional study of oxygen vacancy formation and spin density distribution in octahedral ceria nanoparticles

We report plane wave basis density functional theory (DFT) calculations of the oxygen vacancies formation energy in nanocrystalline CeO _2-x in comparison with corresponding results for bulk and (111) CeO₂ surface. Effects of strong electronic correlation of Ce4f states are taken into account through the use of an effective on-site Coulomb repulsive interaction within DFT+U approach. Different combinations of exchange-correlation functionals and corresponding U values reported in the literature are tested and the obtained results compared with experimental data. We found that both absolute values and trends in oxygen vacancy formation energy depend on the value of U and associated with degree of localization of Ce4f states. Effect of oxygen vacancy and geometry optimization method on spatial spin distribution in model ceria nanoparticles is also discussed.     

QM/MM investigation of structure and spectroscopic properties of a vanadium-containing peroxidase

We present a comprehensive analysis of the most likely ground state configuration of the resting state of vanadium dependent chloroperoxidase (VCPO) based on quantum mechanics/molecular mechanics (QM/MM) evaluations of ground state properties, UV-vis spectra and NMR chemical shifts. Within the QM/MM framework, density functional theory (DFT) calculations are used to characterize the resting state of VCPO via time-dependent density functional theory (TD-DFT) calculations of electronic excitation energies and NMR chemical shifts. Comparison with available experimental data allows us to determine the most likely protonation state of VCPO, a state which results in a doubly protonated axial oxygen, a site largely stabilized by hydrogen bonds. We found that the bulk of the protein that is beyond the immediate layer surrounding the cofactor, has an important electrostatic effect on the absorption maximum. Through examination of frontier orbitals, we analyze the nature of two bound water molecules and the extent to which relevant residues in the active site influence the spectroscopy calculations.     

Molecular interactions of alcohols with zeolite BEA and MOR frameworks

Zeolites can adsorb small organic molecules such as alcohols from a fermentation broth. Also in the zeolite-catalyzed conversion of alcohols to biofuels, biochemicals, or gasoline, adsorption is the first step. Several studies have investigated the adsorption of alcohols in different zeolites experimentally, but computational investigations in this field have mostly been restricted to zeolite MFI. In this study, the adsorption of C1–C4 alcohols in BEA and MOR was investigated using density functional theory (DFT). Calculated adsorption geometries and the corresponding energies of the designed cluster models were comparable to periodic calculations, and the adsorption energies were in the same range as the corresponding computational and experimental values reported in the literature for zeolite MFI. Thus, BEA and MOR may be good adsorption materials for alcohols in the field of downstream processing and catalysis. Aside from the DFT calculations, adsorption isotherms were determined experimentally in this study from aqueous solutions. For BEA, the adsorption of significant amounts of alcohol from aqueous solution was observed experimentally. In contrast, MOR was loaded with only a very small amount of alcohol. Although differences were found between the affinities obtained from gas-phase DFT calculations and those observed experimentally in aqueous solution, the computational data presented here represent molecular level information on the geometries and energies of C1–C4 alcohols adsorbed in zeolites BEA and MOR. This knowledge should prove very useful in the design of zeolite materials intended for use in adsorption and catalytic processes, as it allows adsorption behavior to be predicted via judiciously designed computational models.     

Density functional theory calculations on entire proteins for free energies of binding: Application to a model polar binding site

In drug optimization calculations, the molecular mechanics Poisson‐Boltzmann surface area (MM‐PBSA) method can be used to compute free energies of binding of ligands to proteins. The method involves the evaluation of the energy of configurations in an implicit solvent model. One source of errors is the force field used, which can potentially lead to large errors due to the restrictions in accuracy imposed by its empirical nature. To assess the effect of the force field on the calculation of binding energies, in this article we use large‐scale density functional theory (DFT) calculations as an alternative method to evaluate the energies of the configurations in a “QM‐PBSA” approach. Our DFT calculations are performed with a near‐complete basis set and a minimal parameter implicit solvent model, within the self‐consistent calculation, using the ONETEP program on protein–ligand complexes containing more than 2600 atoms. We apply this approach to the T4‐lysozyme double mutant L99A/M102Q protein, which is a well‐studied model of a polar binding site, using a set of eight small aromatic ligands. We observe that there is very good correlation between the MM and QM binding energies in vacuum but less so in the solvent. The relative binding free energies from DFT are more accurate than the ones from the MM calculations, and give markedly better agreement with experiment for six of the eight ligands. Furthermore, in contrast to MM‐PBSA, QM‐PBSA is able to correctly predict a nonbinder. Proteins 2014; 82:3335–3346. © 2014 Wiley Periodicals, Inc.     

d(A)<Subscript>3</Subscript>d(T)<Subscript>3</Subscript> and d(G)<Subscript>3</Subscript>d(C)<Subscript>3</Subscript> B-DNA mini-helixes: the DFT/M06-2x and DFT/B97-D3 comparison of geometrical and energetic characteristics

We report the comprehensive DFT based comparison of geometrical and energetic parameters of the d(A)₃·d(T)₃ and d(G)₃·d(C)₃ nucleic acid mini-helixes performed at B97-D3 and M06-2× levels of theory. We studied the ability of mini-helixes to retain the conformation of B-DNA in the gas phase and under the influence of water bulk, uncompensated charges, and counter-ions. The def2-SV(P) and 6-31G(d,p) basis sets have been used for B97-D3 and M06-2× calculations, correspondently. To estimate basis set superposition error, the recently developed semi-empirical procedure that calls geometrical counterpoise type correction for inter- and intra—molecular basis set superposition error (gcp) has been used in the case of def2-SV(P) basis set. We found that both considered DFT functionals predict very similar results for geometrical ad energetic characteristics. We also found that in contrast to average classical molecular dynamics and data of simple geometrical models, both considered DFT functionals predict the existence of duplex specific geometries. A prediction of interaction energies of d(A)₃d(T)₃ and d(G)₃d(C)₃ duplexes accomplished in this study also verifies the applied models and confirms reliability of the new computational gcp technique.     

Density functional theory calculations on the thermodynamic properties of polynitrosoprismanes

A series of polynitrosoprismanes, C(6)H(6 - n )(NO)( n ) (n?=?1-6), considered as high energy density compounds (HEDCs), have been designed computationally. We calculated the electronic structures, the heats of formation, the specific enthalpies of combustion, the bond dissociation energies, and the strain energies of the title compounds using density functional theory (DFT) with the 6-311G** basis set. It was found that the ΔE (LUMO-HOMO) values of the title compounds decrease as the number of nitroso groups increase, and the energy gaps of the prismane derivatives are much lower than that of TATB. Their high positive heats of formation indicate that polynitrosoprismanes can store a great deal of energy. Furthermore, the HOFs for the nitrosoprismane series were observed to decrease until three nitroso groups were connected to the prismane skeleton. For the polynitrosoprismanes, the trigger bond was confirmed to be the C-C bond in the skeleton. According to our calculations, all nitrosoprismanes appear to have large strain energies, and these calculations can provide basic information that may prove useful for the molecular design of novel high energy density materials.     

Linkage isomerism of pentaammine(dimethylsulfoxide)ruthenium(II/III) complexes: A theoretical study

Density functional theory (DFT) calculations have been performed for understanding the linkage isomerism of [Ru^II/III(NH₃)₅(dmso)]^2+/3+ (dmso = dimethylsulfoxide) from a theoretical point of view. In particular, we focus on the interchange between O-bonded and S-bonded structures of the dmso ligand by oxidation/reduction. We have examined five different exchange-correlation functionals (SVWN, BP86, mPWPW91, B3PW91, and B3LYP) in our DFT calculations and found that the relative stabilities of the O-bonded and S-bonded structures are largely dependent on the functional employed. From detailed analyses of atomic charge distributions, it has been found that the calculated atomic charges on the central metal ions are strongly correlated with the relative energies. We also studied the effect of solvation on the linkage isomerism using continuum solvation models.     

Relating Frontier Orbital Energies from Voltammetry and Photoelectron Spectroscopy to the Open‐Circuit Voltage of Organic Solar Cells

For 19 diketopyrrolopyrrole polymers, the highest occupied molecular orbital (HOMO) energies are determined from i) the oxidation potential with square‐wave voltammetry (SWV), ii) the ionization potential using ultraviolet photoelectron spectroscopy (UPS), and iii) density functional theory (DFT) calculations. The SWV HOMO energies show an excellent linear correlation with the open‐circuit voltage (V_oc) of optimized solar cells in which the polymers form blends with a fullerene acceptor ([6,6]‐phenyl‐C₆₁‐butyl acid methyl ester or [6,6]‐phenyl‐C₇₁‐butyl acid methyl ester). Remarkably, the slope of the best linear fit is 0.75 ± 0.04, i.e., significantly less than unity. A weaker correlation with V_oc is found for the HOMO energies obtained from UPS and DFT. Within the experimental error, the SWV and UPS data are correlated with a slope close to unity. The results show that electrochemically determined oxidation potentials provide an excellent method for predicting the V_oc of bulk heterojunction solar cells, with absolute deviations less than 0.1 V.     

“Nanostructures in Thin-film Epitaxy: Exploring and Exploiting Substrate-mediated Interactions”

We review recent advances made in understanding the ramifications of substrate-mediated interactions for thin-film growth. Experimental studies and first-principles calculations with density-functional theory (DFT) indicate that substrate-mediated interactions can significantly influence thin-film growth. We review the findings from our kinetic Monte Carlo simulations used to model the growth of thin films, both with and without substrate-mediated interactions. For Ag heteroepitaxy on Pt(1?1?1), the pair interaction energies and adsorbate diffusion barriers were obtained from DFT calculations. Island densities for this system show significant deviations from what is predicted by classical nucleation theory. The electronic interactions created by the adsorbed atoms lead to the formation of repulsive barriers surrounding small islands and, as a result, sharp island-size distributions are produced. The island-size distributions can be manipulated by changing the growth conditions to yield desirable island sizes and shapes.     

Atomistic simulations of materials for optical chemical sensors: DFT-D calculations of molecular interactions between gas-phase analyte molecules and simple substrate models

The structures of complexes of some small molecules (formaldehyde, acetaldehyde, ammonia, methylamine, methanol, ethanol, acetone, benzene, acetonitrile, ethyl acetate, chloroform, and tetrahydrofuran, considered as possible analytes) with ethylbenzene and silanol (C₆H₅C₂H₅ and SiH₃OH, considered as models of polystyrene and silica gel substrates) and with acridine (C₁₃H₉N, considered as a model of an indicator dye molecule of the acridine series) and the corresponding interaction energies have been calculated using the DFT-D approximation. The PBE exchange-correlation potential was used in the calculations. The structures of complexes between the analyte and the substrate were determined by optimizing their ground-state geometry using the SVP split-valence double-zeta plus polarization basis set. The complex formation energies were refined by single-point calculations at the calculated equilibrium geometries using the sufficiently large triple-zeta TZVPP basis set. The calculated interaction energies are used to assess the possibility of using dyes of the acridine series adsorbed on a polystyrene or silica substrate for detecting the small molecules listed above.     

An improved generalized AMBER force field (GAFF) for urea

We describe an improved force field parameter set for the generalized AMBER force field (GAFF) for urea. Quantum chemical computations were used to obtain geometrical and energetic parameters of urea dimers and larger oligomers using AM1 semiempirical MO theory, density functional theory at the B3LYP/6-31G(d,p) level, MP2 and CCSD ab initio calculations with the 6-311++G(d,p), aug-cc-pVDZ, aug-cc-pVTZ, and aug-cc-pVQZ basis sets, and with the CBS-QB3 and CBS-APNO complete basis set methods. Seven different urea dimer structures were optimized at the MP2/aug-cc-pVDZ level to obtain accurate interaction energies. Atomic partial charges were calculated at the MP2/aug-cc-pVDZ level with the restrained electrostatic potential (RESP) fitting approach. The interaction energies computed with these new RESP charges in the force field are consistent with those obtained from CCSD and MP2 calculations. The linear dimer structure calculated using the force field with modified geometrical parameters and the new RESP charge set agrees well with available experimental data.     

Structure,energetics and properties of some molecules with potent anti-HIV activity: a theoretical study

Density functional theory (DFT; B3LYP) and Hartree–Fock (HF; 3-21G, 6-31G(d) and 6-311G(d,p)) calculations with complete geometry optimisations are carried out in the ground state on five 6-aminoquinolone derivatives, which have been proved to be highly effective in inhibiting HIV replication, to study their structures, energetics and HOMO–LUMO correlation with physiological action. The gas-phase calculations and single-point polarisable continuum model water-phase calculations show that the molecules are highly effective in inhibiting HIV replication, which is in excellent agreement with the experiment. Structural features, energies, charge densities and HOMO–LUMO correlation have been found to substantiate the experimental findings. Compound 4 (pyrazine) shows some special features in DFT calculations which are not found in HF calculations. In the present series, HF results are more reliable as expected.     

Theoretical Ab initio study of the effects of methylation on the nature of hydrogen bonding in A:T base pair

We report the results of a theoretical ab initio study of methylation in Watson-Crick A:T base pairs. Equilibrium geometries were obtained without symmetry restrictions by the gradient procedure at DFT level of theory with the standard 6-31G(d) basis set. Each local minima was verified by energy second derivative calculations. Single-point calculations for the DFT geometries have been performed at the MP2/6-31G(d,p), MP2/6-31++G(d,p), and MP2/6-311++G(2d,2p) levels of theory. The geometrical parameters, relative stabilities and counterpoise corrected interaction energies are reported. In addition, using a variation-perturbation energy decomposition scheme, we have found the important contributions to the total interaction energy.     

Combined procedure of distance geometry and restrained molecular dynamics techniques for protein structure determination from nuclear magnetic resonance data: application to the DNA binding domain of lac repressor from Escherichia coli

The technique of two-dimensional nuclear magnetic resonance (2D-NMR) has recently assumed an active role in obtaining information on structures of polypeptides, small proteins, sugars, and DNA fragments in solution. In order to generate spatial structures from the atom-atom distance information obtained by the NMR method, different procedures have been developed. Here we introduce a combined procedure of distance geometry (DG) and molecular dynamics (MD) calculations for generating 3D structures that are consistent with the NMR data set and have reasonable internal energies. We report the application of the combined procedure on the lac repressor DNA binding domain (headpiece) using a set of 169 NOE and 17 "hydrogen bond" distance constraints. Eight of ten structures generated by the distance geometry algorithm were refined within 10 ps MD simulation time to structures with low internal energies that satisfied the distance constraints. Although the combination of DG and MD was designed to combine the good sampling properties of the DG algorithm with an efficient method of lowering the internal energy of the molecule, we found that the MD algorithm contributes significantly to the sampling as well.     

A model for oxygen supply in fixed bed reactors with immobilized hybridoma cells

In fixed bed reactors with animal cells immobilized in macroporous carriers sufficient oxygen supply is a critical parameter. For modelling of the oxygen consumption and the oxygen profile in a fixed bed oxygen gradients within the porous carriers and along the length of the fixed bed have to be considered. For the complex geometry of the fixed bed a model structure was assumed, that allows the calculation of the oxygen profile. The model for oxygen supply of the immobilized cells included the transport resistance from the bulk fluid into the carriers and diffusion within the carriers. The model was compared with experimental data obtained with a hybridoma cell line for production of monoclonal antibodies. Model calculations and experimental data agree rather well. The mean volume-specific oxygen uptake rate as an indicator for the cell activity increased with the superficial flow velocity of the bulk liquid flow, and did not depend on the length of the fixed bed in the range tested. This indicates, that the convective transport from the bulk liquid flow between the carriers to the outer surface of the carriers is a dominating transport resistance besides the diffusive oxygen supply within the carriers.     

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

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

An evaluation of Poisson-Boltzmann electrostatic free energy calculations through comparison with experimental mutagenesis data

For systems involving highly and oppositely charged proteins, electrostatic forces dominate association and contribute to biomolecular complex stability. Using experimental or theoretical alanine-scanning mutagenesis, it is possible to elucidate the contribution of individual ionizable amino acids to protein association. We evaluated our electrostatic free energy calculations by comparing calculated and experimental data for alanine mutants of five protein complexes. We calculated Poisson-Boltzmann electrostatic free energies based on a thermodynamic cycle, which incorporates association in a reference (Coulombic) and solvated (solution) state, as well as solvation effects. We observe that Coulombic and solvation free energy values correlate with experimental data in highly and oppositely charged systems, but not in systems comprised of similarly charged proteins. We also observe that correlation between solution and experimental free energies is dependent on dielectric coefficient selection for the protein interior. Free energy correlations improve as protein dielectric coefficient increases, suggesting that the protein interior experiences moderate dielectric screening, despite being shielded from solvent. We propose that higher dielectric coefficients may be necessary to more accurately predict protein-protein association. Additionally, our data suggest that Coulombic potential calculations alone may be sufficient to predict relative binding of protein mutants.     

Effects of local protein environment on the binding of diatomic molecules to heme in myoglobins. DFT and dispersion-corrected DFT studies

The heme-AB binding energies (AB = CO, O₂) in a wild-type myoglobin (Mb) and two mutants (H64L, V68N) of Mb have been investigated in detail with both DFT and dispersion-corrected DFT methods, where H64L and V68N represent two different, opposite situations. Several dispersion correction approaches were tested in the calculations. The effects of the local protein environment were accounted for by including the five nearest surrounding residues in the calculated systems. The specific role of histidine-64 in the distal pocket was examined in more detail in this study than in other studies in the literature. Although the present calculated results do not change the previous conclusion that the hydrogen bonding by the distal histidine-64 residue plays a major role in the O₂/CO discrimination by Mb, more details about the interaction between the protein environment and the bound ligand have been revealed in this study by comparing the binding energies of AB to a porphyrin and the various myoglobins. The changes in the experimental binding energies from one system to another are well reproduced by the calculations. Without constraints on the residues in geometry optimization, the dispersion correction is necessary, since it improves the calculated structures and energetic results significantly.