The effect of solvent on the internal rotation of formamide: an ab initio study using a polarizable continuum model

The effect of solvent on the barrier to internal rotation of formamide has been studied using an ab initio 4-3IG method coupled with a polarizable continuum model. The trend of increasing rotation barrier with increasing solvent dielectric constant is reproduced. Conformational free energies were calculated by adding dispersion and cavitation terms to the electrostatic energy.     

Alkaline hydrolysis of ethylene phosphate: an ab initio study by supermolecule model and polarizable continuum approach

The alkaline hydrolysis reaction of ethylene phosphate (EP) has been investigated using a supermolecule model, in which several explicit water molecules are included. The structures and single-point energies for all of the stationary points are calculated in the gas phase and in solution at the B3LYP/6-31++G(df,p) and MP2/6-311++G(df,2p) levels. The effect of water bulk solvent is introduced by the polarizable continuum model (PCM). Water attack and hydroxide attack pathways are taken into account for the alkaline hydrolysis of EP. An associative mechanism is observed for both of the two pathways with a kinetically insignificant intermediate. The water attack pathway involves a water molecule attacking and a proton transfer from the attacking water to the hydroxide in the first step, followed by an endocyclic bond cleavage to the leaving group. While in the first step of the hydroxide attack pathway the nucleophile is the hydroxide anion. The calculated barriers in aqueous solution for the water attack and hydroxide attack pathways are all about 22 kcal/mol. The excellent agreement between the calculated and observed values demonstrates that both of the two pathways are possible for the alkaline hydrolysis of EP.     

DFT and ab initio quantum chemical studies on p-cyanobenzoic acid

The Fourier transform infrared (FTIR) and FT-Raman spectra of p-cyanobenzoic acid (CBA) have been recorded in the range 4000-400 and 4000-100 cm(-1), respectively. The complete vibrational assignment and analysis of the fundamental modes of the compound were carried out using the observed FTIR and FT-Raman data. The vibrational frequencies determined experimentally were compared with theoretical wavenumbers obtained from ab initio HF and DFT-B3LYP gradient calculations employing 6-31G**, 6-311++G** and cc-pVTZ basis sets for the optimised geometry of the compound. The geometry and normal modes of vibration obtained from the HF and DFT methods are in good agreement with the experimental data. The normal coordinate analysis was also carried out with ab initio force fields utilising Wilson's FG matrix method. The interactions of cyano and carboxylic acid groups with the skeletal vibrational modes were investigated.     

An ab initio quantum chemical investigation of 43Ca NMR interaction parameters for the Ca2+ sites in organic complexes and in metalloproteins

We have carried out an extensive ab initio quantum chemical (QC)43Ca NMR study on a series of Ca-O organic compounds and three different Ca-bound proteins and found that the HF/6-31G* level of function can reliably predict 43Ca NMR interaction parameters (delta(iso) and chi(q)), especially for organic solids. This QC study finds correlations between Ca-O bond environment (mean distance and coordination number) and delta(iso)(43Ca). Although relatively small values of chi(q)(43Ca) are found for Ca-O organic compounds with a coordination number between 6 and 10, the QC shows that chi(q)(43Ca) is sensitive to the Ca-O coordination geometry of the Ca2+ sites in metalloproteins--a potentially important observation. An application of such ab initio QC 43Ca NMR studies is in characterizing the Ca-O bonding environment around target Ca2+ sites. As an example, we propose a new potential analytical approach using the absolute (43)Ca chemical shielding constant to investigate the hydration shell of Ca2+ in a dilute CaCl2 aqueous solution. Furthermore, by adopting a NMR methodology similar to that reported in Wong et al. Chem. Phys. Lett. 2006, 427, 201, natural abundance 43Ca MAS NMR spectra of Ca(L-glutamate)(2) x 4H2O were recorded, and delta(iso)(43Ca) and the quadrupolar parameter (Pq) were estimated to be 6.6 ppm and 0.8 MHz, respectively.     

Development of a polarizable force field for proteins via ab initio quantum chemistry: first generation model and gas phase tests

We present results of developing a methodology suitable for producing molecular mechanics force fields with explicit treatment of electrostatic polarization for proteins and other molecular system of biological interest. The technique allows simulation of realistic-size systems. Employing high-level ab initio data as a target for fitting allows us to avoid the problem of the lack of detailed experimental data. Using the fast and reliable quantum mechanical methods supplies robust fitting data for the resulting parameter sets. As a result, gas-phase many-body effects for dipeptides are captured within the average RMSD of 0.22 kcal/mol from their ab initio values, and conformational energies for the di- and tetrapeptides are reproduced within the average RMSD of 0.43 kcal/mol from their quantum mechanical counterparts. The latter is achieved in part because of application of a novel torsional fitting technique recently developed in our group, which has already been used to greatly improve accuracy of the peptide conformational equilibrium prediction with the OPLS-AA force field.1 Finally, we have employed the newly developed first-generation model in computing gas-phase conformations of real proteins, as well as in molecular dynamics studies of the systems. The results show that, although the overall accuracy is no better than what can be achieved with a fixed-charges model, the methodology produces robust results, permits reasonably low computational cost, and avoids other computational problems typical for polarizable force fields. It can be considered as a solid basis for building a more accurate and complete second-generation model.     

