Molecular docking and 3D QSAR studies on 1-amino-2-phenyl-4-(piperidin-1-yl)-butanes based on the structural modeling of human CCR5 receptor

In the present study, we have used an approach combining protein structure modeling, molecular dynamics (MD) simulation, automated docking, and 3D QSAR analyses to investigate the detailed interactions of CCR5 with their antagonists. Homology modeling and MD simulation were used to build the 3D model of CCR5 receptor based on the high-resolution X-ray structure of bovine rhodopsin. A series of 64 CCR5 antagonists, 1-amino-2-phenyl-4-(piperidin-1-yl)-butanes, were docked into the putative binding site of the 3D model of CCR5 using the docking method, and the probable interaction model between CCR5 and the antagonists were obtained. The predicted binding affinities of the antagonists to CCR5 correlate well with the antagonist activities, and the interaction model could be used to explain many mutagenesis results. All these indicate that the 3D model of antagonist-CCR5 interaction is reliable. Based on the binding conformations and their alignment inside the binding pocket of CCR5, three-dimensional structure-activity relationship (3D QSAR) analyses were performed on these antagonists using comparative molecular field analysis (CoMFA) and comparative molecular similarity analysis (CoMSIA) methods. Both CoMFA and CoMSIA provide statistically valid models with good correlation and predictive power. The q(2)(r(cross)(2)) values are 0.568 and 0.587 for CoMFA and CoMSIA, respectively. The predictive ability of these models was validated by six compounds that were not included in the training set. Mapping these models back to the topology of the active site of CCR5 leads to a better understanding of antagonist-CCR5 interaction. These results suggest that the 3D model of CCR5 can be used in structure-based drug design and the 3D QSAR models provide clear guidelines and accurate activity predictions for novel antagonist design.     

Structure prediction and molecular dynamics simulations of a G-protein coupled receptor: human CCR2 receptor

CC chemokine receptor type-2 (CCR2) is a member of G-protein coupled receptors superfamily, expressed on the cell surface of monocytes and macrophages. It binds to the monocyte chemoattractant protein-1, a CC chemokine, produced at the sites of inflammation and infection. A homology model of human CCR2 receptor based on the recently available C-X-C chemokine recepor-4 crystal structure has been reported. Ligand information was used as an essential element in the homology modeling process. Six known CCR2 antagonists were docked into the model using simple and induced fit docking procedure. Docked complexes were then subjected to visual inspection to check their suitability to explain the experimental data obtained from site directed mutagenesis and structure-activity relationship studies. The homology model was refined, validated, and assessed for its performance in docking-based virtual screening on a set of CCR2 antagonists and decoys. The docked complexes of CCR2 with the known antagonists, TAK779, a dual CCR2/CCR5 antagonist, and Teijin-comp1, a CCR2 specific antagonist were subjected to molecular dynamics (MD) simulations, which further validated the binding modes of these antagonists. B-factor analysis of 20?ns MD simulations demonstrated that Cys190 is helpful in providing structural rigidity to the extracellular loop (EL2). Residues important for CCR2 antagonism were recognized using free energy decomposition studies. The acidic residue Glu291 from TM7, a conserved residue in chemokine receptors, is favorable for the binding of Teijin-comp1 with CCR2 by ΔG of ?11.4?kcal/mol. Its contribution arises more from the side chains than the backbone atoms. In addition, Tyr193 from EL2 contributes ?0.9?kcal/mol towards the binding of the CCR2 specific antagonist with the receptor. Here, the homology modeling and subsequent molecular modeling studies proved successful in probing the structure of human CCR2 chemokine receptor for the structure-based virtual screening and predicting the binding modes of CCR2 antagonists.     

Alkylsulfone-containing trisubstituted cyclohexanes as potent and bioavailable chemokine receptor 2 (CCR2) antagonists

We describe novel alkylsulfones as potent CCR2 antagonists with reduced hERG channel activity and improved pharmacokinetics over our previously described antagonists. Several of these new alkylsulfones have a profile that includes functional antagonism of CCR2, in vitro microsomal stability, and oral bioavailability. With this improved profile, we demonstrate that two of these antagonists, 2 and 12, are orally efficacious in an animal model of inflammatory recruitment.     

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

