Quantum-chemical and empirical calculations on phospholipids. III. Hydration of the dimethylphosphate anion

The binding of water to the dimethylphosphate anion (DMP^?) was calculated using the PCILO method. We found binding energies of 25.95 kcal·mol^?1 in the O1-P-O3 plane bridging the anionic oxygen atoms and 19.3 kcal·mol^?1 for the one-site association of a water molecule to an anionic oxygen atom of DMP^?. In this range one water molecule added to DMP^? in the O1 … O3 bridged configuration has a significantly higher binding energy to DMP^? than water molecules added to other binding sites. The total binding energy of 5 water molecules to DMP^? is 92 kcal·mol^?1, a quantity which is about 10% less than the sum of the binding energies of the corresponding monohydrates.     

Conformational diversity of the quercetin molecule: a quantum-chemical view

Abstract

This paper focuses on the comprehensive conformational analysis of the quercetin molecule with a broad range of the therapeutic and biological actions. All possible conformers of these molecule, corresponding to the local minima on the potential energy hypersurface, have been obtained by the sequential rotation of its five hydroxyl groups and also by the rotation of its (A?+?C) and B rings relatively each other. Altogether, it was established 48 stable conformers, among which 24 conformers possess planar structure and 24 conformers – nonplanar structure. Their structural, symmetrical, energetical and polar characteristics have been investigated in details. Quantum-mechanical calculations indicate that conformers of the quercetin molecule are polar structures with a dipole moment, which varies within the range from 0.35 to 9.87 Debay for different conformers. Relative Gibbs free energies of these conformers are located within the range from 0.0 to 25.3?kcal·mol^?1 in vacuum under normal conditions. Impact of the continuum with ε?=?4 leads to the decreasing of the Gibbs free energies (–0.19–18.15?kcal·mol^?1) and increasing of the dipole moment (0.57–12.48?D). It was shown that conformers of the quercetin molecule differ from each other by the intramolecular specific contacts (two or three), stabilizing all possible conformers of the molecule – H-bonds (both classical ОН…О and so-called unusual С′Н…О and ОН…С′) and attractive van-der-Waals contacts О…О. Obtained conformational analysis for the quercetin molecule enables to provide deeper understanding of the ‘structure-function’ relationship and also to suggest its mechanisms of the therapeutic and biological actions.

Communicated by Ramaswamy H. Sarma     

A theoretical investigation on the conformation and the interaction of CHF2OCF2CHF2 (desflurane II) with one water molecule

The conformation and the interaction of CHF₂OCF₂CHF₂ (desflurane II) with one water molecule is investigated theoretically using the ab initio MP2/aug-cc-pvdz and DFT-based M062X/6-311++G(d,p) methods. The calculations include the optimized geometries, the harmonic frequencies of relevant vibrational modes along with a natural bond orbital (NBO) analysis including the NBO charges, the hybridization of the C atom and the intra- and intermolecular hyperconjugation energies. In the two most stable conformers, the CH bond of the F₂HCO- group occupies the gauche position. The hyperconjugation energies are about the same for both conformers and the conformational preference depends on the interaction between the non-bonded F and H atoms. The deprotonation enthalpies of the CH bonds are about the same for both conformers, the proton affinity of the less stable conformer being 3 kcal mol^?1 higher. Both conformers of desflurane II interact with water forming cyclic complexes characterized by CH…O and OH…F hydrogen bonds. The binding energies are moderate, ranging from ?2.4 to ?3.2 kcal mol^?1 at the MP2 level. The origin of the blue shifts of the ν(CH) vibrations is analyzed. In three of the complexes, the water molecule acts as an electron donor. Interestingly, in these cases a charge transfer is also directed to the non bonded OH group of the water molecule. This effect seems to be a property of polyfluorinated ethers.     

Base Pairing at the Stem-loop Junction in the SL1 Kissing Complex of HIV-1 RNA: A Thermodynamic Study Probed by Molecular Dynamics Simulation

Abstract

