Physicochemical and graph theoretical descriptors in developmental toxicity SAR: a comparative study

Chemical insults to the developing fetus can lead to growth retardation, malformation, death, and functional deficits. The present study seeks to determine if physicochemical and/or graph theoretical parameters can be used to determine a structure-activity relationship (SAR) for developmental toxicity, and if consistency is observed among the selected features. The biological data utilized consists of a diverse series of compounds evaluated within the Chernoff-Kavlock in vivo mouse assay. Physicochemical parameters calculated correspond to electronic, steric, and transport properties. Graph theoretical parameters calculated include the simple, valence, and kappa indices. Both sets of parameters were independently applied to derive SARs in order to compare the quality of the respective models. Multiple random sampling, without replacement, was utilized to obtain ten training/test partitions. Models were built by linear discriminant analysis, decision trees, and neural networks respectively. Comparisons on identical sets of data were carried out to determine if any of the model building procedures had a significant advantage in terms of predictive performance. Furthermore, comparison of the features selected within and across the model building processes led to the determination of model consistency. Our results indicate that consistent features related to developmental toxicity are observed and that both physicochemical and graph theoretical parameters have equal utility.     

Multivariate Association of Graph-Theoretic Variables and Physicochemical Properties

Abstract

We used canonical correlation analysis to examine the multivariate association between two distinct data sets commonly measured or calculated for approximately 600 chemicals: (1) measured or calculated values of select physieochemical properties (i.e., K _ow, boiling point, heat of vaporization, molecular weight, water solubility, molecular volume, hydrogen bonding potential, and vapor pressure) and (2) calculated algorithmically-derived variables (i.e., topological and neighborhood indices derived from graph theory). Canonical correlation analysis identified eight highly significant associations between linear combinations of graph-theoretic variables and linear combinations of physicochemical properties. The set of graph theoretic variables was significantly related to all physieochemical properties, explaining 55% to 99% of the variation in these properties.     

Molecular Similarity and Risk Assessment: Analog Selection and Property Estimation Using Graph Invariants

Abstract

Four molecular similarity measures have been used to select the nearest neighbor of chemicals in two data sets of 139 hydrocarbons and 15 nitrosamines, respectively. The similarity methods are based on calculated graph invariants which include atom pairs, connectivity indices and information theoretic topological indices. The property of the selected nearest neighbor by each method was taken as the estimate of the property under investigation. The results show that for these data sets, all four methods give reasonable estimates of the properties studied.     

Quantitative molecular similarity analysis (QMSA) methods for property estimation: A comparison of property-based,arbitrary, and tailored similarity spaces

Three classes of arbitrary quantitative molecular similarity analysis (QMSA) methods have been computed using atom pairs, topological indices, and physicochemical properties. Tailored QMSA models have been developed using a selected number of TIs chosen by ridge regression. The methods have been applied to the K -nearest neighbor based estimation of log P of two sets of chemicals. Results show that the property-based and tailored QMSA methods are superior to the arbitrary similarity methods in estimating log P of both sets of chemicals     

FT-Raman and infrared spectra,r0 structural parameters,ab initio calculations and vibrational assignment for nitryl chloride

The FT-Raman spectra (2000-30 cm⁻¹) of liquid and solid nitryl chloride, ClNO₂, along with the infrared spectra (2000-80 cm⁻¹) of the gas and solid have been recorded. All six fundamentals are confidently identified and the potential energy distributions determined from the force fields obtained from ab initio calculations. Several different basis sets have been utilized to determine the harmonic frequencies and force constants which are compared to the previously reported valence force constants. Structural parameters have been calculated with these basis sets including electron correlation with MP2, MP3 and MP4 perturbation. The calculated equilibrium structural parameters are compared to the experimental r₀ structural parameters. The spectra of the solid indicate that there are at least two molecules per primitive cell. All of these results are compared to the corresponding quantities for some similar molecules.     

Spectroscopic constants and molecular properties of CN−, SiH−, PO−, SO−, SF−, and SiS−: Density functional study

