Modeling mechanisms of unusual benzene imine n6 adduct formation in carcinogenic reactions of arylnitrenium ions with adenosine

Reaction mechanisms of the unusual benzene imine N6 adduct formation in carcinogenic reactions of arylnitrenium ions with adenosine have been investigated with density functional theory (DFT) and high-level ab initio methods. The DFT calculations indicate that the reaction pathways initiated by attack of adenine at the ortho C site of 4-biphenylylnitrenium ion are favored. However, high-level MP2 and QCISD calculations provide a contrary conclusion, that is the reaction pathways initiated by attack of adenine at the para C site of 4-biphenylylnitrenium ion are more feasible. Comparing with experimental results, the conclusion from high-level ab initio calculations is ultimately supported. The present study makes a theoretical prediction on the final products in the studied reaction, which is in agreement with experimental observations. In addition, this study provides some inspirations to the attacks of arylnitrenium ions at amino group of purines and pyrimidines in similar carcinogenic reactions.     

An experimental and theoretical study of the 13C and 31P chemical shielding tensors in solid O-phosphorylated amino acids

A series of l and dl forms of O-phosphorylated amino acids (serine, threonine, tyrosine) have been studied by using solid-state multinuclear NMR spectroscopy and ab initio calculations. Principal elements of the (13)C and (31)P chemical shielding tensors have been measured and discussed in relation to zwitterionic structures and intermolecular contacts. DFT calculations have been compared with experimental data showing their ability to reproduce experimentally obtained tensor values in this challenging class of compounds. The changes of orientation of (31)P chemical shielding tensor with respect to the molecular frame in the presence of hydrogen bonds have been revealed and discussed on the ground of theoretical calculations. The measurements of internuclear P...P distances, based on Zeeman magnetization exchange between (31)P spins with differing chemical shielding tensor orientations, were exploited for a clear distinction between enantiomers and racemates.     

Density functional theory study on molecular structure and vibrational IR spectra of isocytosine

Density functional theory with the combined Becke3-LYP exchange-correlation energy functional [DFT(B3-LYP) method] using the 6-31G(d, p) basis set is applied to predict molecular parameters (geometries, rotational constants, dipole moments) and vibrational IR spectra (harmonic wavenumbers, absolute intensities) of six tautomers of the isocytosine molecule. The results are compared with the corresponding data calculated at the conventional ab initio Hartree-Fock (HF) level using the same basis set and with available experimental data. Calculations show that (a) three amino tautomers are slightly nonplanar species with, evidently, a distorted amino group, (b) the DFT (B3-LYP)/6-31G(d, p) method predicts better molecular parameters, than do the HF calculations, and (c) the DFT(B3-LYP)-calculated vibrational IR spectra of isocytosine agree well with the available recorded IR spectra, and they show marked improvement over the IR spectra predicted at the HF/6-31G(d, p) level. Tautomeric stabilities of isocytosine are discussed on the basis of computed electronic energies by the DFT(B3-LYP) and ab initio approaches [including the MP2 and MP4(SDQ) calculations of electronic energies] and predicted zero-point vibrational energies by DFT(B3-LYP) and HF methods. This relative energies at 0 K of the tautomeric forms of isocytosine predicted by both conventional ab initio and DFT(B3-LYP) methods correlate well with the experimental data, showing the predominance of the aminohydroxy tautomer of isocytosine for an isolated molecule. © 1997 John Wiley & Sons, Inc.     

On the question of salt bridges of cationized amino acids in the gas phase: glycine and arginine

The geometrical shapes of the sodiated and cesiated amino acids glycine and arginine were probed in the gas phase by using the ion mobility based ion chromatography method. The data were compared to those obtained for alkali cationized methyl esters and for all the protonated species. Molecular mechanics, semiempirical, and ab initio/density functional theory (DFT) calculations were carried out to generate model structures for comparison with experiment and to determine the relative energies of different structures. For alkali cationized glycine the experimental cross sections agreed with charge solvation structures which were found by calculation to be the most stable forms as well. Both experiment and theory indicated that sodium is solvated by both the amino and the carbonyl groups, while cesium is solvated by one or both oxygen(s) of the carboxyl group. Alkali cationized arginine was found to form a salt bridge structure. The carboxylate group is stabilized by both the charged guanidinium group and the alkali ion. High level (6-311++G7 and DZVP) ab initio/DFT calculations carried out on sodiated and rubidiated N amidinoglycine, which contains a guanidino group and which was used as a model for the larger arginine molecule, indicated that the salt bridge structures are ∼10 kcal/mol more stable than the charge solvation forms for both alkali ions. The structure of protonated arginine, i.e. salt bridge or charge solvation, could not be unambiguously determined.     

