A tentative quantitative structure-toxicity relationship study of benzodiazepine drugs.

Benzodiazepines belong to a large family of drugs, being used as hypnotics, anxiolytics, tranquillizers, anticonvulsants, in pre-medication and intravenous sedation. Several quantitative structure-toxicity (lethal oral dose for mouse) relationship (QSTR) models for 54 benzodiazepine derivatives have been developed. The molecular structure of these compounds was energetically optimized by molecular mechanics calculations. To the lowest energy conformations thus obtained, quantum chemical calculations (RM1 approach) were applied to finally optimize the structures. Several structural descriptors, volumes, molecular surface area, hydrophobicities and quantum chemical descriptors were calculated from the minimized structures. Multiple linear regression (MLR) combined with genetic algorithm for variable selection, artificial neural networks (ANNs), support vector machines (SVMs) and partial least squares (PLS) have been employed. Few satisfactory MLR models with predictive power were obtained. Nonlinear modelling methods of ANNs and SVMs gave somewhat better models than those obtained by MLR using same set of descriptors. Additional information on the factors which influence the benzodiazepine toxicity was given by PLS. The obtained models can be used for a rough evaluation of benzodiazepine toxicity.

[1]  J. F. Gentleman,et al.  Detecting Outliers. II. Supplementing the Direct Analysis of Residuals , 1975 .

[2]  D. W. Osten,et al.  Selection of optimal regression models via cross‐validation , 1988 .

[3]  B. Efron.  Better Bootstrap Confidence Intervals , 1987 .

[4]  Hxugo Kubiny.  Variable Selection in QSAR Studies. I. An Evolutionary Algorithm , 1994 .

[5]  Roberto Todeschini,et al.  Handbook of Molecular Descriptors , 2002 .

[6]  S. Wold.  Cross-Validatory Estimation of the Number of Components in Factor and Principal Components Models , 1978 .

[7]  Chih-Jen Lin,et al.  LIBSVM: A library for support vector machines , 2011, TIST.

[8]  S. Gayen,et al.  Quantitative structure–activity relationship study on some benzodiazepine derivatives as anti-Alzheimer agents , 2004, Journal of molecular modeling.

[9]  C. Hansch,et al.  QUANTITATIVE STRUCTURE-ACTIVITY RELATIONSHIPS OF THE BENZODIAZEPINES. A REVIEW AND REEVALUATION , 1994 .

[10]  Vladimir N. Vapnik,et al.  The Nature of Statistical Learning Theory , 2000, Statistics for Engineering and Information Science.

[11]  P. Breggin.  ANALYSIS OF ADVERSE BEHAVIORAL EFFECTS OF BENZODIAZEPINES WITH A DISCUSSION ON DRAWING SCIENTIFIC CONCLUSIONS FROM THE FDA'S SPONTANEOUS REPORTING SYS TEM , 1998 .

[12]  W. Meylan,et al.  Atom/fragment contribution method for estimating octanol-water partition coefficients. , 1995, Journal of pharmaceutical sciences.

[13]  A. Höskuldsson.  PLS regression methods , 1988 .

[14]  Wenjian Wang,et al.  Determination of the spread parameter in the Gaussian kernel for classification and regression , 2003, Neurocomputing.

[15]  Igor V. Tetko,et al.  Application of Associative Neural Networks for Prediction of Lipophilicity in ALOGPS 2.1 Program , 2002, J. Chem. Inf. Comput. Sci..

[16]  W. Krzanowski,et al.  Cross-Validatory Choice of the Number of Components From a Principal Component Analysis , 1982 .

[17]  Lennart Eriksson,et al.  Model validation by permutation tests: Applications to variable selection , 1996 .

[18]  Alexander Golbraikh,et al.  Rational selection of training and test sets for the development of validated QSAR models , 2003, J. Comput. Aided Mol. Des..

[19]  A. Klamt,et al.  COSMO : a new approach to dielectric screening in solvents with explicit expressions for the screening energy and its gradient , 1993 .

[20]  Roberto Todeschini,et al.  MobyDigs: software for regression and classification models by genetic algorithms , 2003 .

