UV photoelectron spectroscopy and ab initio characterization of valence orbital structures and conformations of neutral phosphate esters

The HeI UV photoelectron spectrum of trimethyl phosphate (TMP) has been measured and interpreted with the aid of SCF molecular orbital calculations carried out with STO-3G, STO-3G* and 4-31G basis functions. The photoelectron spectrum of TMP is more accurately reproduced by results from 4-31G calculations than by results from STO-3G or STO-3G* calculations. However, all three basis sets yield results which predict the same assignment of the photoelectron spectrum. Results at the 4-31G level indicate that whether calculations are based on crystallographic bond angles and bond lengths or on STO-3G optimized geometries has little effect on the energetic ordering of the upper occupied orbitals. The energetic ordering of orbitals is also found to be only weakly dependent upon the torsional angle phi, describing rotation of ester groups about P-O bonds and upon the torsional angle psi, describing rotation of methyl groups about C-O bonds. For trimethyl phosphate, with C3 symmetry, the vertical ionization potentials of the upper occupied orbitals are 10.81 eV (8e), 11.4 eV (9a), 11.93 eV (7e), 12.6-12.9 eV (8a and 6e), 14.4 eV (7a) and 15.0-16.0 eV (5e and 6a). Calculations at the 4-31G level indicate that many of the highest occupied orbitals in neutral dimethyl phosphate and methyl phosphate have energies and electron distributions similar to orbitals in TMP. For TMP, a search for optimized values of phi and psi has been carried out at the STO-3G*level. In agreement with previous NMR studies and with classical potential calculations, the STO-3G* results indicate that both the gauche (phi = 53.1 degrees) and anticlinal (phi = 141.9 degrees) conformations are thermally accessible. Also in agreement with the classical potential calculations, the STO-3G* results predict that in the all gauche conformation energy is minimized when the methyl groups assume a staggered geometry (psi = 60 degrees to 80 degrees) and that an energy maximum occurs for an eclipsed geometry (phi = 0 degrees to 20 degrees). A study of the dependence of optimized values of O-P-O ester bond angles on the torsional angles, phi, was carried out at the STO-3G, STO-3G* and 4-31G levels. The results demonstrate that for C3 symmetry, the coupling of O-P-O angles to phi is influence by repulsive steric interactions.     

Conformation of 1,4-dihydropyridine--planar or boat-like?

The geometry of the 1,4-dihydropyridine molecule was completely optimized employing three different ab initio basis sets (6-31 G*, 4-31G, STO-3G). The most reliable 6-31G* basis set provides a very flat boat conformation which may easily undergo defolding to a planar ring arrangement. This result is discussed with respect to enzymatic redox cofactors and the pharmacological activity of dihydropyridine calcium antagonists.     

UV Photoelectron Spectroscopy and Ab Initio Characterization of Valence Orbital Structures and Conformations of Neutral Phosphate Esters

Abstract

The Hel UV photoelectron spectrum of trimethyl phosphate (TMP) has been measured and interpreted with the aid of SCF molecular orbital calculations carried out with STO-3G, STO-3G* and 4–31G basis functions. The photoelectron spectrum of TMP is more accurately reproduced by results from 4–31G calculations than by results from STO-3G or STO-3G* calculations. However, all three basis sets yield results which predict the same assignment of the photoelectron spectrum. Results at the 4–31G level indicate that whether calculations are based on crystallographic bond angles and bond lengths or on STO-3G optimized geometries has little effect on the energetic ordering of the upper occupied orbitals. The energetic ordering of orbitals is also found to be only weakly dependent upon the torsional angle φ, describing rotation of ester groups about P-O bonds and upon the torsional angle ψ, describing rotation of methyl groups about C-O bonds. For trimethyl phosphate, with C₃ symmetry, the vertical ionization potentials of the upper occupied orbitals are 10.81 eV (8e), 11.4 eV (9a), 11.93 eV (7e), 12.6–12.9 eV (8a and 6e), 14.4 eV (7a) and 15.0–16.0 eV(5e and 6a). Calculations at the 4–31G level indicate that many of the highest occupied orbitals in neutral dimethyl phosphate and methyl phosphate have energies and electron distributions similar to orbitals in TMP.