A new extension of the polarizable continuum model: Toward a quantum chemical description of chemical reactions at extreme high pressure

A quantum chemical method for studying potential energy surfaces of reactive molecular systems at extreme high pressures is presented. The method is an extension of the standard Polarizable Continuum Model that is usually used for Quantum Chemical study of chemical reactions at a standard condition of pressure. The physical basis of the method and the corresponding computational protocol are described in necessary detail, and an application of the method to the dimerization of cyclopentadiene (up to 20 GPa) is reported. © 2015 Wiley Periodicals, Inc.     

Intramolecular interaction energies in model alanine and glycine tetrapeptides. Evaluation of anisotropy, polarization, and correlation effects. A parallel ab initio HF/MP2, DFT, and polarizable molecular mechanics study

An extension of the SIBFA polarizable molecular mechanics procedure to flexible oligopeptides is reported. The procedure is evaluated by computing the relative conformational energies, deltaE(conf), of the alanine tetrapeptide in 10 representative conformations, which were originally derived by Beachy et al. (J Am Chem Soc 1997, 119, 5908) to benchmark molecular mechanics procedures with respect to ab initio computations. In the present study, a particular emphasis is on the separable nature of the components of the energy and the particular impact of the polarization energy component on deltaE(conf). We perform comparisons with respect to single-point HF, DFT, LMP2, and MP2 computations done at the SIBFA-derived energy minima. Such comparisons are made first for the 10 conformers derived from phi/psi torsional angle energy-minimization (the rigid rotor approach), and, in a second step, after allowing additional relaxation of the C(alpha) centered valence angles. In both series of energy-minimization, the SIBFA deltaE(conf) compared best with the LMP2 results using the 6-311G** basis set, the rms being 1.3 kcal/mol. In the absence of the polarization component, the rms is 3.5 kcal/mol. In both series of minimizations, the magnitudes of deltaE(conf), computed as differences with respect to the most stable conformer taken as energy zero, decrease along the series: HF > DFT > LMP2 > SIBFA > MP2, indicative of increasing stabilization of the most highly folded conformers.     

Correlated ab initio quantum chemical calculations of di- and trisaccharide conformations

High level correlated quantum chemical calculations, using MP2 and local MP2 theory, have been performed for conformations of the disaccharide, beta-maltose, and the trisaccharide, 3,6-di-O-(alpha-D-mannopyranosyl)-alpha-D-mannopyranose. For beta-maltose, MP2 and local MP2 calculations using the 6-311++G** basis set are in good agreement, predicting a global minimum gas-phase conformation with a counterclockwise hydrogen bond network and the experimentally-observed intersaccharide hydrogen bonding arrangement. For conformations of 3,6-di-O-(alpha-D-mannopyranosyl)-alpha-D-mannopyranose, MP2/6-311++G**, and local MP2/6-311++G** calculations do not provide a consensus prediction of relative energetics, with the MP2 method finding large differences in stability between extended and folded trisaccharide conformations. Local MP2 calculations, less susceptible to intramolecular basis set superposition errors, predict a narrower range of trisaccharide energetics, in line with estimates from Hartree-Fock theory and B3LYP and BP86 density functionals. All levels of theory predict compact, highly hydrogen-bonded conformations as lowest in energy on the in vacuo potential energy surface of the trisaccharide. These high level, correlated local MP2/6-311++G** calculations of di- and trisaccharide energetics constitute potential reference data in the development and testing of improved empirical and semiempirical potentials for modeling of carbohydrates in the condensed phase.     

An interaction potential between an alanine zwitterion and a water molecule based on ab initio calculations

Interaction energies between an alanine zwitterion and a water molecule at 150 different positions and orientations have been calculated using the ab initio method with the minimal basis set and employing the counterpoise method to eliminate the basis set superposition error. Dispersion energies are estimated using the Slater–Kirkwood formula. Out of a total of 150 computed interaction energies, 140 whose SCF interaction energies are below 5 kcal/mol have been fitted with a summation of atom-atom pair potentials in the form of the Lennard–Jones potential plus an electrostatic term. The standard deviation for this fitting is 0.49 kcal/mol. A sampling scheme regarding geometrical configurations is presented. Twenty rays are uniformly drawn from the origin of coordinates, a floatable division with equal ratios is made along each ray, and one of 60 orientations is randomly taken as the orientation of a water molecule. A nonlinear fitting method is used with a restriction on the sign change of fitting coefficients.     

