Discovery of new thymol-3,4-disubstituted thiazole hybrids as dual COX-2/5-LOX inhibitors with in vivo proof.

New thymol-3,4-disubstitutedthiazole hybrids were synthesised as dual COX-2/5-LOX inhibitors. Compounds 6b, 6d, 6e, and 6f displayed in vitro inhibitory activity against COX-2 (IC50= 0.037, 0.042, 0.046, and 0.039 µM) nearly equal to celecoxib (IC50= 0.045 µM). 6b, 6d, and 6f showed SI (379, 341, and 374, respectively) higher than that of celecoxib (327). 6a-l elicited in vitro 5-LOX inhibitory activity higher than quercetin. 6a-f, 6i-l, 7a, and 7c possessed in vivo inhibition of formalin induced paw edoema higher than celecoxib. 6a, 6b, 6f, 6h-l, and 7b showed gastrointestinal safety profile as celecoxib and diclofenac sodium in the population of fasted rats. Induced fit docking and molecular dynamics simulation predicted good fitting of 6b and 6f without changing the packing and globularity of the apo protein. In conclusion, 6b and 6f achieved the target goal as multitarget inhibitors of inflammation.

The Nitrofuran-Warhead-Equipped Spirocyclic Azetidines Show Excellent Activity against Mycobacterium tuberculosis.

A series of 21 new 7'H-spiro[azetidine-3,5'-furo [3,4-d]pyrimidine]s substituted at the pyrimidine ring second position were synthesized. The compounds showed high antibacterial in vitro activity against M. tuberculosis. Two compounds had lower minimum inhibitory concentrations against Mtb (H37Rv strain) compared with isoniazid. The novel spirocyclic scaffold shows excellent properties for anti-tuberculosis drug development.

Structural Features Affecting the Interactions and Transportability of LAT1-Targeted Phenylalanine Drug Conjugates.

L-type amino acid transporter 1 (LAT1) transfers essential amino acids across cell membranes. Owing to its predominant expression in the blood-brain barrier and tumor cells, LAT1 has been exploited for drug delivery and targeting to the central nervous system (CNS) and various cancers. Although the interactions of amino acids and their mimicking compounds with LAT1 have been extensively investigated, the specific structural features for an optimal drug scaffold have not yet been determined. Here, we evaluated a series of LAT1-targeted drug-phenylalanine conjugates (ligands) by determining their uptake rates by in vitro studies and investigating their interaction with LAT1 via induced-fit docking. Combining the experimental and computational data, we concluded that although LAT1 can accommodate various types of structures, smaller compounds are preferred. As the ligand size increased, its flexibility became more crucial in determining the compound's transportability and interactions. Compounds with linear or planar structures exhibited reduced uptake; those with rigid lipophilic structures lacked interactions and likely utilized other transport mechanisms for cellular entry. Introducing polar groups between aromatic structures enhanced interactions. Interestingly, compounds with a carbamate bond in the aromatic ring's para-position displayed very good transport efficiencies for the larger compounds. Compared to the ester bond, the corresponding amide bond had superior hydrogen bond acceptor properties and increased interactions. A reverse amide bond was less favorable than a direct amide bond for interactions with LAT1. The present information can be applied broadly to design appropriate CNS or antineoplastic drug candidates with a prodrug strategy and to discover novel LAT1 inhibitors used either as direct or adjuvant cancer therapy.

A new series of thiazole-hydrazone hybrids for Akt-targeted therapy of non-small cell lung cancer.

