Structures of inclusion complexes of halogenbenzoic acids and alpha-cyclodextrin based on AM1 calculations

Quantitative ROESY analysis of computational models: structural studies of citalopram and β-cyclodextrin complexes by (1) H-NMR and computational methods

Complexation of racemic citalopram with β-cyclodextrin (β-CD) in aqueous medium was investigated to determine atom-accurate structure of the inclusion complexes. (1) H-NMR chemical shift change data of β-CD cavity protons in the presence of citalopram confirmed the formation of 1 : 1 inclusion complexes. ROESY spectrum confirmed the presence of aromatic ring in the β-CD cavity but whether one of the two or both rings was not clear. Molecular mechanics and molecular dynamic calculations showed the entry of fluoro-ring from wider side of β-CD cavity as the most favored mode of inclusion. Minimum energy computational models were analyzed for their accuracy in atomic coordinates by comparison of calculated and experimental intermolecular ROESY peak intensities, which were not found in agreement. Several least energy computational models were refined and analyzed till calculated and experimental intensities were compatible. The results demonstrate that computational models of CD complexes need to be analyzed for atom-accuracy and quantitative ROESY analysis is a promising method. Moreover, the study also validates that the quantitative use of ROESY is feasible even with longer mixing times if peak intensity ratios instead of absolute intensities are used.

Computational Modeling of Inclusion Complexes of β-Cyclodextrin with enantiomers of Salsolinol, N-Methyl-Salsolinol, and 1-Benzyl-Tetrahydroisoquinoline

Capillary electrophoresis with β-CD as a chiral selector has successfully separated the two enantiomers of salsolinol, N-methyl-salsolinol, and 1-benzyl-tetrahydroisoquinoline. The migration times of each enantiomer in capillary electrophoresis reflect the stability of their β-CD inclusion complexes. This paper reports a computational modeling study of the inclusion complexes of β-cyclodextrin (β-CD) with salsolinol, N-methyl-salsolinol, and 1-benzyl-tetrahydroisoquinoline by using PM3 (Parametric Method 3) semi-empirical molecular orbital calculations and the ONIOM hybrid method. Two types of the inclusion complexes, cis- and trans-orientations, are considered for each enantiomer of the guest molecules, salsolinol, N-methyl-salsolinol, and 1-benzyl-tetrahydroisoquinoline. In the cis-orientation, the nitrogen in the salsolinol, N-methyl-salsolinol, and 1-benzyl-tetrahydroisoquinoline points toward the secondary hydroxyls of the β-CD, while in the trans-orientation, the nitrogen in salsolinol, N-methyl-salsolinol, and 1-benzyl-tetrahydroisoquinoline points toward the primary hydroxyls of the β-CD. We found that the stabilization energies of these inclusion complexes from these PM3 and ONIOM different methods correlate very well with the migration order deduced from the study of capillary electrophoretic separation.

Theoretical and experimental study of the inclusion complexes of the 3-carboxy-5,6-benzocoumarinic acid with cyclodextrins

The 3-carboxy-5,6-benzocoumarinic acid (BzCum) presents important emission properties that could make it suitable to be used as a fluorescence marker for proteins. Since the guest-cyclodextrin interaction can be used as a model for the ligand-protein interaction, the host-guest complexes of BzCum with alpha-, beta-, gamma-, 2-hydroxypropyl-beta-cyclodextrin (2-HP-beta-CD) and 2-hydroxypropyl-gamma-cyclodextrin (2-HP-gamma-CD) were investigated by means of steady-state fluorescence measurements. As the guest may exist either in anionic or in neutral form, the experiments were performed at two different pH values, 9.48 and 2.11. The association constants and the stoichiometry of the complexes were estimated by non-linear regression analysis. A theoretical study of the complexes was made using molecular mechanics (MM) calculation methods. The electrostatic and van der Waals contributions to the energy were estimated and a discussion was made on this basis.

On the relative stabilities of the alkali cations 222 cryptates in the gas phase and in water-methanol solution

The relative stabilities of the alkali [M subset 222]+ cryptates (M = Na, K, Rb and Cs) in the gas phase and in solution (80:20 v/v methanol:water mixture) at 298 K, are computed using a combination of ab initio quantum-chemical calculations (HF/6-31G and MP2/6-31+G*//HF/6-31+G*) and explicit-solvent Monte Carlo free-energy simulations. The results suggest that the relative stabilities of the cryptates in solution are due to a combination of steric effects (compression of large ions within the cryptand cavity), electronic effects (delocalization of the ionic charge onto the cryptand atoms) and solvent effects (dominantly the ionic dessolvation penalty). Thus, the relative stabilities in solution cannot be rationalized solely on the basis of a simple match or mismatch between the ionic radius and the cryptand cavity size as has been suggested previously. For example, although the [K subset 222]+ cryptate is found to be the most stable in solution, in agreement with experimental data, it is the [Na subset 222]+ cryptate that is the most stable in the gas phase. The present results provide further support to the notion that the solvent in which supramolecules are dissolved plays a key role in modulating molecular recognition processes.

Quantum chemical study of the mechanism of ethylene elimination in silylative coupling of olefins

Silylative coupling of olefins differs from olefin metathesis. Although in both these reactions ruthenium catalysts play a crucial role and ethylene product is detected, ruthenium-carbene intermediate is formed only in the course of the metathesis reaction. In this study quantum chemical calculations based on the density functional theory (DFT) have been carried out in order to examine the mechanism of the silylative coupling of olefins leading to ethylene elimination. In the first step of the catalytic cycle, a hydrogen atom from the ruthenium catalytic center is transferred preferentially to the carbon atom bound to Si in a vinylsilane. This H transfer is coupled with the formation of Ru-C bond. Next, the rotation around the newly formed C-C single bond occurs so that silicon atom is placed in the vicinity of the ruthenium center. The following step involves the migration of a silyl moiety, and leads to Ru-Si bond formation, coupled with ethylene elimination. The next reaction, that is the insertion of ethylene (alkene) into Ru-Si bond, has an activation barrier almost as high as the reaction of ethylene elimination. However, the posibility of removing gaseous ethylene from the reactive mixture together with the entropic fators suggests that the insertion of alkene that is larger than C(2)H(4) is the rate limiting step in the silylative coupling of olefins. It also suggests that the substituents attached to the silicon atom or the carbon atoms of an alkene by electronic and steric effects may significantly affect silyl migration and thus the effectiveness of the catalytic reaction.