Mobile phase effects in reversed-phase liquid chromatography: a comparison of acetonitrile/water and methanol/water solvents as studied by molecular simulation

Molecular simulations of water/acetonitrile and water/methanol mobile phases in contact with a C(18) stationary phase were carried out to examine the molecular-level effects of mobile phase composition on structure and retention in reversed-phase liquid chromatography. The simulations indicate that increases in the fraction of organic modifier increase the amount of solvent penetration into the stationary phase and that this intercalated solvent increases chain alignment. This effect is slightly more apparent for acetonitrile containing solvents. The retention mechanism of alkane solutes showed contributions from both partitioning and adsorption. Despite changes in chain structure and solvation, the molecular mechanism of retention for alkane solutes was not affected by solvent composition. The mechanism of retention for alcohol solutes was primarily adsorption at the interface between the mobile and stationary phase, but there were also contributions from interactions with surface silanols. The interaction between the solute and surface silanols become very important at high concentrations of acetonitrile.     

Comparison of retention of aromatic hydrocarbons with polar groups in binary reversed-phase high-performance liquid chromatography systems

The retention of aromatic hydrocarbons with polar groups has been correlated as log k1 versus log k2 for reversed-phase high-performance liquid chromatography systems with different binary aqueous mobile phases containing methanol, acetonitrile or tetrahydrofuran as modifiers. Distinct changes in separation selectivity have been observed between tetrahydrofuran and acetonitrile or methanol systems. Methanol and acetonitrile systems show lower diversity of separation selectivity. The changes in retention and selectivity of aromatic hydrocarbons with various polar groups between any two chromatographic systems with binary aqueous eluents (tetrahydrofuran vs. acetonitrile, tetrahydrofuran vs. methanol and methanol vs. acetonitrile) have been interpreted in terms of molecular interactions of the solute with especially one component of the stationary phase region, i.e. extracted modifier, and stationary phase ordering. The ordering of the stationary phase region caused by modifier type influences the chromatographic selectivity of solutes with different molecular shape.     

Absolute free energy calculations by thermodynamic integration in four spatial dimensions

An optimized technique for calculating the excess chemical potential of small molecules in dense liquids and the binding affinity of molecular ligands to biomolecules is reported. In this method, a molecular species is coupled to the system of interest via a nonphysical fourth spatial dimension w through which insertion or extraction can be carried out [R. Pomes, E. Eisenmesser, C. B. Post et al., J. Chem. Phys. 111, 3387 (1999)]. Molecular simulations are used to compute the potential of mean force (PMF) acting on the solute molecule in the fourth dimension. The excess chemical potential of that molecule is obtained as the difference in the PMF between fully coupled and fully decoupled systems. The simplicity, efficiency, and generality of the method are demonstrated for the calculation of the hydration free energies of water and methanol as well as sodium, cesium, and chloride ions. A significant advantage over other methods is that the 4D-PMF approach provides a single effective and general route for decoupling all nonbonded interactions (i.e., both Lennard-Jones and Coulombic) at once for both neutral and charged solutes. Direct calculation of the mean force from thermodynamic integration is shown to be more computationally efficient than calculating the PMF from umbrella sampling. Statistical error analysis suggests a simple strategy for optimizing sampling. The detailed analysis of systematic errors arising from the truncation of Coulombic interactions in a solvent droplet of finite size leads to straightforward corrections to ionic hydration free energies.     

Characterizing hydrophobicity at the nanoscale: a molecular dynamics simulation study

We use molecular dynamics (MD) simulations of water near nanoscopic surfaces to characterize hydrophobic solute-water interfaces. By using nanoscopic paraffin like plates as model solutes, MD simulations in isothermal-isobaric ensemble have been employed to identify characteristic features of such an interface. Enhanced water correlation, density fluctuations, and position dependent compressibility apart from surface specific hydrogen bond distribution and molecular orientations have been identified as characteristic features of such interfaces. Tetrahedral order parameter that quantifies the degree of tetrahedrality in the water structure and an orientational order parameter, which quantifies the orientational preferences of the second solvation shell water around a central water molecule, have also been calculated as a function of distance from the plate surface. In the vicinity of the surface these two order parameters too show considerable sensitivity to the surface hydrophobicity. The potential of mean force (PMF) between water and the surface as a function of the distance from the surface has also been analyzed in terms of direct interaction and induced contribution, which shows unusual effect of plate hydrophobicity on the solvent induced PMF. In order to investigate hydrophobic nature of these plates, we have also investigated interplate dewetting when two such plates are immersed in water.     

Effects of the molecular properties of mixed solvents on the elution of alkyl benzoates in RPLC.

