Conformational analysis of dimethylphosphate with quantum-mechanical and classical methods

In order to shed light on the conformational behavior of polynucleotide chains, and in particular to clarify the origins of the barriers to internal rotation in the phosphodiester linkage, we computed, with a quantum-mechanical ab initio procedure, the energies associated to 86 combinations of the two torsion angles in the dimethylphosphate anion (CH₃O)₂PO₂ ^–, and then we sought for an analytical expression apt to reproduce these energies with the highest possible accuracy. An excellent agreement (standard deviation of the fitted energies from the ab initio energies 0.28 kcal/mole) with the quantum-mechanical calculations was reached with a potential consisting of four terms: 1) a 6–12 Lennard-Jones contribution, in which different parameters are used to describe the interactions of methyls with the ester oxygens and with the anionic oxygens; 2) a contribution with twofold periodicity, accounting for the anomeric effects connected to the interactions between the lone pair electrons and the polar bonds of phosphorus with the anionic oxygens; 3) a contribution with threefold periodicity, representing the usual bond-staggering term; and 4) a Coulombic contribution, arising from electrostatic interactions between partially charged atoms.     

Molecular orbital studies on the conformation of GABA (γ-aminobutyric acid)

In distinction to Extended Hückel Theory which predicts as the most stable conformation of free zwitterionic GABA a totally extended form, PCILO and SCF ab initio studies show that the intrinsically preferred conformation of the isolated molecule is a highly folded one, resulting from strong interactions between the two charged ends. Computations are also carried out for hydrated GABA in the supermolecule approach allowing moreover for the flexibility of binding of some of the water molecules of the first hydration shell. They predict the coexistence in solution of a large number of conformations showing different degrees of folding (or extension), a result confirmed by recent NMR studies. This and a number of similar results show that we have to adapt our thinking on the role of conformations in pharmacological activity to this situation, which was frequently obscured by the more abundant results of X-ray crystallography yielding a single conformation.     

On the sites and mechanisms of alkylation in the nucleic acids

The intrinsic capacity (computed by ab initio molecular orbital techniques) of ethyl and methyl cations to bind to the anionic oxygens of the phosphate groups of the nucleic acids is found to be larger than their affinity for the carbonyl oxygens or for the nitrogens of the bases cytosine and guanine. Analysis of the components of the binding energies indicates that the preference of the ethyl cation for the oxygens of the phosphate or for those of the bases is due to different reasons. The results are discussed in connection with the available experimental observations.     

Ab initio crystal orbital studies on linear chains of H atoms

Anab initio crystal orbital method is used to calculate the energies of an infinite chain of H atoms and of linear arrangements of H₂ molecules with different interatomic distances. The H₂ arrangements are not stable in respect to isolated molecules. The cohesive energy of an optimized arrangement of H atoms chain is 0.0354 a.u.     

Quantum-mechanical exploration of the properties of the sugar rings

Theab initio SCF method is used for computing the main electronic properties of the ribose unit of the nucleic acids. The present study is devoted to the ribose in the C₃-endo,gg conformation. The properties investigated include the distribution of the electronic charges, the electrostatic molecular potential around the four oxygens of the unit, the hydration and the Na⁺ binding schemes studied in the supermolecule approximation. The possibilities of through-water binding of the cation to the sugar are also explored. The predictions of the computation in particular with regard to cation binding to the ribose ring are correlated with recent experimental results.     

Theoretical study of the electronic structure of diazomethane

The least-energy dissociation path of the ground state of CH₂N₂ was determined fromab initio calculations using in a complementary way basis sets of minimal size (STO-3G) and double-zeta (DZ) quality. The results indicate that the least-energy point of attack of the N₂ molecule on CH₂ (¹ A ₁) is roughly perpendicular to the molecular plane (93 °), the C and N atoms being almost co-linear (angle C-N-N203 ° with outermost N atom pointing away from CH₂). The potential barrier of 1.2 eV found previously on theC _2v dissociation path, disappears completely along the least-energy dissociation path (point groupC _s (out-of-plane)). These findings corroborate the Woodward-Hoffman rules for this process since the outermost orbitals of the two intersecting states found in point groupC _2v (...2b ₁ and ...8a ₁) both correlate to the same irreducible representation (10á) in point groupC _s (out-of-plane).Larger basis set calculations (DZ + polarization functions on all centers, 3d _c and 3d _N developed here), were also carried out on CH₂N₂ (¹ A ₁,³ A ₂ and¹ A ₂) at the¹ A ₁ equilibrium geometry and on CH₂ (³ B ₁) and N₂ (¹ _g ⁺ ) at their respective equilibrium geometries. These calculations, together with consideration of correlation energy differences, yieldD ₀ ⁰ (CH₂N₂,¹ A ₁) = 19 kcal/mole and vertical excitation energies of 67 and 73 kcal/mole for the³ A ₂ and¹ A ₂ states respectively. The latter value is in good agreement with the measured experimental value: 72.4 kcal/mole corresponding to the maximum of intensity in the¹ A ₂¹ A ₁ absorption band.     