The thermodynamics of the opening/closure process of a GC base pair located at the stem-loop junction of the SL1 sequence from HIV-1_Lai genomic RNA was investigated in the context of a loop-loop homodimer (or kissing complex) using molecular dynamics simulation. The free energy, enthalpy and entropy changes for the closing reaction are 0 kcal·mol^?1, ?11 kcal·mol^?1and ?0.037 kcal·mol^?1-K^?1 at 300° K respectively. Furthermore it is found that the free energy change is the same for the formation of a 11 nucleotide loop closed with UG and for the formation of a 9 nucleotide loop closed with GC. Our study evidences the high flexibility of the nucleotides at the stem-loop junction explaining the weak stability of this structure.     

Quantum chemical study of the structure and properties of citrinin

Detailed structures and electronic properties of three tautomeric forms of the toxin citrinin were investigated using several quantum calculation methods. Energetic preference of the predominant p- and o-quinone methide tautomeric forms is dependent on the method of calculation. A previously unstudied carboxylic acid enol tautomer was calculated to be surprisingly stable in vacuo, being within 2.5 kcal mol^? 1 at the B3LYP/6-311++G(2d,2p) level of theory. Despite differences in bond nature and connectivity of tautomers, the natural bond orbital analysis revealed that tautomeric forms share similar natural charges and natural electron configurations. Calculated bond lengths corresponded with experimentally observed values and assignments for the calculated infrared vibrational frequencies are reported.     

Theoretical studies on the tautomerism and intramolecular hydrogen shifts of 5-Amino-tetrazole in the gas phase

The tautomerism and intramolecular hydrogen shifts of 5-amino-tetrazole in the gas phase were studied in the present work. The minimum energy path (MEP) information of 5-amino-tetrazole was obtained at the CCSD(T)/6–311G**//MP2/6–311G** level of theory. The six possible tautomers of 1H, 4H-5-imino-tetrazole (a), 1H-5-amino-tetrazole (b), 2H-5-amino-tetrazole (c), 1H, 2H-5-imino-tetrazole (d), the mesoionic form (e) and 2H, 4H-5-imino-tetrazole (f) were investigated. Among these tautomers, there are 2 amino- forms, 3 imino- forms, and 1 mesoionic structure form. In all the tautomers, 2-H form (c) is the energetically preferred one in the gas phase. In the imino- tautomers, the energy value of the compound d is similar as that of the compound f but it is higher than the energy value of the compound a. The potential energetic surface (PES) and kinetics for five reactions have been investigated. Reaction 2 (b→c) was hydrogen shifts only in which the 1-H and 2-H rearrangement. This means that the reaction 2 (b→c) is energetically favorable having an activation barrier of 45.66 kcal·mol⁻¹ and the reaction energies (ΔE) is only 2.67 kcal·mol⁻¹. However, the reaction energy barrier for tautomerism of reaction 1 (b→e) is 54.90 kcal·mol⁻¹. Reaction 1 (b→a), reaction 3 (c→d), and reaction 5 (c→f) were amino- →imino- tautomerism reactions. The energy barriers of amino- →imino- tautomerism reactions required are 59.39, 65.57, 73.61 kcal·mol⁻¹ respectively in the gas phase. The calculated values of rate constants using TST, TST/Eckart, CVT, CVT/SCT and CVT/ZCT methods using the optimized geometries obtained at the MP2/6–311G** level of theory show the variational effects are small over the whole temperature range, while tunneling effects are big in the lower temperature range for all tautomerism reactions. Graphical Abstract Figure (DOC 45.0 KB)     

A theoretical investigation of the characteristics of hydrogen/halogen bonding interactions in dibromo-nitroaniline

In this work, computations of density functional theory (DFT) were carried out to investigate the nature of interactions in solid 2,6-dibromo-4-nitroaniline (DBNA). This system was selected to mimic the hydrogen/halogen bonding found within crystal structures as well as within biological molecules. DFT (M06-2X/6-311++G**) calculations indicated that the binding energies for different of interactions lie in the range between ?1.66 and ?9.77 kcal mol^?1. The quantum theory of atoms in molecules (QTAIM) was applied to provide more insight into the nature of these interactions. Symmetry-adapted perturbation theory (SAPT) analysis indicated that stability of the Br···Br halogen bonds is predicted to be attributable mainly to dispersion, while electrostatic forces, which have been widely believed to be responsible for these types of interactions, play a smaller role. Our results indicate that, for those nuclei participating in hydrogen/halogen bonding interactions, nuclear quadrupole resonance parameters exhibit considerable changes on going from the isolated molecule model to crystalline DBNA.

