Abstract

Pharmacophore modelling-based virtual screening of compound is a ligand-based approach and is useful when the 3D structure of target is not available but a few known active compounds are known. Pharmacophore mapping studies were undertaken for a set of 50โ€‰N3-phenylpyrazinones possessing Corticotropin-releasing Factor 1 (CRF 1) antagonistic activity. Six point pharmacophores with two hydrogen bond acceptors, one hydrogen bond donor, two hydrophobic regions, and one aromatic ring as pharmacophoric features were developed. Amongst them the pharmacophore hypothesis AADHHR.47 yielded a statistically significant 3D-QSAR model with 0.803 as value and was considered to be the best pharmacophore hypothesis. The developed pharmacophore model was externally validated by predicting the activity of test set molecules. The squared predictive correlation coefficient of 0.91 was observed between experimental and predicted activity values of test set molecules. The geometry and features of pharmacophore were expected to be useful for the design of selective CRF 1 receptor antagonists.

1. Introduction

Anxiety and depression are among the most common disorders seen in medical practice. The coexistence of anxiety and depression with medical illness is a topic of considerable clinical and research interest [1]. Depression is a serious mental health problem, with significant consequences in terms of human suffering, lost productivity, and even loss of life [2]. Corticotropin-releasing factor 1 (CRF 1) receptor antagonists have been sought since the stress-secreted peptide (adrenocorticotropin-releasing hypothalamic peptide) was isolated in 1981. Although evidence is mixed concerning the efficacy of CRF 1 receptor antagonist as antidepressants, CRF 1 receptor antagonist might be novel pharmacotherapies for anxiety and addiction [3]. Two well-characterized receptor subtypes, CRF 1 and CRF 2, have been identified. These G-protein-coupled receptors are widely distributed throughout the central and peripheral nervous systems [4]. Clinical evidence supports the hypothesis that overproduction of CRF 1โ€‰โ€‰may underlie the pathology of depression, anxiety, and stress-related disorders and suggests that CRF 1 receptor antagonists could be useful for the treatment of these conditions [5].

To reduce the overall cost associated with the discovery and development of a new drug, the computer-aided molecular design methods have been identified as the most promising candidates to focus on the experimental efforts in modern medicinal chemistry. Pharmacophore mapping is one of the major elements of drug design in the absence of structural data of the target receptor. With the aim of providing useful pharmacophoric information for the future efforts in the development of more potent molecules in the series of phenylpyrazinones and to get insight into the structural and molecular properties, the ligand-based 3D-QSAR study was performed using pharmacophore techniques with PHASE module from Schrodinger, New York [6, 7].

2. Experimental

2.1. Dataset

The in vitro biological data of a series of 50 phenylpyrazinones as CRF 1 receptor antagonists were used for the present studies [6]. The CRF 1 antagonistic activity was expressed as IC50, that is, concentration in ฮผm required for 50% inhibition of enzyme activity. The dataset was divided randomly into training set and test set by considering the 70% of the total molecules in the training set and 30% in the test set. The basic structures of the -phenylpyrazinones are shown in Figure 1, and various substituents are listed in Table 1. Thirty-five molecules forming the training set were used to generate pharmacophore model and prediction of the activity of test set (15 compounds) molecules was used as a method to validate the proposed models.

2.2. PHASE Methodology

The 3D-QSAR studies were carried out using PHASE [8โ€“10] version 3.0 implemented in the Maestro 8.5 molecular modeling package from Schrodinger, Molecular Modeling Interface Inc., LLC, New York, NY USA. Phase is a versatile product for pharmacophore perception, structure alignment, activity prediction, and 3D database searching. Phase provides support for lead discovery, SAR development, lead optimization, and lead expansion. Phase may also be used as a source of molecular alignments for third-party 3D-QSAR programs. Phase is well suited to drug discovery projects for which no receptor structure is available. Phase utilizes fine-grained conformational sampling and a range of scoring techniques to identify common pharmacophore hypotheses, which convey characteristics of 3D chemical structures that are purported to be critical for binding. A given hypothesis may be combined with known activity data to create 3D-QSAR models that identify overall aspects of molecular structure that govern activity. These models may be used in conjunction with the hypothesis to mine a 3D database for molecules that are most likely to exhibit strong activity toward the target [11].

2.3. Preparing Ligands

LigPrep [12] was used to attach hydrogen, converts 2D structures to 3D, generates stereoisomer, and, optionally, neutralizes charged structures or determines the most probable ionization state at user-defined pH. All the structures were ionized at neutral pH 7. Conformers for each ligand were generated using ConfGen by applying OPLS-2005 force field method [13, 14] with implicit GB/SA distance-dependent dielectric solvent model at cutoff root mean square deviation (RMSD) of 1 (MacroModel 9.6 2010) with 1,000 iterations using water as solvent.

