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 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 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.     

Stereochemical reaction mechanism formulations for enzyme-catalyzed pyrophosphate hydrolysis,ATP hydrolysis,and ATP synthesis

Fundamental concepts pertaining to the stereochemistry paths of polar additionelimination (nucleophilic substitution) reactions at phosphate phosphorus centers are reviewed and employed to analyze ¹⁸O exchange reactions catalyzed by inorganic pyrophosphatase and mitochondrial ATP synthetase. The analysis suggests reasonable choices for the stereochemistry path of the ¹⁸O exchanges. This, in turn, permits reasonable choices for the stereochemistry paths of hydrolysis of pyrophosphate catalyzed by pyrophosphatase and of hydrolysis and synthesis of ATP catalyzed by ATP synthetase.     

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 prediction of thermodynamically feasible reaction directions from biochemical network stoichiometry

Analysis of the stoichiometric structure of metabolic networks provides insights into the relationships between structure, function, and regulation of metabolic systems. Based on knowledge of only reaction stoichiometry, certain aspects of network functionality and robustness can be predicted. Current theories focus on breaking a metabolic network down into non-decomposable pathways able to operate in steady state. The physics underlying these theories is based on mass balance and the laws of thermodynamics. However, due to the inherent nonlinearity of the thermodynamic constraints on metabolic fluxes, computational analysis of large-scale biochemical systems can be expensive. In this study, it is shown how the feasible reaction directions may be determined by either computing the allowable ranges under the mass-balance and thermodynamic constraints or by analyzing the stoichiometric structure of the network. The computed reaction directions translate into a set of linear constraints necessary for thermodynamic feasibility. This set of necessary linear constraints is shown to be sufficient to guarantee feasibility in certain cases, thus translating the nonlinear thermodynamic constraints to linear. We show that for a reaction network of 44 internal reactions representing energy metabolism, the computed linear inequality constraints represent necessary and sufficient conditions for thermodynamic feasibility.     

Ab initio studies of the reaction of hydrogen transfer from DNA to the calicheamicinone diradical

BACKGROUND: The biological activity of enediyne chemotherapeutic (anti-cancer) agents is attributed to their ability to cleave duplex DNA. Part of the reaction of cleavage is the abstraction of hydrogens from the deoxyribose moiety of DNA by the biradical formed via a Bergman rearrangement. METHODS: The mechanism of the reaction of abstraction of two hydrogen atoms from two deoxyribophosphate molecules by the calicheamicinone biradical is studied with ab initio calculations at Hartree-Fock and post-Hartree-Fock level. The Titan program is used to perform the calculations. RESULTS: It is found that the reactions are exothermic and thus thermodynamically reasonable. CONCLUSIONS: The mechanism of DNA cleavage by the enediyne-containing drugs is likely to proceed by the abstraction of the hydrogens from deoxyribose by the biradical formed by the drug. Further studies should determine in which way the modification of the drug's structure would make this reaction even more exothermic and, thus, more likely to occur.     

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 multireference study of Hetero-Diels-Alder reaction of buta-1,3-diene with alkyl glyoxylates

Hetero-Diels-Alder (HDA) reaction of methyl glyoxylate with buta-1,3-diene has been investigated using multireference methods (complete active space SCF and multi-reference perturbation theory) and compared with several single-reference methods (including DFT) often used in calculations of catalysed [4+2] cycloadditions. Concerted and stepwise mechanisms, found in the literature, are compared. It is shown, that the stepwise mechanism may be a result of choosing unbalanced active space. Such choice leads to very close singlet and triplet states in the intermediate geometry - an artificial effect, that disappears if properly balanced active space is used (here, we use active space of 12 orbitals and 12 electrons). Conclusions concerning the mechanism and usefulness of the applied methodology are drawn, which might be important for theoretical investigation of stereoselectivity and specificity of catalysts for the HDA reaction. Figure Hetero-Diels-Alder reaction of alkyl glyoxylates with buta-1,3-diene, investigated using multi- and single-reference ab initio methods Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Ab initio reconstruction of metabolic pathways

We propose a new formulation for the problem of ab initio metabolic pathway reconstruction. Given a set of biochemical reactions together with their substrates and products, we consider the reactions as transfers of atoms between the chemical compounds and we look for successions of reactions transferring a maximal (or preset) number of atoms between a given source and sink compound. We state this problem as the one of finding a composition of partial injections that maximizes the image size. First, we study the theoretical complexity of this problem, state some related problems and then give a practical algorithm to solve them. Finally, we present two applications of this approach to the reconstruction of the tryptophan biosynthesis pathway and to the glycolysis.     

Ab initio investigations on the intermediates of RNA cleaving processes.

Ab initio molecular orbital calculations have been carried out on acyclic oxyphosphoranes in order to elucidate the origin of previously calculated energy differences of cyclic counterparts of RNA cleaving reactions. It appears that the orientation of the equatorial methoxyl group has pronounced effects in energy on the transition states than the metastable intermediates.     

Ab initio self-consistent field and potential-dependent partial equalization of orbital electronegativity calculations of hydration properties of N-acetyl-N'-methyl-alanineamide