Gold adatoms and clusters on PPV: An ab initio investigation

We have performed an ab initio investigation of the energetic, structural, electronic, and vibrational properties of Au atoms and clusters adsorbed on poly-p-phenylene vinylene (PPV) chains, Au(n)/PPV (with n =?1, 2, 6, 7, 10, and 12). We find that the Au(n)/PPV systems are energetically stable by 0.5 eV, compared with the isolated systems, viz., PPV chain and Au(n) clusters, thus supporting the formation of Au(n)/PPV nanocomposites. Further support to the formation of Au(n)/PPV has been provided by examining the vibrational properties of pristine PPV and Au(n)/PPV systems. In agreement with experimental measurements, we find a reduction on the in-plane vibrational frequency of C-C bonds of Au(n)/PPV, when compared with the same vibrational modes of pristine PPV. The electronic properties of isolated Au(n) clusters are modified when adsorbed on PPV. The highest occupied states of Au(n)/PPV are mostly concentrated on the Au(n) cluster, while the lowest unoccupied states are mainly localized along the PPV chain. The HOMO-LUMO energy gap of the Au(n)/PPV systems are smaller than the energy gap of the isolated systems, Au(n) clusters, and pristime PPV chains.     

An ab initio investigation of 2-amino-2-imidazoline: a key moiety in chemical and biochemical processes

The 2-amino-2-imidazoline moiety is currently used not only in drugs, but also in insecticides, and fungicides. Ab initio calculations are performed to evaluate the molecular properties of the two tautomeric forms and the protonated form with extended basis sets ranging from 6-31G* to 6-311++G** at Hartree-Fock and density functional (BLYP and B3LYP) levels. M?ller-Plesset perturbation is tested at the MP2/6-31G* level only. Optimized geometry structures, energies and thermochemical properties are generated. Basis set and correlation effects on geometries, tautomer equilibrium constant and protonation enthalpy are carefully analysed. Although observed for the isolated molecule, these results may be extrapolated to chemical and biochemical systems of interest.     

Polarization energy gradients in combined quantum mechanics, effective fragment potential, and polarizable continuum model calculations

A method that combines quantum mechanics (QM), typically a solute, the effective fragment potential (EFP) discrete solvent model, and the polarizable continuum model is described. The EFP induced dipoles and polarizable continuum model (PCM) induced surface charges are determined in a self-consistent fashion. The gradients of these two energies with respect to molecular coordinate changes are derived and implemented. In general, the gradients can be formulated as simple electrostatic forces and torques among the QM nuclei, electrons, EFP static multipoles, induced dipoles, and PCM induced charges. Molecular geometry optimizations can be performed efficiently with these gradients. The formulas derived for EFPPCM can be generally applied to other combined molecular mechanics and continuum methods that employ induced dipoles and charges.     

Hydrolysis of cyclic phosphates by ribonuclease A: a computational study using a simplified ab initio quantum model

The second step in the enzyme-catalyzed hydrolysis of phosphate esters by ribonuclease A (RNase A) was studied using an ab initio quantum-based model of the active site including constrained parts of three critical residues, His-12, His-119, and Lys-41, and a small substrate. The competition between release of the cyclic phosphate intermediate and subsequent hydrolysis following transphosphorylation was explored to determine the electronic factors that contribute to preferential intermediate product release observed experimentally. The structural and energetic results obtained at both the RHF and MP2 levels reveal several contributing factors consistent with experimental observation. Although the intrinsic electronic effects tend to favor hydrolysis slightly with an overall activation free energy of approximately 70 kJ mol(-1), entropic and environmental effects favor release of the cyclic phosphate intermediate over hydrolysis. Exploration of the second, hydrolysis step also revealed interesting similarity with the transphosphorylation step, including the observation of autocatalysis by the substrate. Moreover, both steps of the overall RNase A reaction reveal multiple pathways involving proton transfers to sites of similar proton affinities. The anionic phosphate in both steps can act as a stable proton binding site as protons are moved around the active site throughout the progress of the reaction. These results suggest autocatalysis may be representative of more general behavior in enzymes containing highly charged substrates, especially phosphates.     

Vibrational spectroscopy investigation using ab initio and density functional theory on p-anisaldehyde

The FTIR and FT Raman spectra of p-anisaldehyde has been recorded in the regions 4,000-400 and 3,500-100 cm(-1), respectively. The optimized geometry, frequency and intensity of the vibrational bands of p-anisaldehyde were obtained by ab initio and DFT levels of theory with complete relaxation in the potential energy surface using 6-31G(d,p) basis set. A complete vibrational assignment aided by the theoretical harmonic frequency analysis has been proposed. The harmonic vibrational frequencies calculated have been compared with experimental FTIR and FT Raman spectra. The observed and the calculated frequencies are found to be in good agreement. The experimental spectra also coincide satisfactorily with those of theoretically constructed bar type spectrograms.     