2.4. Pharmacophore Hypothesis Generation

PHASE can identify the spatial arrangements of functional groups that are common and essential for the biological activity of the ligands under investigation [10, 15]. The most dominating features, hydrogen bond acceptor (A), hydrogen bond donor (D), hydrophobic group (H), and negatively charged group (N) were defined by a set of chemical structural patterns with the requirement that all five match. Pharmacophore-matching tolerance was set to 1โ€‰Aยฐ. Hypotheses were generated by a systematic variation of number of sites and the number of matching active compounds. Common pharmacophore hypotheses (CPH) were considered, which indicated at least five sites common to all molecules. Further, the best CPH was selected depending on the survival score, until at least one hypothesis was found and scored successfully. The hypotheses were scored using default parameters for site, vector, volume, selectivity, number of matches, and energy terms. The regression analysis was performed by constructing a series of models with an increasing number of PLS factors. Pharmacophore-based 3D-QSAR models were generated for the hypotheses using the 35 member training set with three PLS factors and a grid spacing of 1โ€‰Aยฐ. The evaluation of generated CPHs was performed by correlating the observed and the estimated activity for the training set of 35 molecules. PLS analyses were performed in which a series of models were constructed with an increasing number of PLS factors. Score hypotheses step was employed to align the actives to the hypotheses and calculate the score for the actives. CPHs of significant statistical values were selected for molecular alignments.

2.5. Validation of Pharmacophore Model

For accurate and reliable predictions of biological activities of new compounds, the main target was to develop QSAR models, which were statistically robust both internally as well as externally. The data set was divided into a training set and a test set as external validation is considered to be a conclusive proof for judging predictability of a model. The training set was used to generate pharmacophore model and prediction of the activity of test set was used as a method to validate the proposed models. The robustness of the developed pharmacophore hypotheses was internally validated by statistical parameters, that is, squared correlation coefficient ( ) and variance ratio ( ). Validation is a crucial aspect of pharmacophore design, particularly when the model is built for the purpose of predicting activities of molecules in external test series [16]. In the present case, the developed pharmacophore model was externally validated by predicting the activity of test set molecules. The correlation between the experimental and predicted activities of the test set molecules was determined.

3. Results and Discussion

CRF 1 receptor antagonists may offer therapeutic potential for the treatment of diseases resulting from elevated levels of CRF 1 such as anxiety and depression. Efforts to identify structurally diverse CRF 1 receptor antagonists led to the discovery of pyrazinone-based compounds and it was found that CRF 1 receptor binding affinity was improved when the branching point was on the carbon atom bonded directly to the pyrazinone ring [6]. Ligand-based drug design relies on knowledge of other molecules that bind to the biological target of interest. These molecules may be used to derive a pharmacophore which defines the minimum necessary structural characteristics a molecule must possess in order to bind to the target [7]. In other words, a model of the biological target is built based on the knowledge of what binds to it and this model in turn may be used to design new molecular entities that interact with the target.

Thirty-five molecules forming the training set were used to develop the pharmacophore models. The pharmacophoric features selected for creating sites were hydrogen bond acceptor (A), hydrogen bond donor (D), hydrophobic region (H), and aromatic ring (R). Pharmacophore models containing four to six features were generated. The four and five featured pharmacophore hypotheses were rejected due to low value of survival score (less than 2.5), as they were unable to define the complete binding space of the selected molecules. Six featured pharmacophore hypotheses were selected and subjected to stringent scoring function analysis.

The results of six featured pharmacophore hypotheses, labeled AADHHR.47, ADHHHR.203, DHHHRR.611, and AAHHHR.19, are presented in Table 2. The first hypothesis AADHHR.47 is the best hypothesis in this study, characterized by highest survival score (2.980) and best regression coefficient (0.803).

The pharmacophore hypothesis AADHHR.47 is presented in Figure 2. The features represented by this hypothesis are two hydrogen bond acceptors (A), one hydrogen bond donor (D), two hydrophobic regions (H), and one aromatic ring (R). The distances and angles between different sites of AADHHR.47 are given in Tables 3 and 4, respectively.

For each ligand, one aligned conformer based on the lowest RMSE of feature atom coordinates from those of the corresponding reference feature was superimposed on AADHHR.47. The fitness scores for all ligands were observed on the best scored pharmacophore model AADHHR.47. The greater the fitness score, the greater the activity prediction of the compound. The fit function does not only check if the feature is mapped or not but also contains a distance term, which measures the distance that separates the feature on the molecule from the centroid of the hypothesis feature.

Table 5 shows the fitness score for all the molecules of training set. Figure 3 shows the AADHHR.47 aligned with ligand having maximum fitness score, that is, molecule 48 of the training set. Beside this survival score analysis, another validation method to characterize the quality of AADHHR.47 is represented by its capacity for correct activity prediction of training set molecules.