The retention behavior and mechanism of methyl, ethyl, propyl, isopropyl, buthyl and isobuthyl benzoates have been studied at different eluent compositions of aqueous mixtures with water-soluble organic solvents (methanol, ethanol, 1-propanol, 2-propanol, acetonitrile (AN), 1,4-dioxane and tetrahydrofuran (THF)) in RPLC. The retention of the solutes is discussed based on the solvent composition, solvent polarity (ETN value), preferential solvation, hydrogen bonding and solvent clusters of the eluents. The smaller ETN values and the larger preferential solvation of the mixed solvent eluted the solutes faster. The IR spectra of HDO suggested that the solvents, except for methanol and ethanol, break the hydrogen bonding between water molecules, resulting in fast elution of the solutes. Based upon the results, we chose an optimum solvent composition for the separation of benzoates and applied it to the determination of the benzoates in clove.     

Molecular dynamics study of free energy of transfer of alcohol and amine from water phase to the micelle by thermodynamic integration method

Free energy of transfer of methylamine, octylamine, methanol, and octanol from water phase to sodium dodecyl sulfate (SDS) micelle has been calculated using thermodynamic integration method combined with molecular dynamics calculations. Together with the results for alkanes obtained in our previous study [K. Fujimoto, N. Yoshii, and S. Okazaki, J. Chem. Phys. 133, 074511 (2010)], the effect of polar group on the partition of hydrophilic solutes between water phase and the micelle has been investigated in detail at a molecular level. The calculations showed that the molecules with octyl group are more stable in the SDS micelle than in the water phase due to their hydrophobicity of long alkyl chain. In contrast, methanol and methylamine are stable in the water phase as well as in the micelle because of their high hydrophilicity. The spatial distribution of methylamine, octylamine, methanol, and octanol has also been evaluated as a function of the distance, R, from the center of mass of SDS micelle to the solutes. The distribution shows that the methylamine molecule is adsorbed on the SDS micelle surface, while the methanol molecule is delocalized among the whole system, i.e., in the water phase, on the surface of the micelle, and in the hydrophobic core of the micelle. The octylamine and octanol molecules are solubilized in the SDS micelle with palisade layer structure and are not found in the water phase.     

Effect of ion-exchange site and eluent modifiers on the anion-exchange of carboxylic acids

The chromatographic behavior of carboxylic acids has been investigated, on three different latex-based anion-exchange columns, in order to define the effect of the ion-exchange site structure on selectivity. The analytical columns produced are characterized by alkyl amines containing zero, one or two hydroxyl groups on the anion-exchange functional site. Divalent carboxylic acids, namely fumaric, maleic, trans-beta-hydromuconic, trans,trans-muconic, oxalic, malonic, succinic, glutaric, adipic, malic, tartaric and mucic acids, have been chosen as test solutes. The performance of the three stationary phases has been studied employing NaOH eluents and has been discussed with respect to the different hydrophilicity of the ion-exchange sites and analytes. Considering on previous results obtained using organic solvents (methanol and acetonitrile) with carbonate eluents on a highly hydrophilic column, the performance of the three exchangers has also been studied using acetonitrile, methanol and n-propanol. The chromatographic behavior was similar for the three columns studied, but the different organic solvents gave variations in selectivity. In order to characterize these differences, particle size measurements of the latices were performed both in pure water and in the presence of each organic solvent studied.     

Molecular origin of anticooperativity in hydrophobic association

The structuring of water molecules in the vicinity of nonpolar solutes is responsible for hydrophobic hydration and association thermodynamics in aqueous solutions. Here, we studied the potential of mean force (PMF) for the formation of a dimer and trimers of methane molecules in three specific configurations in explicit water to explain multibody effects in hydrophobic association on a molecular level. We analyzed the packing and orientation of water molecules in the vicinity of the solute to explain the effect of ordering of the water around nonpolar solutes on many-body interactions. Consistent with previous theoretical studies, we observed cooperativity, manifested as a reduction of the height of the desolvation barrier for the trimer in an isosceles triangle geometry, but for linear trimers, we observed only anticooperativity. A simple mechanistic picture of hydrophobic association is drawn. The free energy of hydrophobic association depends primarily on the difference in the number of water molecules in the first solvation shell of a cluster and that in the monomers of a cluster; this can be approximated by the molecular surface area. However, there are unfavorable electrostatic interactions between the water molecules from different parts of the solvation shell of a trimer because of their increased orientation induced by the nonpolar solute. These electrostatic interactions make an anticooperative contribution to the PMF, which is clearly manifested for the linear trimer where the multibody contribution due to changes in the molecular surface area is equal to zero. The information theory model of hydrophobic interactions of Hummer et al. also explains the anticooperativity of hydrophobic association of the linear trimers; however, it predicts anticooperativity with a qualitatively identical distance dependence for nonlinear trimers, which disagrees with the results of simulations.     

Molecular dynamics simulations of pressure effects on hydrophobic interactions.