Structure of isolated tryptophyl-glycine dipeptide and tryptophyl-glycyl-glycine tripeptide: ab initio SCC-DFTB-D molecular dynamics simulations and high-level correlated ab initio quantum chemical calculations

The tryptophyl-glycine (Trp-Gly) and tryptophyl-glycyl-glycine (Trp-Gly-Gly) peptides have been studied by means of molecular dynamic simulations combined with high-level correlated ab initio quantum chemical and statistical thermodynamic calculations. The lowest energy conformers were localized in the free energy surface. The structures of the different Trp-Gly and Trp-Gly-Gly conformers coexisting in the gas phase have been for the first time reported and their scaled theoretical IR spectra unambiguously assigned and compared with previous gas-phase experimental results. Common geometrical features have been systematically observed for the sequence Trp, Trp-Gly, and Trp-Gly-Gly. In addition, the peptide backbone of Trp-Gly-Gly has been compared with that of the previously studied Phe-Gly-Gly (Reha, D. et. al. Chem. Eur. J. 2005, 11, 6803). From the observed systematic structural behavior between these peptide analogues, it is expected that the gas-phase conformers of other similar aromatic small peptides would present equivalent geometries. The DFT methodology failed to describe the potential energy surface of the studied peptides since the London dispersion energy (not covered in DFT) plays a significant role in the stabilization of most stable conformers.     

Studies of the rates of thermal decomposition of glycine,alanine, and serine

Rates of thermal decomposition of glycine, alanine, and serine are described by the equation of first order reaction in the temperature range 200–300°C. Apparent rate constants and apparent activation energies of decomposition of α-amino acids were evaluated. It was found that the main gaseos reaction product is carbon dioxide.     

MgH^+,Mg2H^+量子化学从头计算的研究

用量子化学自洽场分子轨道从头计算方法,采用STO-3G基组计算了MgH~+,Mg_2H~+离子的位能曲线与位能面;给出电子波函数与电子集居数;得到MgH~+的平衡构型键长γ_(Mg-H)1.60A.Mg_2H~+的平衡构型有C_(∞v)与D_(∞h)两类,前者为[Mg—Mg—H]~+,键长γ_(Mg-Mg)2.41,γ_(Mg-H)1.63A;后者为[Mg—H—Mg]~+,γ_(Mg-H)1.73A.前者的总能量比后者低50kcal/mol.还讨论了它们的稳定性.MgH~+的键长计算结果和实测值较吻合,Mg_2H~+的平衡构型计算结果支持了Porter从热力学研究提出的假设.     

The ab initio quantum chemical calculation of MgH+ and Mg2H+

The potential energy curve of MgH+ and the potential energy surface of Mg₂H+ are calculated by quantum chemical ab initio SCF MO method with STO-3G basis set. The electronic wave functions and populations are obtained. The eqilibrium internuclear distance of MgH+ is 1.60 Å. There are two possible configurations of Mg₂H+: C and D . The former is (Mg-Mg-H)+, with bond length rM_g-M_g = 2.41Å and rM_g-H =1.63Å. The latter is (Mg-H-Mg)+, with bond length rM_g-H = 1.73Å. The cause of the stability of these species is discussed. The result of calculation about the bond length of MgH+ agrees reasonably well with experimental value. The conclusion about the stability of (Mg-Mg-H)+ ion supports the hypothesis proposed by Porter, based on thermodynamical calculation.     

Ab initio quantum chemical investigation of the first steps of the photocycle of phototropin: a model study

Phototropin is a blue light-activated photoreceptor that plays a dominant role in the phototropism of plants. The protein contains two subunits that bind flavin mononucleotide (FMN), which are responsible for the initial steps of the light-induced reaction. It has been proposed that the photoexcited flavin molecule adds a cysteine residue of the protein backbone, thus activating autophosphorylation of the enzyme. In this study, the electronic properties of several FMN-related compounds in different charge and spin states are characterized by means of ab initio quantum mechanical calculations. The model compounds serve as idealized model chromophores for phototropism. Reaction energies are estimated for simple model reactions, roughly representing the addition of a cysteine residue to the flavin molecule. Excitation energies were calculated with the help of time-dependent density functional theory. On the basis of these calculations we propose the following mechanism for the addition reaction: (1) after photoexcitation of FMN out of the singlet ground state S0, excited singlet state(s) are populated; these relax to the lowest excited singlet state S1, and subsequently by intersystem crossing FMN in the lowest triplet state, T1 is formed; (2) the triplet easily removes the neutral hydrogen atom from the H-S group of the cysteine residue; and (3) the resulting thio radical is added.