Spectroscopic constants and molecular properties of selected diatomic anions namely CN^?, SiH^?, PO^?, SO^?, SF^?, and SiS^? in their ground states have been studied in detail using the hybrid HF/DF B3LYP method. The consistency of the calculated values has been verified with four different basis sets, with improved quality. The spectroscopic constants and molecular properties calculated with the aug‐cc‐pVTZ basis set agree very well with the experimental and theoretical values wherever available. Most of the spectroscopic constants and molecular properties of the selected diatomic anions, particularly the spectroscopic constants and molecular properties of SO^? and SiS^? are reported for the first time. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

Performance of revised STO(1M)-3G basis set for prediction of 5-fluorocytosine chemical shifts

Nuclear shieldings and chemical shifts of 5-fluorocytosine (5FC) were predicted in the gas phase and DMSO solution modeled by polarizable continuum model using B3LYP density functional and revised STO(1M)-3G basis set. For comparison, eight arbitrary selected basis sets including STO-3G and medium-size Pople-type and larger dedicated Jensen-type ones were applied. The former basis sets were significantly smaller, but the calculated structural parameters, harmonic vibrational frequencies, were very accurate and close to those obtained with larger, polarization-consistent ones. The predicted ¹³C and ¹H chemical shieldings of 5FC and cytosine, selected as parent molecule, were acceptable (root mean square for ¹³C chemical shifts in DMSO of about 5 ppm and less) though less accurate than those calculated with large basis sets, dedicated for prediction of nuclear magnetic resonance parameters.     

Molecular dynamics simulation of the hydrocarbon region of a biomembrane using a reduced representation model

The development of a coarse‐grained reduced‐representation model of the hydrocarbon region of a biological membrane is reported. The potential is based on the popular Gay–Berne model of liquid crystals, and involves the linking of individual Gay–Berne ellipsoids by harmonic springs to form each hydrocarbon chain. Diffusion coefficients and order parameters have been calculated by molecular dynamics computer simulations for a range of parameter sets. The results clearly demonstrate the presence of a phase transition from an ordered low‐temperature solid phase reminiscent of the L_β′ phase of phospholipids, to a high‐temperature disordered phase reminiscent of the L_α phase. Order parameters calculated for each layer of the model are consistent with the experimental segmental order parameters reported for dipalmitoyl phosphatidylcholine. The application of this model to the study of small molecule diffusion within the membrane core is proposed. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1622–1633, 2001     

Effects of Chemical Reactivity on the Toxicity of Phosphorus Fluoridates

Abstract

Semiempirical quantum calculations were performed on a series of organophosphorus fluoridates to determine the relative reactivity for hydrolysis. This value was determined by subtracting the energy of the metastable intermediate from the energy of the stable molecule. Plotting this relative reactivity for each compound vs. its toxicity resulted in a parabolic curve with nerve agents and other similarly toxic compounds in the center. The more reactive phosphinates and less reactive phosphates were at the edges of the graph in the region of lower toxicity. The results indicate that for compounds meeting minimal structural requirements, chemical reactivity is the principal determinant of cholinesterase inhibition.     

Decision tree SAR models for developmental toxicity based on an FDA/TERIS database

Humans are exposed to thousands of environmental chemicals for which no developmental toxicity information is available. Structure-activity relationships (SARs) are models that could be used to efficiently predict the biological activity of potential developmental toxicants. However, at this time, no adequate SAR models of developmental toxicity are available for risk assessment. In the present study, a new developmental database was compiled by combining toxicity information from the Teratogen Information System (TERIS) and the Food and Drug Administration (FDA) guidelines. We implemented a decision tree modeling procedure, using Classification and Regression Tree software and a model ensemble approach termed bagging. We then assessed the empirical distributions of the prediction accuracy measures of the single and ensemble-based models, achieved by repeating our modeling experiment many times by repeated random partitioning of the working database. The decision tree developmental SAR models exhibited modest prediction accuracy. Bagging tended to enhance the accuracy of prediction. Also, the model ensemble approach reduced the variability of prediction measures compared to the single model approach. Further research with data derived from animal species- and endpoint-specific components of an extended and refined FDA/TERIS database has the potential to derive SAR models that would be useful in the developmental risk assessment of the thousands of untested chemicals.