Three-Dimensional-QSAR and Relative Binding Affinity Estimation of Focal Adhesion Kinase Inhibitors.

Precise binding affinity predictions are essential for structure-based drug discovery (SBDD). Focal adhesion kinase (FAK) is a member of the tyrosine kinase protein family and is overexpressed in a variety of human malignancies. Inhibition of FAK using small molecules is a promising therapeutic option for several types of cancer. Here, we conducted computational modeling of FAK-targeting inhibitors using three-dimensional structure-activity relationship (3D-QSAR), molecular dynamics (MD), and hybrid topology-based free energy perturbation (FEP) methods. The structure-activity relationship (SAR) studies between the physicochemical descriptors and inhibitory activities of the chemical compounds were performed with reasonable statistical accuracy using CoMFA and CoMSIA. These are two well-known 3D-QSAR methods based on the principle of supervised machine learning (ML). Essential information regarding residue-specific binding interactions was determined using MD and MM-PB/GBSA methods. Finally, physics-based relative binding free energy (ΔΔGRBFEA→B) terms of analogous ligands were estimated using alchemical FEP simulation. An acceptable agreement was observed between the experimental and computed relative binding free energies. Overall, the results suggested that using ML and physics-based hybrid approaches could be useful in synergy for the rational optimization of accessible lead compounds with similar scaffolds targeting the FAK receptor.

Binding Studies and Lead Generation of Pteridin-7(8H)-one Derivatives Targeting FLT3.

Ligand modification by substituting chemical groups within the binding pocket is a popular strategy for kinase drug development. In this study, a series of pteridin-7(8H)-one derivatives targeting wild-type FMS-like tyrosine kinase-3 (FLT3) and its D835Y mutant (FL3D835Y) were studied using a combination of molecular modeling techniques, such as docking, molecular dynamics (MD), binding energy calculation, and three-dimensional quantitative structure-activity relationship (3D-QSAR) studies. We determined the protein-ligand binding affinity by employing molecular mechanics Poisson-Boltzmann/generalized Born surface area (MM-PB/GBSA), fast pulling ligand (FPL) simulation, linear interaction energy (LIE), umbrella sampling (US), and free energy perturbation (FEP) scoring functions. The structure-activity relationship (SAR) study was conducted using comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA), and the results were emphasized as a SAR scheme. In both the CoMFA and CoMSIA models, satisfactory correlation statistics were obtained between the observed and predicted inhibitory activity. The MD and SAR models were co-utilized to design several new compounds, and their inhibitory activities were anticipated using the CoMSIA model. The designed compounds with higher predicted pIC50 values than the most active compound were carried out for binding free energy evaluation to wild-type and mutant receptors using MM-PB/GBSA, LIE, and FEP methods.

Structural Insights from Molecular Modeling of Isoindolin-1-One Derivatives as PI3Kγ Inhibitors against Gastric Carcinoma.

The upregulation of phosphoinositol-3-kinase γ (PI3Kγ) is deemed to be positively correlated with tumor-associated-macrophage (TAM)-mediated gastric carcinoma (GC). PI3Kγ suppresses tumor necrosis factor-alpha (TNF-α) and interleukin-12 (IL-12) through activation of the AKT/mTOR pathway, which promotes the immunosuppressant phenotype of TAM. Unlike α and ß isoforms, δ and γ isoforms are primarily distributed in leucocytes and macrophages. Dual inhibitors against PI3Kδ and PI3Kγ have been proven to have merits in targeting solid tumors. Furthermore, it has been found that PI3Kδ is activated by cytokines, while PI3Kγ is activated by G-protein-coupled receptors (GPCRs). This facilitates determining the functional difference between these two isoforms. For this goal, selective inhibitors would be immensely helpful. In the current manuscript, we conducted various molecular modeling studies with a series of isoindolin-1-one derivatives as potent PI3Kγ inhibitors by combining molecular docking, molecular dynamics (MD), molecular mechanics, Poisson-Boltzmann/generalized Born surface area (MM-PB/GBSA) binding free energy calculation, and three-dimensional structure-activity relationship (3D-QSAR) study. To evaluate the selectivity of γ isoform over δ, the molecular modeling studies of idelalisib analogs reported as PI3Kδ inhibitors were also investigated. The contour polyhedrons were generated from the comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA) around the ligand-bound active site for both isoforms, which could emphasize plausible explanations for the physicochemical factors that affect selective ligand recognition. The binding modalities of the two isoforms using CoMFA and MD models were compared, which suggested some key differences in the molecular interactions with the ligands and could be summarized as three subsites (one affinity subsite near the C-helix and DFG and two hydrophobic subsites). In the context of the structure-activity relationship (SAR), several new compounds were designed using a fragment-substitution strategy with the aim of selectively targeting PI3Kγ. The pIC50 values of the designed compounds were predicted by the 3D-QSAR models, followed by the MM-PB/GBSA binding energy estimation. The overall findings suggest that the designed compounds have the potential to be used as PI3Kγ inhibitors with a higher binding affinity and selectivity.

