3D-QSAR CoMFA study of some heteroarylpyrroles as possible anticandida agents

*Corresponding Author:

P. C. Sharma
University institute of pharmaceutical sciences, Kurukshetra university, Kurukshetra-136 119, India
E-mail: sharma_prabodh@rediffmail.com

Abstract

A three dimensional quantitative structure-activity relationship study using the comparative molecular field analysis method was performed on a series of 3-aryl-4-[α-(1H-imidazol-1-yl) aryl methyl] pyrroles for their anticandida activity. This study was performed using 40 compounds, for which comparative molecular field analysis models were developed using a training set of 33 compounds. Database alignment of all 33 compounds was carried out by root-mean-square fit of atoms and field fit of the steric and electrostatic molecular fields. The resulting database was analyzed by partial least squares analysis with cross-validation; leave one out and no validation to extract optimum number of components. The analysis was then repeated with bootstrapping to generate the quantitative structure-activity relationship models. The predictive ability of comparative molecular field analysis model was evaluated by using a test set of 7 compounds. The 3D- quantitative structure-activity relationship model demonstrated a good fit, having r 2 value of 0.964 and a cross validated coefficient r 2 value as 0.598. Further comparison of the coefficient contour maps with the steric and electrostatic properties of the receptor has shown a high level of compatibility and good predictive capability.

Keywords

CoMFA, 3D-QSAR, Pyrroles, Anticandida activity

Fungal infections in the human, range from superficial and cutaneous (such as dermatomycosis) to deeply invasive and disseminated (such as candidiasis and cryptococcosis) infections. In the past 20 years, fungal infections have increased dramatically-paradoxically, as risk advances. Fungal infections occur more frequently in people whose immune system is suppressed (because of organ transplantation, cancer chemotherapy, or the acquired immune deficiency syndrome) or who have been subject to invasive procedures (catheters, prosthetic device) [1]. Fungal infections are now important causes of morbidity and mortality among immuno-compromised hospitalized patients. The frequency of invasive candidiasis has increased ten-fold during the past decade [2].

Candida albicans (CA) has been identified as the major opportunistic pathogen in the etiology of fungal infections; however, the frequency of other Candida species is increasing [3]. The current standard of therapies is the fungicidal (but toxic) polyene antibiotic, amphotericin B, and the safer (but fungistatic) azoles. In particular, the latter class of drugs is an important antifungal class widely used for AIDS-related mycotic pathologies [4].

Quantitative structure activity relationship (QSAR) enables the investigators to establish a reliable quantitative structure-activity and structure-property relationships to derive in silico QSAR models to predict the activity of novel molecules prior to their synthesis. The overall process of QSAR model development can be divided into three stages namely, data preparation, data analysis, and model validation, representing a standard practice of any QSAR modeling. Successful application of 3D-QSAR methodologies have been used to generate models for various chemotherapeutic agents [5,6].

We have carried out 3D-QSAR studies employing comparative molecular field analysis [5] (CoMFA) techniques in order to study and gain further insight to deduce a correlation between structure and biological activity of 3-aryl-4-[α-(1H-imidazol-1-yl) aryl methyl] pyrroles as potent anticandida agents [7].

In the CoMFA method, introduced by Crammer [8,9], a relationship is established between the biological activities of a set of compounds and their steric and electrostatic properties. An advantage of CoMFA is its ability to predict the biological activity of molecules and represent the relationship between steric and electrostatic properties and biological activity in the form of contour maps [10]. An ‘active conformation’ of the ligands is generated and superimposed as per the predefined rules. These molecules are then placed in a box of predefined grid size. The steric and electrostatic interaction energy between each structure and a probe atom of defi ned size and charge are calculated at each grid point using the molecular mechanics force fields. A multivariate data analysis technique like partial least squares (PLS) [11-13] is used to derive a linear equation from the resulting matrices. PLS is used in combination with cross validation to obtain the optimum number of components. This ensures that the QSAR models are selected on their ability to predict the data rather than to fit the data. The advantages of CoMFA studies are in the ability to predict the target properties of the compounds and to graphically present the QSAR in the form of coefÞ cient contour maps [14].

We present here 3D-QSAR studies using CoMFA method on a series of 3-aryl-4-[α-(1H-imidazol-1-yl) aryl methyl] pyrroles and the contour maps derived reveal the significance of steric and electrostatic fields. The structural variations in the molecular fields at particular regions in the space provide underlying structural requirements and 3D-QSAR models generated give good predictive ability and aid in the design of potent anticandida agents.

Materials and Methods

Biological activity data