In an attempt to identify potent antitumor agents for the fight against non-small cell lung cancer, new thiazolyl hydrazones (2a-n) were synthesized and examined for their in vitro cytotoxic effects on A549 human lung adenocarcinoma and L929 mouse embryonic fibroblast cells by means of the MTT assay. Furthermore, the effects of the most potent anticancer agents on apoptosis and Akt inhibition were investigated. 2-[2-((Isoquinolin-5-yl)methylene)hydrazinyl]-4-(4-methylsulfonylphenyl)thiazole (2k) (IC50 = 1.43 ± 0.12 µM) and 2-[2-((isoquinolin-5-yl)methylene)hydrazinyl]-4-(1,3-benzodioxol-5-yl)thiazole (2l) (IC50 = 1.75 ± 0.07 µM) displayed more pronounced anticancer activity than cisplatin (IC50 = 3.90 ± 0.10 µM) on A549 cell lines; 2-[2-((isoquinolin-5-yl)methylene)hydrazinyl]-4-(4-methoxyphenyl)thiazole (2j) (IC50 = 3.93 ± 0.06 µM) showed anticancer activity close to cisplatin. These compounds were found to induce apoptosis in A549 cells. Compound 2j (IC50 = 3.55 ± 0.64 µM) showed stronger Akt inhibitory activity than GSK690693 (IC50 = 4.93 ± 0.06 µM), while compounds 2k and 2l did not cause Akt inhibition at IC50 concentrations (1.43 and 1.75 µM, respectively). To comprehensively elucidate the binding pose of compound 2j and to provide a detailed understanding on the ligand' binding mechanism, induced-fit docking calculations were also conducted. Both in vitro and in silico studies suggest that compound 2j shows its cytotoxic and apoptotic effects on A549 cell lines via Akt inhibition. However, it is understood that compounds 2k and 2l exert their strong anticancer effects on A549 cells through different pathways.

Virtual Screening Strategy and In Vitro Tests to Identify New Inhibitors of the Immunoproteasome.

Immunoproteasome inhibition is a promising strategy for the treatment of hematological malignancies, autoimmune diseases, and inflammatory diseases. The design of non-covalent inhibitors of the immunoproteasome ß1i/ß5i catalytic subunits could be a novel approach to avoid the drawbacks of the known covalent inhibitors, such as toxicity due to off-target binding. In this work, we report the biological evaluation of thirty-four compounds selected from a commercially available collection. These hit compounds are the outcomes of a virtual screening strategy including a dynamic pharmacophore modeling approach onto the ß1i subunit and a pharmacophore/docking approach onto the ß5i subunit. The computational studies were first followed by in vitro enzymatic assays at 100 µM. Only compounds capable of inhibiting the enzymatic activity by more than 50% were characterized in detail using Tian continuous assays, determining the dissociation constant (Ki) of the non-covalent complex where Ki is also the measure of the binding affinity. Seven out of thirty-four hits showed to inhibit ß1i and/or ß5i subunit. Compound 3 is the most active on the ß1i subunit with Ki = 11.84 ± 1.63 µM, and compound 17 showed Ki = 12.50 ± 0.77 µM on the ß5i subunit. Compound 2 showed inhibitory activity on both subunits (Ki = 12.53 ± 0.18 and Ki = 31.95 ± 0.81 on the ß1i subunit and ß5i subunit, respectively). The induced fit docking analysis revealed interactions with Thr1 and Phe31 of ß1i subunit and that represent new key residues as reported in our previous work. Onto ß5i subunit, it interacts with the key residues Thr1, Thr21, and Tyr169. This last hit compound identified represents an interesting starting point for further optimization of ß1i/ß5i dual inhibitors of the immunoproteasome.

Isomeric Activity Cliffs-A Case Study for Fluorine Substitution of Aminergic G Protein-Coupled Receptor Ligands.

Currently, G protein-coupled receptors (GPCRs) constitute a significant group of membrane-bound receptors representing more than 30% of therapeutic targets. Fluorine is commonly used in designing highly active biological compounds, as evidenced by the steadily increasing number of drugs by the Food and Drug Administration (FDA). Herein, we identified and analyzed 898 target-based F-containing isomeric analog sets for SAR analysis in the ChEMBL database-FiSAR sets active against 33 different aminergic GPCRs comprising a total of 2163 fluorinated (1201 unique) compounds. We found 30 FiSAR sets contain activity cliffs (ACs), defined as pairs of structurally similar compounds showing significant differences in affinity (≥50-fold change), where the change of fluorine position may lead up to a 1300-fold change in potency. The analysis of matched molecular pair (MMP) networks indicated that the fluorination of aromatic rings showed no clear trend toward a positive or negative effect on affinity. Additionally, we propose an in silico workflow (including induced-fit docking, molecular dynamics, quantum polarized ligand docking, and binding free energy calculations based on the Generalized-Born Surface-Area (GBSA) model) to score the fluorine positions in the molecule.