Monocyte chemoattracant-1 (MCP-1) stimulates leukocyte chemotaxis to inflammatory sites, such as rheumatoid arthritis, atherosclerosis, and asthma, by use of the MCP-1 receptor, CCR2, a member of the G-protein-coupled seven-transmembrane receptor superfamily. These studies identified a family of antagonists, spiropiperidines. One of the more potent compounds blocks MCP-1 binding to CCR2 with a K(d) of 60 nm, but it is unable to block binding to CXCR1, CCR1, or CCR3. These compounds were effective inhibitors of chemotaxis toward MCP-1 but were very poor inhibitors of CCR1-mediated chemotaxis. The compounds are effective blockers of MCP-1-driven inhibition of adenylate cyclase and MCP-1- and MCP-3-driven cytosolic calcium influx; the compounds are not agonists for these pathways. We showed that glutamate 291 (Glu(291)) of CCR2 is a critical residue for high affinity binding and that this residue contributes little to MCP-1 binding to CCR2. The basic nitrogen present in the spiropiperidine compounds may be the interaction partner for Glu(291), because the basicity of this nitrogen was essential for affinity; furthermore, a different class of antagonists, a class that does not have a basic nitrogen (2-carboxypyrroles), were not affected by mutations of Glu(291). In addition to the CCR2 receptor, spiropiperidine compounds have affinity for several biogenic amine receptors. Receptor models indicate that the acidic residue, Glu(291), from transmembrane-7 of CCR2 is in a position similar to the acidic residue contributed from transmembrane-3 of biogenic amine receptors, which may account for the shared affinity of spiropiperidines for these two receptor classes. The models suggest that the acid-base pair, Glu(291) to piperidine nitrogen, anchors the spiropiperidine compound within the transmembrane ovoid bundle. This binding site may overlap with the space required by MCP-1 during binding and signaling; thus the small molecule ligands act as antagonists. An acidic residue in transmembrane region 7 is found in most chemokine receptors and is rare in other serpentine receptors. The model of the binding site may suggest ways to make new small molecule chemokine receptor antagonists, and it may rationalize the design of more potent and selective antagonists.     

QSAR studies on CCR2 antagonists with chiral sensitive hologram descriptors

Chemokines are small molecular weight water-soluble proteins playing a key role in immunomodulation and host-defense mechanisms. CCR2 receptor is targeted for diseases like arthritis, multiple sclerosis, vascular disease, obesity, and type 2 diabetes. Reported, herein are the QSAR studies performed on a diverse set of enantiopure analogues reported as CCR2 antagonists by hologram analysis. The best model highlights the importance of chirality feature in comparison with the other models developed without the chirality. The validated model showed high internal and external predictive power. The robustness of the model was achieved with good statistical r(2) of 0.945 and cross-validated r(cv)(2) of 0.837. The challenging test predictivity of the model was confirmed with r(pred)(2) of 0.807. The fragment fingerprints help in understanding essential pharmacophoric features for CCR2 antagonism and provide basis for SAR of the molecules. The 2D contribution maps with fragment information will be useful for the design of novel CCR2 antagonists having improved efficacy.     

Ligand supported homology modeling and docking evaluation of CCR2: docked pose selection by consensus scoring

Chemokine receptor 2 (CCR2) is a G-protein coupled receptor (GPCR) and a crucial target for various inflammatory and autoimmune diseases. The structure based antagonists design for many GPCRs, including CCR2, is restricted by the lack of an experimental three dimensional structure. Homology modeling is widely used for the study of GPCR-ligand binding. Since there is substantial diversity for the ligand binding pocket and binding modes among GPCRs, the receptor-ligand binding mode predictions should be derived from homology modeling with supported ligand information. Thus, we modeled the binding of our proprietary CCR2 antagonist using ligand supported homology modeling followed by consensus scoring the docking evaluation based on all modeled binding sites. The protein-ligand model was then validated by visual inspection of receptor-ligand interaction for consistency of published site-directed mutagenesis data and virtual screening a decoy compound database. This model was able to successfully identify active compounds within the decoy database. Finally, additional hit compounds were identified through a docking-based virtual screening of a commercial database, followed by a biological assay to validate CCR2 inhibitory activity. Thus, this procedure can be employed to screen a large database of compounds to identify new CCR2 antagonists.     