Computer aided designing of novel pyrrolopyridine derivatives as JAK1 inhibitors.

Janus kinases (JAKs) are a family of non-receptor kinases that play a key role in cytokine signaling and their aberrant activities are associated with the pathogenesis of various immune diseases. The JAK1 isoform plays an essential role in the types 1 and II interferon signaling and elicits signals from the interleukin-2, interleukin-4, gp130, and class 2 receptor families. It is ubiquitously expressed in humans and its overexpression has been linked with autoimmune diseases such as myeloproliferative neoplasm. Although JAK1 inhibitors such as Tofacitinib have been approved for medical use, the low potency and off-target effects of these inhibitors have limited their use and calls for the development of novel JAK1 inhibitors. In this study, we used computational methods on a series of pyrrolopyridine derivatives to design new JAK1 inhibitors. Molecular docking and molecular dynamics simulation methods were used to study the protein-inhibitor interactions. 3D-quantitative structure-activity relationship models were developed and were used to predict the activity of newly designed compounds. Free energy calculation methods were used to study the binding affinity of the inhibitors with JAK1. Of the designed compounds, seventeen of the compounds showed a higher binding energy value than the most active compound in the dataset and at least six of the compounds showed higher binding energy value than the pan JAK inhibitor Tofacitinib. The findings made in this study could be utilized for the further development of JAK1 inhibitors.

Molecular Modeling Studies of N-phenylpyrimidine-4-amine Derivatives for Inhibiting FMS-like Tyrosine Kinase-3.

Overexpression and frequent mutations in FMS-like tyrosine kinase-3 (FLT3) are considered risk factors for severe acute myeloid leukemia (AML). Hyperactive FLT3 induces premature activation of multiple intracellular signaling pathways, resulting in cell proliferation and anti-apoptosis. We conducted the computational modeling studies of 40 pyrimidine-4,6-diamine-based compounds by integrating docking, molecular dynamics, and three-dimensional structure-activity relationship (3D-QSAR). Molecular docking showed that K644, C694, F691, E692, N701, D829, and F830 are critical residues for the binding of ligands at the hydrophobic active site. Molecular dynamics (MD), together with Molecular Mechanics Poison-Boltzmann/Generalized Born Surface Area, i.e., MM-PB(GB)SA, and linear interaction energy (LIE) estimation, provided critical information on the stability and binding affinity of the selected docked compounds. The MD study suggested that the mutation in the gatekeeper residue F691 exhibited a lower binding affinity to the ligand. Although, the mutation in D835 in the activation loop did not exhibit any significant change in the binding energy to the most active compound. We developed the ligand-based comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA) models. CoMFA (q2 = 0.802, r2 = 0.983, and QF32 = 0.698) and CoMSIA (q2 = 0.725, r2 = 0.965 and QF32 = 0.668) established the structure-activity relationship (SAR) and showed a reasonable external predictive power. The contour maps from the CoMFA and CoMSIA models could explain valuable information about the favorable and unfavorable positions for chemical group substitution, which can increase or decrease the inhibitory activity of the compounds. In addition, we designed 30 novel compounds, and their predicted pIC50 values were assessed with the CoMSIA model, followed by the assessment of their physicochemical properties, bioavailability, and free energy calculation. The overall outcome could provide valuable information for designing and synthesizing more potent FLT3 inhibitors.

