Ab initio self-consistent field study of deoxyguanosine

The rotational energy of deoxyguanosine around the C₁–N₉ bond is calculated using the Hartree-Fock method with an STO-3G basis set.     

Ab initio calculations on the thermodynamic properties of azaborospiropentanes

Following our recent studies of the thermodynamic properties of azaspiropentane and borospiropentane, in consideration of their usefulness as new potential high energy materials, we follow up with ab initio calculations on the thermodynamic properties of azaborospiropentanes. Properties reported in this study include optimized structural parameters, vibrational frequencies, enthalpies of formation, specific enthalpies of combustion, proton affinities, and hydride affinities. Our results indicate that azatriborospiropentane gives off most energy when combusted, as evidenced by its specific enthalpy of combustion of about −52 kJ per gram. Figure Optimized geometry for R-azatriborospiropentane (10) Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Ab initio molecular orbital studies of closed shell flavins

Ab initio calculations with a (7 s 3 p) basis set are performed on uracil, lumazine, alloxazine and various isoalloxazines. The results as total energies and charge distributions are discussed in relation to the biochemical behaviour of the flavins. The calculations correctly predict equilibrium situations in the alloxazine-isoalloxazine system and explain the high affinity for nucleophilic addition at N₅ in the flavins. The reduction of flavins and their reoxidation by oxygen are discussed.     

Ab initio calculations of the pyrophosphate hydrolysis reaction.

Ab initio quantum mechanical calculations were used to study the hydrolysis reaction H4P2O7 + H2O in equilibrium with 2H3PO4, as well as some molecular properties of the reactants and products. SCF calculations with several basis sets ranging from minimal to extended with polarization functions were used to look at the basis dependency of the reaction enthalpies and optimized geometries. Although the minimal basis sets yield erratic predictions of the enthalpy, when a more extended basis (3-21G*) was used for the geometry optimization, and the total energies of the reactants and products were computed with this and larger basis sets, we obtained more consistent predictions of the structural properties of the P-O-P bridge and of the heat of the hydrolysis reaction (delta E = -7.39 kcal/mol at the SCF/6-31G** level). A comparison is made with previous estimates performed with smaller basis sets and without taking into account the electron correlation effects, which are calculated in the present work. The inclusion of the zero point energy calculated using the harmonic approximation, and of the electronic correlation energy determined at the MBPT(2) level, raised the computed heat of the reaction to -3.83 kcal/mol, and when an estimate for the thermal energy was added, the value obtained was of -3.38 kcal/mol. In conclusion, we found that the hydrolysis of pyrophosphate should be exothermic in the gas phase. The implications of this result in relation to some recent theories about enzyme catalysis are discussed.     

Ab initio calculations of electron distributions in heme-CO models

We have, by the use of ab initio calculations, found a back-bonding state of pi symmetry close to the Fermi level for CO bound to FeN5C14. We thus find it likely that small shifts of the redox potential magnitude of EF - EV magnitude of will cause relatively large changes of the CO vibrational frequency. The separation of Fe 3d orbitals in our heme model is found to agree with what is predicted by ligand field theory for Oh symmetry. This paper presents nonrelativistic Hartree-Fock-Slater calculations of the 5 sigma bonding and 2 pi back-bonding between CO and Fe. The effects of up to 19 additional atoms are discussed for models of heme (COFe to COFeN5C14). The filled back-bonding state is found to be strongly influenced by second nearest neighbor atoms. By use of symmetry orbitals we have resolved the Fe 3d orbitals into the T2g and Eg representations of the Oh point group and find the former states to be occupied whereas the latter are unoccupied. The difference in occupancy is reduced when the CO ligand is removed which also causes an increased density of states at the Fermi level, i.e., the highest occupied and lowest unoccupied orbitals. Possible correlations between our data and experimental results are discussed for heme proteins as well as for metal surfaces.     

Ab initio molecular orbital calculations on furanose sugars: a study with the 6-31G basis set