Prototropic tautomerism and basic molecular principles of hypoxanthine mutagenicity: an exhaustive quantum-chemical analysis

The molecular structures, relative stability order, and dipole moments of a complete family of 21 planar hypoxanthine (Hyp) prototropic molecular–zwitterionic tautomers including ylidic forms were computationally investigated at the MP2/6–311++G(2df,pd)//B3LYP/6–311++G(d,p) level of theory in vacuum and in three different surrounding environments: continuum with a low dielectric constant (??=?4) corresponding to a hydrophobic interface of protein–nucleic acid interactions, dimethylsulfoxide (DMSO), and water. The keto-N1HN7H tautomer was established to be the global minimum in vacuum and in continuum with ??=?4, while Hyp molecule exists as a mixture of the keto-N1HN9H and keto-N1HN7H tautomers in approximately equal amounts in DMSO and in water at T?=?298.15?K. We found out that neither intramolecular tautomerization by single proton transfer in the Hyp base, nor intermolecular tautomerization by double proton transfer in the most energetically favorable Hyp·Hyp homodimer (symmetry C _2h ), stabilized by two equivalent N1H…O6 H-bonds, induces the formation of the enol tautomer (marked with an asterisk) of Hyp with cis-oriented O6H hydroxyl group relative to neighboring N1C6 bond. We first discovered a new scenario of the keto–enol tautomerization of Hyp?·?Hyp homodimer (C _2h ) via zwitterionic near-orthogonal transition state (TS), stabilized by N1⁺H…N1^? and O6⁺H…N1^? H-bonds, to heterodimer Hyp^??·?Hyp (C _s ), stabilized by O6H…O6 and N1H…N1 H-bonds. We first showed that Hyp^??·?Thy mispair (C _s ), stabilized by O6H…O4, N3H…N1, and C2H…O2 H-bonds, mimicking Watson–Crick base pairing, converts to the wobble Hyp?·?Thy base pair (C _s ), stabilized by N3H…O6 and N1H…O2 H-bonds, via high- and low-energy TSs and intermediate Hyp?·?Thy^?, stabilized by O4H…O6, N1H…N3, and C2H…O2 H-bonds. The most energetically favorable TS is the zwitterionic pair Hyp⁺?·?Thy^? (C _s ), stabilized by O6⁺H…O4^?, O6⁺H…N3^?, N1⁺H…N3^?, and N1⁺H…O2^? H-bonds. The authors expressed and substantiated the hypothesis, that the keto tautomer of Hyp is a mutagenic compound, while enol tautomer Hyp^? does not possess mutagenic properties. The lifetime of the nonmutagenic tautomer Hyp^? exceeds by many orders the time needed to complete a round of DNA replication in the cell. For the first time purine–purine planar H-bonded mispairs containing Hyp in the anti-orientation with respect to the sugar moiety – Hyp?·?Ade _syn , Hyp?·?Gua^? _syn , and Hyp?·?Gua _syn , that closely resembles the geometry of the Watson–Crick base pairs, have been suggested as the source of transversions. An influence of the surrounding environment (??=?4) on the stability of studied complexes and corresponding TSs was estimated by means of the conductor-like polarizable continuum model. Electron-topological, structural, vibrational, and energetic characterictics of all conventional and nonconventional H-bonds in the investigated structures are presented. Presented data are key to understanding elementary molecular mechanisms of mutagenic action of Hyp as a product of the adenine deamination in DNA.     

Zwitterionic conformers of pyrrolysine and their interactions with metal ions—a theoretical study