Conformational study of 3β,16β-dihydroxyfriedelane by NMR and ab initio calculations

The complete chemical shift assignments of 3β,16β-dihydroxyfriedelane, a natural compound isolated from the hexane extract of Maytenus acanthophylla leaves was achieved through 1D/2D NMR spectral data. Combining ab initio Hartree-Fock [HF/6-31G(d)] and Density Functional Theory [DFT/B3LYP/6-31G(d,p)] calculations with NMR spectral data, it was possible to establish the predominance of the chair-chair-chair-boat-boat conformation for the five six-membered ring systems in this pentacyclic triterpene.     

Hydrogen-bonding interactions in cinchonidine-2-methyl-2-hexenoic acid complexes: a combined spectroscopic and theoretical study

Molecular interactions between cinchonidine (CD) and 2-methyl-2-hexenoic acid (MHA) have been studied by means of NMR, ATR-IR MES, DFT, and ab initio molecular dynamics. These interactions are of particular interest due to their pivotal role in the chiral induction occurring in the heterogeneous catalytic asymmetric hydrogenation of alpha,beta-unsaturated acids. The population density of the Open(3) conformer of CD, the most populated one at room temperature in apolar solvents, considerably increased to a maximum by addition of MHA to CD in toluene. The CD-MHA complex showed prominent symmetric and asymmetric carboxylate stretching vibrations in the regions of 1350-1410 and 1520-1580 cm(-1), respectively. DFT calculations revealed that these vibrational frequencies are expected to significantly shift depending on the chemical surrounding of MHA, that is, the hydrogen bond network. Earlier postulated 1:1 binding between CD and MHA was considered unlikely; instead, a dynamic equilibrium involving the MHA monomer and dimer, the 1:3 and possibly 1:2 CD-MHA complexes, were rationalized. Stable CD-MHA structures suggested by DFT calculations are the "1:3, halfN, cyclic" and the "1:3, halfN, cyclic tilted" complexes, where three MHA molecules are connected in wire by hydrogen bonding, two having direct interaction with CD. The confinement of CD's torsional motions in the complexes, leading to a slightly distorted Open(3) conformer via specific hydrogen-bonding interactions, was clearly reproduced by ab initio molecular dynamics, and the stable and flexible nature of the interaction was verified. Theoretical IR spectra of the complexes reproduced the characteristic vibrational frequencies of the complexes observed experimentally, supporting the stability of the 1:3 and implying the possibility of even higher molecular weight CD-MHA complexes.     

Intramolecular SN2 reaction caused by photoionization of benzene chloride-NH3 complex: direct ab initio molecular dynamics study

Ionization processes of chlorobenzene-ammonia 1:1 complex (PhCl-NH3) have been investigated by means of full dimensional direct ab initio molecular dynamics (MD) method, static ab initio calculations, and density functional theory (DFT) calculations. The static ab initio and DFT calculations of neutral PhCl-NH3 complex showed that one of the hydrogen atoms of NH3 orients toward a carbon atom in the para-position of PhCl. The dynamics calculation for ionization of PhCl-NH3 indicated that two reaction channels are competitive with each other as product channels: one is an intramolecular SN2 reaction expressed by a reaction scheme [PhCl-NH3]+-->SN2 intermediate complex-->PhNH3++Cl, and the other is ortho-NH3 addition complex (ortho complex) in which NH3 attacks the ortho-carbon of PhCl+ and the trajectory leads to a bound complex expressed by (PhCl-NH3)+. The mechanism of the ionization of PhCl-NH3 is discussed on the basis of the theoretical results.     

Nitrogen 1s near-edge X-ray absorption fine structure spectroscopy of amino acids: Resolving zwitterionic effects