Discovery and pharmacological characterization of a novel rodent-active CCR2 antagonist, INCB3344

This report describes the characterization of INCB3344, a novel, potent and selective small molecule antagonist of the mouse CCR2 receptor. The lack of rodent cross-reactivity inherent in the small molecule CCR2 antagonists discovered to date has precluded pharmacological studies of antagonists of this receptor and its therapeutic relevance. In vitro, INCB3344 inhibits the binding of CCL2 to mouse monocytes with nanomolar potency (IC(50) = 10 nM) and displays dose-dependent inhibition of CCL2-mediated functional responses such as ERK phosphorylation and chemotaxis with similar potency. Against a panel of G protein-coupled receptors that includes other CC chemokine receptors, INCB3344 is at least 100-fold selective for CCR2. INCB3344 possesses good oral bioavailability and systemic exposure in rodents that allows in vivo pharmacological studies. INCB3344 treatment results in a dose-dependent inhibition of macrophage influx in a mouse model of delayed-type hypersensitivity. The histopathological analysis of tissues from the delayed-type hypersensitivity model demonstrates that inhibition of CCR2 leads to a substantial reduction in tissue inflammation, suggesting that macrophages play an orchestrating role in immune-based inflammatory reactions. These results led to the investigation of INCB3344 in inflammatory disease models. We demonstrate that therapeutic dosing of INCB3344 significantly reduces disease in mice subjected to experimental autoimmune encephalomyelitis, a model of multiple sclerosis, as well as a rat model of inflammatory arthritis. In summary, we present the first report on the pharmacological characterization of a selective, potent and rodent-active small molecule CCR2 antagonist. These data support targeting this receptor for the treatment of chronic inflammatory diseases.     

The human specific CCR1 antagonist CP-481,715 inhibits cell infiltration and inflammatory responses in human CCR1 transgenic mice

We previously described the in vitro characteristics of the potent and selective CCR1 antagonist, CP-481,715. In addition to being selective for CCR1 vs other chemokine receptors, CP-481,715 is also specific for human CCR1 (hCCR1), preventing its evaluation in classical animal models. To address this, we generated mice whereby murine CCR1 was replaced by hCCR1 (knockin) and used these animals to assess the anti-inflammatory properties of CP-481,715. Cells isolated from hCCR1 knockin mice were shown to express hCCR1 and migrate in response to both murine CCR1 and hCCR1 ligands. Furthermore, this migration is inhibited by CP-481,715 at dose levels comparable to those obtained with human cells. In animal models of cell infiltration, CP-481,715 inhibited CCL3-induced neutrophil infiltration into skin or into an air pouch with an ED50 of 0.2 mg/kg. CP-481,715 did not inhibit cell infiltration in wild-type animals expressing murine CCR1. In a more generalized model of inflammation, delayed-type hypersensitivity, CP-481,715 significantly inhibited footpad swelling and decreased the amount of IFN-gamma and IL-2 produced by isolated spleen cells from sensitized animals. It did not, however, induce tolerance to a subsequent challenge. These studies illustrate the utility of hCCR1 knockin animals to assess the activity of human specific CCR1 antagonists; demonstrate the ability of the CCR1 antagonist CP-481,715 to inhibit cell infiltration, inflammation, and Th1 cytokine responses in these animals; and suggest that CP-481,715 may be useful to modulate inflammatory responses in human disease.     

Overcoming hERG activity in the discovery of a series of 4-azetidinyl-1-aryl-cyclohexanes as CCR2 antagonists

A series of 4-azetidinyl-1-aryl-cyclohexanes as potent CCR2 antagonists with high selectivity over activity for the hERG potassium channel is discovered through divergent SARs of CCR2 and hERG.     

Molecular modeling and QSAR studies of a set of indole and benzimidazole derivatives as H<Subscript>4</Subscript> receptor antagonists