Designing of the N-ethyl-4-(pyridin-4-yl)benzamide based potent ROCK1 inhibitors using docking, molecular dynamics, and 3D-QSAR.

Rho-associated kinase-1 (ROCK1) has been recognized for its pivotal role in heart diseases, different types of malignancy, and many neurological disorders. Hyperactivity of ROCK phosphorylates the protein kinase-C (PKC), which ultimately induces smooth muscle cell contraction in the vascular system. Inhibition of ROCK1 has been shown to be a promising therapy for patients with cardiovascular disease. In this study, we have conducted molecular modeling techniques such as docking, molecular dynamics (MD), and 3-Dimensional structure-activity relationship (3D-QSAR) on a series of N-ethyl-4-(pyridin-4-yl)benzamide-based compounds. Docking and MD showed critical interactions and binding affinities between ROCK1 and its inhibitors. To establish the structure-activity relationship (SAR) of the compounds, 3D-QSAR techniques such as Comparative Molecular Field Analysis (CoMFA) and Comparative Molecular Similarity Indices Analysis (CoMSIA) were used. The CoMFA (q 2 = 0.774, r 2 = 0.965, ONC = 6, and r p r e d 2 = 0.703) and CoMSIA (q 2 = 0.676, r 2 = 0.949, ONC = 6, and r p r e d 2 = 0.548) both models have shown reasonable external predictive activity, and contour maps revealed favorable and unfavorable substitutions for chemical group modifications. Based on the contour maps, we have designed forty new compounds, among which, seven compounds exhibited higher predictive activity (pIC50). Further, we conducted the MD study, ADME/Tox, and SA score prediction using the seven newly designed compounds. The combination of docking, MD, and 3D-QSAR studies helps to understand the coherence modification of existing molecules. Our study may provide valuable insight into the development of more potent ROCK1 inhibitors.

Molecular Modelling Studies on Pyrazole Derivatives for the Design of Potent Rearranged during Transfection Kinase Inhibitors.

RET (rearranged during transfection) kinase, one of the receptor tyrosine kinases, plays a crucial role in the development of the human nervous system. It is also involved in various cell signaling networks responsible for the normal cell division, growth, migration, and survival. Previously reported clinical studies revealed that deregulation or aberrant activation of RET signaling can cause several types of human cancer. For example, medullary thyroid carcinoma (MTC) and multiple endocrine neoplasia (MEN2A, MEN2B) occur due to sporadic mutation or germline RET mutation. A number of RET kinase inhibitors have been approved by the FDA for the treatment of cancer, such as cabozantinib, vandetanib, lenvatinib, and sorafenib. However, each of these drugs is a multikinase inhibitor. Hence, RET is an important therapeutic target for cancer drug design. In this work, we have performed various molecular modelling studies, such as molecular docking and dynamics simulation for the most active compound of the pyrazole series as RET kinase inhibitors. Furthermore, molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) free energy calculation and 3-dimensional quantitative structure-activity relationship (3D-QSAR) were performed using g_mmpbsa and SYBYL-X 2.1 package. The results of this study revealed the crucial binding site residues at the active site of RET kinase and contour map analysis showed important structural characteristics for the design of new highly active inhibitors. Therefore, we have designed ten RET kinase inhibitors, which showed higher inhibitory activity than the most active compound of the series. The results of our study provide insights to design more potent and selective RET kinase inhibitors.

Molecular modeling studies of pyrrolo[2,3-d]pyrimidin-4-amine derivatives as JAK1 inhibitors based on 3D-QSAR, molecular docking, molecular dynamics (MD) and MM-PBSA calculations.