We report results on the pressure effects on hydrophobic interactions obtained from molecular dynamics simulations of aqueous solutions of methanes in water. A wide range of pressures that is relevant to pressure denaturation of proteins is investigated. The characteristic features of water-mediated interactions between hydrophobic solutes are found to be pressure-dependent. In particular, with increasing pressure we find that (1) the solvent-separated configurations in the solute-solute potential of mean force (PMF) are stabilized with respect to the contact configurations; (2) the desolvation barrier increases monotonically with respect to both contact and solvent-separated configurations; (3) the locations of the minima and the barrier move toward shorter separations; and (4) pressure effects are considerably amplified for larger hydrophobic solutes. Together, these observations lend strong support to the picture of the pressure denaturation process proposed previously by Hummer et al. (Proc. Natl. Acad. Sci. U.S.A. 1998, 95, 1552): with increasing pressure, the transfer of water into protein interior becomes key to the pressure denaturation process, leading to the dissociation of close hydrophobic contacts and subsequent swelling of the hydrophobic protein interior through insertions of water molecules. The pressure dependence of the PMF between larger hydrophobic solutes shows that pressure effects on the interaction between hydrophobic amino acids may be considerably amplified compared to those on the methane-methane PMF.     

Relevance of pi-pi and dipole-dipole interactions for retention on cyano and phenyl columns in reversed-phase liquid chromatography

Previous work suggests that pi-pi interactions between certain solutes and both phenyl and cyano columns can contribute to sample retention and the selectivity of these two column types versus alkylsilica columns. Recent studies also suggest that dipole-dipole interactions are generally unimportant for retention on cyano columns. The present study presents data for 44 solutes, three columns and two different mobile phases that were selected to further test these conclusions. We find that pi-pi interactions can contribute to retention on both cyano and phenyl columns, while dipole-dipole interactions are likely to be significant for the retention of polar aliphatic solutes on cyano columns. When acetonitrile/water mobile phases are used, both pi-pi and dipole-dipole interactions are suppressed, compared to the use of methanol/water.     

Mobile phase effects on retention on a new butylimidazolium-based high-performance liquid chromatographic stationary phase

A new HPLC stationary phase based on n-butylimidazolium bromide has been characterized by a linear solvation energy relationship (LSER) approach in the binary acetonitrile/water mobile phases. The retention properties of the stationary phase were systematically evaluated in terms of intermolecular interactions between 28 test solutes and the stationary phase. The results and further comparisons with conventional reversed phase system confirm that retention properties are similar to phenyl phases in acetonitrile/water mixtures. The results obtained with acetonitrile/water mixtures are also compared with results obtained using methanol/water mixtures.     

Influence of solvent effects on retention in reversed-phase liquid chromatography and solid-phase extraction using a cyanopropylsiloxane-bonded,silica-based sorbent

Summary The solvation parameter model is used to characterize the retention properties of a cyanopropylsiloxanebonded, silica-based sorbent with methanol, acetonitrile, tetrahydrofuran, and isopropanol in water as mobile phases. The system constants over the composition range 1 to 50 % (v/v) organic solvent indicate that retention occurs because of the relative ease of cavity formation in the solvated stationary phase compared to the same process in the predominantly aqueous mobile phase as well as from more favorable stationary phase interactions with solutes containing π- and n-electrons. The capacity of the solute for dipole-type interactions is not important whereas all hydrogen-bond-type interactions result in reduced retention. Graphing the system constants as a function of mobile phase composition provides a simple mechanism for interpreting the change in capacity of the chromatographic system for retention in terms of changes in the relative weighting of fundamental intermolecular interactions. A comparison is also made with the retention properties of an octadecylsiloxane-bonded, silica-based sorbent with 30 % (v/v) methanol in water as the mobile phase and the extraction characteristics of a porous polymer sorbent with 1 % (v/v) methanol, acetonitrile, tetrahydrofuran, and isopropanol in water as the sample processing solvent. Changes in sorbent selectivity due to selective uptake of the processing solvent are much smaller for the cyanopropylsiloxane-bonded sorbent than the results found for a porous polymer sorbent.     

Enthalpies of solvation of alkylated uracils in water,nonaqueous and mixed solvents

The enthalpies of solution of uracil and its alkylated derivatives in water, methanol, N,N-dimethylformamide (DMF) and water+DMF mixtures were measured at 25°C. The enthalpies of solvation were determined. The enthalpies of cavity formation, corresponding to the enthalpies of solvent-solvent interactions were calculated and the enthalpies of solute-solvent interactions were obtained. The presence of the alkyl groups was found to have different effects on the enthalpy of interaction depending on the position and size of the substitution. The effect of alkylation at the nonpolar side of the uracil ring was found to arise mostly from the enhancement of the van der Waals interactions. The alkyl substitutions at the polar side resulted also in the removal of the solvent molecules interacting specifically with the polar groups of uracil. The enthalpy of those specific interactions was determined and found to be stronger in methanol and DMF than in water. Enthalpies of solvation in the binary water+DMF solvent were found to depend in a nonlinear way on the solvent composition. The nonlinearities in the water-rich region were found to arise from the decay of the hydrophobic hydration of the solutes with the increasing DMF content. The substitution of two methyl groups caused the uracil molecule to bahave as a predominantly hydrophobic solute. The nonlinearities in the DMF-rich region were found only for those solutes which can form hydrogen bonds with DMF.     

