Prediction of Tumoricidal Activity and Accumulation of Photosensitizers in Photodynamic Therapy Using Multiple Linear Regression and Artificial Neural Networks¶

Abstract The biological activities of a congeneric series of pyropheophorbides used as sensitizers in photodynamic therapy have been predicted on the basis of their molecular structures, using multiple linear regression and artificial neural network (ANN) computations. Theoretical descriptors (a total of 81) were calculated by the 3DNET program based on the three-dimensional structure (3D) of the geometry-optimized molecules. These input descriptors were tested as independent variables and used for model building. Systematic descriptor selections yielded models with one, two or three descriptors with good cross-validation results. The predictive abilities of the best fitting models were checked by shuffling and cross-validation procedures. ANN was suitable for building models for both linear and nonlinear relationships. Lipophilicity was sufficient to predict the accumulation of the sensitizers in the target tissue. Weighted holistic invariant molecular descriptors weighted by atomic mass, Van der Waals volume or electronegativity were also needed to predict photodynamic activity properly. Our models were able to predict the biological activities of 13 pyropheophorbide derivatives solely on the basis of their 3D molecular structures. Moreover, linear and nonlinear variable selection methods were compared in models built linearly and nonlinearly. It is expedient to use the same method (linear or nonlinear) for variable selection as for parameter estimation.

[1]  Paola Gramatica,et al.  New 3D Molecular Descriptors: The WHIM theory and QSAR Applications , 2002 .

[2]  N. Oleinick,et al.  The Synthesis, Photophysical and Photobiological Properties and in vitro Structure‐Activity Relationships of a Set of Silicon Phthalocyanine PDT Photosensitizers , 1997, Photochemistry and photobiology.

[3]  M. Randic Characterization of molecular branching , 1975 .

[4]  A. Y. Meyer The size of molecules , 1987 .

[5]  Johann Gasteiger,et al.  Electronegativity equalization: application and parametrization , 1985 .

[6]  Gregg D. Wilensky,et al.  Neural Network Studies , 1993 .

[7]  T J Dougherty,et al.  STUDIES ON THE STRUCTURE OF PORPHYRINS CONTAINED IN PHOTOFRIN® II , 1987, Photochemistry and photobiology.

[8]  George M. Whitesides,et al.  FEED-FORWARD NEURAL NETWORKS IN CHEMISTRY : MATHEMATICAL SYSTEMS FOR CLASSIFICATION AND PATTERN RECOGNITION , 1993 .

[9]  Alan R. Katritzky,et al.  Normal Boiling Points for Organic Compounds: Correlation and Prediction by a Quantitative Structure-Property Relationship , 1998, J. Chem. Inf. Comput. Sci..

[10]  Bernd Beck,et al.  Prediction of the n-Octanol/Water Partition Coefficient, logP, Using a Combination of Semiempirical MO-Calculations and a Neural Network , 1997 .

[11]  W Freyer,et al.  Photophysical properties and photodynamic activity in vivo of some tetrapyrroles. , 1990, Biophysical chemistry.

[12]  N. Bodor,et al.  Neural network studies: Part 3. Prediction of partition coefficients , 1994 .

[13]  Michael L. Connolly,et al.  Computation of molecular volume , 1985 .

[14]  Shuichi Hirono,et al.  Estimation of hydrophobicity based on the solvent-accessible surface area of molecules , 1985 .

[15]  W R Greco,et al.  An in vivo quantitative structure-activity relationship for a congeneric series of pyropheophorbide derivatives as photosensitizers for photodynamic therapy. , 1997, Cancer research.

[16]  G. R. Famini,et al.  Using theoretical descriptors in quantitative structure–property relationships: some distribution equilibria , 1998 .

[17]  H. Kubinyi QSAR : Hansch analysis and related approaches , 1993 .

[18]  R. Boyle,et al.  Structure and Biodistribution Relationships of Photodynamic Sensitizers * , 1996, Photochemistry and photobiology.

[19]  K. Héberger,et al.  Comparison of chemometric methods for prediction of rate constants and activation energies of radical addition reactions , 1999 .

[20]  Gerald T. Ankley,et al.  A QSAR analysis of substituent effects on the photoinduced acute toxicity of PAHs , 1995 .

[21]  Joop L. M. Hermens,et al.  Relationships Between Octanol/Water Partition Coefficients and Total Molecular Surface Area and Total Molecular Volume of Hydrophobic Organic Chemicals , 1990 .

[22]  H. Wiener Correlation of Heats of Isomerization, and Differences in Heats of Vaporization of Isomers, Among the Paraffin Hydrocarbons , 1947 .

[23]  Patrick Gaillard,et al.  Molecular Lipophilicity Potential, a tool in 3D QSAR: Method and applications , 1994, J. Comput. Aided Mol. Des..

[24]  R. K. Pandey,et al.  Parabolic Quantitative Structure‐Activity Relationships and Photodynamic Therapy: Application of a Three‐Compartment Model with Clearance to the In Vivo Quantitative Structure‐Activity Relationships of a Congeneric Series of Pyropheophorbide Derivatives Used as Photosensitizers for Photodynamic Ther , 1999, Photochemistry and photobiology.

[25]  András P Borosy,et al.  Quantitative composition–property modelling of rubber mixtures by utilising artificial neural networks , 1999 .

[26]  David J. Craik,et al.  FUNCTIONAL GROUP CONTRIBUTIONS TO DRUG-RECEPTOR INTERACTIONS , 1985 .

[27]  T Takemura,et al.  The tumour-localizing properties of porphyrin derivatives. , 1990, Journal of photochemistry and photobiology. B, Biology.

[28]  R. Pandey Recent advances in photodynamic therapy , 2000 .

[29]  P Mátyus,et al.  Application of neural networks in structure–activity relationships , 1999, Medicinal research reviews.

[30]  A. Debnath,et al.  Three-dimensional structure-activity analysis of a series of porphyrin derivatives with anti-HIV-1 activity targeted to the V3 loop of the gp120 envelope glycoprotein of the human immunodeficiency virus type 1. , 1994, Journal of medicinal chemistry.

[31]  Kenneth J. Miller,et al.  Additivity methods in molecular polarizability , 1990 .

[32]  P. Broto,et al.  Molecular structures: perception, autocorrelation descriptor and sar studies: system of atomic contributions for the calculation of the n-octanol/water partition coefficients , 1984 .