A total of 16 pyrrolysine conformers in their zwitterionic forms are studied in gas and simulated aqueous phase using a polarizable continuum model (PCM). These conformers are selected on the basis of our study on the intrinsic conformational properties of non-ionic pyrrolysine molecule in gas phase [Das and Mandal (2013) J Mol Model 19:1695?1704]. In aqueous phase, the stable zwitterionic pyrrolysine conformers are characterized by full geometry optimization and vibrational frequency calculations using B3LYP/6-311++G(d,p) level of theory. Single point calculations are also carried out at MP2/6-311++G(d,p) level. Characteristic intramolecular hydrogen bonds present in each conformer, their relative energies, theoretically predicted vibrational spectra, rotational constants and dipole moments are systematically reported. The calculated relative energy range of the conformers at B3LYP/6-311++G(d,p) level is 5.19 kcal mol^?1 whereas the same obtained by single point calculations at MP2/6-311++G(d,p) level is 4.58 kcal mol^?1. A thorough analysis reveals that four types of intramolecular H-bonds are present in the conformers; all of which play key roles in determining the energetics and in imparting the observed conformations to the conformers. The vibrational frequencies are found to shift invariably toward the lower side of frequency scale corresponding to the presence of the H-bonds. This study also points out that conformers with diverse structural motifs may differ in their thermodynamical stability by a narrow range of relative energy. The effects of metal coordination on the relative stability order and structural features of the conformers are examined by complexing five zwitterionic conformers of pyrrolysine with Cu⁺² through their carboxylate groups. The interaction enthalpies and Gibbs energies, rotational constants, vibrational frequencies and dipole moments of the metal complexes calculated at B3LYP level are also reported. The zwitterionic conformers of pyrrolysine are not stable in gas phase; after geometry optimization they are converted to the non-ionic forms.     

DFT comparison of the OH-initiated degradation mechanisms for five chlorophenoxy herbicides

To compare the OH-initiated reaction mechanisms of five chlorophenoxy herbicides, density functional theory (DFT) calculations of reactions in which ·OH attacks one of three active positions on each herbicide were carried out at the MPWB1K/6-311 + G(3df,2p)//MPWB1K/6-31 + G(d,p) level. For each herbicide, the calculation results show that ·OH addition to the C1 atom, which is the nexus between the benzene ring and the side group, possesses the lowest energy barrier among the three kinds of reactions, indicating that ·OH addition–substitution of the side chain is the most energetically and kinetically favorable reaction mechanism. Comparisons among the herbicides show that the mechanisms are affected by the steric hindrance and the electronegativities of the –CH₃ and –Cl groups. When comparing the addition of ·OH to the C1 site among the five herbicides, the activation energy for the reaction of ·OH with DCPP reaction is the lowest (3.61 kcal mol^?1), while that for the ·OH and 4-CPA reaction was the highest (5.91 kcal mol^?1). ·OH addition to the C4 site presents the highest energy barriers among the three kinds of reactions, indicating that the para Cl is difficult to break down. When comparing the H-atom abstraction reactions of the five herbicides, the H atoms in the –CH₂– group of 2,4-D are the easiest for ·OH to abstract, whereas those of DCPP and MCPP are more difficult to abstract, due to the steric hindrance of the –CH₃ group. Additionally, the results obtained from the PCM calculations reveal that most of the reactions occur more easily in water than in gas, though the mechanisms involved are the same as those discussed above.     

DFT investigation on the decarboxylation mechanism of <Emphasis Type="Italic">ortho</Emphasis> hydroxy benzoic acids with acid catalysis

A density functional theory (DFT) study was performed to explore the mechanisms of the acid-catalyzed decarboxylation reaction of salicylic acids using the B3LYP method with 6-31++G(d,p) basis set in both gas phase and aqueous environment. The α-protonated cation of carboxylate acid was formed during the decarboxylation process in acidic conditions, and the presence of hydrogen ions promotes decarboxylation greatly by significantly decreasing the overall reaction energy barriers to 20.98 kcal mol^?1 in gas phase and 20.93 kcal mol^?1 in water, respectively. The hydrogen in the α-carbon came directly from the acid rather than from the carboxyl group in neutral state. Compared with the reaction in gas phase, water in aqueous state causes the reaction to occur more easily. Substituents of methyl group, chlorine and fluorine at the ortho-position to the carboxyl of salicylic acid could further lower the decarboxylation energy barriers and facilitate the reaction.     

Kinetic and thermodynamic studies on nitrobenzylthioinosine binding to the nucleoside transporter of chinese hamster ovary cells