The antifungal activity data against Candida albicans for a series of 3-aryl-4-[α-(1H-imidazol-1-yl) aryl methyl] pyrroles containing 40 compounds as anticandida agents was used in this analysis. General structure of the compounds is shown in (fig. 1). Training set was formed by selecting 33 compounds from the original series. Test set compounds were no. 11, 12, 33, 34, 35, 37 and 42 (total 7 compounds), selected randomly. These compounds were not included in the analysis to generate the CoMFA model. The robustness and predictive ability of models were evaluated by selecting biological activity with chemical class similar to training set. CoMFA techniques were used to derive 3D-QSAR models for 3-aryl-4-[α-(1Himidazol- 1-yl)aryl methyl)pyrroles. The MIC data were used for the QSAR analysis as a dependent parameter, after converting to the reciprocal of the logarithm of MIC (pMIC) expressed in μM/ml (Table 1).

Figure 1: Heteroaryl pyrroles used for CoMFA study. General structure of compounds used for CoMFA study. Rotatable bonds selected for conformational analysis are shown in bold face.

Table 1: Experimental activities of molecules used in training set and test set

Molecular modeling

A database of 33 compounds forming the training set was generated by molecular modeling. All molecular modeling and 3D-QSAR studies were performed using SYBYL 6.7 [15] with TRIPOS Force Field [16] on a Silicon Graphics O₂ workstation with IRIX operating system. The crystallographic data for these ligand complexes was not available hence all the molecules were constructed using a grid having a spacing of 1.54 A0 between grid points. This is the default spacing, which represents sp³ carboncarbon bond length. The molecules were cleaned up and quick minimized after sketching. Because no experimental data on the biologically relevant conformations of the selected compounds were available (for example, atomic coordinates derived from X-ray crystallographic studies of their complexes with the putative receptor), we resorted to a general molecular mechanics approach (AM1) to build the conformational models to be used for generation of CoMFA models. A chirality check was performed to identify chiral atoms, after adding hydrogens, it was important to consider all possible enantiomers as the activity was reported for racemic mixtures. Then the molecules were subjected for energy minimization (geometry optimization) at a gradient of 1.0 kcal/mol with delta energy change of 0.001 kcal/mol with the TRIPOS standard force Þ eld. Structures were drawn by using default setting of SYBYL.

After conformational analysis (CA) by adopting AM1 Hamiltonian approach, the least energy conformation was selected, saved, and used for the charge calculation, assuming that it was the active one. We have used two different types of charges, calculated using the Gasteiger-Marsili method and the semi-empirical MOPAC method [17].

Partial Least Squares analysis (PLS)

The PLS analyses were done by following standard protocols [18]. In order to speed up the analysis and reduce the amount of noise, column filter was used by excluding the columns with a variance smaller than 2.0. Equal weights for the steric and electrostatic descriptors were assigned using the CoMFA scaling option.

CoMFA Results

Two CoMFA models were generated by using different types of partial atomic charges, results of which are shown in Table 2. Model A was derived using charges calculated according to Gasteiger-Marsili method, while Model B was obtained using MOPAC charges. From the results it can be observed that both the models are significant in term of their statistical acceptance, however model B was considered to be better due to higher correlation coefficient and Fischer’s statistical value.

Table 2: Summary of comfa results

Prediction of Activity

The 3-D QSAR analysis obtained as Model B was used for predicting the activity of the 33 compounds in the training set. The results are shown in Table 3. From the table, it can be seen that the predicted activities are very close to the experimental activities with minimum residual activity.

Table 3: Predicted and experimental activities of training set

CoMFA contour maps

The QSAR produced by CoMFA were represented as a 3-D coefficient contour map. To visualize the CoMFA steric and electrostatic fields from PLS analysis, contour maps of the product of the standard deviation associated with CoMFA column and coefficient (SD Coeff) at each lattice point were generated. The contour maps were plotted as percentage contribution to the QSAR equation and were associated with difference in biological activity.

The CoMFA contour maps generated for model B were used to explain the structure activity relationship of antifungal drugs.

In CoMFA contour maps, the regions of high and low steric tolerance are shown in green and yellow polyhedra, respectively. CoMFA electrostatic field are shown as blue and red polyhedra in fig. 2. A low electron density within the inhibitors near blue and red polyhedra, respectively increase or decrease the activity.

Figure 2: CoMFA steric and electrostatic contours field plot. Green contours in mesh view indicate regions where bulky groups increase activity, whereas yellow contours indicate regions where bulky groups decrease activity. Blue contours indicate where positive electrostatic groups increase activity, whereas red contours indicate where negative electrostatic groups increase activity. Compound 1 from training set is shown in fig. 2A and database alignment of 33 molecules of training set is shown in fig. 2B.

Results and Discussion

The validity of Model B was further enhanced by bootstrapping process. Bootstrapping of 10 runs gave r² of 0.986±0.005 with a very low standard error of 0.046±0.007 which added to the high confidence limit to this analysis. It can be seen that, both steric and electrostatic fi elds contributed to the QSAR equation by 42.8 and 57.2%, respectively. This suggested that variation in the antifungal activity is predominantly determined by electrostatic properties. Thus the results suggested that there is a good internal consistency in the data set generated in model B.