Identifying Dopamine D3 Receptor Ligands through Virtual Screening and Exploring the Binding Modes of Hit Compounds.

The dopamine D3 receptor (D3R) is an important central nervous system target for treating various neurological diseases. D3R antagonists modulate the improvement of psychostimulant addiction and relapse, while D3R agonists can enhance the response to dopaminergic stimulation and have potential applications in treating Parkinson's disease, which highlights the importance of identifying novel D3R ligands. Therefore, we performed auto dock Vina-based virtual screening and D3R-binding-affinity assays to identify human D3R ligands with diverse structures. All molecules in the ChemDiv library (>1,500,000) were narrowed down to a final set of 37 molecules for the binding assays. Twenty-seven compounds exhibited over 50% inhibition of D3R at a concentration of 10 µM, and 23 compounds exhibited over 70% D3R inhibition at a concentration of 10 µM. Thirteen compounds exhibited over 80% inhibition of D3R at a concentration of 10 µM and the IC50 values were measured. The IC50 values of the five compounds with the highest D3R-inhibition rates ranged from 0.97 µM to 1.49 µM. These hit compounds exhibited good structural diversity, which prompted us to investigate their D3R-binding modes. After trial and error, we combined unbiased molecular dynamics simulation (MD) and molecular mechanics generalized Born surface area (MM/GBSA) binding free-energy calculations with the reported protein−ligand-binding pose prediction method using induced-fit docking (IFD) and binding pose metadynamics (BPMD) simulations into a self-consistent and computationally efficient method for predicting and verifying the binding poses of the hit ligands to D3R. Using this IFD-BPMD-MD-MM/GBSA method, we obtained more accurate and reliable D3R−ligand-binding poses than were obtained using the reported IFD-BPMD method. This IFD-BPMD-MD-MM/GBSA method provides a novel paradigm and reference for predicting and validating other protein−ligand binding poses.

Tuning the Biological Activity of PI3Kδ Inhibitor by the Introduction of a Fluorine Atom Using the Computational Workflow.

As a member of the class I PI3K family, phosphoinositide 3-kinase δ (PI3Kδ) is an important signaling biomolecule that controls immune cell differentiation, proliferation, migration, and survival. It also represents a potential and promising therapeutic approach for the management of numerous inflammatory and autoimmune diseases. We designed and assessed the biological activity of new fluorinated analogues of CPL302415, taking into account the therapeutic potential of our selective PI3K inhibitor and fluorine introduction as one of the most frequently used modifications of a lead compound to further improve its biological activity. In this paper, we compare and evaluate the accuracy of our previously described and validated in silico workflow with that of the standard (rigid) molecular docking approach. The findings demonstrated that a properly fitted catalytic (binding) pocket for our chemical cores at the induced-fit docking (IFD) and molecular dynamics (MD) stages, along with QM-derived atomic charges, can be used for activity prediction to better distinguish between active and inactive molecules. Moreover, the standard approach seems to be insufficient to score the halogenated derivatives due to the fixed atomic charges, which do not consider the response and indictive effects caused by fluorine. The proposed computational workflow provides a computational tool for the rational design of novel halogenated drugs.

Design, synthesis, molecular docking study and molecular dynamics simulation of new coumarin-pyrimidine hybrid compounds having anticancer and antidiabetic activity.

Coumarin-pyrimidine hybrid compounds were synthesized by condensation reaction of α,ß-unsaturated ketones of 6-acetyl-5-hydroxy-4-methylcoumarin with guanidine. The reaction yields were of 42-62%. The antidiabetic and anticancer activities of these compounds were examined. These compounds displayed low toxicity to two cancer cell lines (including KB and HepG2 ones), but exhibited remarkably active against α-amylase with IC50 values of 102.32 ± 1.15 µM to 249.52 ± 1.14 µM and against α-glucosidase with IC50 values of 52.16 ± 1.12 µM to 184.52 ± 1.15 µM. Amongst these compounds, 6c was the best inhibitory activity against α-amylase, and 6f had the highest activity against α-glucosidase. The kinetics of inhibitor 6f was competitive α-glucosidase inhibitor property. ADMET predictions showed that almost all synthesized compounds exhibited drug-like activity. IFD and MD simulations were carried out on enzymes 4W93 and 5NN8 to elucidate inhibitory potential of 6c and 6f against tested enzymes. The binding free energy calculation by MM-GBSA approach showed that Coulomb, lipophilic and van der Waals energy terms are major contributors for the inhibitor binding. Molecular dynamics simulations in water solvent system were carried out for the 6f/5NN8 complex to elucidate the variability of active interactions between ligand 6f and active pockets of this enzyme.