Rheumatoid Arthritis (RA) is an autoimmune disease caused by overproduction of pro-inflammatory cytokines. Janus Kinases (JAKs) mediate cytokines signaling through the Janus Kinase (JAK)/signal transducer and activator of transcription (STAT) signaling pathways. Clinical studies have shown that Janus kinase 1 (JAK1) mediated signaling plays a key role in synovial response in rheumatoid arthritis. Hence, the inhibition JAK1 is considered as an important therapeutic route for treatment of rheumatoid arthritis. In this study, we have performed three-dimensional quantitative structure-activity relationship (3 D-QSAR), molecular docking, molecular dynamics (MD) and free energy calculations on a series of pyrrolo[2,3-d]pyrimidin-4-amine JAK1 inhibitors. Molecular docking studies of the compounds 03, 13, 36 and 49 with JAK1 were performed to study the binding interactions. The binding conformations of the compounds from docking studies were selected based on binding energy and H-bond interactions and were used as initial structure for MD simulations. Using 3 D-QSAR techniques, a ligand-based comparative molecular field analysis (CoMFA) model (q2 = 0.5, r2 = 0.96) and a receptor-based CoMFA model (q2 = 0.78, r2 = 0.98) were developed. Analysis of the MD results of the most active compound (compound 49) with JAK1 showed the formation of H-bond interactions with residues Glu957, Leu959 and Gly887 and water-mediated H-bond interaction with Gly887 and His885. Based on the contour map analyses of the receptor-based CoMFA, a design strategy was proposed and was used for designing new JAK1 inhibitors. Four of the designed compounds (D57, D58, D98 and D99) showed predicted activity values (pIC50> 8.8) greater than the most active compound for JAK1. MM-PBSA based free energy calculations indicated that the designed compounds were able to form stable binding with JAK1 primarily through electrostatic interactions and van der Waal interactions. Collectively, the outcome of this study can be used to further the progress of JAK1 inhibition for the treatment of rheumatoid arthritis. Communicated by Ramaswamy H. Sarma.

Molecular Modeling Study of c-KIT/PDGFRα Dual Inhibitors for the Treatment of Gastrointestinal Stromal Tumors.

Gastrointestinal stromal tumors (GISTs) are the most common Mesenchymal Neoplasm of the gastrointestinal tract. The tumorigenesis of GISTs has been associated with the gain-of-function mutation and abnormal activation of the stem cell factor receptor (c-KIT) and platelet-derived growth factor receptor alpha (PDGFRα) kinases. Hence, inhibitors that target c-KIT and PDGFRα could be a therapeutic option for the treatment of GISTs. The available approved c-KIT/PDGFRα inhibitors possessed low efficacy with off-target effects, which necessitated the development of potent inhibitors. We performed computational studies of 48 pyrazolopyridine derivatives that showed inhibitory activity against c-KIT and PDGFRα to study the structural properties important for inhibition of both the kinases. The derivative of phenylurea, which has high activities for both c-KIT (pIC50 = 8.6) and PDGFRα (pIC50 = 8.1), was used as the representative compound for the dataset. Molecular docking and molecular dynamics simulation (100 ns) of compound 14 was performed. Compound 14 showed the formation of hydrogen bonding with Cys673, Glu640, and Asp810 in c-KIT, and Cys677, Glu644, and Asp836 in PDGFRα. The results also suggested that Thr670/T674 substitution in c-KIT/PDGFRα induced conformational changes at the binding site of the receptors. Three-dimensional quantitative structure-activity relationship (3D-QSAR) models were developed based on the inhibitors. Contour map analysis showed that electropositive and bulky substituents at the para-position and the meta-position of the benzyl ring of compound 14 was favorable and may increase the inhibitory activity against both c-KIT and PDGFRα. Analysis of the results suggested that having bulky and hydrophobic substituents that extend into the hydrophobic pocket of the binding site increases the activity for both c-KIT and PDGFRα. Based on the contour map analysis, 50 compounds were designed, and the activities were predicted. An evaluation of binding free energy showed that eight of the designed compounds have potential binding affinity with c-KIT/PDGFRα. Absorption, distribution, metabolism, excretion and toxicity (ADMET) and synthetic feasibility tests showed that the designed compounds have reasonable pharmaceutical properties and synthetic feasibility. Further experimental study of the designed compounds is recommended. The structural information from this study could provide useful insight into the future development of c-KIT and PDGFRα inhibitors.