For TMP, a search for optimized values of φ and ψ has been carried out at the STO-3G* level. In agreement with previous NMR studies and with classical potential calculations, the STO- 3G* results indicate that both the gauche φ= 53.1 °) and anticlinal (φ = 141.9°) conformations are thermally accessible. Also in agreement with the classical potential calculations, the STO-3G* results predict that in the all gauche conformation energy is minimized when the methyl groups assume a staggered geometry (ψ= 60° to 80°) and that an energy maximum occurs for an eclipsed geometry (ψ = 0° to 20°). A study of the dependence of optimized values of O-P-O ester bond angles on the torsional angles, φ, was carried out at the STO-3G, STO-3G* and 4–31G levels. The results demonstrate that for C₃ symmetry, the coupling of O-P-O angles to φ is influenced by repulsive steric interactions.     

RNA hydrolysis via an oxyphosphorane intermediate.

From calculations of a model reaction scheme for base-catalyzed RNA hydrolysis, a pentacoodinate dianionic intermediate 2a (Storer, et al., J. Am. Chem. Soc., 1991, 113, 5216-5219) as well as two transition states, TS1 and TS2, to the intermediate have been located by ab initio calculations at the 3-21G* level. Although the intermediate, which has the well depth on the order of kBT, is unlikely to be kinetically significant, the overall rate-limiting transition state structure TS2 obtained at 3-21G* level is very close to the corresponding structure at the STO-3G level; it has an extended P-O(5') bond breaking character. These gas-phase calculation results are used to qualitatively interpret mutagenesis results of Barnase and RNase T1 where water molecules are absent from the active site.     

Conformational preferences in glycosylamines. implications for the exo-anomeric effect

The conformational preferences about the C-N bond in N-(4-methoxyphenyl)-2,3,4,6-tetra-O-acetyl-alpha (1) and beta-D-glucopyranosylamine (2), in the solid state and in solution, have been investigated. The crystal structure of the axially substituted alpha anomer (1) indicates a conformational preference about the C-1-N bond in which nN-->sigma*C-O exo-anomeric interactions may be expressed, although this conformational preference is not displayed in solution. The solution conformation relieves steric interactions that result from expression of the exo-anomeric effect in the solid-state conformation. The conformational preference in the equatorially substituted beta anomer (2) both in solution and in the solid state is similar and permits expression of nN-->sigma*C-O exo-anomeric interactions. The structural data for 1 and 2 indicate significant differences in O-5-C-1-N-1 bond angles but insignificant differences in each of the O-5-C-1 or C-1-N-1 bond lengths. The J(C-1-H-1 coupling constants in 1 and 2 indicate a greater coupling constant for the alpha anomer that is consistent with a dominant nO-->sigma*C-H orbital interaction in the beta anomer that weakens the C-1-H-1 bond.     

Lanthanide complex coordination polyhedron geometry prediction accuracies of ab initio effective core potential calculations

The ability to efficiently and accurately predict solid-state geometries of lanthanide coordination compounds efficiently and accurately is central for the design of new ligands capable of forming stable and highly luminescent complexes. Accordingly, we present in this paper a report on the capability of various ab initio effective core potential calculations in reproducing the coordination polyhedron geometries of lanthanide complexes. Starting with all combinations of HF, B3LYP and MP2(Full) with STO-3G, 3-21G, 6-31G, 6-31G* and 6-31+G basis sets for [Eu(H₂O)₉]³⁺ and closing with more manageable calculations for the larger complexes, we computed the fully predicted ab initio geometries for a total of 80 calculations on 52 complexes of Sm(III), Eu(III), Gd(III), Tb(III), Dy(III), Ho(III), Er(III) and Tm(III), the largest containing 164 atoms. Our results indicate that RHF/STO-3G/ECP appears to be the most efficient model chemistry in terms of coordination polyhedron crystallographic geometry predictions from isolated lanthanide complex ion calculations. Moreover, both augmenting the basis set and/or including electron correlation generally enlarged the deviations and aggravated the quality of the predicted coordination polyhedron crystallographic geometry. Our results further indicate that Cosentino et al.’s suggestion of using RHF/3-21G/ECP geometries appears to be indeed a more robust, but not necessarily, more accurate recommendation to be adopted for the general lanthanide complex case. Figure Graphical visualization of unsigned mean errors, UME(Eu-L)s, involving only the interatomic distances between the europium central ion and the oxygen atoms of the coordination polyhedron of the cation nona-aqua-europium(III) for various model chemistries, all compared to the “Cambridge Structural Database 2004” crystallographic geometry     

Basis set validation for polyatomic cation-water interactions

A study of the effects of counterpoise (CP) corrections in polyatomic cation-water interactions is added to the systematic analysis performed in the past of the basis set superposition error (BSSE) in neutral and anionic adducts. The interaction with water of the ammonium cation and its methyl and -CHO derivatives is considered due to the need to model accurately this functional group, which is common in biological molecules. The basis sets employed are the STO-3G and MINI-1 minimal basis sets and the 3-21G, 4-31G and 6-31G** extended ones. In addition, the 6-311++G (2d,p) and 6-311++G (3d, 2p) basis sets have been used for the smallest system, at the SCF and MP2 levels, both without and with CP correction. These basis sets give an equilibrium distance slightly larger and an interaction energy less favourable than the 6-31G** basis set at the corresponding level, while the inclusion of correlation corrections produces a stronger H-bond at a shorter distance. The results confirm the previous hint of a lower incidence of BSSEs in medium size cationic systems, at a different extent for the various basis sets. While the STO-3G basis set is sharply affected by BSSEs in both the equilibrium distance and the interaction energy, the MINI-1 basis set shows a small BSSE, though its trend is not completely satisfactory because the charge transfer component has an anomalous behaviour with respect to our reference basis set. The 4-31G basis set is the only one able to hold comparison with the 6-31G**, even though the interaction energy produced is slightly overestimated. The 3-21G basis set, when corrected, almost parallels the 4-31G one in this set of compounds. The reliability of the CPED correction is checked and discussed.     

Conformation and charge distribution of bicyclic beta-lactams: structure-activity relationships.

The structures of 7-oxo-1-azabicyclo[3.2.0]heptane and its 4-oxa, 3-ethylene-4-oxa, and 3-ethylene-6-methyl-4-oxa derivatives, and of 8-oxo-1-azabicyclo[4.2.0]octane and its 5-oxa derivative, were studied by ab initio methods. Conformations were refined without constraints using the 4-21G and the 4-21G* basis sets, and energies and charge distributions were improved by single-point 6-31G*/4-21G* calculations. The results are are interpreted in terms of structural trends related to beta-lactamase inhibitor capability.     

An AB initio investigation of molecules with a disulfide bond: (HS)2, (CH3S)2 and cystine

ab initio Calculations at the Gaussian-70 STO-3G and 4-31G basis levels have been carried out for (HS)2 and (CH3S)2. Cystine was investigated at the STO-3G level. The STO-3G energy minimized geometry agrees well with experiments for (HS)2 and (CH3S)2. The barriers to internal rotation are predicted to be (at the 4.31G level): (HS)2, cis 8.5 kcal, trans 3.03 kcal; (CH3S)2, cis 18.47 kcal, trans 6.04 kcal.     

Theoretical investigation of the triphosphate forms of azidothymidine and thymidine

In this paper we investigate (using AM1 semi-empirical as well as HF methods at the STO-3G, 3-21G, 6-31G, 6-31G* and 6-31+G** level) the conformations, geometrical parameters, Mulliken charges, and solvation effects of the triphosphate form of AZT (AZTTP), as well as the thymidine nucleotide (dTTP) structure. Our calculated geometrical parameters and Mulliken charges, with and without solvation effects, are correlated with recent experimental results.     

An Investigation of Structural Stability and Internal Rotation in 3-Cyclopropenecarboxaldehyde and 3-Cyclopropenecarboxylic Acid Fluoride by ab initio Calculations

The structural stability and internal rotation in 3-cyclopropenecarboxaldehyde and 3-cyclo-propenecarboxylic acid fluoride were investigated by ab initio calculations with a 6-31G* atomic basis in the latter and a 6-311G* atomic basis set in the former case. For the sake of comparison also results obtained with a 3-21G basis are given in the paper. As expected, it turned out that this basis set is not large enough for three-membered rings. The calculations were carried out both at the Restricted Hartree-Fock (HF) and the second order Moller-Plesset (MP2) levels. The trans-form is predicted to be the lower energy conformer for both molecules. However, in case of the fluoride the two conformers have nearly the same energy. Full optimization was performed at the transition states and the asymmetric potential function for the CXO internal rotations was predicted for both molecules.     

The role of stacking interactions in clonidine binding

Stacking interactions of the clonidine aromatic ring with that of Phe or Tyr in the α₂-adrenoreceptor and Tyr in the tetrodotoxin-resistant sodium channel pore have been studied. Ab initio quantum-chemical calculations for a model system of two parallel aromatic rings have been performed with the GAMESS software using the 6-31G* basis set in the framework of the second-order Muller-Plesset perturbation theory with full geometry optimization without symmetry constraints. The parallel shifted conformation of two aromatic rings was found to be energetically most favorable. The 2′,6′-chlorination of one of the benzene rings enhances the stacking interaction, somewhat increases the shift of these rings, and possibly improves the hypotensive and analgesic functions of clonidine owing to an increase in the binding energy. The 4-fluorination of the clonidine ring can increase its analgesic effect but practically excludes its hypotensive activity.     

A computer-aided quantum chemical study of the N<Stack><Subscript>15</Subscript><Superscript>−</Superscript></Stack> cluster

Ab initio (RHF, MP2) and Density Functional Theory (DFT) methods have been used to examine six isomers of the N15m cluster with the 6-31+G* basis set. Different from the known odd-numbered anionic N7m, N9m, and N11m clusters, in which the open-chain structures are the most stable species, the most stable N15m isomer is structure 1 (C1), which may be considered as a complex between the fragments cyclic N5m (D5h) and staggered N10 (D2d). The decomposition pathways of structure 2 (CS), containing two aromatic N5 rings connected by a N5 chain, and the open-chain structure 3 (C2v) were studied at the B3LYP/6-31+G* level of theory. Relative energies were refined at the level of B3LYP/6-311+G(3df,2p)//B3LYP/6-31+G*+ZPE (B3LYP/6-31+G*). The barriers for N2 and N5m (D5h) fission reactions for structure 2 are predicted to be 18.2 and 14.2 kcal x mol(-1), respectively. The corresponding N2+N3m fission barrier for structure 3 is predicted to be 11.2 kcal x mol(-1). Supplementary material is available for this article if you access the article at http://dx.doi.org/10.1007/s00894-003-0118-0. A link in the frame on the left on that page takes you directly to the supplementary material. Figure Structure 1 of the N15m cluster, showing bond distances in A and bond angles in degrees     

嘌呤受羟基自由基损伤机制的量子化学研究

用量子化学从头算方法 ,在 3— 2 1 G基组水平上对腺嘌呤 ( A)和鸟嘌呤 ( G)受羟基自由基(·OH)进攻形成的各种可能产物自由基进行了几何全优化 ,然后在优化构型下 ,用 UHF/6— 31 G基组做单点计算 .根据计算结果 ,由能量、自旋密度和键长分析了羟基自由基造成 DNA的损伤及其修复机制 .结果表明羟基自由基易进攻腺嘌呤的 C— 4、C— 5位 ,形成的产物自由基 A4OH·比A5OH·易与 N— 1 0位 H脱水 ,且脱水后的产物自由基进行异构化 ,这样可通过加入氧化剂带走未成对电子从而使嘌呤体系得到修复 .鸟嘌呤也易形成 G4OH·和 G5OH· ,但 G5OH·比 G4OH·更易脱水 .另外 ,羟基自由基进攻腺嘌呤和鸟嘌呤的 C— 8位在能量上最为有利 .A8OH·与G8OH·会发生开环反应 ,一旦开环 ,DNA便不易修复     

Vibrational analysis of nucleic acids. V. Force field and conformation-dependent modes of the phosphodiester backbone modeled by diethyl phosphate.

A generalized valence force field is derived for the diethyl phosphate anion [(CH3CH2O)2PO2-] and its deuterium [(CH3CD2O)2PO2-, (CD3CH2O)2PO2- and (CD3CD2O)2PO2-] and carbon-13 [(CH3 13CH2O)2PO2-] derivatives in the stable trans-gauche-gauche-trans conformation. Normal coordinate analysis of the trans-gauche-gauche-trans conformer, which serves as a structural analog of the nucleic acid phosphodiester group, is based on comprehensive infrared and Raman spectroscopic data and vibrational assignments obtained for the diethyl phosphate anion. The generalized valence force field is in good agreement with the scaled ab initio force field of diethyl phosphate and represents significant improvement over earlier modeling of the phosphodiester moiety with dimethyl phosphate. The conformational dependence of skeletal C-C-O-P(O2-)-O-C-C stretching vibrations is also explored. Starting with the trans-gauche-gauche-trans conformation, the frequency dependence of skeletal stretching modes has been obtained by stepwise rotation of the torsion angles of the P-O and C-O bonds corresponding to nucleic acid torsions alpha (P-O5'), beta (O5'-C5'), epsilon (C3'-O3'), and zeta (O3'-P). Both symmetric and antisymmetric phosphoester stretching modes are highly sensitive to P-O and C-O torsions, whereas symmetric and antisymmetric phosphodioxy (PO2-) stretching modes are less sensitive. The present results provide an improved structural basis for understanding previously developed empirical correlations between vibrational marker bands and nucleic acid backbone conformation.     

Molecular structure and effects of intermolecular hydrogen bonding on the vibrational spectrum of trifluorothymine, an antitumor and antiviral agent

In the present work, the experimental and the theoretical vibrational spectra of trifluorothymine were investigated. The FT-IR (400-4000?cm(-1)) and μ-Raman spectra (100-4000?cm(-1)) of trifluorothymine in the solid phase were recorded. The geometric parameters (bond lengths and bond angles) and vibrational frequencies of the title molecule in the ground state were calculated using ab initio Hartree-Fock (HF) method and density functional theory (B3LYP) method with the 6-31++G(d,p) and 6-311++G(d,p) basis sets for the first time. The optimized geometric parameters and the theoretical vibrational frequencies were found to be in good agreement with the corresponding experimental data and with results found in the literature. Vibrational frequencies were assigned based on the potential energy distribution using the VEDA 4 program. The dimeric form of trifluorothymine was also simulated to evaluate the effect of intermolecular hydrogen bonding on the vibrational frequencies. It was observed that the stretching modes shifted to lower frequencies, while the in-plane and out-of-plane bending modes shifted to higher frequencies due to the intermolecular N-H?O hydrogen bonds.     

Coordination and thermodynamics of stable Zn(II) complexes in the gas phase

Ab initio molecular orbital methods in combination with DFT calculations were used to study the structural and thermodynamic properties of 17 complexes containing zinc cation and four first-shell ligands as models of active site of metalloenzymes (e.g. angiotensin converting enzyme, thermolysin). The geometry of the complexes was relaxed by complete optimization by ab initio molecular orbital methods at Hertree-Fock level with 3-21G* basis set. Following single point calculation with tight SCF criteria at the B3LYP level with 6-311+G(2d,p) basis set was used to calculate accurate interaction enthalpies. The structure and thermodynamics of optimized complexes are discussed from the point of view of their biological importance.     

QM/MM study of the active site of free papain and of the NMA-papain complex.

Hybrid quantum mechanical/molecular mechanical (QM/MM) calculations using restricted and unrestricted Hartree-Fock and B3LYP ab initio (QM) and Amber force field (MM), respectively, have been applied to study the catalytic site of papain in both free and substrate bonded forms. Ab initio geometry optimizations have been performed for the active site of papain and the N-methyl-acetamide (NMA)-papain complex within the molecular mechanical treatment of the protein environment. A covalent tetrahedral intermediate structure could be obtained only when the amide N atom of the substrate molecule was protonated through a proton transfer from the His-159 in the catalytic site. Our results support the previous assumption that a proton transfer from His-159 to the amide N atom of the substrate occurs prior to or concerted with the nucleophilic attack of the Cys-25 sulfur atom to the carbonyl group of the substrate. The electron correlation effect will reduce the proton transfer barrier. Therefore, this proton transfer can be easily observed in the B3LYP/6-31G* calculations. The HF/6-31G* method overestimates the reaction barrier against this proton transfer. The sulfur atom of Cys-25 and the imidazole ring of His-159 are found to be coplanar in the free form of the enzyme. However, the rotation of the imidazole ring of His-159 was observed during the formation of the tetrahedral intermediate. Without the papain environment, the coplanar thiolate-imidazolium ion pair RS-...ImH+ is much less stable than the neutral form of RSH....Im. Within the protein environment, however, the thiolate-imidazolium ion pair becomes more stable than its neutral form by 4.1 and 0.4 kcal/mol in HF/6-31G* and B3LYP/6-31G* calculations, respectively. The barrier of proton transfer from S-H group of Cys-25 to the imidazole ring of His-159 was reduced from 22.0 kcal/mol to 15.2 kcal/mol by the protein environment in HF/6-31G* calculations. This barrier is found to be much smaller (2.5 kcal/mol) in B3LYP/6-31G* calculations.     

Theoretical Investigation of the Molecular Structure of the πκ DNA Base Pair

Abstract

The structure of the nonclassical πκ base pair (7–methyl-oxoformycin … 2,4-diaminopyrimidine) was studied at the ab initio Hartree-Fock (HF) and MP2 levels using the 6–31G* and 6–31G** basis sets. The πκ base pair is bound by three parallel hydrogen bonds with the donor-acceptor-donor recognition pattern. Recently, these bases were proposed as an extension of the genetic alphabet from four to six letters (Piccirilli et al. Nature 343, 33(1990)). By the HF/6- 31G* method with full geometry optimization we calculated the 12 degree propeller twist for the minimum energy structure of this complex. The linearity of hydrogen bonds is preserved in the twisted structure by virtue of the pyramidal arrangement of the κ-base amino groups. The rings of both the π and κ molecules remain nearly planar. This nonplanar structure of the πκ base pair is only 0.1 kcal/mol more stable than the planar (Cs) conformation. The HF/6- 31G* level gas-phase interaction energy of πκ (—13.5 kcal/mol) calculated by us turned out to be nearly the same as the interaction energy obtained previously for the adenine-thymine base pair (—13.4 kcal/mol) at the same computational level. The inclusion of p-polarization functions on hydrogens, electron correlation effects (MP2/6–31G** level), and the correction for the basis set superposition error (BSSE) increase this energy to -14.0 kcal/mol.     

Study of the hydrogen bond in different orientations of adenine-thymine base pairs: An <Emphasis Type="Italic">ab initio</Emphasis> study

In order to gain deeper insight into structure, charge distribution, and energies of A-T base pairs, we have performed quantum chemical ab initio and density functional calculations at the HF (Hartree-Fock) and B3LYP levels with 3-21G*, 6-31G*, 6-31G**, and 6-31++G** basis sets. The calculated donor-acceptor atom distances in the Watson-Crick A-T base pair are in good agreement with the experimental mean values obtained from an analysis of 21 high resolution DNA structures. In addition, for further correction of interaction energies between adenine and thymine, the basis set superposition error (BSSE) associated with the hydrogen bond energy has been computed via the counterpoise method using the individual bases as fragments. In the Watson-Crick A-T base pair there is a good agreement between theory and experimental results. The distances for (N2...H23-N19), (N8-H13...O24), and (C1...O18) are 2.84, 2.94, and 3.63 A, respectively, at B3LYP/6-31G** level, which is in good agreement with experimental results (2.82, 2.98, and 3.52 A). Interaction energy of the Watson-Crick A-T base pair is -13.90 and -10.24 kcal/mol at B3LYP/6-31G** and HF/6-31G** levels, respectively. The interaction energy of model (9) A-T base pair is larger than others, -18.28 and -17.26 kcal/mol, and for model (2) is the smallest value, -13.53 and -13.03 kcal/mol, at B3LYP/6-31G** and B3LYP/6-31++G** levels, respectively. The computed B3LYP/6-31G** bond enthalpies for Watson-Crick A-T pairs of -14.4 kcal/mol agree well with the experimental results of -12.1 kcal/mol deviating by as little as -2.3 kcal/mol. The BSSE of some cases is large (9.85 kcal/mol) and some is quite small (0.6 kcal/mol).