The binding of [G-³H]nitrobenzylthioinosine to intact Chinese hamster ovary cells has been studied kinetically and thermodynamically. The association of nitrobenzylthioinosine with cells is a second-order process which proceeds at 24°C with a rate constant of 2·10⁷ M^?1·s^?1. Dissociation of the complex was characterized as a simple first-order process with rate constant on the order of 7·10^?3 s^?1. The quotient of these is comparable to the dissociation constant as measured in equilibrium binding studies, 2.2·10^?10 M. The temperature dependence of the rate of association indicated an Arrhenius activation energy of 8.4 kcal·mol^?1, while that of the equilibrium constant for dissociation indicated a standard enthalpy change of 8.8 kcal·mol^?1. The large increase in affinity of nitrobenzylthioinosine as compared to natural nucleosides is attributable to an entropy-driven interaction with the binding site. Thymidine, dipyridamole and papaverine each decrease the apparent dissociation constant for the nitrobenzylthioinosine-cell complex; the latter, inhibitors of nucleoside transport, decrease the rate of dissociation of the complex.     

Anesthetic-protein interaction Random versus helix polylysine monolayers and interaction with 1-alkanols

Penetration of 1-alkanols into monolayers of hydrophobic polypeptides, poly(ε-benzyloxycarbonyl-l-lysine) and poly(ε-benzyloxycarbonyl-dl-lysine), was compared with their adsorption on the air/water interface in the absence of monolayers. The polypeptide prepared from l-lysine is generally considered to be in the α-helical form whereas dl-copolymer polypeptide contains random-coiled portions due to the structural incompatibility between the two isomers. The free energy of adsorption of 1-alkanols on the air/water interface at dilute concentrations was ?0.68 kcal·mol^?1 per methylene group and 0.15 kcal·mol^?1 for the hydroxyl group at 25°C. In the close-packed state, the surface area occupied by each molecule of 1-alkanols of varying carbon chain-lengths showed nearly a constant value of about 27.2 Ų, indicating perpendicular orientation of the alkanol molecules at the interface. About 75% of the water surface was covered by 1-butanol in this close-packed state. The mode of adsorption of 1-alkanols on the vacant air/water interface followed the Gibbs surface excess while the mode on the polypeptide membranes followed the Langmuir adsorption isotherm, indicating that the latter is characterized by the presence of a finite number of binding sites. The free energies of adsorption of 1-alkanols on the l-polymer monolayers were more negative than those on the vacant air/water interface and less negative than those on the dl-copolymer monolayers. Thus, the affinity of 1-alkanols to the interface was in the order of vacant air/water interface <l-polymer <dl-copolymer. The difference between the air/water interface and l-polymer was about 0.54 kcal·mol^?1 and that between l-polymer and dl-copolymer was 0.17 kcal·mol^?1 at 25°C: the adsorption of 1-alkanols to the dl-copolymer was favored compared to the l-polymer. The polar moieties of the backbone of the dl-copolymer may be exposed to the aqueous phase at the disordered portion. Dipole interaction between this portion and 1-alkanol molecules may account for the enhanced adsorption of the alkanols to the dl-copolymer.     

DPT tautomerization of the long A∙A* Watson-Crick base pair formed by the amino and imino tautomers of adenine: combined QM and QTAIM investigation

Combining quantum-mechanical (QM) calculations with quantum theory of atoms in molecules (QTAIM) and using the methodology of sweeps of the energetic, electron-topological, geometric and polar parameters, which describe the course of the tautomerization along the intrinsic reaction coordinate (IRC), we showed for the first time that the biologically important A?A* base pair (C_s symmetry) formed by the amino and imino tautomers of adenine (A) tautomerizes via asynchronous concerted double proton transfer (DPT) through a transition state (TS), which is the A⁺?A^? zwitterion with the separated charge, with C_s symmetry. The nine key points, which can be considered as electron-topological “fingerprints” of the asynchronous concerted A?A*?A*?A tautomerization process via the DPT, were detected and completely investigated along the IRC of the A?A*?A*?A tautomerization. Based on the sweeps of the H-bond energies, it was found that intermolecular antiparallel N6Н?N6 (7.01 kcal mol^?1) and N1H?N1 (6.88 kcal mol^?1) H-bonds are significantly cooperative and mutually reinforce each other. It was shown for the first time that the A?A*?A*?A tautomerization is assisted by the third C2H?HC2 dihydrogen bond (DHB), which, in contrast to the two others N6H?N6 and N1H?N1 H-bonds, exists within the IRC range from ?2.92 to 2.92 Å. The DHB cooperatively strengthens, reaching its maximum energy 0.42 kcal mol^?1 at IRC?=??0.52 Å and minimum energy 0.25 kcal mol^?1 at IRC?=??2.92 Å, and is accompanied by strengthening of the two other aforementioned classical H-bonds. We established that the C2H?HC2 DHB completely satisfies the electron-topological criteria for H-bonding, in particular Bader’s and all eight “two-molecule” Koch and Popelier’s criteria. The positive value of the Grunenberg’s compliance constant (5.203 Å/mdyn) at the TS_A?A*?A*?A proves that the C2H?HC2 DHB is a stabilizing interaction. NBO analysis predicts transfer of charge from σ(C2–H) bonding orbital to σ*(H–C2) anti-bonding orbital; at this point, the stabilization energy E⁽²⁾ is equal to 0.19 kcal mol^?1 at the TS_A?A*?A*?A.     