Separation selectivity of some phenolic acids in RP HPLC systems with binary mobile phase comprised various modifiers

Retention of phenolic acid has been correlated for reversed-phase high-performance liquid chromatography systems with different binary mobile phases containing methanol, acetonitrile or tetrahydrofuran as modifiers and buffer at pH 3.0, 4.6 and 5.0. The changes of separation selectivity of solutes, when one modifier is replaced by another in the eluent, has been explained taking into consideration molecular interactions of the solutes with components of the stationary phase region, i.e. extracted modifier, and ordering of the stationary phase by the modifier.     

Potential-dependent vibrational spectroscopy of solvent molecules at the Pt(111) electrode in a water/acetonitrile mixture studied by sum frequency generation

Sum frequency generation (SFG) vibrational spectra of D(2)O and/or acetonitrile (CH(3)CN) on a Pt(111) single-crystal electrode were obtained as a function of applied potential in a 5 mol % water/acetonitrile mixed solvent with different 0.1 molar MSO(3)CF(3) salts (M = H(+), Li(+), Na(+), K(+), and Cs(+)). The results provide a very specific model for the composition of the inner Helmholtz layer as a function of potential and surface charge. Acetonitrile dominates the inner layer with the CN group directed toward the metal at potentials where the metal has a positive charge. As the surface becomes negatively charged, the acetonitrile orientation flips 180 degrees, with the CH(3) group pointing toward the surface. At even more negative surface charge, D(2)O displaces acetonitrile from the inner layer and is the predominant molecule on the surface. Here water is present as an oriented molecule with the oxygen end pointing toward the metal. The potential (and surface charge) where water is the dominant molecule in the inner Helmholtz layer is determined by the solvation energy of the cation.     

Thermodynamic and kinetic characterization of a bridged-ethylene hybrid C18 stationary phase

In this study, a series of polycyclic aromatic hydrocarbons (PAHs) and nitrogen-containing polycyclic aromatic hydrocarbons (NPAHs) is separated on a hybrid stationary phase using methanol and acetonitrile mobile phases. Temperature is varied from 283 to 313 K in order to determine thermodynamic and kinetic parameters of the separation. Thermodynamic behavior is characterized by the retention factor and associated changes in molar enthalpy, whereas kinetic behavior is characterized by the rate constants and associated activation energies. In this study, the retention factors for the NPAHs are smaller than those for the parent PAHs in methanol, while they are more similar to the parent PAHs in acetonitrile. The changes in molar enthalpy are very similar for all solutes, yet are more negative in acetonitrile than in methanol. The rate constants for the NPAHs are smaller than those for their parent PAHs in both mobile phases. Moreover, the rate constants in acetonitrile are one to four orders of magnitude smaller than those in methanol. Based on these thermodynamic and kinetic results, the hybrid stationary phase is compared to traditional silica stationary phases. In addition, the relative contributions from the partition and adsorption mechanisms are discussed.     

Dynamics of water trapped between hydrophobic solutes

We describe the model dynamical behavior of the solvent between two nanoscopic hydrophobic solutes. The dynamics of the vicinal water in various sized traps is found to be significantly different from bulk behavior. We consider the dynamics at normal temperature and pressure at three intersolute distances corresponding to the three solvent separated minima in the free energy profile between the solutes with attractions. These three states correspond to one, two, and three intervening layers of water molecules. Results are obtained from a molecular dynamics simulation at constant temperature and pressure (NPT) ensemble. Translational diffusion of water molecules trapped between the two solutes has been analyzed from the velocity correlation function as well as from the mean square displacement of the water molecules. The rotational behavior has been analyzed through the reorientational dynamics of the dipole moment vector of the water molecule by calculating both first and second rank dipole-dipole correlation functions. Both the translational and reorientational mobilities of water are found to be much slower at the smaller separation and increases as the separation between solutes becomes larger. The occupation time distribution functions calculated from the trajectories also show that the relaxation is much slower for the smallest intersolute separation as compared to other wider separations. The sublinear trend in mean square displacement and the stretched exponential decay of the relaxation of dipolar correlation and occupation distribution function indicate that the dynamical behavior of water in the confined region between two large hydrophobic solutes departs from usual Brownian behavior. This behavior is reminiscent of the behavior of water in the vicinity of protein surface clefts or trapped between two domains of a protein.