Structure-based pharmacophore modeling, virtual screening approaches to identifying the potent hepatitis C viral protease and polymerase novel inhibitors.

Hepatitis C is an infectious disease that leads to acute and chronic liver illnesses. Currently, there are no effective vaccines against this deadly virus. Direct acting antiviral (DAA) drugs are given in the combination with ribavirin and pegylated interferon which lead to adverse effects. Through in silico analysis, the structure-based docking study was performed against NS3/4A protease and NS5B polymerase proteins of HCV. In the current study, multiple e-pharmacophore-based virtual screening methods such as HTVS, SP, and XP were carried out to screen natural compounds and enamine databases. Our result outcomes revealed that CID AE-848/13196185 and CID AE-848/36959205 compounds show good binding interactions with protease protein. In addition, CID 15081408 and CID 173568 show better binding interactions with the polymerase protein. Further to validate the docking results, we performed molecular dynamics simulation for the top hit compounds bound with protease and polymerase proteins to illustrate conformational differences in the stability compared with the active site of the cocrystal inhibitor. Thus, the current study emphasizes these compounds could be an effective drug to treat HCV.

CIFDock: A novel CHARMM-based flexible receptor-flexible ligand docking protocol.

Docking studies play a critical role in the current workflow of drug discovery. However, limitations may often arise through factors including inadequate ligand sampling, a lack of protein flexibility, scoring function inadequacies (e.g., due to metals, co-factors, etc.), and difficulty in retaining explicit water molecules. Herein, we present a novel CHARMM-based induced fit docking (CIFDock) workflow that can circumvent these limitations by employing all-atom force fields coupled to enhanced sampling molecular dynamics procedures. Self-guided Langevin dynamics simulations are used to effectively sample relevant ligand conformations, side chain orientations, crystal water positions, and active site residue motion. Protein flexibility is further enhanced by dynamic sampling of side chain orientations using an expandable rotamer library. Steps in the procedure consisting of fixing individual components (e.g., the ligand) while sampling the other components (e.g., the residues in the active site of the protein) allow for the complex to adapt to conformational changes. Ultimately, all components of the complex-the protein, ligand, and waters-are sampled simultaneously and unrestrained with SGLD to capture any induced fit effects. This modular flexible docking procedure is automated using CHARMM scripting, interfaced with SLURM array processing, and parallelized to use the desired number of processors. We validated the CIFDock procedure by performing cross-docking studies using a data set comprised of 21 pharmaceutically relevant proteins. Five variants of the CHARMM-based SWISSDOCK scoring functions were created to quantify the results of the final generated poses. Results obtained were comparable to, or in some cases improved upon, commercial docking program data.

Discovery, synthesis and in combo studies of Schiff's bases as promising dipeptidyl peptidase-IV inhibitors.

Diabetes mellitus is a main global health apprehension. Macrovascular illnesses, neuropathy, retinopathy, and nephropathy are considered some of its severe hitches. Gliptins are a group of hypoglycemic agents that inhibit dipeptidyl peptidase-IV (DPP-IV) enzyme and support blood glucose-lowering effect of incretins. In the current research, synthesis, characterization, docking, and biological evaluation of fourteen Schiff's bases 5a-f and 9a-h were carried out. Compound 9f revealed the best in vitro anti-DPP-IV activity of 35.7% inhibition at a concentration of 100 µM. Compounds 9c and 9f with the highest in vitro DPP-IV inhibition were subjected to the in vivo glucose-lowering test using vildagliptin as a positive inhibitor. Vildagliptin, 9c, and 9f showed significant reduction in the blood glucose levels of the treated mice after 30 min of glucose administration. Moreover, induced fit docking showed that these derivatives accommodated the enzyme binding site with comparable docking scores. Schiff's bases can serve as promising lead for the development of new DPP-IV inhibitors.

