Structure simulations for the 0.22 and 1 molar aqueous dimethylammonium chloride solutions

Solvation structure and dynamics of the dimethylammonium cation diluted in liquid water: A molecular dynamics approach

Classical molecular dynamics simulation techniques were employed to investigate the local solvation structure and related dynamics of the dimethylammonium cation diluted in liquid water at ambient conditions. The translational and orientational order around the dimethylammonium cation was investigated in terms of the corresponding radial and angular distribution functions. The results obtained revealed that the first solvation shell of the dimethylammonium consists mainly of two and, less frequently, three water molecules. The two nearest water neighbors form hydrogen bonds with the ammonium hydrogen atoms of the cation, whereas the third neighbor interacts with the methyl hydrogen atoms as well. The distribution of the trigonal order parameter exhibits a bimodal behavior, signifying the existence of local orientational heterogeneities in the solvation shell of the dimethylammonium cation. The calculated continuous and intermittent residence and hydrogen bond lifetimes for the cation-water pairs have also been found to be longer in comparison with the water-water ones. The very similar self-diffusion coefficients of the dimethylammonium cation and the water molecules in the bulk dilute solution indicate that the translational motions of the cation are mainly controlled by the translational mobility of the surrounding water molecules.

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.

Theoretical studies of the solvent effect on the conformation of the HO-C-C-X (X = F, NH2, NO2) moiety with competing intra- and intermolecular hydrogen bonds

Theoretical calculations up to the ab initio IEF-PCM/CCSD(T)/CBS//IEF-PCM/B3LYP/6-311++G** and IEF-PCM/B97D/aug-cc-pvtz levels have been performed for 2X-ethanol and 2X-phenol systems with X = F, NH(2), NO(2) in chloroform and aqueous solution. The calculated relative free energies by means of the IEF-PCM continuum dielectric method do not differ very much at the DFT and ab initio levels. Application of explicit solvent models and the FEP/MC method for determining relative solvation free energies causes, however, large deviations in the predicted equilibrium compositions, although the predominant conformation for the solute is generally in agreement with that from the corresponding IEF-PCM calculations. Existence of an intramolecular hydrogen bond (HB structure) for species with the HO-C-C-X moiety is preferred compared with a conformation when the hydrogen bond is disrupted (NoHB) for the considered F- and NO(2)-substituted molecules both in chloroform and aqueous solution. For 2NH(2)-ethanol, the HB structure is predominant in chloroform, whereas the 93:7 ratio for the OCCN trans/gauche species was obtained in aqueous solution. 2NH(2)-phenol exhibits a subtle equilibrium of the HB and NoHB conformations in both solvents. Potential of mean force calculations predict about a 10% solute association for the trans 2NH(2)-ethanol solute even in the fairly dilute 0.22 molar solution, whereas direct MC simulations do not support the maintenance of a doubly hydrogen-bonded dimer. Aqueous solution characteristics, as coordination numbers and numbers of strongly bound water molecules to the solute at T = 298 K and p = 1 atm, correspond reasonably to the derived molecular structures.

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.