|Year : 2012 | Volume
| Issue : 2 | Page : 123-133
3D-QSAR studies on CCR2B receptor antagonists: Insight into the structural requirements of (R)-3-aminopyrrolidine series of molecules based on CoMFA/CoMSIA models
Swetha Gade1, Shaik Mahmood2
1 Department of Botany, Bioinformatics Division, Osmania University, Hyderabad, Andhra Pradesh, India
2 Environmental Microbiology Lab, Osmania University, Hyderabad, Andhra Pradesh, India
|Date of Submission||19-Jan-2011|
|Date of Decision||14-May-2011|
|Date of Acceptance||17-Jul-2011|
|Date of Web Publication||10-Apr-2012|
Department of Botany, Bioinformatics Division, Osmania University, Hyderabad, Andhra Pradesh
Source of Support: None, Conflict of Interest: None
| Abstract|| |
Objective: Monocyte chemo attractant protein-1 (MCP-1) is a member of the CC-chemokine family and it selectively recruits leukocytes from the circulation to the site of inflammation through binding with the chemotactic cytokine receptor 2B (CCR2B). The recruitment and activation of selected populations of leukocytes is a key feature in a variety of inflammatory conditions. Thus MCP-1 receptor antagonist represents an attractive target for drug discovery. To understand the structural requirements that will lead to enhanced inhibitory potencies, we have carried out 3D-QSAR (quantitative structure-activity relationship) studies on (R)-3-aminopyrrolidine series of molecules as CCR2B receptor antagonists. Materials and Methods: Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were performed on a series of (R)-3-aminopyrrolidine derivatives as antagonists of CCR2B receptor with Sybyl 6.7v. Results: We have derived statistically significant model from 37 molecules and validated it against an external test set of 13 compounds. The CoMFA model yielded a leave one out r 2 (r 2 loo ) of 0.847, non-cross-validated r 2 (r 2 ncv ) of 0.977, F value of 267.930, and bootstrapped r 2 (r 2 bs ) of 0.988. We have derived the standard error of prediction value of 0.367, standard error of estimate 0.141, and a reliable external predictivity, with a predictive r 2 (r 2 pred ) of 0.673. While the CoMSIA model yielded an r 2 loo of 0.719, r 2 ncv of 0.964,F value of 135.666, r 2 bs of 0.975, standard error of prediction of 0.512, standard error of estimate of 0.180, and an external predictivity with an r 2 pred of 0.611. These validation tests not only revealed the robustness of the models but also demonstrated that for our models r 2 pred, based on the mean activity of test set compounds can accurately estimate external predictivity. Conclusion: The QSAR model gave satisfactory statistical results in terms of q 2 and r 2 values. We analyzed the contour maps obtained, to study the activity trends of the molecules. We have tried to demonstrate structural features of compounds to account for the activity in terms of positively contributing physicochemical properties such as steric, electrostatic, hydrophobic, hydrogen bond donor, and acceptor fields. These contour plots identified several key features, which explain the wide range of activities. The results obtained from models offer important structural insight into designing novel CCR2B antagonists before their synthesis.
Keywords: (R)-3-aminopyrrolidine series, 3D-QSAR, CCR2B, CoMFA, CoMSIA
|How to cite this article:|
Gade S, Mahmood S. 3D-QSAR studies on CCR2B receptor antagonists: Insight into the structural requirements of (R)-3-aminopyrrolidine series of molecules based on CoMFA/CoMSIA models. J Pharm Bioall Sci 2012;4:123-33
|How to cite this URL:|
Gade S, Mahmood S. 3D-QSAR studies on CCR2B receptor antagonists: Insight into the structural requirements of (R)-3-aminopyrrolidine series of molecules based on CoMFA/CoMSIA models. J Pharm Bioall Sci [serial online] 2012 [cited 2020 Jan 25];4:123-33. Available from: http://www.jpbsonline.org/text.asp?2012/4/2/123/94813
Chemokines or chemotactic cytokines are small molecular weight (6-15 KD) proteins that modulate inflammatory and immune responses through promotion of cell migration, cell adhesion, and transmigration. ,,,,,, They are divided into four classes dependent on the arrangement of conserved cysteines in the N-terminal region. Monocyte chemo attractant protein-1 (MCP-1) is a member of the CC-chemokine family and it selectively recruits leukocytes from the circulation to the site of inflammation through binding with the chemotactic cytokine receptor 2B (CCR2b). , This receptor is a member of the seven transmembrane receptor families (7-TM) and is expressed on the surface of monocytes and macrophages.  The recruitment and activation of selected populations of leukocytes is a key feature of a variety of inflammatory conditions. This response is crucial for host defense during inflammation, but the secretory products of white blood cells may increase injury by damaging surrounding healthy tissue. ,,, These effects are mediated principally through activation of intracellular signaling pathways following binding of MCP-1 to the CCR2b. MCP-1 is a potent chemotactic and activating factor for monocytes and memory T-cells. It regulates adhesion molecule expression and cytokine production.  It also induces superoxide anion and lysosomal enzyme release from human monocytes.  The role of MCP-1 is seen in pathophysiology of a wide range of acute and chronic inflammatory diseases such as rheumatoid arthritis, ,, atherosclerosis, ,,, asthma, ,, psoriasis, , and transplant rejection. ,,, This involvement is evidenced by studies showing elevated MCP-1 expression correlated with leukocyte infiltration in vivo, ,, the use of neutralizing antibodies, , and through both animal receptor  and ligand  knockout studies.
An MCP-1 receptor antagonist thus represents an attractive target for drug discovery, and this has prompted an intense period of pharmaceutical research. Several companies have reported the discovery of potent small molecule antagonists of the CCR2b, showing varying degrees of selectivity over closely homologous receptors. ,,,,,
Nowadays, the use of three-dimensional quantitative structure-activity relationship (3D-QSAR) techniques, such as comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA), ,, is routine in modern drug design to help understand drug-receptor interaction. It has been shown in the literature that these computational techniques can strongly support and help the design of novel, more potent inhibitors by revealing the mechanism of drug-receptor interaction. ,, In this study, we have developed predictive 3D-QSAR models using (R)-3-aminopyrrolidine derivatives as antagonists of the CCR2b. 
| Materials and Methods|| |
We took the in vitro inhibitory activity data (IC50, nM) of a series of (R)-3-aminopyrrolidine derivatives, reported by Moree et al,  for the study. We have omitted 12 molecules whose IC50 values were not quantitatively reported from the given 62 antagonists.
We have considered the remaining 50 compounds for QSAR study. The IC50 values used in this study span a range of three orders of magnitude ranging from 3.2 to 4330 nM. Although major portion of compounds weighted toward the high-potency end of the spectrum, activities are acceptably distributed across the range of values. Thus, these molecules provided a broad and homogenous dataset for the 3D-QSAR study. The generation of reliable models is dependent on the creation of appropriate training and test sets. We have carefully chosen the training set molecules considering the following rules: composition of the QSAR training and test sets individually is necessarily representative of the whole data set in terms of structural complexity. Thus, these sets should show the appropriate chemical diversity and distribution of biological property across the range of IC50 values. The most active and the least active molecules should always be a part of training set, as these are the representative compounds for diversity in biological property. Thus, data set was divided into training (37 compounds) and test (13 compounds) sets [Table 1] in the ratio of 2.85:1.
|Table 1: Structures and biological activities of the training and test set molecules|
Click here to view
We have used the partial least squares (PLS) method for all QSAR analyses. We have converted CCR2b inhibitory activities into the corresponding PIC50 (-log IC50) values and used them as dependent variables, whereas CoMFA and CoMSIA descriptors are used as independent variables in the PLS regression analyses to derive QSAR models. The predictive ability of the models was assessed by their q2 values.
We performed the three-dimensional structure building and all modeling using the Sybyl program package, version 6.7 on a Silicon Graphics Fuel workstation.  CoMFA and CoMSIA studies require that the 3D structures of the molecules to be analyzed should be aligned according to a suitable conformational template, which is assumed a bioactive conformation.  As the bioactive conformations of these inhibitors are not known, the lowest energy conformations were reasonable initial structures to perform 3D-QSAR calculations. We performed energy minimizations using the Tripos force field  and the Gasteiger-Hückel  charge with a distance-dependent dielectric and Powell conjugate gradient algorithm. The criterion of convergence was 0.05 kcal/mol. Subsequently, the lowest energy conformation of each structure was submitted to optimization with the semi-empirical program MOPAC 6.0 and applying the AM1 Hamiltonian. 
One of the most important adjustable parameters in 3D-QSAR is the relative alignment of all the molecules to one another so that they have a comparable conformation and a similar orientation in space.  Usually the bound conformation of a ligand from X-ray crystallographic or NMR studies makes a good starting point for the alignment. However, due to the unavailability of experimentally determined structures, we have opted to take the most potent inhibitor of the data set, compound 71  as a template for superimposition, assuming that its conformation represents the most bioactive conformation of the (R)-3-aminopyrrolidine derivatives. Each analog was aligned to the template by rotation and translation to minimize the RMSD between atoms in the template and the corresponding atoms in the analog using the DATABASE-ALIGN option in SYBYL. This method involves the alignment of molecules in a structurally and pharmacologically reasonable manner based on the assumption that each compound acts via a common macromolecular target-binding site. The most active compound is shown in [Figure 1], and the aligned compounds are shown in [Figure 2].
To derive the CoMFA descriptor fields, a 3D cubic lattice with grid spacing of 2Å in x, y, and z directions was created to encompass the aligned molecules. In CoMFA, Lenard-Jones (6-12 interactions), the steric interaction field, and Columbic electrostatic potentials (1/r) were calculated at each lattice intersection. The grid box dimensions were determined automatically in such a way that region boundaries were extended beyond four Å in each direction from coordinates of each molecule. The Vander Waals potentials and Columbic terms, which represent steric and electrostatic fields, respectively, were calculated using Tripos force field. A sp3 hybridized carbon atom with radius 1.52 Å bearing +1 charge served as probe atom to calculate steric and electrostatic fields.
The CoMFA steric and electrostatic fields generated were scaled by the CoMFA standard option available in SYBYL. A PLS approach, an extension of multiple regression analysis, was used to derive 3D-QSAR, in which the CoMFA descriptors were used as independent variables and PIC50 values as dependent variables .,, We performed cross-validation analysis for selecting the model that is most likely to have predictive values. The intensity of the cross-validation process is controlled by selecting the number of 10 groups. After the optimum number of components was determined, we performed a non-cross-validated analysis. We have computed the r2 cv, PRESS (the root mean predictive error sum of squares), r2ncv , F value, and standard error of estimate values according to the definition in the SYBYL. The cross-validated coefficient was calculated using the following equation
where γpred , γactual , and γmean are predicted, actual, and mean values of the target property (PIC50), respectively. We have used the following formula to calculate lowest standard error of prediction
The non-cross-validated PLS analyses were performed with column filtering value of 2.0, to reduce analysis time with small effect on the q 2 values. We further assessed the robustness and statistical confidence of the derived models, through bootstrapping analysis for 100 runs. We have examined the predictive power of 3D-QSAR models, derived by using the training set with an external test set of 13 molecules. The predictive ability of the models is expressed by the predictive r 2 (r 2 pred) value, which is analogous to cross-validated r 2 (r 2 cv ) and is calculated using the following formula:
where SD (standard deviation) is the sum of the squared deviations between the biological activities of the test set and the mean activity of the training set molecules and PRESS is the sum of the squared deviations between predicted and actual activities for every compound in the test set. The activity of the test set was predicted by the CoMFA model using the predict Command. CoMFA coefficient maps were generated by interpolation of the pair wise products between the PLS coefficients and the standard deviations of the corresponding CoMFA descriptor values.
CoMSIA was performed to evaluate steric, electrostatic, hydrophobic, hydrogen bond donor, and hydrogen bond acceptor properties of molecules by using the standard options in SYBYL. The steric, electrostatic, hydrophobic, H-bond donor, and H-bond acceptor fields were calculated separately using the sp3 carbon atom probe with a charge of 1 provided in SYBYL7.0. Similar to CoMFA, a data table has been constructed from similarity indices calculated at the intersections of regularly spaced lattice (2Å spacing). We have calculated similarity indices AF, K between the compounds of interest, and a probe atom according to the following equation:
where q is the grid point for molecule j; ωik is the actual value of the physic chemical property k of atom i; ωprobe, k indicates probe atom with charge 1, radius1Å, hydrophobicity 1, H-bond donor, and acceptor property 1; α is an attenuation factor; and r 2iq is the mutual distance between the probe atom and grid point q and atom i of the test molecule. The default value of α is 0.3.
The predictive power of CoMFA and CoMSIA models was further validated by using an external test set (inhibitors marked with "d" in [Table 1]). The inhibitors in the test set were given exactly the same pretreatment as the inhibitors in the corresponding training set. The correlation between the experimental and predicted activity for models was calculated as r 2pred value. We have also performed a cross-validation that is based on Fischer randomization test method.
| Results and Discussion|| |
We have used CoMFA and CoMSIA techniques to derive 3D-QSAR models on novel series of (R)-3-aminopyrrolidine-based compounds acting as CCR2b antagonists. The biological activity of negative logarithm PIC50 was used as a dependent variable.
We have used the low-energy conformer obtained from the AM1 optimization as template and aligned all compounds using DATABASE ALIGNMENT method. We generated various 3D-QSAR models and selected the best one based on statistically significant parameters obtained. We obtained the final model with 37 and 13 molecules in the training and test sets, respectively. The predictive power of the 3D-QSAR models, derived using the training set, was assessed by predicting biological activities of the test set molecules. In 3D-QSAR studies r 2 of 0.3 is considered statistically significant.  In view of it, we can consider the models having r 2 > 0.5 much better and statistically significant. Final predicted versus experimental PIC50 values for both CoMFA and CoMSIA models and their residuals (for training and test set compounds) are given in [Table 2].
|Table 2: Predicted and residual activities of the training and test set molecules from CoMFA and CoMSIA analyses|
Click here to view
The CoMFA analysis with contribution of the steric and electrostatic fields [Table 2] yielded a cross-validated, r2cv = 0.847 with five components, non-cross-validated r2ncv of 0.977, a conventional r2 (leave one out), r2loo of 0.856, an F value 267.930, and a predictive r2 , r2pr of 0.673. The results of CoMFA study are given in [Table 3]. The graphs of actual versus predicted activities for the training and test sets of molecules are depicted in [Figure 3]. CoMFA contours were generated using this model. To further assess the robustness of the model, bootstrapping analysis (100 runs) was performed and an r2bs of 0.988 (S.Dbs 0.005) was obtained, further establishing the strength of the model. [Figure 4] shows the histogram of residual values obtained from CoMFA analysis. The steric and electrostatic contributions were found to be 54.6% and 45.4%, respectively. We have further used data set and alignment of CoMFA for CoMSIA analysis.
|Figure 3: Graph of actual versus predicted activity of training and test set molecules from CoMFA analysis|
Click here to view
CoMSIA is similar to CoMFA but uses a Gaussian function rather than Columbic and Lennard-Jones potentials to assess the contribution from different fields. CoMSIA was performed using steric, electrostatic, hydrophobic, hydrogen bond donor, and hydrogen bond acceptor fields. 3D-QSAR models were generated using all the above fields, and the results of study are summarized in [Table 4].
The CoMSIA model yielded the cross-validated q2 = 0.715 with six components, non-cross-validated R2 of 0.964, F value of 135.666, bootstrapped R2 of 0.975, and a predictive r2 of 0.611. The steric, electrostatic, hydrophobic, donor, and acceptor field's contributions were 15.8%, 27.2%, 8.2%, 29.3%, and 19.5%, respectively. The graph of the actual versus predicted biological activities for the molecules is shown in [Figure 5]. Histogram of residual values obtained from CoMSIA analysis is depicted in [Figure 6].
The field values generated at each grid point were calculated as the scalar product of the associated QSAR coefficient and the standard deviation of all values in the corresponding column of the data table (STDDEV*COEFF), plotted as the percentage contributions to QSAR equation. The CoMFA and CoMSIA contour maps developed are shown in [Figure 7] and [Figure 8], respectively.
|Figure 7a: CoMFA STDEV COEFF contour maps: steric fields; green contours indicate regions where bulky groups increase activity, while yellow contours indicate regions where bulky groups decrease activity. Active compound 71 is displayed in the background for reference |
Figure 7b: CoMFA STDEV COEFF contour maps: electrostatic fields; blue contours indicate regions where electropositive groups increase activity, while red contours indicate regions where electronegative groups increase activity. Active compound 71 is displayed in the background for reference
Click here to view
|Figure 8a: CoMSIA STDEV*COEFF contour maps: steric fields; green contours indicate regions where bulky groups increase activity, while yellow contours indicate regions where bulky groups decrease activity Active compound 71 is displayed in the background for reference|
Figure 8b: CoMSIA STDEV*COEFF contour maps: electrostatic fields; blue contours indicate regions where electropositive groups increase activity, while red contours indicate regions where electronegative groups increase activity. Active compound 71 is displayed in the background for reference
Figure 8c: CoMSIA STDEV*COEFF contour maps: Acceptor fields; the magenta contour for H-bond acceptor group increase activity, red indicates the disfavor region. Active compound 71 is displayed in the background for reference
Figure 8d: CoMSIA STDEV*COEFF contour maps: Donor fields; the cyan contour for H-bond donor favor region, purple indicates the disfavor region. Active compound 71 is displayed in the background for reference
Figure 8e: CoMSIA STDEV*COEFF contour maps: Hydrophobic fields; the yellow contour for hydrophobic favor region, white indicates the hydrophilic favored region. Active compound 71 is displayed in the background for reference
Click here to view
In model validation, we examined the internal predictive power of the models and their ability to reproduce biological activities of the compounds for the training set. The computed CCR2b inhibitory activity from the CoMFA and CoMSIA showed a good correlation with experimental inhibitory activity.
We have performed external validation of QSAR models to verify the excellent statistical parameters that were obtained and to investigate whether the activity of (R)-3-aminopyrrolidine derivatives from external data could be predicted well with this model. We have selected a set of 13 compounds, as a test set, from the 50 compounds for the validation experiments. The ultimate test for predictability of QSAR analysis in the drug design process is to predict the biological activity of compounds that have not been included in the training set. The r2pred was calculated and we have obtained values of 0.673 and 0.611 for CoMFA and CoMSIA, respectively. Thus, the CoMFA model displays higher predictivity both in regular cross-validation and in the prediction of the test compounds.
Fischer statistics (F test)
Fischer statistics (F) is the ratio between explained and unexplained variance for a given number of degrees of freedom. F test is a variance-related statistics that compares two models differing by one or more variables to see if the more complex model is more reliable than the less complex one. The model is supposed to be good if the F test is above a threshold value, i.e., tabulated value. The larger the value of F, the greater is the probability that the QSAR equation is significant. The F values for the CoMFA and CoMSIA models were 267.930 and 135.666, respectively.
Visualization and analysis of contour maps
Steric analysis: The CoMFA contour maps have permitted an understanding of the steric and electrostatic requirements that represent the QSAR result. [Figure 7]a shows the contour map derived from the CoMFA PLS model. The contours are mapped on to the most active compound 71. The contour plots may help to identify important regions where any change may affect the activity of molecule. Furthermore, they are helpful in identifying important features contributing to increased activity of molecules. The steric interactions are represented by green, and region where steric bulk decreases the activity is depicted in yellow color. The yellow contour map shown extending from the tri-fluoro carbon to the left lower region indicates that the longer chain substituents toward this spatial distribution decrease the activity. The CoMFA steric contour map was less informative than the electrostatic map as the green contours were missing.
Electrostatic analysis: The electrostatic contour maps obtained from CoMFA analyses were mapped on to the compound 71 and are shown in the [Figure 7]b. The presence of red contour map at the C13 position of the methyl benzamide scaffold indicates that the substituents having electronegative group attached at this position show improved activity as observed in the compounds 57-72 against the unsubstituted compounds 12-34. The compound 71 is showing highest activity as its electropositive substituent is mapped to the blue contour appeared at this position. Substitution of electronegative group in this position has detrimental effect on activity as observed in compounds 39, 40, and 45.
Steric contour analysis: [Figure 8]a shows sterically favored (green) and disfavored (yellow) regions. The yellow contour map shown toward the left lower region of this position indicates that the longer chain substituents toward this spatial distribution decrease the activity. The contour at this position is also observed in CoMFA map [Figure 7]. The yellow contour around the pyrrolidin ring signifies that the steric substituents are disfavored (in case of compound 34). Yellow blocks appeared at fourth and fifth positions of phenyl ring indicate that the bulky groups with longer chain substituents decrease the activity.
Electrostatic contour analysis: The electrostatic contour maps shown by the CoMSIA model [Figure 8]b are more informative than that of CoMFA model [Figure 7]b, as they give more detailed picture in the case of substitutions. The blue contour map positioned near CF3 of benzene ring illustrates that electropositive groups are favored at this position.
As depicted from the [Figure 8]b, blue contour map of CoMSIA model shown at the methyl benzamide near the red contour map denotes that electropositive groups are favored at this region. At this position, the longer carbon linker chains show increased activity as noted in compounds 31, 21, and 61. The red contour map observed at methyl group near the blue contour specifies the region for electronegative groups. The electronegative group substituted compounds 58, 59, and 72 are observed to have comparatively higher activity. It was also observed in other compounds having phenyl ring with substituted halogens showing higher activity (56 and 57). The compounds 12, 17, and 15 with electropositive substituted phenyl ring at this position are showing decreased activity.
Hydrogen bond acceptor and donor contour analysis: [Figure 8]c and [Figure 8]d depict the hydrogen bond acceptor and donor contour maps of the CoMSIA models. Magenta contours encompass regions where a hydrogen bond acceptor will lead to improved biological activity, while an acceptor located near the red regions will result in impaired biological activity. In donor contour map, 6d, cyan color indicates the regions where hydrogen bond donor acts as favored and orange color refers to the disfavored regions. There are two magenta contour maps in [Figure 8]c, surrounding the two benzene rings, which supports the requirement of H-bond acceptor. This can be seen in compounds of table 3, 31, and 21. The small red contours in the same figure indicate that any hydrogen bond acceptors are not favored in these areas. This offers an explanation for the worse biological activity of compounds 45, 39, 40, and 15 as their substituents having acceptor group at this position are mapped to this region. The cyan contours in [Figure 8]d indicate that presence of H-bond donors increases the biological activity. The presence of strong H-bond donor in compounds with first scaffold is responsible for their high activity.
Hydrophobic contour analysis: Yellow and white contours enclose regions favorable for hydrophobic and hydrophilic groups, respectively. The white contours in [Figure 8]e support the importance of hydrophilic substitutions. This is responsible for the higher biological activities of molecules (64, 68, 69, and 70).This hydrophilic interaction might be very important for binding affinity, since this feature was also observed in CoMFA and CoMSIA steric contour maps. A small yellow contour mapped very near to methyl group of benzene ring illustrates that hydrophobic group at this position is important for increased biological activity.
| Conclusions|| |
In this study, we have investigated the CoMFA and CoMSIA models based on a training set of 37 structurally diverse (R)-3-aminopyrrolidine series, followed by validation of the results by an external test set of 13 analogues. Despite the lack of structural information on ccr0 2b, the design of potent inhibitors can be attempted by means of CoMFA and comparative molecular similarity indices, well-established 3D-QSAR techniques. These models demonstrated excellent internal and external predictive ability, which was shown by several strategies including cross-validation, predictive r2 , and test set predictions.
Overall, the CoMFA model gave good results. The CoMSIA model was more valuable for the three fields that contributed significantly (hydrophobic, hydrogen bond acceptor, and hydrogen bond donor). The CoMSIA analysis indicated that variations in the activity are dominated by hydrophobic interactions. The excellent correlation with several experimental studies suggests that these 3D-QSAR models are reliable, helping us to understand the binding interaction of these compounds and providing a helpful guideline for further lead optimization. The features derived from the above models bear a close correlation with the structural variations inherent in the training set, so other structurally distinct data results in diverse features causing different conclusions. In summary, our preliminary findings may aid in identifying potent and specific (R)-3-aminopyrrolidine series that may be used as potent antagonists of the CCR2b and offer more significant insights into the overall pharmacology of this system.
| References|| |
|1.||Rollins BJ. Blood. Chemokines 1997;90:909-28. |
|2.||Saunders J, Tarby CM. Opportunities for novel therapeutic agents acting at chemokine receptors. Drug Discov Today 1999;4:80-92. |
|3.||Murphy PM, Baggiolini M, Charo IF, Hébert CA, Horuk R, Matsushima K, et al. International union of pharmacology. XXII. Nomenclature for chemokine receptors. Pharmacol Rev 2000;52:145-76. |
|4.||Horuk R. Chemokine receptors. Cytokine Growth Factor Rev 2001;12:313-35. |
|5.||Carter PH. Chemokine receptor antagonism as an approach to anti-inflammatory therapy: 'just right' or plain wrong? Curr Opin Chem Biol 2002;6:510-25. |
|6.||Schwarz MK, Wells TN. New therapeutics that modulate chemokine networks. Nat Rev Drug Discov 2002;1:347-58. |
|7.||Charo IF, Ransohoff RM. The many roles of chemokines and chemokine receptors in inflammation. N Engl J Med 2006;354:610-21. |
|8.||Rollins BJ, Morrison ED, Stiles CD. Cloning and expression of JE, a gene inducible by platelet-derived growth factor and whose product has cytokine-like properties. Proc Natl Acad Sci U S A 1988;85:3738-42. |
|9.||Charo IF, Myers SJ, Herman A, Franci C, Connolly AJ, Coughlin SR. Molecular cloning and functional expression of two monocyte chemoattractant protein 1 receptors reveals alternative splicing of the carboxyl-terminal tails. Proc Natl Acad Sci U S A 1994;91:2752-6. |
|10.||Moree WJ, Kataoka K, Ramirez-Weinhouse MM, Shiota T, Imai M, Tsutsumi T, et al. Potent antagonists of the CCR2b receptor. Part 3: SAR of the (R)-3-aminopyrrolidine series. Bioorg Med Chem Lett 2008;18:1869-73. |
|11.||Oppenheim JJ, Zachariae CO, Mukaida N, Matsushima K. Properties of the novel proinflammatory supergene "intercrine" cytokine family. Annu Rev Immunol 1991;9:617-48. |
|12.||Rollins BJ. JE/MCP-1: an early-response gene encodes a monocyte-specific cytokine. Cancer Cells 1991;3:517-24. |
|13.||Schall TJ. Biology of the RANTES/SIS cytokine family. Cytokine 1991;3:165-83. |
|14.||Murphy PM. The molecular biology of leukocyte chemoattractant receptors. Annu Rev Immunol 1994;12:593-633. |
|15.||Jiang Y, Beller DI, Frendl G, Graves DT. Monocyte chemoattractant protein-1 regulates adhesion molecule expression and cytokine production in human monocytes. J Immunol 1992;148:2423-8. |
|16.||Zachariae CO, Anderson AO, Thompson HL, Appella E, Mantovani A, Oppenheim JJ, et al. Properties of monocyte chemotactic and activating factor (MCAF) purified from a human fibrosarcoma cell line. J Exp Med 1990;171:2177-82. |
|17.||Kunkel SL, Lukacs N, Kasama T, Strieter RM. The role of chemokines in inflammatory joint disease. J Leukoc Biol 1996;59:6-12. |
|18.||Akahoshi T, Wada C, Endo H, Hirota K, Hosaka S, Takagishi K, et al. Expression of monocyte chemotactic and activating factor in rheumatoid arthritis. Regulation of its production in synovial cells by interleukin-1 and tumor necrosis factor. Arthritis Rheum 1993;36:762-71. |
|19.||Koch AE, Kunkel SL, Harlow LA, Johnson B, Evanoff HL, Haines GK, et al. Enhanced production of monocyte chemoattractant protein-1 in rheumatoid arthritis. J Clin Invest 1992;90:772-9. |
|20.||Valente AJ, Graves DT, Vialle-Valentin CE, Delgado R, Schwartz CJ. Purification of a monocyte chemotactic factor secreted by nonhuman primate vascular cells in culture. Biochemistry 1988;27:4162-8. |
|21.||Navab M, Imes SS, Hama SY, Hough GP, Ross LA, Bork RW, et al. Monocyte transmigration induced by modification of low density lipoprotein in cocultures of human aortic wall cells is due to induction of monocyte chemotactic protein 1 synthesis and is abolished by high density lipoprotein. J Clin Invest 1991;88:2039-46. |
|22.||Liao F, Berliner JA, Mehrabian M, Navab M, Demer LL, Lusis AJ, et al. Minimally modified low density lipoprotein is biologically active in vivo in mice. J Clin Invest 1991;87:2253-7. |
|23.||Taubman MB, Rollins BJ, Poon M, Marmur J, Green RS, Berk BC, et al. JE mRNA accumulates rapidly in aortic injury and in platelet-derived growth factor-stimulated vascular smooth muscle cells. Circ Res 1992;70:314-25. |
|24.||Kurashima K, Mukaida N, Fujimura M, Schröder JM, Matsuda T, Matsushima K. Increase of chemokine levels in sputum precedes exacerbation of acute asthma attacks. J Leukoc Biol 1996;59:313-6. |
|25.||Hsieh KH, Chou CC, Chiang BL. Immunotherapy suppresses the production of monocyte chemotactic and activating factor and augments the production of IL-8 in children with asthma. J Allergy Clin Immunol 1996;98:580-7. |
|26.||Sousa AR, Lane SJ, Nakhosteen JA, Yoshimura T, Lee TH, Poston RN. Increased expression of the monocyte chemoattractant protein-1 in bronchial tissue from asthmatic subjects. Am J Respir Cell Mol Biol 1994;10:142-7. |
|27.||Gillitzer R, Wolff K, Tong D, Müller C, Yoshimura T, Hartmann AA, et al. MCP-1 mRNA expression in basal keratinocytes of psoriatic lesions. J Invest Dermatol 1993;101:127-31. |
|28.||Kristensen MS, Deleuran BW, Larsen CG, Thestrup-Pedersen K, Paludan K. Expression of monocyte chemotactic and activating factor (MCAF) in skin related cells. A comparative study. Cytokine 1993;5:520-4. |
|29.||Grandaliano G, Gesualdo L, Ranieri E, Monno R, Stallone G, Schena FP. Monocyte chemotactic peptide-1 expression and monocyte infiltration in acute renal transplant rejection. Transplantation 1997;63:414-20. |
|30.||Nadeau KC, Azuma H, Tilney NL. Sequential cytokine dynamics in chronic rejection of rat renal allografts: Roles for cytokines RANTES and MCP-1. Proc Natl Acad Sci U S A 1995;92:8729-33. |
|31.||Russell ME, Adams DH, Wyner LR, Yamashita Y, Halnon NJ, Karnovsky MJ. Early and persistent induction of monocyte chemoattractant protein 1 in rat cardiac allografts. Proc Natl Acad Sci U S A 1993;90:6086-90. |
|32.||Russell ME, Wallace AF, Hancock WW, Sayegh MH, Adams DH, Sibinga NE, et al. Upregulation of cytokines associated with macrophage activation in the Lewis-to-F344 rat transplantation model of chronic cardiac rejection. Transplantation 1995;59:572-8. |
|33.||Villiger PM, Terkeltaub R, Lotz M. Production of monocyte chemoattractant protein-1 by inflamed synovial tissue and cultured synoviocytes. J Immunol 1992;149:722-7. |
|34.||Cushing SD, Berliner JA, Valente AJ, Territo MC, Navab M, Parhami F, et al. Minimally modified low density lipoprotein induces monocyte chemotactic protein 1 in human endothelial cells and smooth muscle cells. Proc Natl Acad Sci U S A 1990;87:5134-8. |
|35.||Alam R, York J, Boyars M, Stafford S, Grant JA, Lee J, et al. Increased MCP-1, RANTES, and MIP-1alpha in bronchoalveolar lavage fluid of allergic asthmatic patients. Am J Respir Crit Care Med 1996;153:1398-404. |
|36.||Ogata H, Takeya M, Yoshimura T, Takagi K, Takahashi K. The role of monocyte chemoattractant protein-1 (MCP-1) in the pathogenesis of collagen-induced arthritis in rats. J Pathol 1997;182:106-14. |
|37.||Jones ML, Mulligan MS, Flory CM, Ward PA, Warren JS. Potential role of monocyte chemoattractant protein 1/JE in monocyte/macrophage-dependent IgA immune complex alveolitis in the rat. J Immunol 1992;149:2147-54. |
|38.||Boring L, Gosling J, Cleary M, Charo IF. Decreased lesion formation in CCR2-/- mice reveals a role for chemokines in the initiation of atherosclerosis. Nature 1998;394:894-7. |
|39.||Gosling J, Slaymaker S, Gu L, Tseng S, Zlot CH, Young SG, et al. MCP-1 deficiency reduces susceptibility to atherosclerosis in mice that overexpress human apolipoprotein B. J Clin Invest 1999;103:773-8. |
|40.||Forbes IT, Cooper DG, Dodds EK, Hickey DM, Ife RJ, Meeson M, et al. CCR2B receptor antagonists: Conversion of a weak HTS hit to a potent lead compound. Bioorg Med Chem Lett 2000;10:1803-6. |
|41.||Witherington J, Bordas V, Cooper DG, Forbes IT, Gribble AD, Ife RJ, et al. Conformationally restricted indolopiperidine derivatives as potent CCR2B receptor antagonists. Bioorg Med Chem Lett 2001;11:2177-80. |
|42.||Yang L[US], Butora G[US], Parsons WH[US], Pasternak A[US], inventors; MERCK & CO. [US], assignee; Cyclopentyl modulators of chemokine receptor activity. International Patent Application WO 0213824 (A1). 2002 feb 21. |
|43.||Mirzadegan T, Diehl F, Ebi B, Bhakta S, Polsky I, McCarley D, et al. Identification of the binding site for a novel class of CCR2b chemokine receptor antagonists: Binding to a common chemokine receptor motif within the helical bundle. J Biol Chem 2000;275:25562-71. |
|44.||Horino H, Ishiyama T, Tatematsu T, Makino C, Ochiai Y, Kito K, et al., inventors; Daiichi Pharmaceutical Co., Ltd, assignee; Benzoxazine derivative and medicament containing the same. Japanese patent publication JP 2001-278886 (P2001-278886A). 2001 oct 10.. |
|45.||Shiota T (Japan), Teijin Limited (Japan), inventors; Teijin Limited, COMBICHEM, INC. (US), assignee; CYCLIC AMINE DERIVATIVES AND THEIR USE AS DRUGS. Internatinal patent application WO 99/25686. 1999 may 27.. |
|46.||Cramer RD, Patterson DE, Jeffrey D. Comparative Molecular Field Analysis (CoMFA). 1. Effect of Shape on Binding of Steroids to Carrier Proteins. J Am Chem Soc 1988;110:5959-67. |
|47.||Klebe G, Abraham U, Mietzner T. Molecular similarity indices in a comparative analysis (CoMSIA) of drug molecules to correlate and predict their biological activity. J Med Chem 1994;37:4130-46. |
|48.||Böhm M, St rzebecher J, Klebe G. Three-dimensional quantitative structure-activity relationship analyses using comparative molecular field analysis and comparative molecular similarity indices analysis to elucidate selectivity differences of inhibitors binding to trypsin, thrombin, and factor Xa. J Med Chem 1999;42:458-77. |
|49.||Bursi R, Sawa M, Hiramatsu Y, Kondo H. A three-dimensional quantitative structure-activity relationship study of heparin-binding epidermal growth factor shedding inhibitors using comparative molecular field analysis. J Med Chem 2002;45:781-8. |
|50.||Sun WS, Park YS, Yoo J, Park KD, Kim SH, Kim JH, et al. Rational design of an indolebutanoic acid derivative as a novel aldose reductase inhibitor based on docking and 3D QSAR studies of phenethylamine derivatives. J Med Chem 2003;46:5619-27. |
|51.||Kuo CL, Assefa H, Kamath S, Brzozowski Z, Slawinski J, Saczewski F, et al. Application of CoMFA and CoMSIA 3D-QSAR and docking studies in optimization of mercaptobenzenesulfonamides as HIV-1 integrase inhibitors. J Med Chem 2004;47:385-99. |
|52.||The program SYBYL 6.7 is available from Tripos Inc.-1699 South Hanley Rd., St.Louis, MO 63144-2913, USA. . |
|53.||Matthew C, Cramer RD, Opdenbosch NV. Validation of the general purpose tripos 5.2 force field. J Comput Chem 1989;10:982-1012. |
|54.||Gasteiger J, Marsili M. Iterative partial equalization of orbital electronegativity-a rapid access to atomic charges. Tetrahedron 1980;36:3219-28. |
|55.||Dewar M J S, Zoebisch E G,, Healy EF, Stewart JJ. Development and use of quantum mechanical molecular models. 76. AM1: A new general purpose quantum mechanical molecular model. J Am Chem Soc 1985;107:3902-9. |
|56.||Zou XJ, Lai LH, Jin GY, Zhang ZX. Synthesis, fungicidal activity, and 3D-QSAR of pyridazinone-substituted 1,3,4-oxadiazoles and 1,3,4-thiadiazoles. J Agric Food Chem 2002;50:3757-60. |
|57.||Rutan, S. C. Kalman filtering approaches for solving problems in analytical chemistry. J Chemom 1987;1:7-18. |
|58.||Geladi P. Notes on the history and nature of partial least squares (PLS) modeling. J Chemom 1988;2:231-46. |
|59.||Wold S. Cross-validatory estimation of the number of components in factor and principal components models. Technometrics 1978;20:397-405. |
|60.||Böhm M, St rzebecher J, Klebe G. Three-dimensional quantitative structure-activity relationship analyses using comparative molecular field analysis and comparative molecular similarity indices analysis to elucidate selectivity differences of inhibitors binding to trypsin, thrombin, and factor Xa. J Med Chem 1999;42:458-77. |
[Figure 1], [Figure 2], [Figure 3], [Figure 4], [Figure 5], [Figure 6], [Figure 7], [Figure 8]
[Table 1], [Table 2], [Table 3], [Table 4]