Design, Synthesis, α-Amylase/α-Glucosidase Inhibition Assay, Induced Fit Docking Study of New Hybrid Compounds Containing 4H-Pyrano[2,3-d]pyrimidine, 1H-1,2,3-Triazole and D-Glucose Components.

In this study, the click chemistry between N-propargyl derivatives of substituted 4H-pyrano[2,3-d]pyrimidines and tetra-O-acetyl-α-d-glucopyranosyl azide carried out under catalytic conditions using catalyst CuI@Montmorillonite and additive N,N-diisopropylethylamine (DIPEA). The yields of obtained hybrid compounds having 4H-pyrano[2,3-d]pyrimidine connected to 1H-1,2,3-triazole rings were about 85-94 %. All these synthesized hybrid compounds were examined for in vitro α-amylase (with IC50 values in the range of 103.63±1.13 µM to 295.45±1.11 µM) and α-glucosidase (with IC50 values in the range of 45.63±1.14 µM to 184.52±1.15) inhibitory activity. Amongst this series, ethyl ester 8m showed the best inhibitory activity against α-amylase with IC50 of 103.63±1.13 µM, while ethyl ester 8t exhibited the highest activity against α-glucosidase with IC50 of 45.63±1.14 µM. The kinetics of the inhibition of compound 8t showed the competitive α-glucosidase inhibitor property of this compound. Furthermore, the most potent compounds had any cytotoxicity against human normal cells. Induced fit docking and molecular dynamics simulation calculations indicated that the inhibition potential compounds 8m and 8t had the active interactions with the residues in receptors of corresponding tested enzymes. The calculated binding free energy from MM-GBSA approach showed that the major energy components contributed to the active binding of these studied inhibitors, including Coulomb, lipophilic and van der Waals energy. Further, 300 ns MD simulation showed that studied ligand-protein complexes were stable and indicated the structural observations into mode of binding in these complexes.

The Development of Pharmacophore Models for the Search of New Natural Inhibitors of SARS-CoV-2 Spike RBD-ACE2 Binding Interface.

To date, some succeeding variants of SARS-CoV-2 have become more contagious. This virus is known to enter human cells by binding the receptor-binding domain (RBD) of spike protein with the angiotensin-converting enzyme 2 (ACE2), the latter being a membrane protein that regulates the renin-angiotensin system. Since the host cell receptor plays a critical role in viral entry, inhibition of the RBD-ACE2 complex is a promising strategy for preventing COVID-19 infection. In the present communication, we propose and utilize an approach based on the generation of a complex of pharmacophore models and subsequent Induced Fit Docking (IFD) to identify potential inhibitors of the main binding sites of the Omicron SARS-CoV-2 RBD(S1)-ACE2 complex (PDB ID: 7T9L) among a number of natural products of various types and origins. Several natural compounds have been found to provide a high affinity for the receptor of interest. It is expected that the present results will stimulate further research aimed at the development of specialized drugs against this virus.

Computational-Based Discovery of the Anti-Cancer Activities of Pyrrole-Based Compounds Targeting the Colchicine-Binding Site of Tubulin.

Despite the tubulin-binding agents (TBAs) that are widely used in the clinic for cancer therapy, tumor resistance to TBAs (both inherited and acquired) significantly impairs their effectiveness, thereby decreasing overall survival (OS) and progression-free survival (PFS) rates, especially for the patients with metastatic, recurrent, and unresectable forms of the disease. Therefore, the development of novel effective drugs interfering with the microtubules' dynamic state remains a big challenge in current oncology. We report here about the novel ethyl 2-amino-1-(furan-2-carboxamido)-5-(2-aryl/tert-butyl-2-oxoethylidene)-4-oxo-4,5-dihydro-1H-pyrrole-3-carboxylates (EAPCs) exhibiting potent anti-cancer activities against the breast and lung cancer cell lines in vitro. This was due to their ability to inhibit tubulin polymerization and induce cell cycle arrest in M-phase. As an outcome, the EAPC-treated cancer cells exhibited a significant increase in apoptosis, which was evidenced by the expression of cleaved forms of PARP, caspase-3, and increased numbers of Annexin-V-positive cells. By using the in silico molecular modeling methods (e.g., induced-fit docking, binding metadynamics, and unbiased molecular dynamics), we found that EAPC-67 and -70 preferentially bind to the colchicine-binding site of tubulin. Lastly, we have shown that the EAPCs indicated above and colchicine utilizes a similar molecular mechanism to inhibit tubulin polymerization via targeting the T7 loop in the ß-chain of tubulin, thereby preventing the conformational changes in the tubulin dimers required for their polymerization. Collectively, we identified the novel and potent TBAs that bind to the colchicine-binding site and disrupt the microtubule network. As a result of these events, the compounds induced a robust cell cycle arrest in M-phase and exhibited potent pro-apoptotic activities against the epithelial cancer cell lines in vitro.