[21]  Riccardo Leardi,et al.  Nature-Inspired Methods in Chemometrics: Genetic Algorithms and Artificial Neural Networks , 2005 .

[22]  H. Frank,et al.  Statistics: concepts and applications , 1996 .

[23]  D. Nutt,et al.  Radioligands for PET studies of central benzodiazepine receptors and PK (peripheral benzodiazepine) binding sites--current status. , 1993, Nuclear medicine and biology.

[24]  R. Todeschini,et al.  Detecting bad regression models: multicriteria fitness functions in regression analysis , 2004 .

[25]  Zbigniew Michalewicz,et al.  Evolutionary algorithms , 1997, Emerging Evolutionary Algorithms for Antennas and Wireless Communications.

[26]  O. Kempthorne,et al.  STATISTICS AND MATHEMATICS IN BIOLOGY , 1955 .

[27]  M. Cronin,et al.  Pitfalls in QSAR , 2003 .

[28]  Andreas Klamt,et al.  COSMOfrag: A Novel Tool for High-Throughput ADME Property Prediction and Similarity Screening Based on Quantum Chemistry , 2005, J. Chem. Inf. Model..

[29]  S. Geisser.  A predictive approach to the random effect model , 1974 .

[30]  Svante Wold,et al.  Multivariate quantitative structure-activity relationships (QSAR): conditions for their applicability , 1983, J. Chem. Inf. Comput. Sci..

[31]  J. Stewart.  Optimization of parameters for semiempirical methods V: Modification of NDDO approximations and application to 70 elements , 2007, Journal of molecular modeling.

[32]  L. Eriksson.  Multi- and megavariate data analysis , 2006 .

[33]  G. Webb,et al.  Electronic factors in the structure-activity relationship of some 1,4-benzodiazepin-2-ones. , 1977, Journal of medicinal chemistry.

[34]  Brain receptor binding and lipophilic character of benzodiazepines. , 1983, Biochemical pharmacology.

[35]  Johann Gasteiger,et al.  The Anesthetic Activity and Toxicity of Halogenated Ethyl Methyl Ethers, a Multivariate QSAR Modelled by PLS , 1985 .

[36]  Igor V. Tetko,et al.  Prediction of n-Octanol/Water Partition Coefficients from PHYSPROP Database Using Artificial Neural Networks and E-State Indices , 2001, J. Chem. Inf. Comput. Sci..

[37]  S. Hirono,et al.  Simple Method of Calculating Octanol/Water Partition Coefficient. , 1992 .

[38]  M L Connolly,et al.  The molecular surface package. , 1993, Journal of molecular graphics.

[39]  Alfredo Mayall Simas,et al.  RM1: A reparameterization of AM1 for H, C, N, O, P, S, F, Cl, Br, and I , 2006, J. Comput. Chem..

[40]  Svante Wold,et al.  Pattern recognition by means of disjoint principal components models , 1976, Pattern Recognit..

[41]  D. Maddalena,et al.  Prediction of receptor properties and binding affinity of ligands to benzodiazepine/GABAA receptors using artificial neural networks. , 1995, Journal of medicinal chemistry.

[42]  M. Stone,et al.  Cross‐Validatory Choice and Assessment of Statistical Predictions , 1976 .

[43]  G H Loew,et al.  Theoretical structure-activity studies of benzodiazepine analogues. Requirements for receptor affinity and activity. , 1984, Molecular pharmacology.

[44]  G M Crippen,et al.  Modeling the benzodiazepine receptor binding site by the general three-dimensional structure-directed quantitative structure-activity relationship method REMOTEDISC. , 1990, Molecular pharmacology.

[45]  Roberto Todeschini,et al.  The K correlation index: theory development and its application in chemometrics , 1999 .

[46]  B. Efron.  Estimating the Error Rate of a Prediction Rule: Improvement on Cross-Validation , 1983 .

[47]  Ettore Novellino,et al.  Study of Benzodiazepines Receptor Sites Using a Combined QSAR‐CoMFA Approach , 1992 .