Model B performed exceptionally well in predicting the activity of most compound used in the test set. However, it must be emphasized that molecular alignment and conformations used in this study were selected in the absence of X-ray crystallographic coordinates of these molecules; still, the CoMFA model generated in the study showed very good prediction capability. From the Table 4, it can be observed that the predictions made using CoMFA model were satisfactory in most cases. In general, the percentage difference in the predicted activity of the synthesized compound ranges from 1.2 to 28.5%. The relative difference in the prediction is not unexpected and is within the acceptable limits.

Table 4: Comparative variation in the experimental and predicted anticandida activity of compounds in test set

From these results, it is inferred that the 3-D QSAR model generated in this study has a potential to predict the activity of diverse compounds belonging to similar structural class. The investigations concerning the design of new chemical entities based on the proposed CoMFA models, predicting their antifungal activity prior to the synthesis would be part of our forthcoming communication.

Acknowledgements

The authors wish to thank His Holiness Jagadguru Sri Sri Shivarathri Deshikendra Mahaswamigalavaru of Suttur Mutt, Mysore, for providing the facilities to carry out the research work and All India Council of Technical Education (AICTE), New Delhi for financial support through MODROBS scheme.

References

1. Diamond RD. The growing problem of mycoses in patients infected with the human immunodeficiency virus. Rev Infect Dis 1991;3:480-6.
2. Georgopapadakou NH, Walsh TJ. Human mycoses: Drugs and targets for emerging pathogens. Science 1994;15:371-3.
3. Beck-Sague CM, Jarvis WR. Secular trends in the epidemiology of nosocomial fungal infections in the United States, 1980-1990. J Inf Dis 1993;167:1247-51.
4. Koltin Y. Targets for antifungal drug discovery. Annu Rep Med Chem 1990;25:141-8.
5. Gokhale VM, Kulkarni VM. Understanding the Antifungal Activity of Terbinafine Analogues Using Quantitative Structure-Activity Relationship (QSAR) Models. Bioorg Med Chem 2000;8:2487-99.
6. Karki RG, Kulkarni VM. Three-dimensional Quantitative Structure-Activity Relationship (3D-QSAR) of 3-Aryloxazolidin-2-one Antibacterials. Bioorg Med Chem 2001;9:3153-60.
7. Artico M, Santo RD, Costi R, Massa S, Retico A, Artico M, et al. Antifungal Agents. 3-Aryl-4-[a-1(H -imidazol-l-yl)arylmethyl]pyrroles:A New Class of Potent Anticandida Agents. J Med Chem 1995;38:4223-33.
8. Cramer RD, Patterson DE, Bunce JD. Comparative molecular field analysis (CoMFA): Effect of shape on binding of steroids to carrier proteins. J Am ChemSoc 1988;110:5959-67.
9. Clark M, Cramer RD 3rd, Jones DM, Patterson DE, Simeroth PE. Comparative molecular field analysis (CoMFA) Toward its use with 3D structural databases. Tetrahedron Comput Met 1990;3:47-59.
10. Bhongade BA, Gadad AK. 3D-QSAR CoMFA/CoMSIA studies on Urokinase plasminogen activator (uPA) inhibitors: A strategic design in novel anticancer agents. Bioorg Med Chem 2004;12:2797-805.
11. Wold S, Ruhe A, Wold H, Dunn III. The collinearity problem in linear regression: The partial least square approach to generalized inverse. J Sci Stat Comput 1984;5:735-43.
12. Cramer RD, Bunce JD, Patterson DE, Crossvalidation, Bootstrapping, and Partial Least Squares Compared with Multiple Regression in Conventional QSAR Studies. Quant Str Act Relat 1988;7:18-25.
13. Wold S, Albano C, Dunn WJ, Edlund U, Esbensen K, Geladi P, et al. Multivariate data analysis in chemistry. In: Kowalski BR, editor. Chemometrics - Mathematics and Statistics in Chemistry, Reidel, Dordrecht; 1984. p. 17-95.
14. Puntambekar D, Giridhar R, Yadav MR. 3D-QSAR studies of farnesyltransferase inhibitors: A comparative molecular field analysis approach. Bioorg Med ChemLett 2006;16:1821-7.
15. Sybyl 6.7, Tripos Associates, 1699 South Hanley Road, St. Louis, USA, MO 63144.
16. Matthew C, Cramer RD, Opdenbosch NV. Validation of the general purpose tripos 5.2 forceÞeld. J ComputChem 1989;10:982-1012.
17. Gasteiger J, Marsili M. Iterative partial equalization of orbital electronegativity-a rapid access to atomic charges. Tetrahedron 1980;36:3219-28.
18. GeladiP, Notes on the history and nature of partial least squares (PLS) modeling. J Chemom 1988;2:231-46.