Structure-based virtual screening and computational study towards identification of novel inhibitors of hypoxanthine-guanine phosphoribosyltransferase of Trypanosoma cruzi.

Hypoxanthine-guanine phosphoribosyltransferase (HGPRT) is the key regulatory enzyme of the purine salvage pathway present in the members of trypanosomatids. The parasite solely depends on this pathway for the synthesis of nucleotides due to the absence of the de novo pathway. This study intends to identify putative inhibitors towards Trypanosoma cruzi HGPRT (TcHGPRT). Initial virtual screening was performed with substructures of phosphoribosyl pyrophosphate (PRPP), an original substrate of HGPRT. Twenty compounds that had greater binding energy than the substrate was treated as hits and was further screened and narrowed down through induced fit docking which resulted in top five compounds which was distinguished into two groups based on the ligand occupancy within the PRPP binding site of TcHGPRT. Group-I compounds (PubChem CID 130316561 and 134978234) are analogous to PRPP structure with greater occupancy, were preferred over Group-II compounds which had lesser occupancy than the substrate. However, one compound (22404820) among Group II was chosen for further analysis considering its significant electrostatic interactions. Molecular docking studies revealed the requirement of an electronegative moiety like phosphate group to be present in the ligand due to the presence of metal ions in the substrate binding site. The three chosen compounds along with PRPP were subjected to molecular dynamics analysis, which indicated a strong presence of electrostatic interaction. Considering the dynamic stability of interactions as well as pharmacological properties of ligands based on absorption, distribution, metabolism, excretion prediction, Group-I compounds were selected as lead compounds and were subjected to molecular electrostatic potential analysis to determine the charge distribution of the compound. The overall analysis thus suggests both 130316561 and 134978234 can be used as TcHGPRT inhibitors. Furthermore, these computational results emphasize the requirement of phosphorylated ligands which are essential in mediating electrostatic interactions and to compete with the binding affinity of the original substrate.

Exploring aromatic cage flexibility of the histone methyllysine reader protein Spindlin1 and its impact on binding mode prediction: an in silico study.

Some of the main challenges faced in drug discovery are pocket flexibility and binding mode prediction. In this work, we explored the aromatic cage flexibility of the histone methyllysine reader protein Spindlin1 and its impact on binding mode prediction by means of in silico approaches. We first investigated the Spindlin1 aromatic cage plasticity by analyzing the available crystal structures and through molecular dynamic simulations. Then we assessed the ability of rigid docking and flexible docking to rightly reproduce the binding mode of a known ligand into Spindlin1, as an example of a reader protein displaying flexibility in the binding pocket. The ability of induced fit docking was further probed to test if the right ligand binding mode could be obtained through flexible docking regardless of the initial protein conformation. Finally, the stability of generated docking poses was verified by molecular dynamic simulations. Accurate binding mode prediction was obtained showing that the herein reported approach is a highly promising combination of in silico methods able to rightly predict the binding mode of small molecule ligands in flexible binding pockets, such as those observed in some reader proteins.

Ligand-Receptor Interactions and Machine Learning in GCGR and GLP-1R Drug Discovery.