Racemization of the gastrointestinal antisecretory chiral drug esomeprazole magnesium via the pyramidal inversion mechanism: A theoretical study

The pyramidal inversion mechanisms of the 6‐methoxy and the 5‐methoxy tautomers of (S)‐omeprazole were studied, employing ab initio and DFT methods. The conformational space of the model molecule (S)‐2‐[(3‐methyl‐2‐pyridinyl)methyl]sulfinyl‐1H‐benzimidazole was calculated, with respect to rotations around single bonds, at the B3LYP/6‐311G(d,p) level. All of the resulting conformations were used as starting points for full optimizations of (S)‐omeprazole, at B3LYP/6‐31G(d), B3LYP/6‐311G(d,p), B3LYP/6‐311++G(d,p), B3LYP/6‐311G(2df,2pd), MP2/6‐31G(d), and MP2/6‐311G(d,p) levels. Four distinct pathways were found for enantiomerization via the pyramidal inversion mechanism for each of the tautomers of (S)‐omeprazole. Each transition state, in which the sulfur, the oxygen and the two carbon atoms connected directly to the sulfur are in one plane, connects two diastereomeric minima. The enantiomerization is completed by free rotation around the sulfur–methylene bond, and around the methylene–pyridine ring bond. The effective Gibbs' free energy barrier for racemization ΔG of the two tautomers of (S)‐omeprazole are 39.8 kcal/mol (5‐methoxy tautomer) and 40.0 kcal/mol (6‐methoxy tautomer), indicating that the enantiomers of omeprazole are stable at room temperature (in the gas phase). The 5‐methoxy tautomer of (S)‐omeprazole was found to be slightly more stable than the 6‐methoxy tautomer, in the gas phase. The energy barrier (ΔG^?) for the(S,M) (S,P) diastereomerization of (S)‐omeprazole due to the rotation around the pyridine chiral axis was very low, 5.8 kcal/mole at B3LYP/6‐311G(d,p). Chirality 2010. © 2010 Wiley‐Liss, Inc.     

Hybrid MC/QC simulations of water-assisted proton transfer in nucleosides. Guanosine and its analog acyclovir

To provide an in-depth insight into the molecular basis of spontaneous tautomerism in DNA and RNA base pairs, a hybrid Monte Carlo (MC)–quantum chemical (QC) methodology is implemented to map two-dimensional potential energy surfaces along the reaction coordinates of solvent-assisted proton transfer processes in guanosine and its analog acyclovir in aqueous solution. The solvent effects were simulated by explicit inclusion of water molecules that model the relevant part of the first hydration shell around the solute. The position of these water molecules was estimated by carrying out a classical Metropolis Monte Carlo simulation of dilute water solutions of the guanosine (Gs) and acyclovir (ACV) and subsequently analyzing solute–solvent intermolecular interactions in the statistically-independent MC-generated configurations. The solvent-assisted proton transfer processes were further investigated using two different ab initio MP2 quantum chemical approaches. In the first one, potential energy surfaces of the ‘bare’ finite solute–solvent clusters containing Gs/ACV and four water molecules (MP2/6-31+G(d,p) level) were explored, while within the second approach, these clusters were embedded in ‘bulk’ solvent treated as polarizable continuum (C-PCM/MP2/6-31+G(d,p) level of theory). It was found that in the gas phase and in water solution, the most stable tautomer for guanosine and acyclovir is the 1H-2-amino-6-oxo form followed by the 2-amino-6-(sZ)-hydroxy form. The energy barriers of the water-assisted proton transfer reaction in guanosine and in acyclovir are found to be very similar – 11.74 kcal mol^?1 for guanosine and 11.16 kcal mol^?1 for acyclovir, and the respective rate constants (k = 1.5?×?10¹ s^?1, guanosine and k = 4.09?×?10¹ s^?1, acyclovir), are sufficiently large to generate the 2-amino-6-(sZ)-hydroxy tautomer. The analysis of the reaction profiles in both compounds shows that the proton transfer processes occur through the asynchronous concerted mechanism.     