Ab initio molecular orbital calculations were performed on 2-deoxy-beta-D-glycero-tetrofuranose (1) using the 6-31G* basis set to evaluate the effect of ring conformation on the molecular parameters (bond lengths, angles, and torsions). Geometric optimizations were conducted on the planar and ten envelope conformers of 1, and these data were compared to those obtained from previous calculations using the STO-3G and 3-21G basis sets. Conformational energy profiles derived from 3-21G and 6-31G* data were found to be qualitatively comparable. The effect of furanose ring conformation on key bond lengths (e.g., C-H, C-O), bond angles (e.g., COC), and bond torsions (e.g., the exoanomeric C-1-O-1 torsion) was examined, and a qualitative agreement was observed between the 3-21G and 6-31G* analyses. The results indicate that, for semi-quantitative ab initio studies of intact carbohydrates, the 3-21G basis set is sufficient, and that the STO-3G basis set should not be employed unless crude structural approximations are desired. The observed concerted behavior of C-O bond lengths in the vicinity of the anomeric carbon of the aldofuranose ring has suggested a possible role of C-1-O-1 bond orientation in affecting the mechanism of glycoside bond hydrolysis.     

Ab initio molecular dynamics study of the hydration of the formohydroxamate anion

We apply ab initio molecular dynamics (AIMD) to study the hydration structures and electronic properties of the formohydroxamate anion in liquid water. We consider the cis- nitrogen-deprotonated, cis- oxygen-deprotonated, and trans- oxygen-deprotonated formohydroxamate tautomers. They form an average of 6.3, 6.9, and 6.0 hydrogen bonds with water molecules, respectively. The predicted pair correlation functions and time dependence of the hydration numbers suggest that water is highly structured around the nominally negatively charged oxime oxygen in O-deprotonated tautomers but significantly less so around the nitrogen atom in the N-deprotonated species. Wannier function analysis suggests that, in the O-deprotonated anions, the negative charge is concentrated on the oxime oxygen, while in the N-deprotonated case, it is partially delocalized between the nitrogen and the adjoining oxime oxygen atom.     

Ab initio calculations on hidden modulators of theta class glutathione transferase activity

The glutathione transferases decrease the pKa of glutathione, allowing its deprotonation and the formation of the more reactive thiolate anion. The thiolate is maintained in the active site through a weak conventional hydrogen bond first sphere interaction donated by a Tyr hydroxyl in the Alpha, Mu, Pi, and Sigma glutathione transferase classes that can be modified by other second sphere or indirect thiolate contacts. However, the Theta and Delta class isoforms use a Ser hydroxyl for stabilizing the GSH thiolate, and as such, have a different chemical system compared with that of the Tyr possessed by other classes. We have used high level ab initio methods to investigate this interaction by using a simple methanol methanethiol system as a model. The hydrogen bond strength of this initial first sphere interaction was calculated to be less than that of the Tyr interaction. A putative second sphere interaction exists in the Theta and Delta class structures between Cys or Ser-14 and Ser-11 in the mammalian Theta subclass 1 and 2, respectively. The effect of this interaction on the first sphere interaction has also been investigated and found to significantly increase the energy of the bond.     

Effects of non-aqueous solvents on CO2 absorption in monoethanolamine: Ab initio calculations

Monoethanolamine (MEA) is the most typical alkanolamine and its aqueous solutions are widely used for CO₂ absorption with mature technology, but the regeneration process is energy consuming. To reduce the energy demand, non-aqueous solvents, such as methanol and ethanol are proposed to substitute water in amine solutions. To understand the influence of the aqueous and non-aqueous solvents on CO₂ capture process, the chemical reactions of MEA absorbing CO₂ were conducted via ab initio calculations. The non-aqueous solvents discussed in this paper are methanol, ethanol, 1-propanol and 2-propanol. The reaction patterns were investigated and energy barriers were observed. The results show that zwitterion formation and the followed intermolecular hydrogen transfer are proven to be the most possible reaction pattern in both aqueous and non-aqueous solvents. The energy analysis shows that the forward reaction energy barriers increase while the backward barriers decrease as the solvent changes from water to methanol, ethanol, 1-propanol and 2-propanol in turn. The decreases of the energy barriers for backward processes are much higher than the corresponding increases for forward processes. These results indicate that lower energies are required in non-aqueous solvents than in water during the desorption reactions and the non-aqueous solvents are very promising to reduce the regeneration energy consumption in MEA capturing CO₂ process. Moreover, the reaction energy gaps between different solvation effects were found to have linear relationship with the logarithm of the dielectric constant difference, which could provide an easy way to theoretically predict the reaction energies of monoethanolamine absorbing CO₂ in other solvation effect and can be used to screen appropriate CO₂ capture solvent.     

