Identification of Potent and Selective JAK1 Lead Compounds Through Ligand-Based Drug Design Approaches.

JAK1 plays a significant role in the intracellular signaling by interacting with cytokine receptors in different types of cells and is linked to the pathogenesis of various cancers and in the pathology of the immune system. In this study, ligand-based pharmacophore modeling combined with virtual screening and molecular docking methods was incorporated to identify the potent and selective lead compounds for JAK1. Initially, the ligand-based pharmacophore models were generated using a set of 52 JAK1 inhibitors named C-2 methyl/hydroxyethyl imidazopyrrolopyridines derivatives. Twenty-seven pharmacophore models with five and six pharmacophore features were generated and validated using potency and selectivity validation methods. During potency validation, the Guner-Henry score was calculated to check the accuracy of the generated models, whereas in selectivity validation, the pharmacophore models that are capable of identifying selective JAK1 inhibitors were evaluated. Based on the validation results, the best pharmacophore models ADHRRR, DDHRRR, DDRRR, DPRRR, DHRRR, ADRRR, DDHRR, and ADPRR were selected and taken for virtual screening against the Maybridge, Asinex, Chemdiv, Enamine, Lifechemicals, and Zinc database to identify the new molecules with novel scaffold that can bind to JAK1. A total of 4,265 hits were identified from screening and checked for acceptable drug-like properties. A total of 2,856 hits were selected after ADME predictions and taken for Glide molecular docking to assess the accurate binding modes of the lead candidates. Ninety molecules were shortlisted based on binding energy and H-bond interactions with the important residues of JAK1. The docking results were authenticated by calculating binding free energy for protein-ligand complexes using the MM-GBSA calculation and induced fit docking methods. Subsequently, the cross-docking approach was carried out to recognize the selective JAK1 lead compounds. Finally, top five lead compounds that were potent and selective against JAK1 were selected and validated using molecular dynamics simulation. Besides, the density functional theory study was also carried out for the selected leads. Through various computational studies, we observed good potency and selectivity of these lead compounds when compared with the drug ruxolitinib. Compounds such as T5923555 and T5923531 were found to be the best and can be further validated using in vitro and in vivo methods.

Structure and dynamics of the somatostatin receptor 3-ligand binding in the presence of lipids examined using computational structural biology methods.

In the past two decades, the structural biology studies on G-protein coupled receptors (GPCRs) are on the rise. Understanding the relation between the structure and function of GPCRs is important as they play a huge role in various signaling mechanisms in a eukaryotic cell. Somatostatin receptor 3 (SSTR3), one of the GPCRs, is one such important receptor which oversees different cellular processes including cell-to-cell signaling. However, the information available regarding the structural features of SSTR3 responsible for their bioactivity is scarce. In this study, we report a structural understanding of SSTR3-ligand binding that could be helpful in demystifying the structural complexities related to functioning of the receptor. An integrated protocol consisting of different computational structural biology tools including protein structure prediction via comparative modeling, binding site characterization, three-dimensional quantitative structure-activity relationship based on comparative molecular field analysis and comparative molecular similarity indices analysis, density functional theory, and molecular dynamics simulations were performed. Different understandings from the simulation of SSTR3-ligand complexes, mainly the conditions that are favorable for the formation of lowest bioactive state of SSTR3 ligands are reported. In addition to that, we report the important physicochemical descriptors of SSTR3 ligands that could significantly influence their bioactivity. The results of the study could be helpful in developing novel SSTR3 ligands (both agonists and antagonists) with high potency and receptor selectivity.

Understanding the influence of lipid bilayers and ligand molecules in determining the conformational dynamics of somatostatin receptor 2.

Somatostatin receptor 2 (SSTR2) is a G-protein coupled receptor (GPCR) that controls numerous cellular processes including cell-to-cell signaling. In this study, we report how the lipid and ligand molecules influence the conformational dynamics of the membrane-bound SSTR2. Molecular simulations of different holo and apoenzyme complexes of SSTR2 in the presence and absence of a lipid bilayer were performed, observed, and correlated with previously reported studies. We identified the important SSTR2 residues that take part in the formation of the SSTR2-ligand complex. On analyzing the molecular simulation trajectories, we identified that the residue D3.32 is crucial in determining the bioactive conformation of SSTR2 ligands in the binding site. Based on the results, we suggest that designing a novel SSTR2 ligand with an H-bond donor group at the R1 position, and hydrophobic groups at R2 and R3 might have higher activity and SSTR2-selectivity. We analyzed the simulated systems to identify other important structural features involved in SSTR2-ligand binding and to observe the different conformational changes that occur in the protein after the ligand binding. Additionally, we studied the conformational dynamics of N- and C-terminal regions of SSTR2 in the presence and absence of the lipid bilayer. Both the systems were compared to understand the influence of lipid molecules in the formation of secondary structural domains by these extracellular regions. The comparative study revealed that the secondary structural elements formed by C-terminal residues in presence of lipid molecules is crucial for the functioning of SSTR2. Our study results highlight the structural complexities involved in the functioning of SSTR upon binding with the ligands in the presence and absence of lipid bilayer, which is essential for designing novel drug targets.