Considerable variation is observed in the near-edge X-ray absorption fine structure (NEXAFS) spectra of amino acids. To unambiguously characterize the chemical origin of this variation, we have acquired the nitrogen 1s NEXAFS spectra of several amino acids and other model compounds and complemented these experimental measurements with ab initio calculations of isolated molecules and molecular clusters. The systematic differences observed between the zwitterionic and un-ionized forms of amino acids arise directly from the structural difference (-NH2 vs -NH3+), which leads to a change in the degree of Rydberg-valence mixing. Further change arises from quenching of this Rydberg character in the spectra of condensed amino acids. Ab initio calculations are used to explore the degree of Rydberg-valence mixing in the solid state.     

Exploring the Molecular Structure of Imidazolium–Silica‐Based Nanoparticle Networks by Combining Solid‐State NMR Spectroscopy and First‐Principles Calculations

A DFT‐based molecular model for imidazolium–silica‐based nanoparticle networks (INNs) is presented. The INNs were synthesized and characterized by using small‐angle X‐ray scattering (SAXS), NMR spectroscopy, and theoretical ab initio calculations. ¹¹B and ³¹P HETCOR CP MAS experiments were recorded. Calculated ¹⁹F NMR spectroscopy results, combined with the calculated anion–imidazolium (IM) distances, predicted the IM chain density in the INN, which was also confirmed from thermogravimetric analysis/mass spectrometry results. The presence of water molecules trapped between the nanoparticles is also suggested. First considerations on possible π–π stacking between the IM rings are presented. The predicted electronic properties confirm the photoluminescence emissions in the correct spectral domain.     

Synthesis and guest recognition ability of 2,3-dimethoxy-1,4-phenylene-containing porphyrinoids

The acid-catalyzed condensation of the bispyrrolylbenzene derivative and benzaldehyde yielded macrocycles 1 and 2 bearing three and four dipyrrin units which are connected by 2,3-dimethoxy-1,4-phenylene rings. The cationic guest recognition ability of 1 was investigated by UV-vis absorption, (1)H NMR spectroscopic techniques, and ab initio calculations (HF/3-21G(*)). The tris-dipyrrin macrocycle 1 was found to recognize alkali metals in the O6 binding cavity.     

Benchmark database of barrier heights for heavy atom transfer, nucleophilic substitution, association, and unimolecular reactions and its use to test theoretical methods

A benchmark database of forward and reverse barrier heights for 19 non-hydrogen-transfer reactions has been developed by using Weizmann 1 calculations, and 29 DFT methods and 6 ab initio wave-function theory (WFT) methods have been tested against the new database as well as against an older database for hydrogen atom transfer reactions. Among the tested hybrid DFT methods without kinetic energy density, MPW1K is the most accurate model for calculations of barrier heights. Among the tested hybrid meta DFT methods, BB1K and MPWB1K are the two most accurate models for the calculations of barrier heights. Overall, the results show that BB1K and MPWB1K are the two best DFT methods for calculating barrier heights, followed in order by MPW1K, MPWKCIS1K, B1B95, MPW1B95, BH and HLYP, B97-2, mPW1PW91, and B98. The popular B3LYP method has a mean unsigned error four times larger than that of BB1K. Of the methods tested, QCISD(T) is the best ab initio WFT method for barrier height calculations, and QCISD is second best, but QCISD is outperformed by the BB1K, MPWB1K, MPWKCIS1K, and MPW1K methods.     

Artarborol, a nor-caryophyllane sesquiterpene alcohol from Artemisia arborescens. stereostructure assignment through concurrence of NMR data and computational analysis

A nor-caryophyllane derivative, artarborol, has been isolated from wormwood (Artemisia arborescens) and its stereostructure established by using a combination of chemical derivatization, NMR data, molecular modeling, and quantum-mechanical calculations. In particular, comparison of experimental 13C NMR data with a Boltzmann-weighed average of 13C NMR chemical shifts, calculated by ab initio DFT method, supported the stereochemical assignment.     

Experimental and theoretical structural study of 2,9-bis-(4-fluorophenyl)-azaadamantan-4-one ethylene ketal

The title ethylene ketal 1, C₂₃H₂₃F₂NO₂, has a tricyclic azaadamantane cage with aryl rings (4-fluorophenyl) substituted at the 2 and 9 positions. Rings A (C11-C16) and B (C17-C22) are attached to atoms C2 and C9 in axial and equatorial positions, respectively; and the presence of the five-membered 1,3-dioxolane ring attached to atom C4. The structural analysis of the title was carried out by experimental (single crystal X-ray diffraction and NMR) and confirmed by theoretical calculations (ab initio RHF, DFT and NMR-DFT-GIAO).     

