2 D - QSAR studies on CYP26A1 inhibitory activity of 1-[benzofuran-2-yl-(4-alkyl/aryl-phenyl)-methyl]- 1 H-triazoles.

The Quantitative Structure Activity Relationship (QSAR) study is performed over a set of 15, 4-alkyl/aryl-substituted 1- [benzofuran-2-yl-phenylmethyl]-1 H-triazoles derivatives. This study is based on the application of physicochemical parameters in QSAR. The parameters include (MR (molar refractivity), MW (molecular weight), Pc (parachor), St (surface tension), D (density), Ir (index of refraction) and log P (partition coefficient). The parameters describing physiochemical properties are used as independent variables and the biological activity (IC(50)) is considered as dependent variable in multiple regression analysis. Different models were generated with high co-efficient of determination (R(2)). The 2D-QSAR study identified compounds capable of inhibiting the metabolic breakdown of the retinoid (trans-retinoic acid (ATRA)) involved in the activation of specific nuclear Retinoic acid receptors (RARs). This study identifies R115866 as a potential inhibitor of the cytochrome P450 (CYP) mediated metabolism with increased RA levels for retinoid actions.


Background:
In human, CYP26A1 mapped to chromosomes 10q23-q24 is expressed in the liver, heart, pituitary gland, adrenal gland, testis, brain and placenta. It is thought that the principle role of CYP26A1 is homeostatic for the regulation of intercellular ATRA (Anti Trans Retinoic Acid) steady state levels via a negative feedback loop (similar to CYP26A1 cholecalciferol metabolism). The enzyme is an important regulator of differentiation and a possible modulator of disease states in indirectly controlling ATRA and other retinoid concentrations.
Retinoic acid has been used in a number of clinical situations, especially oncology and dermatology. In Oncology ATRA has shown spectacular success in the treatment of acute promyelocytic leukaemia. Although other CYPs are also induced by ATRA, it is thought that CYP26A1 is likely to be the most important enzyme involved in its degradation [1]. ATRA regulates epithelial differentiation and growth through activation of specific nuclear retinoic acid receptors (RARs). We identified compounds capable of inhibiting the metabolic breakdown of the retinoid because high rate metabolism largely impairs the biological efficiency of RA.
QSAR studies are of importance in molecular biochemistry. It is essential that appropriate descriptors are employed, whether they are theoretical, empirical or derived from available experimental characteristics of structure to obtain significant correlation. QSAR represents an attempt to correlate structural or property descriptors of compound with activities. These physiochemical descriptors include parameters to account for hydrophobicity, topology, electronic parameters and steric effects. These are determined empirically by computational methods. Quantitative Structure Activity Relationship (QSAR) is currently being applied in many disciplines pertaining to drug design and environmental risk assessments [2]. ISSN  An insight to structure activity relationship is established using molecular descriptors that effectively characterize molecular size, molecular branching or the variations in molecular shapes affecting the structure and its molecular activity. Structural invariants derived from molecular structure are important for the understanding of a particular structure/property activity relationship. Several graph theoretical invariants have been generalized to produce structure dependent descriptors. It should be noted that most biological activities are dominated by molecular size in such generalization.
The activities and properties are related using the general mathematical function, F: Biological activity = f [structure (physicochemical descriptors as structural parameters)].
Biological activity is often defined as a measure of log 1/C, ki, IC₅₀, EC₅₀, log k and Km as described elsewhere [3][4][5]. The relationship is often not a mathematical expression derived by statistical or related techniques. The parameters describe structural and physiochemical properties as independent variables and the biological activities as dependent variable for usage in multiple linear regression (MLR) analysis. We describe the QSAR of 15 triazoles derivatives with relationship between structural and physiochemical parameters to known biological activity log IC50 values.

Methodology:
The physiochemical properties used in the study are molar refractivity, molar volume, parachor, index of refraction, surface tension and polarizibility to co-relate with biological activity (log IC50 inhibitory concentration). The physiochemical properties were estimated using Chemsketch 5.0 (from ACD Lab) and multiple linear regressions are done by SPSS Software [4,6]. The physicochemical parameters used in the study are described below.

Density (D):
Density is a steric parameter calculated by ACD Lab Chem Sketch Software. This parameter is related to the bulk and size of the substituent.

Molecular weight (MW):
Molecular weight descriptor has been used as a descriptor in systems such as transport studies where diffusion is the mode of operation. It is an important variable in QSAR studies pertaining to cross resistance of various drugs in multi-drug resistant cell lines.

Molar refractivity (MR):
Molar refractivity is generally considered to be a measure of overall bulkiness and is related to London dispersion forces using MR=4πNα/3; where N is Avogadro number and α is the polarizability of the molecule. It gives no information about shape and is generally scaled 0.1 and has been extensively used in QSAR.

Parachor (Pc):
This parameter is calculated by ACD Lab Chem Sketch Software (Chem Sketch). The Parachor may be defined as the molar volume of a liquid at a temperature at which its surface tension is unity. Parachor is obtained experimentally using Pc= (MW/d) γ ¼ where, Pc is a molar volume term that is corrected by surface tension (γ) rise to the power 1/4.

Indicator variables:
These are not QSAR parameters but are used to indicate the significance of any particular group or species at a particular substitution site in a given series of drugs. Their activity data and the physicochemical parameters evaluated in the correlation are listed in (Table 3, supplementary material). Linear regression analysis was performed on present series of triazoles derivaties according to physiochemical properties given in (Table 2, See supplementary material). These parameters were found to be useful in QSAR based drug modeling (please see supplementary material for modeling equations 1 to 7). The effect of substituent indicator parameters -I6 (CH2CH2CH2), + I₂₀ (CH-(CH2CH3)2) at R (see Figure 1a), IR2 position was studied (see Table 2, supplementary material) in equations 2 to 7.

Discussion
Analysis through univariate correlation produced very low correlation coefficient (see equation 1 in supplementary material). Bivariate analysis gave a slightly improved value. Tetra and penta variate analysis results are encouraging with high correlation coefficient. The high value of F-test ratio and low value of standard error of estimation support reports described elsewhere [7,8].
We estimated the log IC50 values with the observed ones (

Conclusion:
The study shows that the biological activity log IC₅₀ is structure specific in nature for a particular series of derivatives.