Molecular-Level Understanding of the Somatostatin Receptor 1 (SSTR1)-Ligand Binding: A Structural Biology Study Based on Computational Methods.

Somatostatin receptor 1 (SSTR1), a subtype of somatostatin receptors, is involved in various signaling mechanisms in different parts of the human body. Like most of the G-protein-coupled receptors (GPCRs), the available information on the structural features of SSTR1 responsible for the biological activity is scarce. In this study, we report a molecular-level understanding of SSTR1-ligand binding, which could be helpful in solving the structural complexities involved in SSTR1 functioning. Based on a three-dimensional quantitative structure-activity relationship (3D-QSAR) study using comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA), we have identified that an electronegative, less-bulkier, and hydrophobic atom substitution can substantially increase the biological activity of SSTR1 ligands. A density functional theory (DFT) study has been followed to study the electron-related properties of the SSTR1 ligands and to validate the results obtained via the 3D-QSAR study. 3D models of SSTR1-ligand systems have been embedded in lipid-lipid bilayer membranes to perform molecular dynamics (MD) simulations. Analysis of the MD trajectories reveals important information about the crucial residues involved in SSTR1-ligand binding and various conformational changes in the protein that occur after ligand binding. Additionally, we have identified the probable ligand-binding site of SSTR1 and validated it using MD. We have also studied the favorable conditions that are essential for forming the most stable and lowest-energy bioactive conformation of the ligands inside the binding site. The results of the study could be useful in constructing more potent and novel SSTR1 antagonists and agonists.

Understanding the structural features of JAK2 inhibitors: a combined 3D-QSAR, DFT and molecular dynamics study.

JAK2 plays a critical role in JAK/STAT signaling pathway and in patho-mechanism of myeloproliferative disorders and autoimmune diseases. Thus, effective JAK2 inhibitors provide a promising opportunity for the pharmaceutical intervention of many diseases. In this work, 3D-QSAR study was performed on a series of 1-amino-5H-pyrido-indole-4-carboxamide derivatives as JAK2 inhibitors to obtain reliable comparative molecular field analysis (CoMFA) and comparative molecular similarity analysis (CoMSIA) models with three different alignment methods. Among the different alignment methods, ligand-based (CoMFA: q2 = 0.676, r2 = 0.979; CoMSIA: q2 = 0.700, r2 = 0.953) and pharmacophore-based alignment (CoMFA: q2 = 0.710, r2 = 0.982; CoMSIA: q2 = 0.686, r2 = 0.960) has produced better statistical results when compared to receptor-based alignment (CoMFA: q2 = 0.507, r2 = 0.979; CoMSIA: q2 = 0.544, r2 = 0.917). Statistical parameters indicated that data are well fitted and have high predictive ability. The presence of electrostatic and hydrophobic field is highly desirable for potent inhibitory activity, and the steric field plays a minor role in modulating the activity. The contour analysis indicates ARG980, ASN981, ASP939 and LEU937 have more possibility of interacting with bulky, hydrophobic groups in pyrido and positive and negative groups in pyrazole ring. Based on our findings, we have designed sixteen molecules and predicted its activity and drug-like properties. Subsequently, molecular docking, molecular dynamics and DFT calculations were performed to evaluate its potency.

Toward a better understanding of the interaction between somatostatin receptor 2 and its ligands: a structural characterization study using molecular dynamics and conceptual density functional theory.

This study is a part of the extensive research intending to provide the structural insights on somatostatin and its receptor. Herein, we have studied the structural complexity involved in the binding of somatostatin receptor 2 (SSTR2) with its agonists and antagonist. A 3D QSAR study based on comparative molecular field analysis and comparative molecular similarity analysis (CoMSIA) discerned that a SSTR2 ligand with electronegative, less-bulkier, and hydrogen atom donating/accepting substitutions is important for their biological activity. A conceptual density functional theory (DFT) study was followed to study the chemical behavior of the ligands based on the molecular descriptors derived using the Fukui's molecular orbital theory. We have performed molecular dynamics simulations of receptor-ligand complexes for 100 ns to analyze the dynamic stability of the backbone Cα atoms of the receptor and strength and approachability of the receptor-ligand complex. The findings of this study could be efficacious in the further studies understanding intricate structural features of the somatostatin receptors and in discovering novel subtype-specific ligands with higher affinity. Communicated by Ramaswamy H. Sarma.

3D-QSAR studies on indole and 7-azoindole derivatives as ROCK-2 inhibitors: An integrative computational approach.