Conformations of the thiathiophthenes and related molecules

A conformational ab initio MO study has been carried out for the thiathiophthene molecule (TTP) and two related model compounds, thiomalonaldehyde (TMA) and its conjugate base (TMA(-)). The conformational energy surfaces for TMA, TMA(-) and TTP were generated using a least squares fit to the calculated data and plotted on a CALCOMP plotter. The results of the calculations showed that the cis-cis planar conformation of TTP is the most stable in agreement with experimental findings. For TMA and TMA(-) the cis-cis planar conformation is not the most stable. Contour plots of the five occupied -MO's of TTP show great similarity to those of naphthalene.Less detailed calculations were carried out for 3-hydroxy-prop-2-en-1-thione (HPT) and 3-mercapto-prop-2-en-1-thione (MPT). HPT was shown to be most stable in the cis planar hydrogen bonded conformation in agreement with the experimentally obtained results. For MPT the non-hydrogen bonded planar structure was found to be the most stable.     

Quantum calculation of the intrinsic torsional barrier of C-C and C-O bonds

PCILO and ab initio calculations have been performed to investigate the energies associated to rotation about the central bond in n-butane and methyl ethyl ether. Quantum mechanical energies have been fit to a classical intramolecular force field, containing torsional and nonbonded (Lennard-Jones 6–12 plus Coulomb) contributions, with a standard deviation comprised between 0.03 and 0.09 kcal mol^–1. Two conditions have proved indispensable to reach such level of accuracy: (a) the use of a torsional potential with threefold periodicity, which corrects for the part of the rotation barrier not covered by van der Waals repulsions and may be interpreted as bond-bond repulsion; (b) the introduction in the force field for ethers of terms accounting for orbital interaction effects of different nature than the normal molecular mechanics nonbonded interactions; these terms are represented either by low order rotational potential functions or preferably by interactions of atoms simulating lone-pair orbitals and bonded to oxygen in such a way as to render it sp ³-hybridized. According to ab initio, the height of the threefold torsional potential about C-C and C-O bonds is comparable and is of the order of 3 kcal mol^–1. According to PCILO, it is larger for C-C (ca. 1.5 kcal mol^–1) than for C-O (ca. 0.5 kcal mol^–1).     

Ab initio study of hydrazoic acid

SCF and CI calculations were carried out on the ground^1A state of HN₃. The equilibrium geometry and vibration frequencies were computed. The results point to a planar structure (groupC _s) but to a non-linear (170 °) N-N-N conformation. The calculated vibration frequencies are in fair agreement with experimental assignments.The dissociation path of the molecule to NH and N₂ products was investigated and compared to the isoelectronic reaction of diazomethane. The dissociation energy of hydrazoic acid is estimated to be about –8 kcal/mole, with a potential barrier to dissociation of about 30 kcal/mole.Boursier IRSIA     

Non-empirical SCF MO studies on the protonation of biopolymer constituents

Proton affinities of a large number of tautomeric adenine and guanine structures have been calculated using theab initio Hartree-Fock-Roothaan SCF method. For guanine, and to a lesser extent for adenine, it is likely that several protonated forms may co-exist in acidic solutions. Protonation at the free adenine and guanine positions in DNA may introduce effective acceptor levels between the energy bands of the polymer that may cause fundamental changes in the conduction properties of DNA.     

A theoretical study of electrophilic reactions

