Modeling Enantiomeric Separations as an Interfacial Process Using Amylose Tris(3,5-dimethylphenyl carbamate) (ADMPC) Polymers Coated on Amorphous Silica

Chiral high-performance liquid chromatography (HPLC) is commonly performed to isolate the biologically active enantiomer of a drug from the ineffective or even harmful ones. Understanding the molecular-level recognition that underlies this process is necessary for trimming down the very large number of possible combinations of chiral stationary phases, solvent systems, and other experimental HPLC conditions, a particularly important consideration when only small quantities of the racemate are available. Fully atomistic molecular dynamics (MD) simulation is a useful tool to provide this molecular-level understanding and predict experimental separation factors under a given set of conditions. To predict the chiral separation results for drug enantiomers by amylose tris(3,5-dimethylphenyl carbamate) (ADMPC) chiral stationary phase, we design a model of multiple ADMPC polymer strands coated on an amorphous silica slab. Using various MD techniques, we successfully coat ADMPCs onto the surface without losing the structural character of the backbone in the presence of the solvent system. Not only is this model more representative of the polymer surface on a solid support that is encountered by the enantiomers, but it also provides more opportunities for chiral molecules interacting with ADMPC, provides the possibility for large drug molecules to interact with two polymer strands at the same instant, and provides better agreement with experiment when we use the overall average quantities as the predictive metric. For a better understanding of why some metrics are better predictors than others, we use charts of the distribution of hydrogen-bonding lifetimes for various donor–acceptor pairs that contribute to the interaction events determining the relative retention times for the enantiomers. We also examine the contribution of ring–ring interactions to the molecular recognition process and ultimately to the differential retention of enantiomers. The results are more consistent than previous models and resolve the problematic case of two drugs, thalidomide and valsartan.