Combined ab Initio/DFT and Monte Carlo Calculation of Relative Standard Chemical Potentials in Solution

Diaryliodonium Salt-Mediated Intramolecular C-N Bond Formation Using Boron-Masking N-Hydroxyamides

Intramolecular aromatic C-N bond formation reactions using electron-rich aromatic tethered boron-masking N-hydroxyamide as substrate were realized. These new C-N bond formation reactions involve the in situ generation of a diaryliodonium salt by treatment with hypervalent iodine, deborylation by base treatment, spontaneous N → O acyl migration, cyclization, reductive elimination, elimination of benzoic acid, and tautomerization to indole formation. Hereby, we obtained highly functionalized electron-rich indoles and quinoline in practical yields.

In-solution conformational analysis of the XCYCH3 moiety for small esters and ethers with all combinations of X, Y = O, S

Favorable steric and electrostatic fit of a ligand to a receptor is of central interest in theoretical drug design. This paper considers the effects of non-protic solvents, in comparison with the gas phase, on the preferred conformation of the XCYCH3 moiety of simple aliphatic esters and heterocyclic methyl ethers with all combinations of the X and Y atoms as oxygen and sulfur. An IEF-PCM/B97D/aug-cc-pv(t+d)z continuum dielectric solvent study in chloroform and acetonitrile explores the through-space polarization effect of the environment on the conformational preference, not affected by possible solute-solvent hydrogen bond formation. The inherently favored structure for the present molecules is important, since the hypothetical oxygen and sulfur lone-pairs point approximately in opposite directions in the cis conformation of esters, whereas the trans and gauche conformations for the methyl group in ethers define nearly parallel or perpendicular directionality for the lone pairs of the ring heteroatoms and the O or S atoms connecting to the ring. These different preferences for the studied two families of compounds allow for designing formation of hydrogen bonds with a protein in fairly different regions of the latter still within the ligand-binding cavity. For a fine-tuning of these hydrogen bonds, a replacement of an oxygen atom of the ligand by a sulfur atom could be a straightforward possibility.

Are the intramolecular O-H···F and O-H···Cl hydrogen bonds maintained in solution? A theoretical study

The present case study aims at calculating the equilibrium conformer compositions for 2X-ethanol and 2X-phenol (X = F, Cl) in solution, and exploring the effect of the applied theoretical method and basis set on the obtained results, as well as considering the usefulness of the continuum solvent approach in comparison with the explicit solvent Monte Carlo model utilizing the free energy perturbation method. Gas-phase optimizations at the DFT/B97D/aug-cc-pvtz and ab initio MP2/aug-cc-pvtz levels predicted structures in good agreement with the available experimental data for three test molecules. Because in-solution geometries change only slightly according to the IEF-PCM continuum solvent calculations in carbon tetrachloride and water, the two theoretical levels were applied further on, and complete basis set (CBS) relative internal free energies were estimated for the conformers under study. The predicted OCCF gauche/trans ratio for 2F-ethanol was well reproduced in comparison with available experimental compositions. The predominant gauche structure maintains an intramolecular hydrogen bond in carbon tetrachloride (HB structure), whereas HB and NoHB gauche conformers appear in nearly the same fraction in aqueous solution. The internally hydrogen-bonded conformer is predominant also for 2X-phenol species, as calculated on the basis of relative CBS internal free energies and IEF-PCM and FEP/MC solvation free energies. Use of a trihydrate supermolecule model for 2F-ethanol conformers leads to the prediction of the aqueous-solution composition in contrast to the experiment. Solution structure modeling predicts weak hydrogen-bond formation capacity for both the covalently bound F and Cl atoms, even in conformations where they are fully exposed to hydration.

Theoretical study of the gauche-trans equilibrium with and without an intramolecular hydrogen bond for +H3NCH2CH2X systems (X = OH, NH2, COO-) in solution

The gauche-trans conformational equilibrium has been studied for (+)H(3)NCH(2)CH(2)X systems (X = OH, NH(2), COO(-)) and for the neutral model of the simplest β-amino acid in aqueous solution and chloroform. Each structure exhibits an intramolecular hydrogen bond in the NCCX gauche conformation, which is necessarily disrupted in the local energy minimum trans form. Geometries were optimized at the IEF-PCM/B97D/aug-cc-pvtz level of theory and indicated that the solvent effect of chloroform vs. water is moderate when the geometries for the corresponding gauche and trans conformers are compared. The only remarkable difference was found for β-alanine, which can exist in gauche, zwitterionic form in aqueous solution but not in chloroform. Total relative free energies were obtained as the sum of the relative IEF-PCM/B97D/aug-cc-pvtz internal free energy and ΔG(solv). The latter was calculated both by means of IEF-PCM and the explicit solvent FEP/Monte Carlo methods. The resulting ΔG(tot) values could differ by 1-2 kcal mol(-1) depending on the accepted ΔG(solv) value, but any calculation indicated that the internally bound gauche conformers are far more populated than the corresponding trans species. A difference of 2.6 kcal mol(-1) for ΔG(solv) by the two methods resulted in the preference of the zwitterionic gauche β-alanine structure by FEP/MC in contrast to IEF-PCM favoring the neutral form in aqueous solution. Monte Carlo characterization of the solution structure in the first solvation shell of polar groups indicates that solvation of a N-H(+) bond is sensitive to its involvement in an intramolecular hydrogen bond. For the zwitterionic, gauche β-alanine, almost no water oxygen can reach that N-H(+) bond, which is in strong interaction with the -COO(-) group by forming a six-membered ring and favorable local geometry for the hydrogen bond.

On the interaction of aliphatic amines and ammonium ions with carboxylic acids in solution and in receptor pockets

Association energies of the acetate ion with cationic amines bearing one to three methyl groups were calculated in the range of -14 to -17 kcal/mol in aqueous solution by means of the IEF-PCM method at the CCSD(T)/CBS//MP2/aug-cc-pvdz and DFT/B97D/CBS//B97D/aug-cc-pvtz levels. The main stabilization factor for the association is the possibility for the formation of an ionic intermolecular hydrogen bond between the elements of the complex. For a quaternary ammonium ion, the favorable electrostatic interaction energy is the only driving force, and the stabilization energy for the complex is reduced to -4 kcal/mol. The internal free energies of the ion-pair tautomers of the studied species are higher by 10-15 kcal/mol in water than those for the neutral, hydrogen-bonded forms. Monte Carlo free energy perturbation calculations at T = 298 K and p = 1 atm predict -11 to -16 kcal/mol relative solvation free energy in favor of the corresponding ionic form. As a result, the ion-pair tautomer is the prevailing form in aqueous solution and on the extracellular surface of a receptor. Modeling the complex of a protonated ligand interacting with an Asp/Glu carboxylate side-chain in the binding cavity of a receptor, two strongly bound water molecules were considered so as to form hydrogen-bonded water bridges between the elements of the ion-pair. Nonetheless, the low polarity environment mimicked by a chloroform solvent cannot stabilize the ionic tautomer. A proton jump was predicted, which suggests that acetylcholine, an inherent cation by structure, might have evolved as the natural agonist for muscarinic receptors because a quaternary ammonium system assures the maintenance of the ion-pair form with a carboxylate side-chain in a protein cavity, the latter perhaps then being needed for receptor activation.