The electronic structures of protonated formyl and acetylium cations and their deprotonation paths leading to HCO⁺, COH⁺ and CH₃CO⁺have been studied by means of ab initio calculations. The results support Olah's theory that dipositive species can be the de facto reagents in electrophilic reactions.     

Quantum chemical studies on electrophilic addition

The geometries of the 2-hydroxyethyl and isomeric oxiranium cations have been fully optimized using ab initio molecular orbital calculations employing the split valence shell 4-31G basis set. These species are possible intermediates in both the electrophilic addition of OH to ethylene and in the acid catalysed ring opening of oxirane. The optimized structures were then used to compute more accurate wave functions using Dunning's double-zeta basis set, and with this large basis set the bridged oxiranium ion was found to be the more stable by 7.2 kcal/mole. The barrier to interconversion of these two C₂H₄OH ions was computed to be 25.0 kcal/mole above the oxiranium ion.     

Cation binding to biomolecules

SCF ab initio computations are carried out on the binding of alkali and alkaline-earth cations to the phosphate monoanion. The effect of the binding on the conformational properties of the phosphodiester linkage and of the polar head of phospholipids is investigated. The results indicate that following the nature of the cation and the site of its binding, the interaction may have a profound influence on the conformation of the ligand. The consequences of this situation on the use of lanthanide probes in NMR studies are considered.     

A non-empirical molecular orbital study on the relative stabilities of adenine and guanine tautomers

The relative stabilities of a series of adenine and guanine tautomers have been calculated using anab initio Hartree-Fock-Roothaan SCF MO method. The calculated relative stabilities agree in general with the results of earlier semiempirical studies. According to the present study, tautomeric forms with regular Kekulé structure for the six-membered purine ring are the most stable. The amine-imine tautomerization of purine bases is not likely to be responsible for spontaneous mutations in DNA.     

Theoretical studies on the conformation of saccharides

Conformations of 2-methoxytetrahydropyran as a model for the six-membered ring in aldopyranosides have been calculated by the PCILO method using the algorithm of the conjugated gradient to optimize the geometry. The calculated geometry of the fourteen basic forms of 2-methoxytetrahydropyran was found to be in agreement with the available data obtained by X-ray diffraction of pyranosides. The results indicate differences in the geometry of 2-methoxytetrahydropyran resulting from the change of the axial vs. equatorial position of the methoxyl group. These changes are particularly meaningful in the values of bond angles and they are in agreement with the anomeric and exoanomeric effects. The experimentally found differences in the energies of an axial (⁴ C ₁) and equatorial (¹ C ₄) conformer, G = 2.9–3.7 kJ/mol, and the dipole moment, = 1.20 ± 0.05 D (1D = 3.33 10^–30mAs) agree well with the calculated values E = 3.18 kJ/mol and <> = 1.18 D which, in turn, suggest that the axial conformer is preferred over the equatorial one by a ratio a:e = 78:22.     

Non-empirical SCF MO studies on the protonation of biopolymer constituents

Side-chain proton affinities for a series of peptide-forming amino acids have been calculated using ab initio Hartree-Fock-Roothaan SCF method; the order of proton affinities is Arg > His > Ser Tyr Lys. Protonation of some side-chains in a protein may introduce new energy levels in the band gaps thus fundamentally altering the conduction properties of the proteins.This paper is dedicated to Professor H. Hartmann for his 65th birthday.Guest at the Laboratory of the National Foundation for Cancer Research at the Chair of Theoretical Chemistry, University of Erlangen-Nürnberg.     

Extended Hückel calculations on the chemisorption of acetylene on tungsten

The chemisorption of acetylene on various faces of bcc tungsten, modelled by clusters of various sizes, has been studied by the Extended Hückel method. Relative binding energies on the various sites have been obtained, and the bonding modes of acetylene are discussed, also in comparison with experiment. It is concluded that, whenever possible, acetylene is di- bonded to the surface, and sp ²hybridization is considered likely in all cases.This paper is dedicated to Professor Hermann Hartmann on the occasion of his 65th birthday.     

An ab initio study of formamide

Calculated energy and molecular properties of the ground and low-energy excited states of formamide are presented at the ground state geometry. Satisfactory results are obtained except for the ¹* energy which remains too high by 1 eV (which is nevertheless a large improvement over previous calculations). The predicted triplet energies lie at 5.4 eV (³ n*) and 5.8 eV (³*).