Rational approach toward COVID-19 main protease inhibitors via molecular docking, molecular dynamics simulation and free energy calculation.

In the rapidly evolving coronavirus disease (COVID-19) pandemic, repurposing existing drugs and evaluating commercially available inhibitors against druggable targets of the virus could be an effective strategy to accelerate the drug discovery process. The 3C-Like proteinase (3CLpro) of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been identified as an important drug target due to its role in viral replication. The lack of a potent 3CLpro inhibitor and the availability of the X-ray crystal structure of 3CLpro (PDB-ID 6LU7) motivated us to perform computational studies to identify commercially available potential inhibitors. A combination of modeling studies was performed to identify potential 3CLpro inhibitors from the protease inhibitor database MEROPS ( https://www.ebi.ac.uk/merops/index.shtml ). Binding energy evaluation identified key residues for inhibitor design. We found 15 potential 3CLpro inhibitors with higher binding affinity than that of an α-ketoamide inhibitor determined via X-ray structure. Among them, saquinavir and three other investigational drugs aclarubicin, TMC-310911, and faldaprevir could be suggested as potential 3CLpro inhibitors. We recommend further experimental investigation of these compounds.

3D-QSAR-aided design of potent c-Met inhibitors using molecular dynamics simulation and binding free energy calculation.

Mesenchymal-epithelial transition factor (c-Met) is a member of receptor tyrosine kinase. It involves in various cellular signaling pathways which includes proliferation, motility, migration, and invasion. Over-expression of c-Met has been reported in various cancers. Hence, it is an ideal therapeutic target for cancer. The main objective of the study is to identify crucial residues involved in the inhibition of c-Met kinase and to design a series of potent imidazo [4,5-b] pyrazine derivatives as c-Met inhibitors. Docking was used to identify important active site residues involved in the inhibition of c-Met kinase which was further validated by 100 ns of molecular dynamics simulation and free energy calculation using molecular mechanics generalized born surface area. Furthermore, binding energy decomposition identified that residues Tyr1230, Met1211, Asp1222, Tyr1159, Met1160, Val1092, Ala1108, and Leu1157 contributed favorably to the binding stability of compound 32. Receptor-guided Comparative Molecular Field Analysis (CoMFA) (q2 = 0.751, NOC = 6, r2 = 0.933) and Comparative Molecular Similarity Indices Analysis (COMSIA) (q2 = 0.744, NOC = 6, r2 = 0.950) models were generated based on the docked conformation of the most active compound 32. The robustness of these models was tested using various validation techniques and found to be predictive. The results of CoMFA and CoMSIA contour maps exposed the regions favorable to enhance the activity. Based on this information, 27 novel c-Met inhibitors were designed. These designed compounds exhibited potent activity than the most active compound of the existing dataset. Communicated by Ramaswamy H. Sarma.

Receptor-guided 3D-QSAR studies, molecular dynamics simulation and free energy calculations of Btk kinase inhibitors.

BACKGROUND: Bruton tyrosine kinase (Btk) plays an important role in B-cell development, differentiation, and signaling. It is also found be in involved in male immunodeficiency disease such as X-linked agammaglobulinemia (XLA). Btk is considered as a potential therapeutic target for treating autoimmune diseases and hematological malignancies. RESULTS: In this work, a combined molecular modeling study was performed on a series of thieno [3,2-c] pyridine-4-amine derivatives as Btk inhibitors. Receptor-guided COMFA (q 2 = 0.574, NOC = 3, r 2 = 0.924) and COMSIA (q 2 = 0.646, NOC = 6, r 2 = 0.971) models were generated based on the docked conformation of the most active compound 26. All the developed models were tested for robustness using various validation techniques. Furthermore, a 5-ns molecular dynamics (MD) simulation and binding free energy calculations were carried out to determine the binding modes of the inhibitors and to identify crucial interacting residues. The rationality and stability of molecular docking and 3D-QSAR results were validated by MD simulation. The binding free energies calculated by the MM/PBSA method showed the importance of the van der Waals interaction. CONCLUSIONS: A good correlation between the MD results, docking studies, and the contour map analysis were observed. The study has identified the key amino acid residues in Btk binding pocket. The results from this study can provide some insights into the development of potent, novel Btk inhibitors.

Molecular modeling studies on series of Btk inhibitors using docking, structure-based 3D-QSAR and molecular dynamics simulation: a combined approach.

Bruton tyrosine kinase (Btk) is a non-receptor tyrosine kinase. It is a crucial component in BCR pathway and expressed only in hematopoietic cells except T cells and Natural killer cells. BTK is a promising target because of its involvement in signaling pathways and B cell diseases such as autoimmune disorders and lymphoma. In this work, a combined molecular modeling study of molecular docking, 3D-QSAR and molecular dynamic (MD) simulation were performed on a series of 2,5-diaminopyrimidine compounds as inhibitors targeting Btk kinase to understand the interaction and key residues involved in the inhibition. A structure based CoMFA (q (2) = 0.675, NOC = 5, r (2) = 0.961) and COMSIA (q (2) = 0.704, NOC = 6, r (2) = 0.962) models were developed from the conformation obtained by docking. The developed models were subjected to various validation techniques such as leave-five-out, external test set, bootstrapping, progressive sampling and rm (2) metrics and found to have a good predictive ability in both internal and external validation. Our docking results showed the important residues that interacts in the active site residues in inhibition of Btk kinase. Furthermore, molecular dynamics simulation was employed to study the stability of the docked conformation and to investigate the binding interactions in detail. The MD simulation analyses identified several important hydrogen bonds with Btk, including the gatekeeper residue Thr474 and Met477 at the hinge region. Hydrogen bond with active site residues Leu408 and Arg525 were also recognized. A good correlation between the MD results, docking studies and the contour map analysis are observed. This indicates that the developed models are reliable. Our results from this study can provide insights in the designing and development of more potent Btk kinase inhibitors.

High-temperature in situ crystallographic observation of reversible gas sorption in impermeable organic cages.

Crystallographic observation of adsorbed gas molecules is a highly difficult task due to their rapid motion. Here, we report the in situ single-crystal and synchrotron powder X-ray observations of reversible CO2 sorption processes in an apparently nonporous organic crystal under varying pressures at high temperatures. The host material is formed by hydrogen bond network between 1,3,5-tris-(4-carboxyphenyl)benzene (H3BTB) and N,N-dimethylformamide (DMF) and by π-π stacking between the H3BTB moieties. The material can be viewed as a well-ordered array of cages, which are tight packed with each other so that the cages are inaccessible from outside. Thus, the host is practically nonporous. Despite the absence of permanent pathways connecting the empty cages, they are permeable to CO2 at high temperatures due to thermally activated molecular gating, and the weakly confined CO2 molecules in the cages allow direct detection by in situ single-crystal X-ray diffraction at 323 K. Variable-temperature in situ synchrotron powder X-ray diffraction studies also show that the CO2 sorption is reversible and driven by temperature increase. Solid-state magic angle spinning NMR defines the interactions of CO2 with the organic framework and dynamic motion of CO2 in cages. The reversible sorption is attributed to the dynamic motion of the DMF molecules combined with the axial motions/angular fluctuations of CO2 (a series of transient opening/closing of compartments enabling CO2 molecule passage), as revealed from NMR and simulations. This temperature-driven transient molecular gating can store gaseous molecules in ordered arrays toward unique collective properties and release them for ready use.

3D-QSAR study of tetrahydro-3H-imidazo[4,5-c]pyridine derivatives as VEGFR-2 kinase inhibitors using various charge schemes.

Vascular endothelial growth factor-2 receptor (VEGFR-2) kinase is a promising target for the development of novel anticancer drugs. Three-dimensional quantitative structure-activity relationship (3D-QSAR) study was performed on a series of tetrahydro-3H-imidazo[4,5-c]pyridine derivatives to understand the structural basis for VEGFR-2 inhibitory activity. Several 3D-QSAR models were developed using various partial atomic charge schemes. Comparative molecular field analysis (CoMFA) and Comparative molecular similarity indices analysis (CoMSIA) methods were employed to derive these models. The CoMFA models performed better than the CoMSIA models. The reliable CoMFA model was obtained with the Gasteiger-Marsili charge scheme. The model produced statistically significant results with a cross-validated correlation coefficient (q(2)) of 0.635 and a coefficient of determination (r(2)) of 0.930. The model showed reasonable predictive power with predictive correlation coefficient ([Formula: see text]) of 0.582. Robustness of the model was further checked by leave-five-out cross-validation, bootstrapping and progressive scrambling analysis. The model was found to be statistically robust and expected to assist in the design of novel compounds with enhanced VEGFR-2 inhibitory activity.

In silico characterization of binding mode of CCR8 inhibitor: homology modeling, docking and membrane based MD simulation study.

Human CC-chemokine receptor 8 (CCR8) is a crucial drug target in asthma that belongs to G-protein-coupled receptor superfamily, which is characterized by seven transmembrane helices. To date, there is no X-ray crystal structure available for CCR8; this hampers active research on the target. Molecular basis of interaction mechanism of antagonist with CCR8 remains unclear. In order to provide binding site information and stable binding mode, we performed modeling, docking and molecular dynamics (MD) simulation of CCR8. Docking study of biaryl-ether-piperidine derivative (13C) was performed inside predefined CCR8 binding site to get the representative conformation of 13C. Further, MD simulations of receptor and complex (13C-CCR8) inside dipalmitoylphosphatidylcholine lipid bilayers were performed to explore the effect of lipids. Results analyses showed that the Gln91, Tyr94, Cys106, Val109, Tyr113, Cys183, Tyr184, Ser185, Lys195, Thr198, Asn199, Met202, Phe254, and Glu286 were conserved in both docking and MD simulations. This indicated possible role of these residues in CCR8 antagonism. However, experimental mutational studies on these identified residues could be effective to confirm their importance in CCR8 antagonism. Furthermore, calculated Coulombic interactions represented the crucial roles of Glu286, Lys195, and Tyr113 in CCR8 antagonism. Important residues identified in this study overlap with the previous non-peptide agonist (LMD-009) binding site. Though, the non-peptide agonist and currently studied inhibitor (13C) share common substructure, but they differ in their effects on CCR8. So, to get more insight into their agonist and antagonist effects, further side-by-side experimental studies on both agonist (LMD-009) and antagonist (13C) are suggested.

The nociceptin receptor (NOPR) and its interaction with clinically important agonist molecules: a membrane molecular dynamics simulation study.

The nociceptin receptor (NOPR) is an orphan G protein-coupled receptor that contains seven transmembrane helices. NOPR has a distinct mechanism of activation, though it shares a significant homology with other opioid receptors. Previously there have been reports on homology modeling of NOPR and also molecular dynamics simulation studies for a short period. Recently the crystal structure of NOPR was reported. In this study, we analyzed the time dependent behavior of NOPR docked with clinically important agonist molecules such as NOP (natural agonist) peptide and compound 10 (SCH-221510 derivative) using molecular dynamics simulations (MDS) for 100 ns. Molecular dynamics simulations of NOPR-agonist complexes allowed us to refine the system and to also identify stable structures with better binding modes. Structure activity relationships (SAR) for SCH221510 derivatives were investigated and reasons for the activities of these derivatives were determined. Our molecular dynamics trajectory analysis of NOPR-peptide and NOPR-compound 10 complexes found residues to be crucial for binding. Mutagenesis studies on the residues identified from our analysis could prove useful. Our results could also provide useful information in the structure-based drug design of novel and potent agonists targeting NOPR.

A Protoberberine derivative HWY336 selectively inhibits MKK4 and MKK7 in mammalian cells: the importance of activation loop on selectivity.

A protoberberine derivative library was used to search for selective inhibitors against kinases of the mitogen-activated protein kinase (MAPK) cascades in mammalian cells. Among kinases in mammalian MAPK pathways, we identified a compound (HWY336) that selectively inhibits kinase activity of mitogen-activated protein kinase kinase 4 and 7 (MKK4 and MKK7). The IC50 of HWY336 was 6 µM for MKK4 and 10 µM for MKK7 in vitro. HWY336 bound to both kinases reversibly via noncovalent interactions, and inhibited their activity by interfering with access of a protein substrate to its binding site. The binding affinity of HWY336 to MKK4 was measured by surface plasmon resonance to determine a dissociation constant (Kd) of 3.2 µM. When mammalian cells were treated with HWY336, MKK4 and MKK7 were selectively inhibited, resulting in inhibition of c-Jun NH2-terminal protein kinases in vivo. The structural model of HWY336 bound to either MKK4 or MKK7 predicted that HWY336 was docked to the activation loop, which is adjacent to the substrate binding site. This model suggested the importance of the activation loop of MKKs in HWY336 selectivity. We verified this model by mutating three critical residues within this loop of MKK4 to the corresponding residues in MKK3. The mutant MKK4 displayed similar kinase activity as wild-type kinase, but its activity was not inhibited by HWY336 compared to wild-type MKK4. We propose that the specific association of HWY336 to the activation loop of MKK4/MKK7 is responsible for its selective inhibition.

In silico study of 1-(4-Phenylpiperazin-1-yl)-2-(1H-pyrazol-1-yl) ethanones derivatives as CCR1 antagonist: homology modeling, docking and 3D-QSAR approach.

C-C chemokine receptor type 1 (CCR1) is a chemokine receptor with seven transmembrane helices and it belongs to the G-Protein Coupled receptor (GPCR) family. It plays an important role in rheumatoid arthritis, organ transplant rejection, Alzheimer's disease and also causes inflammation. Because of its role in disease processes, CCR1 is considered to be an important drug target. In the present study, we have performed three dimensional Quantitative Structure activity relationship (3D-QSAR) studies on a series of 1-(4-Phenylpiperazin-1-yl)-2-(1H-pyrazol-1-yl) ethanone derivatives targeting CCR1. Homology modeling of CCR1 was performed based on a template structure (4EA3) which has a high sequence identity and resolution. The highest active molecule was docked into this model. Ligand-based and Receptor-based quantitative structure-activity relationship (QSAR) study was performed and CoMFA models with reasonable statistics was developed for both ligand-based (q(2)=0.606; r(2)=0.968) and receptor-guided (q(2)=0.640; r(2)=0.932) alignment methods. Contour map analyses identified favorable regions for high affinity binding. The docking results highlighted the important active site residues. Tyr113 was found to interact with the ligand through hydrogen bonding. This residue has been considered responsible for anchoring ligands inside the active site. Our results could also be helpful to understand the inhibitory mechanism of 1-(4-Phenylpiperazin-1-yl)-2-(1H-pyrazol-1-yl) ethanone derivatives thereby to design more effective ligands in the future.

Novel ionophores with 2n-crown-n topology: anion sensing via pure aliphatic C-H···anion hydrogen bonding.

A series of novel coronands having a 2n-crown-n topology based on trioxane (6-crown-3) derivatives are designed and characterized. These neutral hosts can sense anions through pure aliphatic C-H hydrogen bonding (HB) in condensed phases due to the unusual topology of 2n-crown-n. C-H bonds are strongly polarized by two adjacent oxygen atoms in this scaffold. These hosts provide a rare opportunity to modulate anion binding strength by changing the electronic nature of aliphatic C-H bonds and offer ease of synthesis.