Ab initio vibrational calculations on ara-T molecule: application to analysis of IR and Raman spectra

The FTIR and FT-Raman spectra are reported for the arabinonucleoside ara-T (1-beta-D-arabinofuranosylthymine), which shows antiviral activity. The accurate knowledge of the vibrational modes is a prerequisite for the elucidation of drug-nucleotide and drug-enzyme interactions. The FTIR and FT-Raman spectra of ara-T were recorded from 4000 to 30 cm(-1). A tetradeuterated derivative (deuteration at N3, and hydroxyl groups O'2, O'3, and O'5) was synthesized and the observed isotopic shifts in its spectra were used for the vibrational analysis of ara-T. The theoretical frequencies and the potential energy distribution (PED) of the vibrational modes of ara-T were calculated using the ab initio Hartree-Fock/3-21G method. An assignment of the vibrational spectra of ara-T is proposed considering the scaled PED and the observed band shifts under deuteration. The scaled ab initio frequencies were in reasonable agreement with the experimental data.     

Ab initio studies of the properties of intracellular thiols ergothioneine and ovothiol

Intracellular naturally occurring aromatic thiols such as ergothioneine and the ovothiols have been shown to play a variety of roles in cellular function. A detailed ab initio electronic structure analysis of these thiols is reported evaluating the thermodynamics of the reactions of these intracellular thiols with alkyl thiols, HO*, H2O2, ascorbate and their disulfides.     

Ab initio molecular orbital calculations on the conformations of some unsaturated fragments of lipid molecules: 1-Butene, cis-2-pentene and 1,4-pentadiene

Conformational energies have been calculated for the title compounds which are basic fragments of physiologically active polyunsaturated fatty acids. For 1-butene extremely good accordance has been obtained with experimental data. In cis-2-pentene a strong influence of the cis methyl group on the potential energies for internal rotation has been found. For 1,4-pentadiene no great influence has been found from a possible coupling of the two rotors on the relative stabilities of the two stable conformations. Reasonable agreement has been observed with similar calculations published recently.     

Ab initio protein structure assembly using continuous structure fragments and optimized knowledge-based force field

Ab initio protein folding is one of the major unsolved problems in computational biology owing to the difficulties in force field design and conformational search. We developed a novel program, QUARK, for template-free protein structure prediction. Query sequences are first broken into fragments of 1-20 residues where multiple fragment structures are retrieved at each position from unrelated experimental structures. Full-length structure models are then assembled from fragments using replica-exchange Monte Carlo simulations, which are guided by a composite knowledge-based force field. A number of novel energy terms and Monte Carlo movements are introduced and the particular contributions to enhancing the efficiency of both force field and search engine are analyzed in detail. QUARK prediction procedure is depicted and tested on the structure modeling of 145 nonhomologous proteins. Although no global templates are used and all fragments from experimental structures with template modeling score >0.5 are excluded, QUARK can successfully construct 3D models of correct folds in one-third cases of short proteins up to 100 residues. In the ninth community-wide Critical Assessment of protein Structure Prediction experiment, QUARK server outperformed the second and third best servers by 18 and 47% based on the cumulative Z-score of global distance test-total scores in the FM category. Although ab initio protein folding remains a significant challenge, these data demonstrate new progress toward the solution of the most important problem in the field.     

Ab initio DFT study of structural and mechanical properties of methane and carbon dioxide hydrates

The structural and mechanical properties of methane and carbon dioxide hydrates were investigated using density functional theory simulations. Well-established equations of state of solids and exchange-correlation functionals were used for fitting the unit lattice total energy as a function of volume, and the full second-order elastic constants of these two gas hydrates were determined by energy–strain analyses. The polycrystalline elastic properties were also calculated from the unit lattice results. The final results for methane hydrate agree well with available experimental data and with other theoretical results. The two gas hydrates were found to be highly elastically isotropic, but they differed significantly in shear properties. The presented results for carbon dioxide hydrates are the first complete set reported so far. The results are a significant contribution to the ab initio material characterisation of gas hydrates required for ongoing fundamental studies and technological applications.     

Ab initio molecular orbital study of the mispairing ability of a nucleotide base analogue, N4-aminocytosine

The intrinsic properties of N4-aminocytosine, a base analogue of cytosine, are analyzed by an ab initio molecular orbital method. Relative stabilities of four possible isomeric structures of N4-aminocytosine are shown. The more stable isomer has the smaller dipole moment, so the relative stabilities of the isomers in solutions are subject to solvent polarity. The mutagenicity of this base analogue must arise because it can behave like either cytosine or thymine. It can form a guanine-cytosine-like base pair more easily than cytosine, and an adenine-thymine-like base pair less easily than thymine.     

Ab initio parameterization of YFF1, a universal force field for drug-design applications

The YFF1 is a new universal molecular mechanic force field designed for drug discovery purposes. The electrostatic part of YFF1 has already been parameterized to reproduce ab initio calculated dipole and quadrupole moments. Now we report a parameterization of the van der Waals interactions (vdW) for the same atom types that were previously defined. The 6–12 Lennard-Jones potential terms were parameterized against homodimerization energies calculated at the MP2/6-31 G* level of theory. The Boys-Bernardi counterpoise correction was employed to account for the basis-set superposition error. As a source of structural information we used about 2,400 neutral compounds from the ZINC2007 database. About 6,600 homodimeric configurations were generated from this dataset. A special “closure” procedure was designed to accelerate the parameters fitting. As a result, dimerization energies of small organic compounds are reproduced with an average unsigned error of 1.1 kcal mol^-1. Although the primary goal of this work was to parameterize nonbonded interactions, bonded parameters were also derived, by fitting to PM6 semiempirically optimized geometries of approximately 20,000 compounds.     

Ab initio structure prediction of two alpha-helical oligomers with a multiple-chain united-residue force field and global search

A hierarchical methodology for ab initio structure prediction is extended to treat oligomeric proteins. Modifications are made to a united-residue (UNRES) force field and a Conformational Space Annealing (CSA) global search method. The computational cost of including additional chains and the increase in speed from symmetry optimizations are evaluated. The native structures of two oligomeric proteins from the CASP3 exercise, the retro-GCN4 leucine zipper and the synthetic domain-swapped dimer, were identified as the lowest-energy families resulting from the search of the proteins when rotational symmetry was imposed. Additional searches in different symmetries and oligomerization states were carried out, and the results indicate some problems in the thoroughness of the search and in the search of packing arrangements if symmetry constraints are not imposed.     

Ab initio predictions of structural and elastic properties of struvite: contribution to urinary stone research

In the present work, we carried out density functional calculations of struvite – the main component of the so-called infectious urinary stones – to study its structural and elastic properties. Using a local density approximation and a generalised gradient approximation, we calculated the equilibrium structural parameters and elastic constants C _ijkl . At present, there is no experimental data for these elastic constants C _ijkl for comparison. Besides the elastic constants, we also present the calculated macroscopic mechanical parameters, namely the bulk modulus (K), the shear modulus (G) and Young's modulus (E). The values of these moduli are found to be in good agreement with available experimental data. Our results imply that the mechanical stability of struvite is limited by the shear modulus, G. The study also explores the energy-band structure to understand the obtained values of the elastic constants.     

Ab initio predictions of structural and elastic properties of struvite: contribution to urinary stone research

In the present work, we carried out density functional calculations of struvite--the main component of the so-called infectious urinary stones--to study its structural and elastic properties. Using a local density approximation and a generalised gradient approximation, we calculated the equilibrium structural parameters and elastic constants C(ijkl). At present, there is no experimental data for these elastic constants C (ijkl) for comparison. Besides the elastic constants, we also present the calculated macroscopic mechanical parameters, namely the bulk modulus (K), the shear modulus (G) and Young's modulus (E). The values of these moduli are found to be in good agreement with available experimental data. Our results imply that the mechanical stability of struvite is limited by the shear modulus, G. The study also explores the energy-band structure to understand the obtained values of the elastic constants.