Rho Kinases (ROCK) has been found to regulate a wide range of fundamental cell functions such as contraction, motility, proliferation, and apoptosis. Recent experiments have defined new functions of ROCKs in cells, including centrosome positioning and cell-size regulation, which might contribute to various physiological and pathological states. In this study, we have performed pharmacophore modeling and 3D QSAR studies on a series of 36 indoles and 7-azoindoles derivatives as ROCK2 inhibitors to elucidate the structural variations with their inhibitory activities. Ligand based CoMFA and CoMSIA models were generated based on three different alignment methods such as systematic search, simulated annealing and pharmacophore. A total of 15 CoMFA models and 27 CoMSIA were generated using different alignments. One model from each alignment is selected based on the statistical values. Contour maps of the selected models were compared, analysed and reported. The 3D QSAR study revealed that electro positive group linked to the methoxy-benzene ring position of the structure will enhance the biological activity and bulkier substitutions are preferred in the methyl dihydroindole region. Also, it is found that the hydrogen bond donor substituted at the R1 position enhances the inhibitory activity. In future, this study would give proper guidelines to further enhance the activity of novel inhibitors for ROCK2.

Theoretical analysis of somatostatin receptor 5 with antagonists and agonists for the treatment of neuroendocrine tumors.

We report on SSTR5 receptor modeling and its interaction with reported antagonist and agonist molecules. Modeling of the SSTR5 receptor was carried out using multiple templates with the aim of improving the precision of the generated models. The selective SSTR5 antagonists, agonists and native somatostatin SRIF-14 were employed to propose the binding site of SSTR5 and to identify the critical residues involved in the interaction of the receptor with other molecules. Residues Q2.63, D3.32, Q3.36, C186, Y7.34 and Y7.42 were found to be highly significant for their strong interaction with the receptor. SSTR5 antagonists were utilized to perform a 3D quantitative structure-activity relationship study. A comparative molecular field analysis (CoMFA) was conducted using two different alignment schemes, namely the ligand-based and receptor-based alignment methods. The best statistical results were obtained for ligand-based ([Formula: see text], [Formula: see text] = 0.988, noc = 4) and receptor-guided methods (docked mode 1:[Formula: see text], [Formula: see text], noc = 5), (docked mode 2:[Formula: see text] = 0.555, [Formula: see text], noc = 5). Based on CoMFA contour maps, an electropositive substitution at [Formula: see text], [Formula: see text] and [Formula: see text] position and bulky group at [Formula: see text] position are important in enhancing molecular activity.

Molecular Modeling Study on Diazine Indole Acetic Acid Derivatives for CRTH2 Inhibitory Activity.

In the present work, molecular modeling studies have been reported on a series of diazine indole acetic acid derivatives to analyze the structure-activity relationship studies of CRTH2 using fragment (Topomer CoMFA and HQSAR) and field (CoMFA and CoMSIA) based QSAR methods. Twenty-six compounds were used as a training set to establish the model, and six compounds were used as a test set to validate the model. The generated models exhibited good statistical results such as correlation coefficient (r2) and the cross-validated correlation coefficient (q2). Topomer CoMFA analysis yielded the q2 of 0.610 and r2 of 0.981. HQSAR model generated using bond and connectivity as fragment distinction and 3-6 as fragment size has the q2 value of 0.707 and conventional r2 value of 0.892 with five components. CoMFA model was assessed by cross-validated q2 value of 0.543 and r2 value of 0.901 with steric and electrostatic fields. CoMSIA model generated using steric, hydrophobic and donor fields with q2 value of 0.550 and r2 value of 0.888 was found to be the optimal model among the various models generated. The contour maps were generated to analyze the important structural features that regulate their inhibitory potency. From the result of contour maps we have suggested the critical sites for chemical modification which will be useful in designing potent compounds with improved activity.

Computational Analysis of CRTh2 receptor antagonist: A Ligand-based CoMFA and CoMSIA approach.

CRTh2 receptor is an important mediator of inflammatory effects and has attracted much attention as a therapeutic target for the treatment of conditions such as asthma, COPD, allergic rhinitis and atopic dermatitis. In pursuit of better CRTh2 receptor antagonist agents, 3D-QSAR studies were performed on a series of 2-(2-(benzylthio)-1H-benzo[d]imidazol-1-yl) acetic acids. There is no crystal structure information available on this protein; hence in this work, ligand-based comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were performed by atom by atom matching alignment using systematic search and simulated annealing methods. The 3D-QSAR models were generated with 10 different combinations of test and training set molecules, since the robustness and predictive ability of the model is very important. We have generated 20 models for CoMFA and 100 models for CoMSIA based on two different alignments. Each model was validated with statistical cut off values such as q(2)>0.4, r(2)>0.5 and r(2)pred>0.5. Based on better q(2) and r(2)pred values, the best predictions were obtained for the CoMFA (model 5 q(2)=0.488, r(2)pred=0.732), and CoMSIA (model 45 q(2)=0.525, r(2)pred=0.883) from systematic search conformation alignment. The high correlation between the cross-validated/predicted and experimental activities of a test set revealed that the CoMFA and CoMSIA models were robust. Statistical parameters from the generated QSAR models indicated the data is well fitted and have high predictive ability. The generated models suggest that steric, electrostatic, hydrophobic, H-bond donor and acceptor parameters are important for activity. Our study serves as a guide for further experimental investigations on the synthesis of new CRTh2 antagonist.