Theoretical study of the geometrical and electronic structures and thermochemistry of spherophanes

A set of supramolecular cage-structures—spherophanes—was studied at the density functional B3LYP level. Full geometrical structure optimisations were made with 6–31G and 6–31G(d) basis sets followed by frequency calculations, and electronic energies were evaluated at B3LYP/6–31++G(d,p). Three different symmetries were considered: C1, Ci, and Oh. It was found that the bonds between the benzene rings are very long to allow π-electron delocalisation between them. These spherophanes show portal openings of 2.596 Å in Spher1, 4.000 Å in Meth2, 3.659 Å in Oxa3, and 4.412 Å in Thia4. From the point of view of potential host–guest interaction studies, it should also be noted that the atoms nearest to the centre of the cavities are carbons bonded to X groups. These supramolecules seem to exhibit relatively large gap HOMO?LUMO: 2.89 eV(Spher1), 5.26 eV(Meth2), 5.73 eV(Oxa3), and 4.82 eV(Thia4). The calculated ΔH°_f (298.15 K) values at B3LYP/6–31G(d) are (in kcal mol^?1) 750.98, 229.78, ?10.97, and 482.49 for Spher1, Meth2, Oxa3, and Thia4, respectively. Using homodesmotic reactions, relative to Spher1, the spherophanes Meth2, Oxa3, and Thia4 are less strained by ?399.13 kcal mol^?1, ?390.40 kcal mol^?1, and ?411.38 kcal mol^?1, respectively. Their infrared and ¹³C NMR calculated spectra are reported.     

A mechanistic investigation on the formation and rearrangement of silaspiropentane: A theoretical study

The formation of silaspiropentane from addition of singlet silacyclopropylidene 1 and silacyclopropylidenoid 8 to ethylene has been investigated separately at the B3LYP, X3LYP, WB97XD, and M05–2X theories using the 6–31+G(d,p) basis set. The silacycloproylidenoid addition follows a stepwise route. In contrast, a concerted mechanism occurs for silacyclopropylidene addition. Moreover, the intramolecular rearrangements of silaspiropentane 9 to methylenesilacyclobutane 11 and 2-silaallene?+?ethylene 12 have been studied extensively. The required energy barrier for the isomerization of 9 to 10 was determined to be 44.0 kcal mol^?1 at the B3LYP/6–31+G(d,p) level. After formation of 10, the rearrangement to methylenesilacyclobutane 12 is highly exergonic by ?15.9 kcal mol^?1, which makes this reaction promising. However, the conversion of 9 to 11 is calculated to be quite endergonic, by 26.5 kcal mol^?1.     

Theoretical studies of the interaction between enflurane and water

Increase of the atmospheric concentration of halogenated organic compounds is partially responsible for a change of the global climate. In this work we have investigated the interaction between halogenated ether and water, which is one of the most important constituent of the atmosphere. The structures of the complexes formed by the two most stable conformers of enflurane (a volatile anaesthetic) with one and two water molecules were calculated by means of the counterpoise CP-corrected gradient optimization at the MP2/6–311++G(d,p) level. In these complexes the CH…O_w hydrogen bonds are formed, with the H…O_w distances varying between 2.23 and 2.32 Å. A small contraction of the CH bonds and the blue shifts of the ν(CH) stretching vibrations are predicted. There is also a weak interaction between one of the F atoms and the H atom of water, with the H_w…F distances between 2.41 and 2.87 Å. The CCSD(T)/CBS calculated stabilization energies in these complexes are between ?5.89 and ?4.66 kcal?mol^?1, while the enthalpies of formation are between ?4.35 and ?3.22 kcal?mol^?1. The Cl halogen bonding between enflurane and water has been found in two complexes. The intermolecular (Cl···O) distance is smaller than the sum of the corresponding van der Waals radii. The CCSD(T)/CBS stabilization energies for these complexes are about ?2 kcal?mol^?1.