The large amount of data that has been collected so far for G protein-coupled receptors requires machine learning (ML) approaches to fully exploit its potential. Our previous ML model based on gradient boosting used for prediction of drug affinity and selectivity for a receptor subtype was compared with explicit information on ligand-receptor interactions from induced-fit docking. Both methods have proved their usefulness in drug response predictions. Yet, their successful combination still requires allosteric/orthosteric assignment of ligands from datasets. Our ligand datasets included activities of two members of the secretin receptor family: GCGR and GLP-1R. Simultaneous activation of two or three receptors of this family by dual or triple agonists is not a typical kind of information included in compound databases. A precise allosteric/orthosteric ligand assignment requires a continuous update based on new structural and biological data. This data incompleteness remains the main obstacle for current ML methods applied to class B GPCR drug discovery. Even so, for these two class B receptors, our ligand-based ML model demonstrated high accuracy (5-fold cross-validation Q2 > 0.63 and Q2 > 0.67 for GLP-1R and GCGR, respectively). In addition, we performed a ligand annotation using recent cryogenic-electron microscopy (cryo-EM) and X-ray crystallographic data on small-molecule complexes of GCGR and GLP-1R. As a result, we assigned GLP-1R and GCGR actives deposited in ChEMBL to four small-molecule binding sites occupied by positive and negative allosteric modulators and a full agonist. Annotated compounds were added to our recently released repository of GPCR data.

Immunoproteasome and Non-Covalent Inhibition: Exploration by Advanced Molecular Dynamics and Docking Methods.

The selective inhibition of immunoproteasome is a valuable strategy to treat autoimmune, inflammatory diseases, and hematologic malignancies. Recently, a new series of amide derivatives as non-covalent inhibitors of the ß1i subunit with Ki values in the low/submicromolar ranges have been identified. Here, we investigated the binding mechanism of the most potent and selective inhibitor, N-benzyl-2-(2-oxopyridin-1(2H)-yl)propanamide (1), to elucidate the steps from the ligand entrance into the binding pocket to the ligand-induced conformational changes. We carried out a total of 400 ns of MD-binding analyses, followed by 200 ns of plain MD. The trajectories clustering allowed identifying three representative poses evidencing new key interactions with Phe31 and Lys33 together in a flipped orientation of a representative pose. Further, Binding Pose MetaDynamics (BPMD) studies were performed to evaluate the binding stability, comparing 1 with four other inhibitors of the ß1i subunit: N-benzyl-2-(2-oxopyridin-1(2H)-yl)acetamide (2), N-cyclohexyl-3-(2-oxopyridin-1(2H)-yl)propenamide (3), N-butyl-3-(2-oxopyridin-1(2H)-yl)propanamide (4), and (S)-2-(2-oxopyridin-1(2H)-yl)-N,4-diphenylbutanamide (5). The obtained results in terms of free binding energy were consistent with the experimental values of inhibition, confirming 1 as a lead compound of this series. The adopted methods provided a full dynamic description of the binding events, and the information obtained could be exploited for the rational design of new and more active inhibitors.

Identification and neuroprotective evaluation of a potential c-Jun N-terminal kinase 3 inhibitor through structure-based virtual screening and in-vitro assay.

The c-Jun N-terminal kinase 3 (JNK3) signaling cascade is activated during cerebral ischemia leading to neuronal damage. The present study was carried out to identify and evaluate novel JNK3 inhibitors using in-silico and in-vitro approach. A total of 380 JNK3 inhibitors belonging to different organic groups was collected from the previously reported literature. These molecules were used to generate a pharmacophore model. This model was used to screen a chemical database (SPECS) to identify newer molecules with similar chemical features. The top 1000 hits molecules were then docked against the JNK3 enzyme coordinate following GLIDE rigid receptor docking (RRD) protocol. Best posed molecules of RRD were used during induced-fit docking (IFD), allowing receptor flexibility. Other computational predictions such as binding free energy, electronic configuration and ADME/tox were also calculated. Inferences from the best pharmacophore model suggested that, in order to have specific JNK3 inhibitory activity, the molecules must possess one H-bond donor, two hydrophobic and two ring features. Docking studies suggested that the main interaction between lead molecules and JNK3 enzyme consisted of hydrogen bond interaction with methionine 149 of the hinge region. It was also observed that the molecule with better MM-GBSA dG binding free energy, had greater correlation with JNK3 inhibition. Lead molecule (AJ-292-42151532) with the highest binding free energy (dG = 106.8 Kcal/mol) showed better efficacy than the SP600125 (reference JNK3 inhibitor) during cell-free JNK3 kinase assay (IC50 = 58.17 nM) and cell-based neuroprotective assay (EC50 = 7.5 µM).