The predicted CRF 1 antagonistic activity of training set molecule exhibited value of 0.803 (RMSD = 0.618) with experimental CRF 1 antagonistic activity using model AADHHR.47 (Figure 4). The validity and predictive character of AADHHR.47 was further assessed by using the test set prediction.

The test set having fifteen molecules was analyzed. All the test set molecules were built and minimized as well as used in conformational analysis like all training set molecules. Then the activities of test set molecules were predicted using AADHHR.47 and compared with the actual activity. Actual and predicted activity values of test set molecules are given in Table 6.

The predicted CRF 1 antagonistic activity of test molecule exhibited value of 0.91 (RMSD = 0.2961) with experimental CRF 1 antagonistic activity using model AADHHR.47 (Figure 4). For a reliable model, the squared predictive correlation coefficient should exceed 0.60 [17, 18]. The results of this study reveal that model AADHHR.47 can be used for the prediction of CRF 1 antagonistic activity.

3.1. 3D QSAR Analysis

Additional insight into the CRF 1 antagonistic activity can be gained by visualizing the 3D QSAR model in the context of one or more ligands in the series with varying activity. This information can then be used to design new or more active analogues. 3D QSAR model based on the molecules of training and test set using various features, that is, hydrogen bond acceptor (A), hydrogen bond donor (D), hydrophobic regions (H), and aromatic ring (R) has been studied.

3.1.1. Hydrogen Bond Donor Field Predictions

The 3D QSAR model based on molecule 48 of the training set using hydrogen bond donor feature is shown in Figure 5(a). Red region near and around the meta- and parahydrogen of pyridine ring substituted at position 3 through NH terminal indicates that the substitutions at these positions by groups having more hydrogen bond donor property favor the CRF 1 antagonistic activity. Green region around the double bond between C5 and C6 on pyridine ring, NH terminal connecting pyrazinone and pyridine ring, methoxy group at C4 on pyridine ring, indicates that substitutions at these positions by groups having hydrogen bond donor property do not favor CRF 1 antagonistic activity.

3.1.2. Hydrogen Bond Acceptor Field Predictions

The 3D QSAR model based on molecule 48 of the training set using hydrogen bond acceptor feature is shown in Figure 5(b). Red region around N1, carbonyl group at C2 and in between C2 and C3, double bond between C5โ€‰โ€‰and C6โ€‰โ€‰(all 3 on pyrazinone ring), indicates that the substitutions at these positions by groups having more hydrogen bond acceptor property favour the CRF 1 antagonistic activity. Green region around C5 on pyrazinone ring, NH terminal and methoxy group on pyridine ring, indicates that the substitutions at these positions by groups having more hydrogen bond acceptor property do not favour the CRF 1 antagonistic activity.

3.1.3. Hydrophobicity Field Prediction

The 3D QSAR model based on molecule 48 of the training set using hydrophobicity feature is shown in Figure 5(c). Green region around ethynyl group substituted at position 5, at C6 (both on pyrazinone ring), double bond between C5 and C6, methoxy group at C4 (both on pyridine ring), indicates that the substitutions at these positions by groups having more hydrophobicity favour CRF 1 antagonistic activity and substitutions at these positions by more hydrophobic groups will result in increase in CRF 1 antagonistic activity. Blue region around methyl substituent of N1, position 4 (both on pyrazinone ring), methoxy group at C4 on pyridine ring, indicates that groups having more hydrophobic property do not favour CRF 1 antagonistic activity.

4. Conclusions

This study shows the generation of a pharmacophore model AADHHR.47 for โ€‰โ€‰phenylprazinones acting as CRF 1 antagonists. Pharmacophore modelling correlates activities with the spatial arrangement of various chemical features. The first hypothesis AADHHR.47 is the best hypothesis in this study, characterized by the best regression coefficient (0.803), degree of freedom (138.5), and highest survival score (2.980). Hypothesis AADHHR.47 represents the best pharmacophore model for determining CRF 1 antagonistic activity. AADHHR.47 consists of two hydrogen bond acceptors, one hydrogen bond donor, two hydrophobic regions, and one aromatic ring features. AADHHR.47 model had strong correlation between experimental and estimated activity of the training ( ) and test ( ) set molecules. Thus, AADHHR.47 pharmacophore model was able to accurately predict CRF 1 antagonistic activity, and the validation results also provide additional confidence in the proposed pharmacophore model. The obtained results suggested that the proposed 3D-QSAR model AADHHR.47 can be useful to rationally design new -phenylprazinones molecules as CRF 1 antagonists and also to identify new promising molecules as CRF 1 antagonists in large 3D database of molecules.