Structure of Au8: planar or nonplanar?

We have studied the structures and stabilities of Au6 and Au8 at the density-functional theory (DFT) and ab initio correlated levels of theory. For Au8, our ab initio calculations predict the lowest Au8 isomer to be planar, in line with the DFT calculations.     

NMR crystallography of p-tert-butylcalix[4]arene host-guest complexes using 1H complexation-induced chemical shifts

(1)H magic-angle spinning (MAS) NMR spectra of p-tert-butylcalix[4]arene inclusion compounds with toluene and pyridine show large complexation-induced shifts of the guest proton resonances arising from additional magnetic shielding caused by the aromatic rings of the cavities of the host calixarene lattice. In combination with ab initio calculations, these observations can be employed for NMR crystallography of host-guest complexes, providing important spatial information about the location of the guest molecules in the host cavities.     

Multireference ab initio quantum mechanics/molecular mechanics study on intermediates in the catalytic cycle of cytochrome P450(cam)

We have investigated the elusive reactive species of cytochrome P450(cam) (Compound I), the hydroxo complex formed during camphor hydroxylation, and the ferric hydroperoxo complex (Compound 0) by combined quantum mechanical/molecular mechanical (QM/MM) calculations, employing both density functional theory (DFT) and correlated ab initio methods. The first two intermediates appear multiconfigurational in character, especially in the doublet state and less so in the quartet state. DFT(B3LYP)/MM calculations reproduce the relative energies from correlated ab initio QM/MM treatments quite well, except for the splitting of the lowest A(1u)-A(2u) radical states. The inclusion of dynamic correlation is crucial for the proper ab initio treatment of these intermediates.     

Electronic structure of the PYP chromophore in its native protein environment

We report on supermolecular ab initio calculations which clarify the role of the local amino acid environment in determining the unique electronic structure properties of the photoactive yellow protein (PYP) chromophore. The extensive ab initio calculations, at the level of the CC2 and EOM-CCSD methods, allow us to explicitly address how the interactions between the deprotonated p-coumaric thio-methyl ester (pCTM-) chromophore and the surrounding amino acids act together to create a specifically stabilized pCTM- species. Particularly noteworthy is the role of the Arg52 amino acid in stabilizing the chromophore against autoionization, and the role of the Tyr42 and Glu46 amino acids in determining the hydrogen-bonding properties that carry the dominant energetic effects.     

Ab initio and DFT studies of the vibrational spectra of benzofuran and some of its derivatives

The vibrational spectra of benzofuran and some of its derivatives have been systematically investigated by ab initio and density functional B3LYP methods. The harmonic vibrational wavenumbers and intensity of vibrational bands were calculated at ab initio and DFT levels invoking different basis sets up to 6-311++g**. Vibrational assignments have been made and it has been found that the calculated DFT frequencies agree well in most cases with the observed frequencies for each molecule. Conformational studies have also been carried out and it is evident from ab initio calculations that 2(3H) benzofuranone is more stable than 3(2H) benzofuranone in support to our earlier semiempirical results.     

Attacking boron nucleophiles: NMR properties of five-membered diazaborole rings

This paper reports computed NMR spectral data for the diazaborole anion (C2H4B1N2-) and the corresponding neutral five-membered rings with B-H (C2H5B1N2, diazaborole) and B-Li (C2H4B1Li1N2, Li-diazaborole) bonds, which are the central moieties of newly synthesized nucleophilic organoboryl five-membered rings, recently reported by Segawa et al. (Science 2006, 314, 113). Both spin-spin coupling constants and chemical shifts were obtained using high-level ab initio calculations. These data are a necessary complement to the very scarce experimental information available.     

Ab initio studies of NMR chemical shifts for calix[4]arene and its derivatives

Ab initio calculations are performed for the calix[4]arene (1) and its derivatives (2 and 3), in this study. ¹H and ¹³C NMR measured spectral data given in our previous work are used to elucidate the structures of the prepared calix[4]arenes (1–3). The molecular geometry and chemical shift are calculated by using ab initio calculations based on the Hartree-Fock (HF) and the density functional theory (DFT) in the ground state. The results obtained from both methods are in agreement with the experimental results. The results of molecular geometry and chemical shifts show that DFT approach is closer to the experimental data than HF method.