QSAR-based toxicity classification and prediction for single and mixed aromatic compounds

Quantitative structure-activity relationships (QSARs) based on the octanol/water partition coefficient were employed to predict acute toxicities of 36 substituted aromatic compounds and their mixtures. In this study, the model developed by Verhaar et al. was modified and used to calculate octanol/water partition coefficients of chemical mixtures. To validate the model, acute toxicities of these chemicals were measured to Vibrio fischeri in terms of EC⁵⁰. The results indicated that the obtained QSAR models could be used to predict toxicities of samples consisting of these substituted aromatic compounds, individually or in combinations. The obtained equations were proved to be robust enough by using the leave-one-out test method. By classifying these chemicals into two groups, polar and non-polar, the toxicities of chemical mixtures within each group can be predicted accurately from their calculated partition coefficients.     

Profiling primaquine metabolites in primary human hepatocytes using UHPLC‐QTOF‐MS with 13C stable isotope labeling

Therapeutic efficiency and hemolytic toxicity of primaquine (PQ), the only drug available for radical cure of relapsing vivax malaria are believed to be mediated by its metabolites. However, identification of these metabolites has remained a major challenge apparently due to low quantities and their reactive nature. Drug candidates labeled with stable isotopes afford convenient tools for tracking drug‐derived metabolites in complex matrices by liquid chromatography‐tandem mass spectrometry (LC‐MS‐MS) and filtering for masses with twin peaks attributable to the label. This study was undertaken to identify metabolites of PQ from an in vitro incubation of a 1:1 w/w mixture of ¹³C₆‐PQ/PQ with primary human hepatocytes. Acquity ultra‐performance LC (UHPLC) was integrated with QTOF‐MS to combine the efficiency of separation with high sensitivity, selectivity of detection and accurate mass determination. UHPLC retention time, twin mass peaks with difference of 6 (originating from ¹³C₆‐PQ/PQ), and MS‐MS fragmentation pattern were used for phenotyping. Besides carboxy‐PQ (cPQ), formed by oxidative deamination of PQ to an aldehyde and subsequent oxidation, several other metabolites were identified: including PQ alcohol, predictably generated by oxidative deamination of PQ to an aldehyde and subsequent reduction, its acetate and the alcohol's glucuronide conjugate. Trace amounts of quinone‐imine metabolites of PQ and cPQ were also detected which may be generated by hydroxylation of the PQ/cPQ quinoline ring at the 5‐position and subsequent oxidation. These findings shed additional light on the human hepatic metabolism of PQ, and the method can be applied for identification of reactive PQ metabolites generated in vivo in preclinical and clinical studies. Copyright © 2013 John Wiley & Sons, Ltd.     

Computational Study of the One‐ and Two‐Dimensional Infrared Spectra of a Proton‐Transfer Mode in a Hydrogen‐Bonded Complex Dissolved in a Polar Nanocluster

The signatures of nanosolvation on the one‐ and two‐dimensional (1D and 2D) IR spectra of a proton‐transfer mode in a hydrogen‐bonded complex dissolved in polar solvent molecule nanoclusters of varying size are elucidated by using mixed quantum–classical molecular dynamics simulations. For this particular system, increasing the number of solvent molecules successively from N=7 to N=9 initiates the transition of the system from a cluster state to a bulk‐like state. Both the 1D and 2D IR spectra reflect this transition through pronounced changes in their peak intensities and numbers, but the time‐resolved 2D IR spectra also manifest spectral features that uniquely identify the onset of the cluster‐to‐bulk transition. In particular, it is observed that in the 1D IR spectra, the relative intensities of the peaks change such that the number of peaks decreases from three to two as the size of the cluster increases from N=7 to N=9. In the 2D IR spectra, off‐diagonal peaks are observed in the N=7 and N=8 cases at zero waiting time, but not in the N=9 case. It is known that there are no off‐diagonal peaks in the 2D IR spectrum of the bulk version of this system at zero waiting time, so the disappearance of these peaks is a unique signature of the onset of bulk‐like behavior. Through an examination of the trajectories of various properties of the complex and solvent, it is possible to relate the emergence of these off‐diagonal peaks to an interplay between the vibrations of the complex and the solvent polarization dynamics.     

Assessing Solvation Effects on Chemical Reactions with Globally Optimized Solvent Clusters

We have examined the influence of water solvent on the Menshutkin reaction of methyl chloride with ammonia by performing static, quantum chemical calculations. We have employed large, explicit, and globally structure‐optimized water clusters around the reaction center, in a mixed explicit/implicit solvent model. This approach deliberately deviates from attempts to capture the most likely solvent‐molecule distribution around a reaction center. Instead, it explores extremes on the scale of rearrangement speed in terms of the surrounding solvent cluster, relative to the reaction progress itself. A comparison to traditional theoretical and experimental results enables us to quantify the energy penalty that is induced by the inability of the water cluster to instantaneously and completely follow the reaction progress. In addition, the influence of water clusters on the reaction energy profile can be much larger than merely changing it somewhat. Certain clusters can completely annihilate the sizeable activation barrier of 23.5 kcal mol^?1.