Histamine is an important biogenic amine, which acts with a group of four G-protein coupled receptors (GPCRs), namely H₁ to H₄ (H₁R – H₄R) receptors. The actions of histamine at H₄R are related to immunological and inflammatory processes, particularly in pathophysiology of asthma, and H₄R ligands having antagonistic properties could be helpful as antiinflammatory agents. In this work, molecular modeling and QSAR studies of a set of 30 compounds, indole and benzimidazole derivatives, as H₄R antagonists were performed. The QSAR models were built and optimized using a genetic algorithm function and partial least squares regression (WOLF 5.5 program). The best QSAR model constructed with training set (N = 25) presented the following statistical measures: r ² = 0.76, q ² = 0.62, LOF = 0.15, and LSE = 0.07, and was validated using the LNO and y-randomization techniques. Four of five compounds of test set were well predicted by the selected QSAR model, which presented an external prediction power of 80%. These findings can be quite useful to aid the designing of new anti-H₄ compounds with improved biological response.     

In vitro selection of RNA aptamers that block CCL1 chemokine function

CCL1, the CCR8 ligand, is a CC chemokine secreted by activated monocytes and lymphocytes and is a potent chemoattractant for these cell types. The in vivo role of the CCL1/CCR8 axis in Th2-mediated inflammation is far from clear. Ligand neutralisation studies reported discrepancies in the effect of CCL1/CCR8 and CCR8 knockout studies showed very different insights into the functional role of the CCR8. To further study the biological function of CCL1, we focused on the generation and characterisation of RNA aptamers. We report here the in vitro isolation of the first nuclease resistant and selective RNA aptamer (T48) with high-binding affinity for human and mouse CCL1. The T48 aptamer but not a random control aptamer antagonises CCL1 function in a dose-dependent fashion in both heparin binding and chemotaxis assays. To our knowledge, the T48 aptamer constitutes one of the most potent CCL1 antagonists reported to date and is an excellent tool to dissect CCL1-specific function in vivo. The T48 aptamer may also have potential as new generation of therapeutic tools.     

Functional expression and characterization of macaque C-C chemokine receptor 3 (CCR3) and generation of potent antagonistic anti-macaque CCR3 monoclonal antibodies

Eosinophils are major effector cells implicated in a number of chronic inflammatory diseases in humans, particularly bronchial asthma and allergic rhinitis. The beta-chemokine receptor C-C chemokine receptor 3 (CCR3) provides a mechanism for the selective recruitment of eosinophils into tissue and thus has recently become an attractive biological target for therapeutic intervention. In order to develop in vivo models of inflammatory diseases, it is essential to identify and characterize the homologues of human eotaxin (C-C chemokine ligand 11) and CCR3 from other species, such as non-human primates. Accordingly, we cloned the macaque eotaxin and CCR3 genes and revealed that they were 91 and 92% identical at the amino acid level to their human homologues, respectively. Macaque CCR3 expressed in the murine pre-B L1-2 cell line bound macaque eotaxin with high affinity (K(d) = 0.1 nm) and exhibited a robust eotaxin-induced Ca(2+) flux and chemotaxis. Characterization of beta-chemokines on native macaque CCR3 on eosinophils was performed by means of eotaxin-induced shape change in whole blood using a novel signaling assay known as gated autofluorescence forward scatter. Additionally, mAbs were raised against macaque CCR3 using two different immunogens: a 30-amino acid synthetic peptide derived from the predicted NH(2) terminus of macaque CCR3 and intact macaque CCR3-transfected cells. These anti-macaque CCR3 monoclonal antibodies exhibited potent antagonist activity in receptor binding and functional assays. The characterization of the macaque eotaxin/CCR3 axis and development of antagonistic anti-macaque CCR3 monoclonal antibodies will facilitate the development of CCR3 small molecule antagonists with the hope of ameliorating chronic inflammatory diseases in humans.     

Site-directed mutagenesis of CC chemokine receptor 1 reveals the mechanism of action of UCB 35625, a small molecule chemokine receptor antagonist