Using a recently developed parallel computation algorithm, ab initio self-consistent field (SCF) calculations were carried out to estimate the relative hydration energies for 12 low-energy conformations of N-acetyl-N'-methyl-alanineamide. The requisite SCF calculations were carried out using 6-31G and 6-31G* basis sets, both in the absence and presence of a perturbing potential arising from a model solvent. The alpha R, alpha L, polyproline II (PII), and pi helical conformations were preferentially stabilized by the solvent potential, whereas conformations with intramolecular hydrogen-bonding C5 and C7 were preferred in the gas phase. Average vicinal nmr coupling constants (JNH-C alpha H), calculated using the total energies of the various solvated conformations, were consistent with observed coupling constants for this peptide in aqueous solution. Substantial alteration of the solute charge density occurred upon equilibration with the reaction field, as was exemplified in changes both in the molecular dipole moments and in atom-centered multipoles, when the molecule was transferred from a medium of low dielectric constant to one of high dielectric constant. In order to model these changes in charge density with an empirical scheme, we have implemented a novel monopolar representation of the solute charge density based on a potential-dependent form of partial equalization of orbital electronegativities (PDPEOE). In the atom-centered point charge PDPEOE representation, charge flows from one region of the solute to another in response to external fields. Hydration energies calculated using the PDPEOE representation are similar to those calculated by the SCF procedure. Also, the PDPEOE calculations yielded changes in molecular dipole moments upon solvation that agreed closely with the changes in the calculated ab initio SCF dipole moments.     

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 construction of all-atom loop conformations

In this study, a new ab initio method named CLOOP has been developed to build all-atom loop conformations. In this method, a loop main-chain conformation is generated by sampling main-chain dihedral angles from a restrained varphi/psi set, and the side-chain conformations are built randomly. The CHARMM all-atom force field was used to evaluate the loop conformations. Soft core potentials were used to treat the non-bond interactions, and a designed energy-minimization technique was used to close and optimize the loop conformations. It is shown that the two strategies improve the computational efficiency and the loop-closure rate substantially compared to normal minimization methods. CLOOP was used to construct the conformations of 4-, 8-, and 12-residue loops in Fiser's test set. The average main-chain root-mean-square deviations obtained in 1,000 trials for the 10 different loops of each size are 0.33, 1.27, and 2.77 A, respectively. CLOOP can build all-atom loop conformations with a sampling accuracy comparable with previous loop main-chain construction algorithms. [Figure: see text].     

Ab initio multiconfiguration reference perturbation theory calculations on the energetics of low-energy spin states of iron(III) porphyrins

Although a major goal of inorganic spectroscopy is to determine the energetics of the low-lying spin states of transition metal complexes, surprisingly little has been accomplished in this respect by means of accurate ab initio calculations. Against this context, we present ab initio multiconfiguration reference perturbation theory (CASPT2) calculations with large basis sets on the low-lying spin states of Fe(III)(P)Cl and [Fe(P)Cl](+) (P(2-)=porphinato). The CASPT2 results on the energetics of various low-lying spin states studied differ significantly, sometimes even dramatically, from those obtained from density functional theory calculations.     

Ab initio protein structure prediction using chunk-TASSER

We have developed an ab initio protein structure prediction method called chunk-TASSER that uses ab initio folded supersecondary structure chunks of a given target as well as threading templates for obtaining contact potentials and distance restraints. The predicted chunks, selected on the basis of a new fragment comparison method, are folded by a fragment insertion method. Full-length models are built and refined by the TASSER methodology, which searches conformational space via parallel hyperbolic Monte Carlo. We employ an optimized reduced force field that includes knowledge-based statistical potentials and restraints derived from the chunks as well as threading templates. The method is tested on a dataset of 425 hard target proteins < or =250 amino acids in length. The average TM-scores of the best of top five models per target are 0.266, 0.336, and 0.362 by the threading algorithm SP(3), original TASSER and chunk-TASSER, respectively. For a subset of 80 proteins with predicted alpha-helix content > or =50%, these averages are 0.284, 0.356, and 0.403, respectively. The percentages of proteins with the best of top five models having TM-score > or =0.4 (a statistically significant threshold for structural similarity) are 3.76, 20.94, and 28.94% by SP(3), TASSER, and chunk-TASSER, respectively, overall, while for the subset of 80 predominantly helical proteins, these percentages are 2.50, 23.75, and 41.25%. Thus, chunk-TASSER shows a significant improvement over TASSER for modeling hard targets where no good template can be identified. We also tested chunk-TASSER on 21 medium/hard targets <200 amino-acids-long from CASP7. Chunk-TASSER is approximately 11% (10%) better than TASSER for the total TM-score of the first (best of top five) models. Chunk-TASSER is fully automated and can be used in proteome scale protein structure prediction.     

Ab initio gene identification in metagenomic sequences

We describe an algorithm for gene identification in DNA sequences derived from shotgun sequencing of microbial communities. Accurate ab initio gene prediction in a short nucleotide sequence of anonymous origin is hampered by uncertainty in model parameters. While several machine learning approaches could be proposed to bypass this difficulty, one effective method is to estimate parameters from dependencies, formed in evolution, between frequencies of oligonucleotides in protein-coding regions and genome nucleotide composition. Original version of the method was proposed in 1999 and has been used since for (i) reconstructing codon frequency vector needed for gene finding in viral genomes and (ii) initializing parameters of self-training gene finding algorithms. With advent of new prokaryotic genomes en masse it became possible to enhance the original approach by using direct polynomial and logistic approximations of oligonucleotide frequencies, as well as by separating models for bacteria and archaea. These advances have increased the accuracy of model reconstruction and, subsequently, gene prediction. We describe the refined method and assess its accuracy on known prokaryotic genomes split into short sequences. Also, we show that as a result of application of the new method, several thousands of new genes could be added to existing annotations of several human and mouse gut metagenomes.