The chemokine receptor CCR1 and its principal ligand, CCL3/MIP-1alpha, have been implicated in the pathology of several inflammatory diseases including rheumatoid arthritis, multiple sclerosis, and asthma. As such, these molecules are the focus of much research with the ultimate aim of developing novel therapies. We have described previously a non-competitive small molecule antagonist of CCR1 (UCB 35625), which we hypothesized interacted with amino acids located within the receptor transmembrane (TM) helices (Sabroe, I., Peck, M. J., Jan Van Keulen, B., Jorritsma, A., Simmons, G., Clapham, P. R., Williams, T. J., and Pease, J. E. (2000) J. Biol. Chem. 275, 25985-25992). Here we describe an approach to identifying the mechanism by which the molecule antagonizes CCR1. Thirty-three point mutants of CCR1 were expressed transiently in L1.2 cells, and the cells were assessed for their capacity to migrate in response to CCL3 in the presence or absence of UCB 35625. Cells expressing the mutant constructs Y41A (TM helix 1, or TM1), Y113A (TM3), and E287A (TM7) were responsive to CCL3 but resistant to the antagonist, consistent with a role for the TM helices in CCR1 interactions with UCB 35625. Subsequent molecular modeling successfully docked the compound with CCR1 and suggests that the antagonist ligates TM1, 2, and 7 of CCR1 and severely impedes access to TM2 and TM3, a region thought to be perturbed by the chemokine amino terminus during the process of receptor activation. Insights into the mechanism of action of these compounds may facilitate the development of more potent antagonists that show promise as future therapeutic agents in the treatment of inflammatory disease.     

Ligand-based molecular design of 4-benzylpiperidinealkylureas and amides as CCR3 antagonists

Asthma is an inflammatory disease of the lungs. Clinical studies suggest that eotaxin and chemokine receptor-3 (CCR3) play a primary role in the recruitment of eosinophils in allergic asthma. Development of novel and potent CCR3 antagonists could provide a novel mechanism for inhibition of this recruitment process, thereby preventing asthma. With the intention of designing new ligands with enhanced inhibitor potencies against CCR3, a 3D-QSAR CoMFA study was carried out on 41 4-benzylpiperidinealkylureas and amide derivatives. The best statistics of the developed CoMFA model were r ² = 0.960, r_cv² = 0.589 r_{cv}^2 = 0.589 , n = 32 for the training set and r_pred² = 0.619 r_{pred}^2 = 0.619 , n = 9 for the test set. The generated 3D-QSAR contribution maps shed some light on the effects of the substitution pattern related to CCR3 antagonist activity.     

Three-dimensional quantitative structure-activity relationship analyses of piperidine-based CCR5 receptor antagonists

The CCR5 chemokine receptor has recently been found to play a crucial role in the viral entry stage of HIV infection and has therefore become an attractive potential target for anti-HIV therapeutics. On the other hand, the lack of CCR5 crystal structure data has impeded the development of structure-based CCR5 antagonist design. In this paper, we compare two three-dimensional Quantitative Structure-Activity Relationship (3D-QSAR) methods: Comparative Molecular Field Analysis (CoMFA) and Comparative Molecular Similarity Indices Analysis (CoMSIA) on a series of piperidine-based CCR5 antagonists as an alternative approach to investigate the interaction between CCR5 antagonists and their receptor. Superimposition of antagonist structures was performed using two alignment rules: atomic/centroid rms fit and rigid body field fit techniques. The 3D QSAR models were derived from a training set of 72 compounds, and were found to have predictive capability for a set of 19 holdout test compounds. The resulting contour maps produced by the best CoMFA and CoMSIA models were used to identify the structural features relevant to biological activity in this series of compounds. Further analyses of these interaction-field contour maps also showed a high level of internal consistency.     

Structural exploration of tetrahydroisoquinoline derivatives as HDAC8 inhibitors through multi-QSAR modeling study

Abstract

Histone deacetylase 8 (HDAC8) is one of the crucial HDACs responsible for influencing the epigenetic functions of the body. Overexpression of HDAC8 is found to be involved in numerous disease conditions such as tumorigenesis, cell proliferation, cancer, viral infections, neuronal disorders and other epigenetic diseases. Therefore, inhibition of HDAC8 is a primary method to combat these diseases. In this article, a multi-QSAR modeling study on tetrahydroisoquinoline derivatives was conducted to identify important contributions of the structural features of these compounds toward HDAC8 inhibition. All these QSAR modeling techniques were individually validated and justified the observations of each other. The results implied that the tetrahydroisoquinoline moiety may be effective as a cap group than as a linker moiety for HDAC8 inhibition. Different substitutions at the tetrahydroisoquinoline scaffold were also found to be crucial in modulating HDAC8 inhibition.

Communicated by Ramaswamy H. Sarma